Acessibilidade / Reportar erro

Monitoring studies should consider temporal variability to reveal relations between cyanobacterial abundance and environmental variables

Abstracts

One of the main goals of monitoring cyanobacteria blooms in aquatic environments is to reveal the relationship between cyanobacterial abundance and environmental variables. Studies typically correlate data that were simultaneously sampled. However, samplings occur sparsely over time and may not reveal the short-term responses of cyanobacterial abundance to environmental changes. In this study, we tested the hypothesis that stronger cyanobacteria x environment relationships in monitoring are found when the temporal variability of sampling points is incorporated in the statistical analyses. To this end, we investigated relationships between cyanobacteria and seven environmental variables that were sampled twice yearly for three years across 11 reservoirs, and data from an intensive monitoring in one of these reservoirs. Poor correlations were obtained when correlating data simultaneously sampled. In fact, the 'highly recurrent' role of phosphorus in cyanobacteria blooms is not properly observed in all sampling periods. On the other hand, the strongest correlation values for the total phosphorus x cyanobacteria relationship were observed when we used the variation of sampling points. We have also shown that environment variables better explain cyanobacteria when a time lag is considered. We conclude that, in cyanobacteria monitoring, the best approach to reveal determinants of cyanobacteria blooms is to consider environmental variability.

Harmful algae; reservoir; management; limnological variables; temporal variation; Neotropics


Um dos principais objetivos de monitorar as florações de cianobactérias em ambientes aquáticos é revelar a relação entre a abundância de cianobactérias e as variáveis ​​ambientais. Esses estudos tipicamente correlacionam dados coletados simultaneamente. No entanto, as amostragens ocorrem esparsamente ao longo do tempo e não são capazes de revelar respostas de curto prazo na abundância de cianobactérias frente às mudanças ambientais. Neste estudo, testamos a hipótese de que relações mais fortes entre cianobactérias e variáveis ambientais no monitoramento de cianobactérias são encontradas quando a variabilidade temporal dos pontos de amostragem é incorporada nas análises estatísticas. Para isso, nós investigamos as relações entre cianobactérias e sete variáveis ​​ambientais que foram amostradas duas vezes por ano, durante três anos, em 11 reservatórios, e dados de um monitoramento intensivo em um desses reservatórios. Relações fracas foram encontradas correlacionando dados coletados simultaneamente. De fato, o papel "altamente recorrente" do fósforo nas florações de cianobactérias não foi devidamente observado em todos os períodos de amostragem. Por outro lado, coeficientes de correlação mais fortes entre fósforo total e cianobactérias foram observados quando usamos a variabilidade dos pontos de amostragem. Nós também mostramos que as variáveis ambientais explicam melhor a densidade de cianobactérias quando um atraso é considerado. Nós concluímos que, no monitoramento de cianobactérias, a melhor abordagem para revelar os determinantes das florações de cianobactérias é considerar a variabilidade ambiental.

algas nocivas; reservatório; gestão; variáveis limnológicas; variação temporal; Neotropical


INTRODUCTION

It is well established that increasing temperatures and nutrient concentrations in aquatic habitats are causing cyanobacterial blooms (Paerl and Huisman 2008PAERL HW AND HUISMAN J. 2008. Blooms Like It Hot. Science 320: 57-58., 2009PAERL HW AND HUISMAN J. 2009. Climate change: a catalyst for global expansion of harmful cyanobacterial blooms. Environ Microbiol Rep 1: 27-37.). For this reason, cyanobacterial toxin release in freshwater is a major environmental issue, causing serious concern for sanitation companies (Chorus and Bartram 1999CHORUS I AND BARTRAM J. 1999. Toxic cyanobacteria in water: a guide to their public health consequences, monitoring and management. London: E & FN SPON, 416 p.). In fact, cyanotoxins may cause the death of wildlife, affect livestock and result in negative effects to human health (Yoo et al. 1995YOO RS, CARMICHEL WW, HOEHN RC AND HRUDEY SE. 1995. Cyanobacterial (Blue-Green Algal) Toxins: A Resource Guide. AWWA Research Foundation and American Water Works Association, 229 p., Azevedo et al. 2002AZEVEDO SMFO, CARMICHAEL WW, JOCHIMSEN EM, RINEHART KL, LAU S, SAHW GR AND EAGLESHAM GK. 2002. Human intoxications by microcytins during renal dialysis treatment in Caruaru - Brazil. Toxicology 181: 441-446.).

