Is thermal stability a factor that influences environmental heterogeneity and phytoplankton distribution in tropical lakes?

Acta Limnologica Brasiliensia, May 2019

Luciana Gomes Barbosa, Francisco Antônio Rodrigues Barbosa, Carlos Eduardo de Mattos Bicudo

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

Is thermal stability a factor that influences environmental heterogeneity and phytoplankton distribution in tropical lakes?

Thematic Section: Homage to Carlos Eduardo de Mattos Bicudo Is thermal stability a factor that influences environmental heterogeneity and phytoplankton distribution in tropical lakes? A estabilidade térmica é um fator que influencia a heterogeneidade ambiental e a distribuição do fitoplâncton em lagos tropicais? Luciana Gomes Barbosa1  *  Francisco Antônio Rodrigues Barbosa2  Carlos Eduardo de Mattos Bicudo3  1 Laboratório de Limnologia, Departamento de Fitotecnia e Ciências Ambientais, Universidade Federal da Paraíba – UFPB, Campus II, Rodovia BR 079, Km 12, CEP 58397-000, Areia, PB, Brasil 2 Laboratório de Limnologia, Ecotoxicologia e Ecologia Aquática – LIMNEA, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais – UFMG, Av. Presidente Antônio Carlos, 6627, Pampulha, CEP 31270-901, Belo Horizonte, MG, Brasil 3 Núcleo de Pesquisas em Ecologia, Instituto de Botânica, Av. Miguel Stéfano, 3687, Água Funda, CP 68041, CEP 04045-972, São Paulo, SP, Brasil Abstract Aim To evaluate the effects of environmental heterogeneity promoted by thermal stratification on the distribution of limnologic variables and phytoplankton functional groups (FGs) in two natural lakes. Methods Monthly measurements were performed over a five-year period in the vertical profile of a warm monomitic shallow lake (Lago Carioca) and in a deep and meromictic (Lake Dom Helvécio). Results The vertical zonation generated by the high thermal stability during the stratification period promoted an increase in the spatial heterogeneity and, consequently, in the richness of functional groups of the two lakes. In the epilimnion, the dominance of small chlorophytes and desmids (NA, A and X1, ≤20μm) and larger dinophyceans (Lo, ≥20μm) were associated with high turbulence and light availability and soluble phosphorus limitation. In the metalimnion, the presence of filamentous cyanobacteria (R) and colonial chlorophyceans (F) of larger size (≥20µm) were associated with stable habitats with high concentrations of N-NH4 and P-PO4-3 and light limiting conditions. Comparatively, Lake Dom Helvécio presented a higher richness of FGs in the meta-hipolimnetic layers (SN, P) as well as a higher number of species per functional group. Seasonal changes in the climatic conditions (e.g. the decrease in air temperature with the consequent heat loss) caused the break of the water column stability, which promoted the redistribution of the dissolved nutrient forms and the increase of light limitation in the two lakes during the mixing period. Therefore, there was a drastic reduction in the richness and population biomass of FGs (≤80%). Conclusions Thermal stability and atelomixis were the main driving forces of vertical heterogeneity during the stratification, favoring the coexistence of FGs and, consequently, their increase in richness and biomass. Keywords:  atelomixis; environmental changes; functional groups; thermal stratification Resumo Objetivo Avaliar os efeitos da heterogeneidade ambiental promovida pela estratificação térmica sobre a distribuição de variáveis limnológicas e de grupos funcionais fitoplanctônicos (GF´s) em dois lagos naturais. Métodos Foram realizadas medidas mensais no perfil vertical de um lago raso monomítico quente (Lago Carioca) e em um profundo e meromítico (Lago Dom Helvécio), ao longo de cinco anos. Resultados A zonação vertical evidente gerada pela elevada estabilidade térmica durante o período de estratificação favoreceram o aumento da heterogeneidade espacial e, consequentemente, da riqueza de grupos funcionais nos dois lagos. No epilímnio, ocorreu dominância de clorofíceas e desmídias de pequeno tamanho (NA, A e X1, ≤20µm) e dinofíceas de maior porte (Lo, ≥20µm), associadas a alta turbulência, disponibilidade de luz e limitação por fósforo solúvel. No metalímnio, foram observadas cianobactérias filamentosas (R) e clorofíceas coloniais (F) de maior tamanho (≥20µm), associadas a hábitats estáveis com elevadas concentrações de N-NH4 e P-PO4-3 e condições limitantes de luz. Comparativamente, o lago Dom Helvécio registrou maior riqueza de GF´s nas camadas meta-hipolimnéticas (SN, P) bem como maior de número de espécies acomodadas por grupo. Sazonalmente, mudanças nas condições climáticas (ex. redução da temperatura do ar provocando perda de calor) causaram a quebra da estabilidade da coluna d’água, o que promoveu a redistribuição das formas dissolvidas de nutrientes e o aumento da limitação por luz nos dois lagos durante o período de mistura. Como consequência, foi identificada uma drástica redução de riqueza de GF´s e da biomassa das populações (≤80%). Conclusões A estabilidade térmica e atelomixia foram funções de força fundamentais na geração de heterogeneidade vertical durante a estratificação, favorecendo a coexistência e, consequentemente, o aumento da riqueza e biomassa de GF´s. Palavras-chave:  atelomixia; mudanças ambientais; grupos funcionais; estratificação térmica 1. Introduction Phytoplankton periodicity has been associated with seasonal stratification and turbulence cycles for a long time ( Margalef, 1978 ). In lakes, water physical factors control the mechanisms of adjustment and the strategies of phytoplankton for adaptation to environmental changes ( Margalef, 1978 ; Estrada & Berdalet, 1997; Moreno-Ostos et al., 2006 ). Aquatic ecosystems have a complex spatial structure, being the phytoplankton patchiness physically controlled by light conditions and turbulence mixing ( Moreno-Ostos et al., 2006 ). Studies based on growth simulations indicated that phytoplankton responses are depending on the position of the layer hosting algal cells, the depth of the water column, the light availability and the shape and size of the water body ( Rigosi & Rueda, 2012 ). Changes in the phytoplankton dynamics in stratified lakes can be explained by the synergy between turbulence and the light conditions underwater ( Moreno-Ostos et al., 2006 ), since light limitation strongly affects the succession and regulation of the plankton ( Zohary et al., 2010 ). Lewis Junior (1973) described the complete atelomixis as the mixing of the water column where a permanent thermocline can be established during summer in tropical monomitic lakes. Atelomixis has been frequently reported affecting phytoplankton in tropical lakes, however, since the thermocline is developed during the summer and contributes to epilimnetic mixing, it has been considered as a partial atelomixis (Barbosa & Pádisak, 2002). Partial atelomixis possibly drives the selection of algal groups, favoring the maintenance of desmids ( Barbosa & Padisák, 2002 ; Padisák et al., 2009 ; Tavera & Martinez-Almeida, 2005) and the suspension of diatoms throughout the layer ( Lopes et al., 2005 ). The characteristics of the phytoplankton in lakes are selected according to the spatial complexity and temporal changes in environmental factors. MacArthur & MacArthur (1961) considered that environments that are more heterogeneous turn resources more available and increase the number of niches, and consequently, support higher species diversity than less heterogeneous environments. Hutchinson (1961) stated that a large number of species might be able to coexist in a relatively isotropic environment. Therefore, environmental heterogeneity enhanced by vertical gradients of temperature can have great influence in the vertical distribution of phytoplankton ( Mellard et al., 2011 ). Currently, 40 polyphyletic functional groups of phytoplankton are recognized ( Padisák et al., 2009 ) according to the classical and widely used system of phytoplankton classification, which is highly predictive for tropical and subtropical water bodies ( Salmaso et al., 2015 ). In this paper, we characterized the phytoplankton functional groups and analyzed their responses to seasonal variation in the thermal structure of two lakes in the state of Minas Gerais, Brazil. In addition, we sought to answer the following question: Does the apparent increase in heterogeneity during stratification increase the number of functional groups? 2. Material and Methods 2.1. Study area The typology of the lake system of the Rio Doce middle section (Minas Gerais state, Brazil) includes rare and deep dendritic lakes with slow circulation (e.g. Lake Dom Helvécio), and shallow circular lakes with more effective circulation, which are most common in this landscape (e.g. Lake Carioca) ( De Meis & Tundisi, 1997 ). Two lakes of the Rio Doce middle section were sampled ( Figure 1 ): lake Carioca, which is small, mesotrophic and shallow (maximum depth of 11.8 m) and lake Dom Helvécio, which is oligotrophic, dendritic and deep (maximum depth of 39.2 m) ( Bezerra-Neto & Pinto-Coelho, 2008 ). Figure 1 Carioca and Dom Helvécio lakes, southeast Brazil.  2.2. Sampling and analyses Sampling was carried out monthly from January 2002 to December 2006. Water samples at different layers (subsurface, 10% incident light, 1% incident light and the aphotic zone) of the limnetic zone of both lakes were collected using van Dorn bottles. Thermocline depth and thermal stability were used for characterizing the thermal structure. Water temperature, dissolved oxygen, electric conductivity, and pH were monthly measured ‘in situ’, using a multiparameter probe (Horiba sensor model U 22), at 0.50 m intervals from the surface to the bottom of the lakes. Water transparency was estimated using a Secchi disk ( Cole, 1994 ). The euphotic zone (Zeu) was calculated as 2.7 times the Secchi disk depth ( Cole, 1994 ). Vertical strata (epilimnion, metalimnion and hypolimnion) were defined according to the temperature profile ( Reynolds, 1984 ), and considering density gradients greater than 0.02 kg m-3 m-1 , from water density tables ( Hutchinson, 1957 ). The euphotic zone (Zeu): mixing zone (Zmix) ratio was used as an index for light availability in the mixing zone ( Jensen et al., 1994 ). Thermal stability (g-cm cm-2) was calculated according to Idso (1973) . Total phosphorus (TP), soluble reactive phosphorus (SRP), total nitrogen (TN = particulate N and total dissolved N) and dissolved inorganic nitrogen (NO3- , NO2- and NH4+ nitrogen) were measured according to Golterman et al. (1978) . The Carlson´s index adapted by Toledo Junior et al. (1983) for tropical regions was used for characterization of the trophic state. Phytoplankton quantification was performed according to Utermöhl (1958) and sedimentation time was calculated according to Lund et al. (1958) . The sedimentation time was 4 h for each centimeter of the chamber height. Counting was carried out under a Zeiss inverted microscope at 400 × magnification, considering a minimum of 400 individuals of the most frequent species ( Lund et al., 1958 ). The biovolume (μm3 ml-1 → mm3 L -1) of each species was estimated from the geometric solid shape that best represented the cell, taken either isolated or combined ( Sun & Liu, 2003 ; Hillebrand et al., 1999 ). Biovolume was expressed in dry weight units, where 1 mm3 L-1 = 1 mg L-1 ( Wetzel & Likens, 2000 ). Species were classified into functional groups according to Reynolds et al. (2002) and Padisák et al. (2009) . 2.3. Statistical data Vertical heterogeneity was estimated according to vertical zonation (epilimnion, metalimnion, and hypolimnion). Biotic and abiotic variables were compared with the Kruskal-Wallis test (Statsoft Inc. v. 4.2). Multivariate analyses were performed considering a single monthly value of Zmix , Secchi disk (as an indicator of water transparency) and light attenuation coefficient. Phytoplankton biomass and dissolved and total nutrients of the epilimnion, metalimnion, and hypolimnion were measured and integrated arithmetically. A Principal Component Analysis (PCA) of the environmental variables of the five-year series was performed from the covariance matrices, using the abiotic data transformed by Log (x + 1). A multivariate descriptive analysis of the five year abiotic and biotic variables was performed with canonical correspondence analysis (CCA). The significance of environmental variables (p<0.05) was tested using the Monte Carlo´s test with 999 unrestricted permutations. The Pearson correlation coefficient (r) resulting from the relationship between the ordination values was used, as well as the individual variables for the ordination ( McCune & Mefford, 1997 ). The functional groups data was transformed by ranging. FITOPAC ( Shepherd, 1996 ) was used for data transformation, whereas PC-ORD v. 3.0 for Windows ( McCune & Mefford, 1997 ) was employed for all data analyses. 3. Results 3.1. Thermal stability The stratification period was characterized by high physical stability values (maximum 408 g cm cm-2 and 98 g cm cm-2 at lakes Dom Helvécio and Carioca, respectively) related to a high heterogeneity enhanced by the temperature distribution ( Figure 2 ). The epilimnion layer (Zmix) was about 6.5 m at Lake Carioca and 12 m at Dom Helvécio. At Lake Carioca, maximum metalimnion and hypolimnion depths were 8.5 m and 2-3 m, respectively. The beginning of the metalimnion at Lake Dom Helvécio was between 13 and 16 m, and its maximum extension reached 26 m. The mixing zone (Zmix) progressively increased during the turnover period (May-September), mainly in Lake Carioca (Zmix=Z max). Figure 2 Depth-time diagram of water temperature (ºC) and thermal stability (filled line; g cm cm-2) over 5 years (2002-2006) in lakes Carioca (above) and Dom Helvécio (below), Rio Doce State Park, southeast Brazil (M: mixing period).   Long periods of thermal stratification (an average of 200 days per year) resulted in the establishment of chemical and nutrient stratification, with nutrient accumulation associated with anoxia at the bottom of both lakes. Despite the total phosphorus (TP) low concentrations, both lakes were classified as mesotrophic (Carioca) and oligotrophic (Dom Helvecio) and ( Table 1 and Table 2 ). The total phosphorus (p = 0.000) and electrical conductivity (p = 0.000) were the only factors that differed between the lakes, with higher values in Lake Carioca ( Table 1 ). Table 1 Water physical and chemical variables (averages and standard deviations of the five years of Lake Carioca, Rio Doce State Park, Minas Gerais State, southeast Brazil.   Carioca Stratification Period Mixing Period Transparency (m) 1.6 (±0.4) 1.2 (±0.3) Light attenuation vertical coefficient (k) 1.2 (±0.3) 1.5 (±0.5) Zmix 4 (±1.8) 10 (±1.9) Zeu/ Zmix 1.5 (±0.9) 0.5 (±0.4) Trophic state Mesotrophic Mesotrophic Epiliminion Metalimnion Hypolimnion Sub-surface Bottom pH 6.07 (±0.8) 6.11 (±0.8) 6.7 (±1.1) 6.1 (±1.5) 5.9 (±0.8) Electrical Conductivity (μS cm-1) 29.6 (±3.7) 47.3 (±30.6) 104.8 (±34.8) 32.9 (±2.9) 51.2 (±28) Dissolved Oxygen (mg L-1) 9.3 (2.7) 6.4 (±3) 4.1 (±2.7) 10.5 (±4.2) 5.7 (±3.4) Temperature (º C) 30.3 (±1.7) 26.3 (±1.8) 23.8 (±1.6) 25.2 (±1.9) 22.8 (±1.4) Total Dissolved Solids (TDS) (mg L-1) 20.3 (±1.7) 32.6 (±20.5) 64.1 (±24.9) 20 (±0) 29.6 (±19.5) Total Phosphorus (μg L-1) 20.5 (±12.6) 33.4 (±15.8) 54.2 (±132) 24.7 (±17.9) 36.6 (±32) P-PO4 -3 (μg L-1) 3.2 (±2.1) 5.7 (±4.8) 8 (±8.6) 4.0 (±3.3) 3 (±2.5) Total Nitrogen (μg L-1) 302.5 (±150) 593.6 (±435.3) 1159 (±734.5) 411 (±242) 650.7 (±397) N-NH4 (μg L-1) 15.2 (±10.5) 206.2 (±349) 766.8 (±593) 65.4 (±49.9) 232.3 (±240) N-NO3 (μg L-1) 6.5 (±7) 7.1 (±5) 8.5 (±6) 13.1 (±11.6) 11.7 (±9.1) N-NO2 (μg L-1) 1.6 (±1.1) 3.3 (±4.3) 6.2 (±5.8) 2.2 (±1.6) 2.5 (±1.6) Table 2 Water physical and chemical variables (averages and standard deviations of the five years of Lake Dom Helvécio, Rio Doce State Park, Minas Gerais State, southeast Brazil.   Dom Helvécio Stratification Period Mixing Period Transparency (m) 3.6 (±0.8) 2.5 (±0.5) Light attenuation vertical coefficient (k) 0.58 (±0.5) 0.78 (±1.1) Zmix 12.0 (±0.7) 32.0 (±0.3) Zeu/Zmix 1.5 (±0.5) 0.7 (±1.7) Trophic state Oligotrophic Mesotrophic Epiliminion Metalimnion Hypolimnion Sub-surface Bottom pH 6.2 (±1.3) 6.0 (±1.2) 6.4 (± 1.5) 6.3 (±1.1) 6.3 (±1.15) Electrical Conductivity (μS cm-1) 38.4 (±8.7) 41.6 (± 13.8) 72.5 (± 19.3) 44 (±5.7) 65.9 (±22.5) Dissolved Oxygen (mg L-1) 9.2 (±2.8) 6.4 (±3.2) 4.1 (± 3) 10 (±4.3) 4.2 (±3.3) Temperature (º C) 29.8 (±1.3) 26.5 (±1.8) 23.8 (± 0.4) 25.9 (±1.35) 23.7 (±0.7) Total Dissolved Solids (mg L-1) 23.2 (±6.8) 25.7 (±10) 46.3 (± 12.8) 28.8 (± 5.3) 40.8 (±15.8) Total Phosphorus (μg L-1) 14.6 (±9.3) 27.8 (±32.7) 28.7 (±35.7) 26.5 (±23) 34.2 (34.4) P-PO4 -3 (μg L-1) 3.2 (±2.8) 3.7 (±2.7) 4.9 (±3.3) 3.4 (±3.9) 4.2 (±3.9) Total Nitrogen (μg L-1) 337.8 (±163.8) 448.3 (±248) 807.4 (±456) 447.6 (± 167) 805.7 (333.1) N-NH4 (μg L-1) 47.5 (±63) 102.3 (±153.2) 489.5 (±428.1) 136.1 (± 125.2) 409.1 (±326.8) N-NO3 (μg L-1) 9.4 (±8.2) 7.7 (±5.7) 10.6 (±10) 9.9 (± 7) 7.97 (±5.9) N-NO2 (μg L-1) 1.7 (±1.4) 2.2 (±2.1) 2.6 (±2.8) 1.6 (± 0.9) 3.1 (±2.96) Regarding the vertical profile of light, two distinct phases associated with thermal stability were observed: one of high light availability (Zeu/Zmix≥1) in September-April and another of low availability in May-August (Zeu/Zmix ≤1). Principal Component Analysis showed that the first two axes accounted for 66% of total data variability ( Figure 3 ). Four variables contributed to the Axis 1: mixing zone (ZMIX; r = 0.6); thermal stability (S; r = 0.6); transparency (SD; r = 0.4); and the coefficient of light attenuation (K0; r = –0.4). In the case of Axis 2, thermal stability (S; r = 0.55), P-PO4–3 (r = 0.12), mixing zone (ZMIX; r = –0.8) and the coefficient of light attenuation (K0; r = –0.15) showed the highest contribution. Figure 3 PCA biplots of water physical and chemical variables from lakes Carioca and Dom Helvécio (southeast Brazil), over the study period (2000-2006).   Samples from Dom Helvécio were positively associated with Axis 1 and evidenced differences between thermal stratification and mixing periods, whereas samples from Lake Carioca were associated with the coefficient of light attenuation and P-PO4–3 , and showed no differences between those periods. 3.2. Phytoplankton functional groups Seven functional groups were identified at Lake Carioca and 12 at Lake Dom Helvécio, which corresponded to c. 80% of the total biovolume ( Table 3 ). In Lake Carioca, total biomass ranged from 0.0003 to 47.2 mg L-1, and the maximum values (>10 mg L-1) were found at the metalimnion during the stratification period. In Lake Dom Helvécio, the biomass ranged from 0.1 to 19 mg L-1, and the maximum values were identified between the epi and metalimnion. Table 3 Phytoplankton functional groups for Carioca and Dom Helvécio lakes between 2002-2006.   Period Layer Species Functional Group Carioca Stratification Epilimnion Cosmarium asphaerosporum Nordst. var. strigosum Nordst.; S. tetracerum (Kütz.) Ralfs ex Ralfs var. tetracerum f. tetracerum ; Staurodesmus crassus (West) Flor; Teilingia granulata (Roy & Biss.) Bourr. NA Gymnodinium sp.; Peridinium volzii Lemmermann Lo Chlorella vulgaris Beyerinck X1 Metalimnion Botryococcus braunii Kützing; Coelastrum sphaericum Nägeli F Pseudanabaena galeata Böcher R Hypolimnion Planktolyngbya limnetica (Lemmermann) Komárková-Legnerová & Cronberg; Pseudanabaena galeata Böcher; Pseudanabaenaceae NI R Staurastrum rotula Nordstedt P Mixing Period Lyngbya sp. S1 (Zmix=Zmax) Dom Helvécio Stratification Epilimnion Cosmarium asphaerosporum Nordst. var. strigosum Nordst.; C. contractum O.Kirchner; Staurastrum laeve Ralfs, S. tetracerum (Kütz.) Ralfs ex Ralfs var. tetracerum f. tetracerum; S. taylorii Grönblad; S. smithii; Staurodesmus crassus (West) Flor; Teilingia granulata (Roy & Biss.) Bourr. NA Chlorella vulgaris Beyerinck X1 Closterium aciculare T.West P Planktolyngbya limnetica (Lemmermann) Komárková-Legnerová & Cronberg; Pseudanabaena galeata Böcher R Urosolenia longiseta (O.Zacharias) Edlund & Stoermer A Metalimnion Planktolyngbya limnetica (Lemmermann) Komárková-Legnerová & Cronberg; Pseudanabaena galeata Böcher R Cylindrospermopsis raciborskii (Woloszynska) Seenayya & Subba Raju SN Hypolimnion Staurastrum rotula Nordstedt P Limnothrix sp.; Oscillatoria limosa C.Agardh ex Gomont, Planktolyngbya limnetica (Lemmermann) Komárková-Legnerová & Cronberg; Pseudanabaena galeata Böcher R Mixing Period (Zmix=Zmax) Cryptomonas brasiliensis A. Castro Y Lyngbya sp. S1 Peridinium inconspicuum Lemmermann; P. volzii Lemmermann Lo 3.3. Lake carioca In the epilimnion, group NA accounted for 15-30% of the total biomass. Larger specimens (>20 µm), as well as the most morphologically complex ones (e.g. Staurastrum rotula), were more abundant in the hypolimnion ( Table 3 ). Colonial green algae (group F) and small green algae (group X1 ) exhibited a remarkable contribution between epi and metalimnion. Dinoflagellates (group Lo) also constituted a typical epilimnetic assemblage. The group R had populations with maximum biomass (30.5 mg L-1) at the metalimnion and hypolimnion. During the mixing period, group S1 showed the largest biomass contribution. 3.4. Lake Dom Helvécio In the epilimnion, maximum values of desmid biomass (13.2 mg L-1) were recorded in 2002. In the metalimnion and hypolimnion, only S. rotula contributed substantially to total biomass. In addition, in the epi-metalimnion interface, group A and bloom-forming cyanobacteria (SN) were registered. Other groups, such as group Y represented by Cryptomonas brasiliensis , exhibited the highest biomass in the epilimnion during the stratification period and in the aphotic zone during mixing. The group R was detected during stratification in the epilimnion-metalimnion as well as in the hypolimnion, with a high biomass of Planktolyngbya limnetica , Limnothrix sp., Pseudanabaena galeata and Oscilatoria limosa (maximum of 3 mg L-1). During the mixing period, their contributions were above 50%, with Limnothrix sp. being replaced by Lyngbya sp. and group S1. 3.5. Canonical analyses CCA resulted in eigenvalues of 0.50 e 0.073 for axes 1 and 2, respectively. Pearson correlation for axes 1 (0.88) and 2 (0.67) indicated a relationship between abiotic factors and functional groups at a monthly scale. Monte Carlo´s test showed that the first two canonic axes were statistically significant (p < 0.01; Figure 4 ). Figure 4 CCA ordination diagram of the sampling units in relation to the functional groups (mg L-1) at lakes Dom Helvécio and Carioca and main abiotic variables during the period 2002-2006.   According to the canonic coefficients, the thermal stability (S; r = 0.88), mixing zone (ZMIX; r = 0.84) and coefficient of vertical light attenuation (K0; r = -0.20) were associated with Axis 1. Intra-set correlations confirmed thermal stability (r = 0.70), mixing zone (r = 0.67), and the coefficient of vertical light attenuation (r = -0.58) as determining variables for Axis 1. Axis 2 presented higher canonic coefficients for pH (r = 0.89), electrical conductivity (EC; r = 0.25), coefficient vertical light attenuation (K0; r = 0.21), N-NH4 (r = 0.15), P-PO4 (r = -0.19) and transparency (DS; r = -0.28) than Axis 1. Intra-set correlations indicated pH (r = 0.72), electrical conductivity (r = 0.38), coefficient of vertical light attenuation (r = 0.35), N-NH4 (r = 0.26), P-PO4 (r = -0.14) and transparency (r = -0.35) as the ordination determining factors. Along Axis 1, the coefficient of vertical light attenuation (r = -0.20), and the highest concentrations of nutrients (N-NH4 and P-PO4) were associated with the most representative functional groups of Lake Carioca, especially groups X1 (r = -0.2), F (r = -1.26), S1 (r = -1.27) and R (r = -1.86). Along this axis, functional groups S1 (r = 1.5), SN (r = 1.4), A (r = 1.22), Y (r = 0.65), NA (r = 0.40) and LO (r = 0.18) were associated with the mixing zone (ZMIX; r = 0.67) and the transparency (SD; r = -0.28) in Lake Dom Helvécio. 4. Discussion The vertical distribution of functional groups appears to be controlled by a set of factors, such as thermal stability, nutrient limitations, underwater light conditions and phytoplankton characteristics (buoyance capacity) (e.g. Reynolds, 1984 ). These factors are highly influenced by the compartmentalization resulting from thermal stability during the stratification period. The drastic decrease in the biomass of functional groups during the mixing period in both lakes (≥80%) was an indicative of the seasonal break of thermal stability. Some models of phytoplankton vertical distribution include stratification, spatially varying mixing conditions and limitation of different resources within the mixed layer. Multiple nutrients sources can create multi-modal phytoplankton distributions ( Mellard et al., 2011 ). Our results showed that the vertical distribution of phytoplankton was influenced by light conditions and nutrient availability during stratification. The biomass of cyanobacteria increased at the deepest layers while of desmids increased at the upper layers of lakes (meta e hypolimnion), which was related to a lower concentration of nutrients and a higher light availability in Zmix during stratification, as previously reported for other lakes ( Moreti et al., 2013 ). Mellard et al. (2011) suggested that phytoplankton biomass might sustain source populations in both mixed and deep layers simultaneously, with favorable conditions and resource availability for growth at both depths, which represents a bimodal distribution. Thermal stability and partial atelomixis synergistically supported the presence of groups NA and A, and of unicellular Chlorococcales (group X1) in the epilimnion, while of cenobial coccoid Chlorococcales (group F) in the metalimnion. Previous studies have reported that the increase in desmid biomass in tropical lakes is related to the epilimnetic mixing resulting from diurnal variations in the air temperature (Barbosa & Pádisak 2002; Tavera & Martínez-Almeida, 2005 ; Souza et al., 2008 ). Therefore, the mixing regime (turbulence) and the physical variables (Zeu:Zmix) are the most important factors for the maintenance and growth of those functional groups within the epilimnion, especially of non-motile species with a sinking tendency. The persistence of cyanobacteria may be directly associated with the existence of a thermal stratification in the water column ( Maier et al., 2001 ). The ability of members of this group to take advantage of intensified and prolonged periods of thermal stratification is well-known ( Huisman et al., 2005 ). We found that thermal stability increased the cyanobacteria biomass in the metalimnion and in the hypolimnion of both lakes ( Reynolds et al., 1983 ; Barbosa & Padisák, 2002 ; Padisák et al., 2003), evidencing an inter-annual recurrent pattern since the 1980s ( Reynolds et al., 1983 ). At the thermocline depth, cyanobacteria find the appropriate conditions for growth, such as suitable light levels, high stability, and nutrients from the hypolimnion that frequently entrained into the epilimnion ( Moreno-Ostos et al., 2006 ). However, the biomass of group R was higher in Lake Carioca than Dom Helvécio, probably due to the interaction between the high trophic state, physical factors and low transparency associated with an increase in the concentration of pigmented organic substances ( Bezerra-Neto et al., 2006 ). Singularities related to high water temperature and the relative hydrological isolation found in the two lakes of the middle section of the Rio Doce system indicate the almost absence of physical disturbances ( Reynolds, 1997 ). Some conditions, such as the strength of gradients (e.g. resulting from the relative density of the different layers) and the depth of light penetration, may represent effective mechanisms of niche segregation ( Naselli-Flores et al., 2003 ). In the studied lakes, such small-scale disturbances occurred due to breakage of superficial secondary thermoclines by allochthonous factors of disturbance (e.g. wind and rain), resulting in partial atelomixis ( Barbosa & Padisák, 2002 ). The seasonal disruption of stability may contribute to reducing the number of habitats and consequently the availability of niches that are suitable for colonization. The decrease in habitat heterogeneity by disruption of thermal stability also modifies the resource availability (light and nutrients), intensifying some processes such as competitive exclusion. Seasonal climatic changes induced shifts in the thermal structure of lakes, which caused the loss of the stability in the water column. Consequently, changes in the phytoplankton composition of the lakes promoted the collapse of populations during the stratified period, and cyanobacteria (group S1) dominated during the mixing. Despite adaptations to stratification, cyanobacteria are known to dominate or co-dominate during mixing ( Soares et al., 2013 ), probably due to their ability to tolerate shading and to succeed in turbid and eutrophic lakes, which are typical conditions of lakes during mixing (Padisak et al., 2009). In conclusion, thermal stability was the main factor generating vertical heterogeneity during the stratification, favoring the coexistence of FGs. In addition, the highest biomass of FGs supported by a higher thermal stability, the larger vertical extension of the euphotic zone and of the epi, meta and hypolimnetic layers were found in the deeper and more stable Lake Dom Helvécio. Our results corroborated previous studies in which phytoplankton variations in deep lakes are predictable and less susceptible to stochastic events ( Salmaso, 2003 ). Cite as: Barbosa, L.G., Barbosa, F.A.R. and Bicudo, C.E.M. Is thermal stability a factor that influences environmental heterogeneity and phytoplankton distribution in tropical lakes? Acta Limnologica Brasiliensia, 2018, vol. 30, e207. References BARBOSA, F.A.R. and PADISÁK, J. The forgotten lake stratification pattern: atelomixis and its ecological importance. Verhandlungen der Internationalen Vereinigung fur Theoretische und Angewandte Limnologie, 2002, 28, 1385-1395. [ Links ] BEZERRA-NETO, J.F. and PINTO-COELHO, R.M. Morphometric study of lake Dom Helvécio, Parque Estadual do Rio Doce (PERD), Minas Gerais, Brazil: a reevaluation. Acta Limnologica Brasiliensia, 2008, 20(2), 161-167. [ Links ] BEZERRA-NETO, J.F., BARBOSA, P.M. and BARBOSA, F.A.R. Temporal variation of attenuation of the visible radiation in four lakes of Parque Estadual do Rio Doce (PERD), Minas Gerais, Brazil. Acta Limnologica Brasiliensia, 2006, 18, 39-46. [ Links ] COLE, G.A. Textbook of limnology. 4th ed. Illinois: Waveland Press, Prospect Heights, 1994. [ Links ] DE MEIS, M.R.M. and TUNDISI, J.G. Geomorphological and limnological process as a basis for lake tipology: the middle Rio Doce lake system. In J.G. TUNDISI and Y. SAIJO, eds. Limnological studies in the Rio Doce Valley Lakes. São Carlos: Brazilian Academy of Sciences, 1997, pp. 25-50. [ Links ] ESTRADA, M. and BERDALET, E. Phytoplankton in a turbulent world. Scientia Marina , 1997, 61(1), 125-140. [ Links ] GOLTERMAN, H.L., CLYMO, R.S. and OHMSTAD, M.A.M. Methods for chemical analyses of freshwaters. 2nd ed. Oxford & Edinburg: Blackwell Scientific Publications, 1978. [ Links ] HILLEBRAND, H., DÜRSELEN, C.D., KIRSCHTEL, D., POLLINGHER, D. and ZOHARY, T. Biovolume calculation for pelagic and benthic microalgae. Journal of Phycology , 1999, 35(2), 403-424. [ Links ] HUISMAN, J., MATTHIJS, H.C.P. and VISSER, P.M. Harmful cyanobacteria. Netherlands: Springer, 2005. [ Links ] HUTCHINSON, G. E. Concluding remarks. Cold Spring Harbor Symposia on Quantitative Biology, 1957, 22: 415-442. [ Links ] HUTCHINSON, G.E. The paradox of the plankton. American Naturalist, 1961, 95(882), 137-145. [ Links ] IDSO, S. On the concept of lake stability. Limnology and Oceanography , 1973, 18(4), 681-683. [ Links ] JENSEN, P., JEPPESEN, E., OLRIK, K. and KRISTENSEN, P. Impact of nutrients and physical factors on the shift from Cyanobacterial to Chlorophyte dominance in shallow Danish lakes. Canadian Journal of Fisheries and Aquatic Sciences, 1994, 51(8), 1692-1699. [ Links ] LEWIS JUNIOR, W.M. The thermal regime of Lake Lanao (Phillipines) and its theoretical implications for tropical lakes. Limnology and Oceanography, 1973, 18(2), 200-217. [ Links ] LOPES, M.R.M., BICUDO, C.E.M. and FERRAGUT, C. Spatial and temporal variation of phytoplankton in a shallow tropical oligotrophic reservoir, southeast Brazil. Hydrobiologia , 2005, 542(1), 235-247. [ Links ] LUND, J.W.G., KIPLING, C. and LECREN, E.D. The inverted microscope method of estimating algal numbers and the statistical basis of estimations by counting. Hydrobiologia , 1958, 11(2), 143-170. [ Links ] MACARTHUR, R.H. and MACARTHUR, J.W. On bird species diversity. Ecology , 1961, 42(3), 594-598. [ Links ] MAIER, H.R., BURCH, M.D. and BORMANS, M. Flow management strategies to control blooms of the cyanobacterium Anabaena circinalis, in the river Murray at Morgan, South Australia. Regulated Rivers: Research and Management, 2001, 17(6), 637-650. [ Links ] MARGALEF, R. Life-forms of phytoplankton as survival alternatives in an unstable environment. Oceanologica Acta, 1978, 1(4), 493-509. [ Links ] MCCUNE, B. and MEFFORD, M.J. PC-ORD. Multivariate analysis of ecological data, Version 3.0. Oregon: MjM Software Design, 1997. [ Links ] MELLARD, J.P., YOSHIYAMA, K., LITCHMAN, E. and KLAUSMEIER, C.A. The vertical distribution of phytoplankton in stratified water columns. Journal of Theoretical Biology , 2011, 269(1), 16-30. PMid:20932846. [ Links ] MORENO-OSTOS, E., CRUZ-PIZARRO, L., BASANTA, A., ESCOT, C. and GEORGE, D.G. Algae in the motion: spatial distribution of phytoplankton in thermally stratified reservoirs. Limnetica , 2006, 25(1-2), 205-217. [ Links ] MORETI, L.O.R., MARTOS, L., BOVO-SCOMPARIN, V.M. and RODRIGUES, L.C. Spatial and temporal fluctuation of phytoplankton functional groups in a tropical reservoir. Acta Scientiarum. Biological Sciences, 2013, 35(3), 359-366. [ Links ] NASELLI-FLORES, L., PADISÁK, J., DOKULIL, M.T. and CHORUS, I. Equilibrium/steady-state concept in phytoplankton ecology. Hydrobiologia, 2003, 502(1-3), 395-403. [ Links ] PÁDISAK, J., BARBOSA, F.A.R., KOSHEL, R. and KRIENITZ, L. Deep layer Cyanoprokariota maxima in temperate and tropical lakes. Archiv für Hydrobiologie Special Issues Advances in Limnology, 2003, 58, 175-199. [ Links ] PADISÁK, J., CROSSETTI, L.O. and NASELLI-FLORES, L. Use and misuse in the application of the phytoplankton functional classification: a critical review with updates. Hydrobiologia, 2009, 621(1), 1-19. [ Links ] REYNOLDS, C.S. The ecology of freshwater phytoplankton. Cambridge: Cambridge University Press, 1984. [ Links ] REYNOLDS, C.S. Vegetation processes in the pelagic: a model for ecosystem theory . Oldendorf/Luhe: Ecology Institute, 1997. [ Links ] REYNOLDS, C.S., HUSZAR, V., KRUK, C., NASELLI-FLORES, L. and MELO, S. Towards a functional classification of the freshwater phytoplankton. Journal of Plankton Research , 2002, 24(5), 417-428. [ Links ] REYNOLDS, C.S., TUNDISI, J.G. and HINO, K. Observations on a metalimneticLyngbya population in a stably stratified tropical lake (Lagoa Carioca, eastern Brasil). Archiv für Hydrobiologie, 1983, 97, 7-17. [ Links ] RIGOSI, A. and RUEDA, F. J. Propagation of uncertainty in ecological models of reservoirs: from physical to population dynamic predictions. Ecological Modelling , 2012, 247:199-209. [ Links ] SALMASO, N. Life strategies, dominance patterns and mechanisms promoting species coexistence in phytoplankton communities along complex environmental gradients. Hydrobiologia , 2003, 502(1-3), 13-36. [ Links ] SALMASO, N., NASELLI-FLORES, L. and PADISÁK, J. Functional classifications and their application in phytoplankton ecology. Freshwater Biology, 2015, 60(4), 603-619. [ Links ] SHEPHERD, G.J. FITOPAC1: manual de usuário. Campinas: Universidade Estadual de Campinas, 1996. [ Links ] SOARES, M.C.S., HUSZAR, V.L., MIRANDA, M.N., MELLO, M.M., ROLAND, F. and LÜRLING, M. Cyanobacterial dominance in Brazil: distribution and environmental preferences. Hydrobiologia, 2013, 717(1), 1-12. [ Links ] SOUZA, M.B.G., BARROS, C., BARBOSA, F.A.R., HAJNAL, E. and PADISÁK, J. Role of atelomixis in replacement of phytoplankton assemblages in Dom Helvécio Lake, Southeast Brazil. Hydrobiologia, 2008, 607(1), 211-224. [ Links ] SUN, J. and LIU, D. Geometric models for calculating cell biovolume and surface area for phytoplankton. Journal of Plankton Research, 2003, 25(11), 1331-1346. [ Links ] TAVERA, R. and MARTÍNEZ-ALMEIDA, V. Atelomixis as a possible driving force in the phytoplankton composition of Zirahuén, a warm-monomictic tropical lake. Hydrobiologia, 2005, 533(1-3), 199-208. [ Links ] TOLEDO JUNIOR, A.P., TALARICO, M., CHINEZ, S.J. and AGUDO, E.G. 1983. The application of simple models for evaluating eutrophication processes in tropical lakes and reservoirs. In: Annals of the Brazilian Congress of Sanitary and Environmental Engineering . Rio de Janeiro: ABES, pp. 1-34. [ Links ] UTERMÖHL, H. Zur Vervolkomnung der quantitative Phytoplankton Methodik. Mitteilungen Internationale Vereinigung Theoretische und Angewandte Limnologie , 1958, 9, 1-38. [ Links ] WETZEL, R.G. and LIKENS, G.E. Limnological analyses. New York: Springer-Street, 2000. [ Links ] ZOHARY, T., PADISÁK, J. and NASELLI-FLORES, L. Phytoplankton in the physical environment: beyond nutrients, at the end, there is some light. Hydrobiologia, 2010, 639(1), 261-269. [ Links ] Received: June 30, 2017; Accepted: June 04, 2018 * e-mail:   This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This is a preview of a remote PDF:

Luciana Gomes Barbosa, Francisco Antônio Rodrigues Barbosa, Carlos Eduardo de Mattos Bicudo. Is thermal stability a factor that influences environmental heterogeneity and phytoplankton distribution in tropical lakes?, Acta Limnologica Brasiliensia, DOI: 10.1590/s2179-975x8817