Stability of the water column (Huszar et al. 2000HUSZAR VLM, SILVA LHS, MARINHO M, DOMINGOS P AND SANT'ANNA CL. 2000. Cyanoprokaryote assemblages in eight productive tropical Brazilian waters. Hydrobiologia 424: 67-77., Dantas et al. 2011DANTAS EW, MOURA AN AND BITTENCOURT-OLIVEIRA MC. 2011. Cyanobacterial blooms in stratified and destratified eutrophic reservoirs in semi-arid region of Brazil. An Acad Bras Cienc 83: 1327-1338.), high nutrient concentrations (mainly nitrogen and phosphorus) (Bouvy et al. 2000BOUVY M, FALCÃO D, MARINHO M, PAGANO M AND MOURA A. 2000. Occurrence of Cylindrospermopsis (cyanobacteria) in 39 Brazilian tropical reservoirs during the 1998 drought. Aquat Microb Ecol 23: 13-27.), high temperature ( Paerl and Huisman 2008PAERL HW AND HUISMAN J. 2008. Blooms Like It Hot. Science 320: 57-58., 2009PAERL HW AND HUISMAN J. 2009. Climate change: a catalyst for global expansion of harmful cyanobacterial blooms. Environ Microbiol Rep 1: 27-37.), and high luminosity (Chorus and Bartram 1999CHORUS I AND BARTRAM J. 1999. Toxic cyanobacteria in water: a guide to their public health consequences, monitoring and management. London: E & FN SPON, 416 p., Huszar et al. 2000HUSZAR VLM, SILVA LHS, MARINHO M, DOMINGOS P AND SANT'ANNA CL. 2000. Cyanoprokaryote assemblages in eight productive tropical Brazilian waters. Hydrobiologia 424: 67-77.) are the main factors that promote cyanobacterial blooms. Such evidence for the causes of cyanobacterial growth has been identified mainly in experiments (Paerl et al. 1985PAERL HW, BLAND PT, BOWLES ND AND HAIBACH ME. 1985. Adaptation to High-Intensity, Low-Wavelength Light among Surface Blooms of the Cyanobacterium Microcystis aeruginosa. Appl Environ Microb 49: 1046-1052., Coles and Jones 2000COLES JF AND JONES RC. 2000. Effect of temperature on photosynthesis-light response and growth of four phytoplankton species isolated from a tidal freshwater river. J Phycol 36: 7-16., Aubriot and Bonilla 2012AUBRIOT L AND BONILLA S. 2012. Rapid regulation of phosphate uptake in freshwater cyanobacterial blooms. Aquat Microb Ecol 67: 251-263.). Observational studies have also tried to relate environmental variables and cyanobacteria densities (Havens et al. 1998HAVENS KE, PHLIPS EJ, CICHRA MF AND LI B. 1998. Light availability as a possible regulator of cyanobacteria species in a shallow subtropical lake. Freshwater Biol 39: 547-556., Burford et al. 2007BURFORD MA, JOHNSON SA, COOK AJ, PACKER TV, TAYLOR BM AND TOWNSLEY ER. 2007. Correlations between watershed and reservoir characteristics, and algal blooms in subtropical reservoirs. Water Res 41: 4105-4114., Gomes et al. 2012GOMES LNL, OLIVEIRA SMAC, GIANI A AND VON SPERLING E. 2012. Association between biotic and abiotic parameters and the occurrence of cyanobacteria in a Brazilian reservoir. Environ Monit Assess 184: 4635-4645.). However, relationships between cyanobacteria densities or cyanotoxin concentrations with the variables described above are not always evident (Wicks and Thiel 1990WICKS RJ AND THIEL PG. 1990. Environmental factors affecting the production of peptide toxins in floating scums of the cyanobacterium Microcystis aeruginosa in a hypertrophic african reservoir. Environ Sci Technol 24: 1413-1418., Quesada and Fernández-Valiente 1996QUESADA A AND FERNÁNDEZ-VALIENTE E. 1996. Relationship between abundance of N2-fixing cyanobacteria and environmental features of Spanish rice fields. Microbial Ecol 32: 59-71., Vezie et al. 1998VEZIE C, BRIENT L, SIVONEN K, BERTRU G, LEFEUVRE JC AND SALKINOJA-SALONEN M. 1998. Variation of microcystin content of cyanobacterial blooms and isolated strains in Lake Grand-Lieu (France). Microbial Ecol 35: 126-135.).

Observational studies are based mainly on simultaneous samplings of cyanobacterial abundance and environmental variables (e.g. Dantas et al. 2011DANTAS EW, MOURA AN AND BITTENCOURT-OLIVEIRA MC. 2011. Cyanobacterial blooms in stratified and destratified eutrophic reservoirs in semi-arid region of Brazil. An Acad Bras Cienc 83: 1327-1338., Moura et al. 2011MOURA AN, DANTAS EW, OLIVEIRA HSB AND BITTENCOURT-OLIVEIRA MC. 2011. Vertical and temporal dynamics of cyanobacteria in the Carpina potable water reservoir in northeastern Brazil. Braz J Biol 71: 451-459.,Bittencourt-Oliveira et al. 2012BITTENCOURT-OLIVEIRA MC, DIAS SN, MOURA NA, CORDEIRO-ARAÚJO MK AND DANTAS EW. 2012. Seasonal dynamics of cyanobacteria in a eutrophic reservoir (Arcoverde) in a semi-arid region of Brazil. Braz J Biol 72: 533-544.). The implicit assumption in these studies is that snapshot samples in a certain sampling location are a good representation of the environmental conditions in this location in comparison to other locations. This assumption can be met if samplings occur continuously over time, given that cyanobacteria respond to short-term environmental variations. Another issue that can affect relationships between cyanobacteria and environmental variables is the fact that cyanobacteria comprise species that may respond differently to the environmental variations. For instance, picocyanobacteria and large cyanobacterial colonies fundamentally have different ecological restrictions (Huisman et al. 2005HUISMAN J, MATTHIJS HCP AND VISSER PM. 2005. Harmful Cyanobacteria. Netherlands: Springer, 241 p.). The presence of heterocyst, which determines Nitrogen fixation, is observed in only some cyanobacteria (Chorus and Bartram 1999CHORUS I AND BARTRAM J. 1999. Toxic cyanobacteria in water: a guide to their public health consequences, monitoring and management. London: E & FN SPON, 416 p.). However, cyanobacteria monitoring by legal agencies evaluate cyanobacteria as one biological group.

Moreover, cyanobacteria monitoring usually carries out simultaneous and long-term periodic samplings (LTPS) of cyanobacteria and environmental variables. For instance, the Environmental Institute of the state of Paraná (Brazil) carries out samplings of cyanobacterial abundance and environmental variables, only twice a year in reservoirs for water supply and hydroelectricity generation, with the main goal of monitoring cyanobacteria blooms and environmental variations (IAP 2009IAP - INSTITUTO AMBIENTAL DO PARANÁ. 2009. Monitoramento da qualidade das águas dos reservatórios do Estado do Paraná, no período de 2005 a 2008. Curitiba. Available in: www.iap.pr.gov.br, 120 p.
www.iap.pr.gov.br...
). Periodic samplings may not reveal the cyanobacterial responses to environmental variables, and correlating such data may not be the best strategy to investigate relationships between cyanobacteria and environmental variables. We argue this given that cyanobacteria may respond rapidly to environmental variations.

A better alternative for LTPS may be to represent how environmental variables vary over time across sampling points. To our knowledge, there have been no studies relating environmental and cyanobacterial variability. Here, we aimed to investigate the relationship between cyanobacteria and environmental variables, with the hypothesis that stronger relationships are found when the variability of sampling points is used as the sampling unit in statistical analyses, instead of raw long-term periodic sampling values. We dug further into the reasons why LTPS may not be a good strategy to correlate cyanobacterial abundance and environmental variables by investigating correlations with environmental variables sampled at the same time, and sampled days before cyanobacterial survey.

MATERIALS AND METHODS

We investigated the relationship between cyanobacteria and environmental variables using data sampled in 11 reservoirs that are used for water supply and hydroelectricity generation in the state of Paraná, Brazil (Table I). These reservoirs vary on the level of degradation (Andreoli and Carneiro 2005ANDREOLI CV AND CARNEIRO C. 2005. Gestão integrada de mananciais de abastecimento eutrofizados. Curitiba: Companhia de Saneamento do Paraná SANEPAR/FINEP, 494 p.), resulting in persistent toxic cyanobacterial blooms in "Irai" and "Alagados" reservoirs (Fernandes et al. 2005FERNANDES LF, WOSIACK AC, DOMINGUES L, PACHECO CV AND LAGOS PED. 2005. Cianobactérias e cianotoxinas. In: Andreoli CV and Carneiro C (Eds), Gestão integrada de mananciais de abastecimento eutrofizados. Curitiba: Sanepar-Finep, p. 369-388., 2011FERNANDES LF, GUTSEIT KC, WOJCIECHOWSKI J, LAGOS PED, XAVIER CF AND WOSIACK AC. 2011. Ecologia do fitoplâncton no reservatório Rio Verde. In: Cunha CLN, Carneiro C, Gobbi EF and Andreoli CV (Eds), Eutrofização em Reservatórios: Gestão Preventiva. Estudo Interdisciplinar na Bacia do Rio Verde, PR. Curitiba: Ed. da UFPR, p. 286-306., IAP 2009IAP - INSTITUTO AMBIENTAL DO PARANÁ. 2009. Monitoramento da qualidade das águas dos reservatórios do Estado do Paraná, no período de 2005 a 2008. Curitiba. Available in: www.iap.pr.gov.br, 120 p.
www.iap.pr.gov.br...
). Data were provided by the program "Monitoring of water quality in the reservoirs of the state of Paraná" coordinated by the Paraná State Environmental Institute (IAP 2009IAP - INSTITUTO AMBIENTAL DO PARANÁ. 2009. Monitoramento da qualidade das águas dos reservatórios do Estado do Paraná, no período de 2005 a 2008. Curitiba. Available in: www.iap.pr.gov.br, 120 p.
www.iap.pr.gov.br...
) and by "SANEPAR" (The Sanitation Company of the state of Paraná, Brazil). Samples in the 11 reservoirs were taken approximately twice a year (such as in LTPS) between 2005 and 2008 (hereafter referred to as "IAP sampling"). In 2008, some reservoirs were not sampled twice, so the number of samplings vary by reservoir (Table I). In addition, the "Alagados" reservoir was monitored from 2008 to 2014 in a higher frequency than the other reservoirs, when cyanobacterial abundance reached 20,000 cell.ml-1 (hereafter referred to as "Alagados sampling"). In this reservoir, environmental variables were measured once a month and cyanobacteria was sampled weekly.

TABLE I
Geographic location of reservoirs and number of samples taken in each reservoir from 2005 to 2008.

Samples of both cyanobacterial abundance (cells.mL-1) and environmental variables were taken in the sub-surface, in the lacustrine portion of the reservoirs. The environmental variables sampled for IAP sampling were: ammonium (NH4 +, mg.L-1), chemical oxygen demand (COD, mg.L-1), dissolved oxygen (DO, %), nitrate (NO3 -, mg.L-1), nitrite (NO2 2-, mg.L-1), transparency (m, using Secchi disk depth), and total phosphorus (TP, mg.L-1).

For Alagados, the sampled environmental variables were: water column depth (m), water temperature (ºC), pH, dissolved oxygen (DO, mg.L-1), ammonium (NH4 +, mg.L-1), nitrate (NO3 -, mg.L-1), nitrite (NO2 2-, mg.L-1), total phosphorus (TP, mg.L-1), soluble reactive phosphorus (SRP, mg.L-1), turbidity (NTU), conductivity (μS.cm-1), chemical oxygen demand (COD, mg.L-1) and biochemical oxygen demand (BOD, mg.L-1).

We used three approaches to correlate environmental variables and abundance of cyanobacteria in IAP sampling. First, we used data from simultaneous samplings of cyanobacteria and environmental variables, hereafter referred to as the simultaneous approach (SA). As for the approaches described below, the reservoir was used as the sampling unit in correlation analyses (N = 11). Given that we have several sampling periods for each reservoir (up to nine, see Table I), we calculated SA correlations for each sampling period in which all reservoirs were sampled (six samples during summer and winter from 2005 to 2008,Table I). Then we considered the mean correlation across all these temporal samples as the estimate of the relationship between cyanobacteria and environmental variables when they were simultaneously sampled.

Second, we used the temporal average of variables for each reservoir as values in the correlation analyses (N = 11), hereafter called the mean approach (MA). We believe that the mean value of each reservoir over sampling periods is a better representation of how the habitat in one location compares to another, and possibly how environmental processes explain cyanobacterial abundance. It could also be understood as a metric that captures the temporal variability of the reservoir as it relates the mean values over a time spam.

Third, we used the temporal variance of environmental variables and cyanobacterial abundance for each reservoir in the correlation analyses (N = 11), hereafter referred to as the variability approach (VA). We used this metric to represent how each sampling point varies over time. The interpretation of correlations in VA is different from SA and MA: a strong correlation in VA means that reservoirs with high environmental variations also have high variations in cyanobacterial abundance. In this case, a possible prediction is related to the cyanobacterial variability (not to the amount of cyanobacteria) explained by the variability of ecological processes. In SA and MA, in contrast, a strong correlation means that reservoirs with high values of a certain environmental variable also have high values of cyanobacterial abundance. We used variance in VA instead of other measures of dispersion because our goal was to represent overall dispersion and not a dispersion measure that is relative to the mean value (e.g. coefficient of variation).

We hypothesized that MA and VA would result in stronger cyanobacteria-environment correlations when compared to SA. Correlations generated by these three approaches can be compared by considering the magnitude of the coefficient. Usually, studies focus on significance values to evaluate if there is or there is not, an association between environmental data and cyanobacteria (e.g. Oliveira et al. 2014OLIVEIRA FHPC, CAPELA E, ARA ALS, MOREIRA CHP, LIRA OO, PADILHA MRF AND SHINOHARA NKS. 2014. Seasonal changes of water quality in a tropical shallow and eutrophic reservoir in the metropolitan region of Recife (Pernambuco-Brazil). An Acad Bras Cienc 86: 1863-1872.). It is clear that significance depends on the sampling effort (McBride 2005MCBRIDE GB. 2005. Using Statistical Methods for Water Quality Management: Issues, Problems and Solutions. New Jersey: J Wiley & Sons, 344 p.), and thus we used the same number of sampling units (i.e. 11 reservoirs) for all approaches. However, the significance of each relationship is not important here, since our goal is to reveal the stronger cyanobacteria-environment correlations, and not the significant ones.

In addition to not representing the ecological process affecting cyanobacteria, since algae may respond to previous environmental variation, the rationale of proposing MA and VA as an alternative for SA was also to unravel a common misinterpretation regarding ecological data. Absolute values of both environmental variables and cyanobacterial abundance vary over time in most ecosystems (Dantas et al. 2011DANTAS EW, MOURA AN AND BITTENCOURT-OLIVEIRA MC. 2011. Cyanobacterial blooms in stratified and destratified eutrophic reservoirs in semi-arid region of Brazil. An Acad Bras Cienc 83: 1327-1338., Gomes et al. 2012GOMES LNL, OLIVEIRA SMAC, GIANI A AND VON SPERLING E. 2012. Association between biotic and abiotic parameters and the occurrence of cyanobacteria in a Brazilian reservoir. Environ Monit Assess 184: 4635-4645.), highlighting the importance of temporally structured ecological processes. However, it is unclear how environmental determinism on cyanobacteria varies. If SA correlations have a large variation - for instance, in some sampling periods the correlation is strong, and in others, the correlation is weak - it also indicates that simultaneous sampling and correlating environmental variables and cyanobacteria is not a suitable approach, so conclusions will be uncertain. Even though the values of MA or VA correlations may be within the range of variation of SA correlations, these approaches indicate that the variability of environmental and cyanobacterial data are related. We argue that if the MA and VA correlations are higher than the mean SA correlations, then the real effect of environmental variables on cyanobacteria should consider the temporal variation of raw data, and not the variation of correlations.

We used both Pearson and Spearman correlations in order to identify whether one was better than the other for each approach used (SA, MA or VA). A higher Pearson correlation indicates linear relationships, and higher Spearman correlation may suggest non-linear relationships between cyanobacterial abundance and environmental variables. We believe that there are other methods to relate cyanobacteria and environmental variables, however, our goal was to compare correlation coefficients commonly used to associate cyanobacteria and environmental variables, but only to compare strategies of finding predictors of cyanobacteria (i.e. the approaches described above).

In the Alagados reservoir, cyanobacterial samplings made by the SANEPAR Company were intensive. Cyanobacterial abundance was available not only for the same sampling periods of environmental variables sampling, but also for one week after environmental variable sampling. Therefore, we correlated environmental variables with cyanobacterial abundance sampled simultaneously (N = 38 sampling periods spread from 2008 to 2014) and the same environmental variables with cyanobacteria sampled after c. one week ("time lag correlations"; N = 36 sampling periods spread from 2008 to 2014) in the Alagados reservoir. The rationale is that "time lag correlations" will be higher than simultaneous correlations, given that algae may respond to a previous environmental condition of the ecosystem. We used as response variables: both, total cyanobacteria abundance and the abundance of the most abundant taxa in blooms (in this case, Aphanizomenon sp.,Cylindrospermopsis sp., and Dolichospermumsp.). All analyses were done in the software STATISTICA 7.1 (StatSoft 2005STATSOFT, INC. 2005. STATISTICA (data analysis software system), version 7.1. www.statsoft.com.
www.statsoft.com...
).

RESULTS

IAP SAMPLING

MA was the approach that resulted in the highest correlations between cyanobacterial abundance and environmental variables (5 out of 7 for Pearson; and 4 out of 7 for Spearman) (Table II). As expected, SA was the worst approach, having the highest correlation values only for 1 out of 7 environmental variables for both Pearson and Spearman correlations (Table II; SA and MA had the highest, and similar correlations considering Nitrate).

TABLE II
Pearson (r) and Spearman (ρ) correlation coefficients between environmental variables and cyanobacterial density (cells mL-1) considering the following approaches (see methods): the simultaneous approach (SA), showing mean correlations (weakest/strongest); the mean approach (MA); and the variability approach (VA). Six samples were simultaneously carried out in all reservoirs, so SA values are based on the mean for six correlations. COD = chemical oxygen demand. The highest correlations for each variable are in bold.

COD showed the strongest Pearson correlation with cyanobacterial abundance in SA (r = 0.65; Table II). This was the only SA Pearson correlation that was stronger than MA or VA (Table II). For Spearman correlations, SA showed the strongest correlation for ammonium, although its absolute value could be considered weak (ρ = 0.13; Table II). On the other hand, correlation values were higher in MA than in SA for almost all variables (except ammonium, Table II). The highest Pearson correlations in MA were for total phosphorus and dissolved oxygen; and the highest Spearman correlations were for nitrate, transparency and COD. Correlations in VA were, in general, also higher than in SA (Table II). In VA, strong coefficients were found between cyanobacterial abundance and dissolved oxygen and total phosphorus by Pearson correlation; and nitrate, nitrite, and transparency by Spearman correlation (Table II).

We used the mean SA correlation over sampling periods to represent the correlation according to this approach. In this case, strong correlations were found in some sampling periods (e.g. r = 0.94 for the Pearson correlation with COD), whereas extremely low values were recorded for other sampling periods - nearly 0 for the Pearson correlation with phosphorus and Spearman correlation with ammonium (see Table II). Relatedly, although strong phosphorus-cyanobacteria relationships can be found in SA for some sampling periods (maximum values: r = 0.77 or ρ = 0.51, Table III), weak relationships can also be found (minimum values: r = 0.05 or ρ = 0.22, Table III), resulting in a low mean correlation. On the other hand, strong correlations were found for almost all variables (including for total phosphorus) when MA was used. For VA, a strong Pearson correlation between total phosphorus and cyanobacteria abundance was also evident.

TABLE III
Pearson and Spearman correlation coefficients between environmental variables and cyanobacterial density (cell.mL-1) considering the simultaneous correlation (SC) and "time lag" correlation (TL) (see methods, highest values are in bold for each comparison). COD = chemical oxygen demand, BOD = biochemical oxygen demand, SRP = soluble reactive phosphorus.

ALAGADOS SAMPLING

Correlations for data simultaneously sampled were only observed between: total cyanobacterial abundance and depth, pH and DO; Aphanizomenonabundance and temperature, pH, total phosphorus and COD;Cylindrospermopsis abundance and depth, pH and DO;Dolichospermum abundance and temperature (Table III).

In turn, stronger and numerous coefficients were observed when using time lag correlations. Correlations were found between: total cyanobacterial abundance and depth, temperature, pH, nitrate, and dissolved oxygen;Aphanizomenon abundance and depth, temperature, pH, conductivity, and nitrate; Cylindrospermopsis abundance and depth, temperature, pH, total phosphorus, nitrate, and COD; andDolichospermum and depth, temperature, pH, and total phosphorus (Table III).

DISCUSSION

IAP SAMPLING

The role of phosphorus in cyanobacterial blooms is particularly clear. Phosphorus concentrations have being considered one of the main drivers of cyanobacterial growth (Moisander et al. 2003MOISANDER PH, STEPPE TF, HALL NS, KUPARINEN J AND PAERL HW. 2003. Variability in nitrogen and phosphorus limitation for Baltic Sea phytoplankton during nitrogen-fixing cyanobacterial blooms. Mar Ecol-Prog Ser 262: 81-95., Tanaka et al. 2004TANAKA T, RASSOULZADEGAN F AND THINGSTAD TF. 2004. Orthophosphate uptake by heterotrophic bacteria, cyanobacteria, and autotrophic nanoflagellates in Villefranche Bay, northwestern Mediterranean: Vertical, seasonal, and short-term variations of the competitive relationship for phosphorus. Limnol Oceanogr 49: 1063-1072., Nausch et al. 2008NAUSCH M, NAUSCH G, WASMUND N AND NAGEL K. 2008. Phosphorus pool variations and their relation to cyanobacteria development in the Baltic Sea: A three-year study. J Marine Syst 71: 99-111., Posselt et al. 2009POSSELT AJ, BURFORD MA AND SHAW G. 2009. Pulses of phosphate promote dominance of the toxic cyanophyte Cylindrospermopsis raciborskii in a subtropical water reservoir. J Phycol 45: 540-546., Rangel et al. 2012RANGEL LM, SILVA LHS, ROSA P, ROLAND F AND HUSZAR VLM. 2012. Phytoplankton biomass is mainly controlled by hydrology and phosphorus concentrations in tropical hydroelectric reservoirs. Hydrobiologia 693: 13-28.). We have reached this 'highly recurrent' conclusion mainly by analyzing the variability of sampling points (i.e. MA and VA). Our results suggest that the influence of abiotic variables on cyanobacteria, which is recurrently demonstrated in experimental studies, is evidenced only if the variation of the environment is represented in monitoring studies. Correlating snapshot samples in LTPS are clearly insufficient to represent how a certain location is compared to others, resulting in weak cyanobacteria-environment relationships.

The negative relationships between cyanobacteria and nitrogen (mainly observed in MA and SA) can be explained by the fact that some cyanobacteria (Nostocales order) are able to exploit molecular nitrogen when sources of inorganic-N are depleted (Herrero et al. 2001HERRERO A, MURO-PASTOR AM AND FLORES E. 2001. Nitrogen control in Cyanobacteria. J Bacteriol 183: 411-425.). In high nitrogen concentrations, nostocales cyanobacteria may be outcompeted by other algae that use this limited nutrient. Among the available nitrogen ions, ammonium is an important determinant of cyanobacterial growth (Domingues et al. 2011DOMINGUES RB, BARBOSA AB, SOMMER U AND GALVÃO HM. 2011. Ammonium, nitrate and phytoplankton interactions in a freshwater tidal estuarine zone: potential effects of cultural eutrophication. Aquat Sci 73: 331-343.). However, ammonium showed the weakest correlation with cyanobacterial density in the studied reservoirs. The low stability of this ion may explain the lack of strong relationships. This can also explain why the relationship between ammonium and cyanobacterial abundance is sometimes positive and sometimes negative considering SA.

The negative correlation between cyanobacterial density and transparency is probably a consequence, not a cause, of cyanobacterial growth. The reservoirs located in Paraná are typically dominated by Microcystis aeruginosa (Kützing) Kützing (IAP 2009IAP - INSTITUTO AMBIENTAL DO PARANÁ. 2009. Monitoramento da qualidade das águas dos reservatórios do Estado do Paraná, no período de 2005 a 2008. Curitiba. Available in: www.iap.pr.gov.br, 120 p.
www.iap.pr.gov.br...
), which form a scum on the water surface, decreasing underwater radiation available to other phytoplankton groups. This is an important competitive advantage of cyanobacteria. Considering the correlations with other environmental variables, COD represents the amount of dissolved organic matter in the water and has been used as an indicator of environmental pollution in many countries (Zhao et al. 2004ZHAO H, JIANG D, ZHANG S, CATTERALL K AND JOHN R. 2004. Development of a direct photoelectrochemical method for determination of chemical oxygen demand. Anal Chem 76: 155-160.). Relatedly, previous studies have also shown that cyanobacteria may vary depending on COD concentrations (Wei et al. 2001WEI B, SUGIURA N AND MAEKAWA T. 2001. Use of artificial neural network in the prediction of algal blooms. Water Res 35: 2022-2028., Amé et al. 2003AMÉ MV, DÍAZ M AND WUNDERLIN DA. 2003. Occurrence of toxic cyanobacterial blooms in San Roque Reservoir (Córdoba, Argentina): A field and chemometric study. Environ Toxicol 18: 192-201.). This was the only variable for which SA showed the strongest Pearson correlation among the three approaches. Dissolved oxygen can vary according to temperature, presence of pollutants, photosynthetic activity, and other reasons (Smith and Bella 1973SMITH SA AND BELLA DA. 1973. Dissolved Oxygen and Temperature in a Stratified Lake. J Water Pollut Control Fed 45: 119-133.). This can also explain the variation in correlation values when using SA. Oxygen concentrations were positively or negatively correlated with cyanobacterial abundance, depending on the sampling period.

It is interesting to note that Pearson and Spearman correlations may generate different outputs. Our focus was not to choose the best analysis for each variable considering its distribution. Both Pearson and Spearman correlations are commonly used to evaluate relationships between cyanobacteria and environmental variables (e.g. Bouvy et al. 1999BOUVY M, MOLICA R, OLIVEIRA S, MARINHO M AND BEKER B. 1999. Dynamics of a toxic cyanobacterial Bloom (Cylindrospermopsis raciborskii) in a shallow reservoir in the semi-arid region of northeast Brazil. Aquat Microb Ecol 20: 285-297., Rinta-Kanto et al. 2009RINTA-KANTO JM, KONOPKO EA, DEBRUYN JM, BOURBONNIERE RA, BOYER GL AND WILHELM SW. 2009. Lake Erie Microcystis: Relationship between microcystin production, dynamics of genotypes and environmental parameters in a large Lake. Harmful Algae 8: 665-673., Leigh et al. 2010LEIGH C, BURFORD MA, ROBERTS DT AND UDY JW. 2010. Predicting the vulnerability of reservoirs to poor water quality and cyanobacterial blooms. Water Res 44: 4487-4496., Conradie and Barnard 2012CONRADIE KR AND BARNARD S. 2012. The dynamics of toxic Microcystis strains and microcystin production in two hypertrofic South African reservoirs. Harmful Algae 20: 1-10., Xie et al. 2012XIE L, REDISKE RR, HONG Y, O'KEEFE J, GILLETT ND, DYBLE J AND STEINMAN AD. 2012. The role of environmental parameters in the structure of phytoplankton assemblages and cyanobacteria toxins in two hypereutrophic lakes. Hydrobiologia 691: 255-268.). Spearman coefficients were stronger for some variables and weaker for others, indicating that some variables have linear relationships with cyanobacteria while others are non-linear.

Furthermore, correlating environmental variables and cyanobacteria abundance give us only a clue for cyanobacteria determinants, considering that a correlation does not mean causality. However, numerous experimental studies have demonstrated the importance of most environmental variables used here as predictors of cyanobacterial growth (Paerl et al. 1985PAERL HW, BLAND PT, BOWLES ND AND HAIBACH ME. 1985. Adaptation to High-Intensity, Low-Wavelength Light among Surface Blooms of the Cyanobacterium Microcystis aeruginosa. Appl Environ Microb 49: 1046-1052., Coles and Jones 2000COLES JF AND JONES RC. 2000. Effect of temperature on photosynthesis-light response and growth of four phytoplankton species isolated from a tidal freshwater river. J Phycol 36: 7-16., Aubriot and Bonilla 2012AUBRIOT L AND BONILLA S. 2012. Rapid regulation of phosphate uptake in freshwater cyanobacterial blooms. Aquat Microb Ecol 67: 251-263.), and conclusions based on correlations are commonly used for cyanobacteria management (Izydorczyk et al. 2008IZYDORCZYK K, JURCZAK T, WOJTAL-FRANKIEWICZ A, SKOWRON A, MANKIEWICZ-BOCZEK J AND TARCZYNSKA M. 2008. Influence of abiotic and biotic factors on microcystin content in Microcystis aeruginosa cells in a eutrophic temperate reservoir. J Plankton Res 30: 393-400., Rinta-Kanto et al. 2009RINTA-KANTO JM, KONOPKO EA, DEBRUYN JM, BOURBONNIERE RA, BOYER GL AND WILHELM SW. 2009. Lake Erie Microcystis: Relationship between microcystin production, dynamics of genotypes and environmental parameters in a large Lake. Harmful Algae 8: 665-673., Lu et al. 2013LU X, TIAN C, PEI H, HU W AND XIE J. 2013. Environmental factors influencing cyanobacteria community structure in Dongping Lake, China. J Environ Sci 25: 2196-2206.). Therefore, identifying the best approach to reveal determinants of cyanobacteria in LTPS monitoring are crucial.

ALAGADOS SAMPLING

Compared to correlations using data simultaneously sampled, stronger and numerous correlations were found in time lag correlations between cyanobacterial abundance and environmental variables. Particularly, the recurrent relationships between cyanobacteria and temperature (Paerl and Huisman 2008PAERL HW AND HUISMAN J. 2008. Blooms Like It Hot. Science 320: 57-58.) and phosphorus (Isvánovics et al. 2000ISVÁNOVICS V, SHAFIK HM, PRÉSING M AND JUHOS S. 2000. Growth and phosphate uptake kinetics of the cyanobacterium, Cylindrospermopsis raciborskii (Cyanophyceae) in throughflow cultures. Freshwater Biol 43: 257-275., Aubriot et al. 2010AUBRIOT L, WAGNER F AND FALKNER G. 2010. The phosphate uptake behaviour of phytoplankton communities in eutrophic lakes reflects alterations in the phosphate supply. Eur J Phycol 35: 255-262.) were better revealed through the time lag correlations. Therefore, our results highlight that cyanobacteria abundance is promoted by environmental conditions that acted previously at the site, and reliable relationships (indicating possible causal relations) are revealed only using time lag samplings.

Surely, the best strategy to identify cyanobacterial determinants in observational studies is to continuously sample environmental variables and cyanobacterial abundance to track for short-term responses. However, such strategy is rarely applied by monitoring agencies, probably due to logistic and financial restrictions. Instead, monitoring studies usually carry out periodic snapshot samplings to reveal cyanobacteria determinants. We argue that, at least for traditional cyanobacteria monitoring, the SA for correlating environmental variables and cyanobacterial abundance may not reveal the determinants of cyanobacteria in reservoirs.

We advocate that investigating the temporal variation of sampling points (i.e. MA and VA) may provide conclusive answers regarding the relationship between environmental variables and cyanobacteria. We conclude that environmental variables sampled synchronously with cyanobacteria are worse predictors of the population density of these microalgae than metrics that represent temporal variation in habitats in long-term monitoring studies.

ACKNOWLEDGMENTS

We are grateful to "Instituto Ambiental do Paraná" and "SANEPAR" for providing the data necessary for this study. We acknowledge Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for continuous financial support. J.W. received a PhD scholarship from CAPES. We would also like to thank Dr. Jeffrey D. Muehlbauer for an extensive English review. We also acknowledge Dr. João C. Nabout and other anonymous reviewers for their valuable suggestions to a previous draft of the manuscript.

REFERENCES

  • AMÉ MV, DÍAZ M AND WUNDERLIN DA. 2003. Occurrence of toxic cyanobacterial blooms in San Roque Reservoir (Córdoba, Argentina): A field and chemometric study. Environ Toxicol 18: 192-201.
  • ANDREOLI CV AND CARNEIRO C. 2005. Gestão integrada de mananciais de abastecimento eutrofizados. Curitiba: Companhia de Saneamento do Paraná SANEPAR/FINEP, 494 p.
  • AUBRIOT L AND BONILLA S. 2012. Rapid regulation of phosphate uptake in freshwater cyanobacterial blooms. Aquat Microb Ecol 67: 251-263.
  • AUBRIOT L, WAGNER F AND FALKNER G. 2010. The phosphate uptake behaviour of phytoplankton communities in eutrophic lakes reflects alterations in the phosphate supply. Eur J Phycol 35: 255-262.
  • AZEVEDO SMFO, CARMICHAEL WW, JOCHIMSEN EM, RINEHART KL, LAU S, SAHW GR AND EAGLESHAM GK. 2002. Human intoxications by microcytins during renal dialysis treatment in Caruaru - Brazil. Toxicology 181: 441-446.
  • BITTENCOURT-OLIVEIRA MC, DIAS SN, MOURA NA, CORDEIRO-ARAÚJO MK AND DANTAS EW. 2012. Seasonal dynamics of cyanobacteria in a eutrophic reservoir (Arcoverde) in a semi-arid region of Brazil. Braz J Biol 72: 533-544.
  • BOUVY M, FALCÃO D, MARINHO M, PAGANO M AND MOURA A. 2000. Occurrence of Cylindrospermopsis (cyanobacteria) in 39 Brazilian tropical reservoirs during the 1998 drought. Aquat Microb Ecol 23: 13-27.
  • BOUVY M, MOLICA R, OLIVEIRA S, MARINHO M AND BEKER B. 1999. Dynamics of a toxic cyanobacterial Bloom (Cylindrospermopsis raciborskii) in a shallow reservoir in the semi-arid region of northeast Brazil. Aquat Microb Ecol 20: 285-297.
  • BURFORD MA, JOHNSON SA, COOK AJ, PACKER TV, TAYLOR BM AND TOWNSLEY ER. 2007. Correlations between watershed and reservoir characteristics, and algal blooms in subtropical reservoirs. Water Res 41: 4105-4114.
  • CHORUS I AND BARTRAM J. 1999. Toxic cyanobacteria in water: a guide to their public health consequences, monitoring and management. London: E & FN SPON, 416 p.
  • COLES JF AND JONES RC. 2000. Effect of temperature on photosynthesis-light response and growth of four phytoplankton species isolated from a tidal freshwater river. J Phycol 36: 7-16.
  • CONRADIE KR AND BARNARD S. 2012. The dynamics of toxic Microcystis strains and microcystin production in two hypertrofic South African reservoirs. Harmful Algae 20: 1-10.
  • DANTAS EW, MOURA AN AND BITTENCOURT-OLIVEIRA MC. 2011. Cyanobacterial blooms in stratified and destratified eutrophic reservoirs in semi-arid region of Brazil. An Acad Bras Cienc 83: 1327-1338.
  • DOMINGUES RB, BARBOSA AB, SOMMER U AND GALVÃO HM. 2011. Ammonium, nitrate and phytoplankton interactions in a freshwater tidal estuarine zone: potential effects of cultural eutrophication. Aquat Sci 73: 331-343.
  • FERNANDES LF, GUTSEIT KC, WOJCIECHOWSKI J, LAGOS PED, XAVIER CF AND WOSIACK AC. 2011. Ecologia do fitoplâncton no reservatório Rio Verde. In: Cunha CLN, Carneiro C, Gobbi EF and Andreoli CV (Eds), Eutrofização em Reservatórios: Gestão Preventiva. Estudo Interdisciplinar na Bacia do Rio Verde, PR. Curitiba: Ed. da UFPR, p. 286-306.
  • FERNANDES LF, WOSIACK AC, DOMINGUES L, PACHECO CV AND LAGOS PED. 2005. Cianobactérias e cianotoxinas. In: Andreoli CV and Carneiro C (Eds), Gestão integrada de mananciais de abastecimento eutrofizados. Curitiba: Sanepar-Finep, p. 369-388.
  • GOMES LNL, OLIVEIRA SMAC, GIANI A AND VON SPERLING E. 2012. Association between biotic and abiotic parameters and the occurrence of cyanobacteria in a Brazilian reservoir. Environ Monit Assess 184: 4635-4645.
  • HAVENS KE, PHLIPS EJ, CICHRA MF AND LI B. 1998. Light availability as a possible regulator of cyanobacteria species in a shallow subtropical lake. Freshwater Biol 39: 547-556.
  • HERRERO A, MURO-PASTOR AM AND FLORES E. 2001. Nitrogen control in Cyanobacteria. J Bacteriol 183: 411-425.
  • HUISMAN J, MATTHIJS HCP AND VISSER PM. 2005. Harmful Cyanobacteria. Netherlands: Springer, 241 p.
  • HUSZAR VLM, SILVA LHS, MARINHO M, DOMINGOS P AND SANT'ANNA CL. 2000. Cyanoprokaryote assemblages in eight productive tropical Brazilian waters. Hydrobiologia 424: 67-77.
  • IAP - INSTITUTO AMBIENTAL DO PARANÁ. 2009. Monitoramento da qualidade das águas dos reservatórios do Estado do Paraná, no período de 2005 a 2008. Curitiba. Available in: www.iap.pr.gov.br, 120 p.
    » www.iap.pr.gov.br
  • ISVÁNOVICS V, SHAFIK HM, PRÉSING M AND JUHOS S. 2000. Growth and phosphate uptake kinetics of the cyanobacterium, Cylindrospermopsis raciborskii (Cyanophyceae) in throughflow cultures. Freshwater Biol 43: 257-275.
  • IZYDORCZYK K, JURCZAK T, WOJTAL-FRANKIEWICZ A, SKOWRON A, MANKIEWICZ-BOCZEK J AND TARCZYNSKA M. 2008. Influence of abiotic and biotic factors on microcystin content in Microcystis aeruginosa cells in a eutrophic temperate reservoir. J Plankton Res 30: 393-400.
  • LEIGH C, BURFORD MA, ROBERTS DT AND UDY JW. 2010. Predicting the vulnerability of reservoirs to poor water quality and cyanobacterial blooms. Water Res 44: 4487-4496.
  • LU X, TIAN C, PEI H, HU W AND XIE J. 2013. Environmental factors influencing cyanobacteria community structure in Dongping Lake, China. J Environ Sci 25: 2196-2206.
  • MCBRIDE GB. 2005. Using Statistical Methods for Water Quality Management: Issues, Problems and Solutions. New Jersey: J Wiley & Sons, 344 p.
  • MOISANDER PH, STEPPE TF, HALL NS, KUPARINEN J AND PAERL HW. 2003. Variability in nitrogen and phosphorus limitation for Baltic Sea phytoplankton during nitrogen-fixing cyanobacterial blooms. Mar Ecol-Prog Ser 262: 81-95.
  • MOURA AN, DANTAS EW, OLIVEIRA HSB AND BITTENCOURT-OLIVEIRA MC. 2011. Vertical and temporal dynamics of cyanobacteria in the Carpina potable water reservoir in northeastern Brazil. Braz J Biol 71: 451-459.
  • NAUSCH M, NAUSCH G, WASMUND N AND NAGEL K. 2008. Phosphorus pool variations and their relation to cyanobacteria development in the Baltic Sea: A three-year study. J Marine Syst 71: 99-111.
  • OLIVEIRA FHPC, CAPELA E, ARA ALS, MOREIRA CHP, LIRA OO, PADILHA MRF AND SHINOHARA NKS. 2014. Seasonal changes of water quality in a tropical shallow and eutrophic reservoir in the metropolitan region of Recife (Pernambuco-Brazil). An Acad Bras Cienc 86: 1863-1872.
  • PAERL HW, BLAND PT, BOWLES ND AND HAIBACH ME. 1985. Adaptation to High-Intensity, Low-Wavelength Light among Surface Blooms of the Cyanobacterium Microcystis aeruginosa. Appl Environ Microb 49: 1046-1052.
  • PAERL HW AND HUISMAN J. 2008. Blooms Like It Hot. Science 320: 57-58.
  • PAERL HW AND HUISMAN J. 2009. Climate change: a catalyst for global expansion of harmful cyanobacterial blooms. Environ Microbiol Rep 1: 27-37.
  • POSSELT AJ, BURFORD MA AND SHAW G. 2009. Pulses of phosphate promote dominance of the toxic cyanophyte Cylindrospermopsis raciborskii in a subtropical water reservoir. J Phycol 45: 540-546.
  • QUESADA A AND FERNÁNDEZ-VALIENTE E. 1996. Relationship between abundance of N2-fixing cyanobacteria and environmental features of Spanish rice fields. Microbial Ecol 32: 59-71.
  • RANGEL LM, SILVA LHS, ROSA P, ROLAND F AND HUSZAR VLM. 2012. Phytoplankton biomass is mainly controlled by hydrology and phosphorus concentrations in tropical hydroelectric reservoirs. Hydrobiologia 693: 13-28.
  • RINTA-KANTO JM, KONOPKO EA, DEBRUYN JM, BOURBONNIERE RA, BOYER GL AND WILHELM SW. 2009. Lake Erie Microcystis: Relationship between microcystin production, dynamics of genotypes and environmental parameters in a large Lake. Harmful Algae 8: 665-673.
  • SMITH SA AND BELLA DA. 1973. Dissolved Oxygen and Temperature in a Stratified Lake. J Water Pollut Control Fed 45: 119-133.
  • STATSOFT, INC. 2005. STATISTICA (data analysis software system), version 7.1. www.statsoft.com.
    » www.statsoft.com
  • TANAKA T, RASSOULZADEGAN F AND THINGSTAD TF. 2004. Orthophosphate uptake by heterotrophic bacteria, cyanobacteria, and autotrophic nanoflagellates in Villefranche Bay, northwestern Mediterranean: Vertical, seasonal, and short-term variations of the competitive relationship for phosphorus. Limnol Oceanogr 49: 1063-1072.
  • VEZIE C, BRIENT L, SIVONEN K, BERTRU G, LEFEUVRE JC AND SALKINOJA-SALONEN M. 1998. Variation of microcystin content of cyanobacterial blooms and isolated strains in Lake Grand-Lieu (France). Microbial Ecol 35: 126-135.
  • WEI B, SUGIURA N AND MAEKAWA T. 2001. Use of artificial neural network in the prediction of algal blooms. Water Res 35: 2022-2028.
  • WICKS RJ AND THIEL PG. 1990. Environmental factors affecting the production of peptide toxins in floating scums of the cyanobacterium Microcystis aeruginosa in a hypertrophic african reservoir. Environ Sci Technol 24: 1413-1418.
  • XIE L, REDISKE RR, HONG Y, O'KEEFE J, GILLETT ND, DYBLE J AND STEINMAN AD. 2012. The role of environmental parameters in the structure of phytoplankton assemblages and cyanobacteria toxins in two hypereutrophic lakes. Hydrobiologia 691: 255-268.
  • YOO RS, CARMICHEL WW, HOEHN RC AND HRUDEY SE. 1995. Cyanobacterial (Blue-Green Algal) Toxins: A Resource Guide. AWWA Research Foundation and American Water Works Association, 229 p.
  • ZHAO H, JIANG D, ZHANG S, CATTERALL K AND JOHN R. 2004. Development of a direct photoelectrochemical method for determination of chemical oxygen demand. Anal Chem 76: 155-160.

Publication Dates

  • Publication in this collection
    28 July 2015
  • Date of issue
    Sept 2015

History

  • Received
    13 Aug 2014
  • Accepted
    18 Dec 2014
Academia Brasileira de Ciências Rua Anfilófio de Carvalho, 29, 3º andar, 20030-060 Rio de Janeiro RJ Brasil, Tel: +55 21 3907-8100 - Rio de Janeiro - RJ - Brazil
E-mail: aabc@abc.org.br