HUANG, LI, WANG, HAN, LI, XIAO, and LIU: Responses of phytoplankton functional groups to environmental factors in the Maixi River, southwest China

Responses of phytoplankton functional groups to environmental factors in the Maixi River, southwest China


The functional groups approach is an efficient way to analyze seasonal changes in phytoplankton biomass as it is based on the physiological, morphological, and ecological attributes of the species. In this study, we identified the functional groups and driving factors behind short-term succession in phytoplankton communities. We analyzed physical, chemical, and biological factors of the Maixi River in Baihua Reservoir (BHR) from August to September, 2013 (summer, phase I) and March to May, 2014 (late spring and early summer, phase II). The 226 samples collected were divided into 23 functional groups. In phase I, phytoplankton biomass ranged from 4.88 to 30.59 mg·L–1, and the group S1 (Pseudanabaena limnetica) had the greatest biomass. In phase II, phytoplankton biomass ranged from 2.22 to 50.61 mg·L–1, and groups Y (Cryptomonas sp.) and S1 (P. limnetica) had the greatest biomass. Dominant functional groups in the Maixi River changed from S1 + D + Y + Lo in phase I to Y + S1 in summer. Changes in the phytoplankton community varied between 0 and 0.144 day–1 in phase I and between 0.008 and 0.389 day–1 in later spring and early summer. This showed a steady-state phytoplankton community during the two phases, in which the functional groups S1 (P. limnetica) and Y (Cryptomonas sp.) were dominant. Pseudanabaena limnetica, Synedra sp., Peridinium sp., and Cryptomonas sp. were dominant during summer, contributing to 70% of the total biomass in the steady-state community, and P. limnetica, Synedra sp., Cryptomonas sp., and Chlamydomonas sp. were dominant during later spring and early summer, contributing to 60% of the total biomass in the community. Groups S1, D, and Y formed easily in the Maixi River, but P. limnetica was the dominant species in the eutrophic conditions of the Maixi River. According to biotic and abiotic factors, we concluded that the Maixi River is hypertrophic, and water resource management should take blooms of P. limnetica occurring in May into account. Temperature and dissolved oxygen were the critical factors affecting the steady-state of the phytoplankton community in late spring and early summer and summer, respectively. Because the Maixi River is an important source in the BHR, its phytoplankton functional groups directly affect the ecological balance of the water environment.


Phytoplankton is usually the primary producer in aquatic habitats, and its community composition and populations respond directly and rapidly to changes in the environment (Reynolds, 2006). Density and specific composition reorganization of the phytoplankton community are responses to the interactive processes of physical, chemical, and biological variables (Crossetti and Bicudo, 2008a). An approach to studying phytoplankton communities is by defining phytoplankton functional groups, which are based on phytoplankton ecological characteristics and growth strategies (Reynolds et al., 2002). Functional groups include species that have similar physiology, morphology, and ecology, and frequently coexist. Moreover, functional groups share similar tolerances and sensitivities (Salmaso et al., 2015) and, at present, 38 functional groups are described using alphanumerical codes (Padisák et al., 2009). Compared to phylogenetic groups, functional groups are more effective in evaluating the phytoplankton responses to changes in environmental conditions (Kruk et al., 2002). The phytoplankton functional groups approach has also been applied to reservoirs (Nixdorf et al., 2003; Naselli-Flores and Barone, 2003; Crossetti and Bicudo, 2005; Fonseca and Bicudo, 2008; Crossetti and Bicudo, 2008b; Becker et al., 2009; Xiao et al., 2011), estuaries (Smayda and Reynolds, 2003; Anderson and Rengefors, 2006; Costa et al., 2009), lakes (Kruk et al., 2002; Mieleitner et al., 2008), and subtropical coastal system (Bonilla et al., 2005).

Human activities have increased the nitrogen and phosphorus loads into water, accelerating eutrophication, and affecting the biological organization (Howarth, 1998). Phytoplankton growth and density can determine the potential productivity of the entire ecosystem (Wissel and Fry, 2005), by responding quickly to changes in the aquatic environment and indicating when physical, chemical, or biological changes occur in the water column (Reynolds, 2006). Succession is the temporal change in the phytoplankton community and ecosystem, and a “steady-state phytoplankton community” is likely to appear in a period of relatively stable environmental factors. The definition of “steady-state” has been disputed: some researchers suggested that total biomass should allow a deviation of ±15% (Mischke and Nixdorf, 2003); others suggested that dominant species should represent at least 50–80% of the total abundance, and that the steady-state period should last at least two weeks (Dokulil and Teubner, 2003; Morabito et al., 2003). The steady-state concept of succession emerged in phytoplankton ecology to explain the diversity-disturbance relationship (Naselli-Flores et al., 2003). Establishing equilibrium requires consideration of many different factors, such as water depth, stratification, spatial heterogeneity, disturbance, stress, and the morphological and physiological plasticity of phytoplankton species (Naselli-Flores et al., 2003). Padisák (2003) studied 80 Hungarian lakes and found equilibrium phases in only 17 lakes.

In this study, we analyzed phytoplankton functional groups, as well as phytoplankton composition, density, biomass, chlorophyll a (Chl a), and physical and chemical variables to assess the driving forces behind two shortterm succession phases in the Maixi River, based on the adaptations of functional groups to environmental changes. An advantage of a short-term study is that the space-and-time substitution problem is avoided (Cardoso and Marques, 2003). Our primary objectives were to understand phytoplankton response mechanisms to environmental changes, by evaluating the phytoplankton functional groups in two succession phases and if there is a steady-state period of phytoplankton in this succession process. We expected variability in phytoplankton succession to be driven by the environment factors, which were steady-state phenomena during the short-term duration of the study and, therefore, functional groups would allow assessing the ecological status of the Maixi River.


Study site

Maixi River is one of the four main tributaries of the Baihua Reservoir (BHR), flowing into its wharf area. It is a shallow (6-8 m depth), and the water quality is affected by the surrounding sewage from residents, agriculture, and industry. The BHR is located within the catchment of the Maotiao River, which is a branch of the Wujiang tributary of the Yangtze River in China, and acts as a second cascade of the hydropower station on the Maotiao River and provides drinking water for Guiyang City (Li et al., 2013). As the linkage between rivers and lakes, Maixi River ecological environment has a direct impact in phytoplankton composition structure and water quality in BHR. Previous studies showed that frequent activities, large amount of pollutants, and complex hydrological conditions affected phytoplankton composition structure of the Maixi River (Li et al. 2011, 2013).

The sampling site (26°65’78.02’’ N, 106°54’28.06’’ E) was at the mouth of the Maixi River (Fig. 1). Sampling was conducted from August to September 2013 (summer, phase I) and from March to May 2014 (late spring and early summer, phase II). Daily water samples were taken at the surface (0.5 m layer) and weekly water samples were taken at depth (from the 0.5 m to the bottom, at 1 m intervals). Phase I and phase II sampling of surface water samples lasted 43 and 63 days, respectively. Stratified water samples were divided into seven layers (0.5-6 m) in phase I and into nine layers (0.5-8 m) in phase II; each phase was sampled eight times. This sampling strategy had two purposes: understanding the succession of phytoplankton from spring to summer and verifying the existence of a steady-state in the process of succession. Thus, a short-term continuous monitoring of phytoplankton was needed to monitor rate changes in phytoplankton communities.

Sample collection, processing, and evaluation

The 105 surface samples and the 128 vertically stratified samples collected in the two phases were fixed with formalin (3-5%) in the field. Firstly, each sample was concentrated by sedimentation (50 mL) after sedimentation for 48h in the laboratory, registering the volume before and after concentrate. After complete mixing, 0.1 mL of concentrated sample was counted in a 0.1 mL counting chamber under microscope to quantify phytoplankton. The phytoplankton were counted on a cell-by-cell basis. Taxa were identified, and the number of individuals of each taxon was counted under a direct microscope (Olympus BX41). Phytoplankton biomass was estimated using phytoplankton density and specific volumes, assuming a phytoplankton- specific density of 1 mg L−1 (Sunjun et al., 1999). Phytoplankton functional groups were used as described by Reynolds et al. (2002) and Padisák et al. (2009).

Physical and chemical parameters

Water temperature (WT), dissolved oxygen (DO), and pH were measured using an YSI meter (6600V2) on site. Transparent depth (SD) was measured using a Secchi disk. Physical and chemical characteristics of water, including total nitrogen (TN), ammonium (NH4-N), nitrite (NO2 –N), nitrate (NO3-N), total phosphorus (TP), dissolved inorganic phosphorus (PO4-P) and chlorophyll a (Chl.a), were determined using Chinese standard methods for water quality analysis GB3838-2002.

Phytoplankton data and statistical analyses

The rate of succession of the community composition (σ) was calculated according to the methods of Cardoso and Marques (2003) and Pannard et al. (2008):


where bi(t) is the density of the n species; B(t) is the sum of individuals constituting the community sampled; and t1 and t2 are phase I and phase II, respectively. The rate value (σ) is low when the composition of the community does not fluctuate much and high values of σ generally represent abrupt changes (Cardoso and Marques, 2003).

The relationship between the vertical distribution of phytoplankton and environmental factors was drawn on an isoline map using Surfer 8.0. The relationships among environmental variables and phytoplankton functional groups were evaluated using Pearson’s correlation in SPSS 19.0. Ordination analysis was performed in CANOCO 4.5. A detrended correspondence analysis (DCA) was first run to test whether the explored gradient fitted a linear or unimodal model, which allowed deciding on using a redundancy analysis (RDA). This RDA, which is a direct gradient analysis technique, was conducted to analyze the relationship between environmental factors and phytoplankton functional groups in the Maixi River, to determine the main environmental factors affecting phytoplankton functional groups succession. Linear response models were used for the ordination analyses conducted within short-term gradients. More than 1% of total biomass was used for the statistical analysis of data. For the RDA, environmental variables (except pH) and biological data were log (x + 1)-transformed to reduce skewness. The statistical significance of environmental variables explaining the variance in phytoplankton species obtained in the RDA was tested using a Monte Carlo permutation test and considering a P-level of 0.05.


Changes in phytoplankton composition and community

Among the 101 species identified in the two phases, Chlorophyta were the most numerous, accounting for 56.43% of the total species count, followed by Bacillariophyta (20.79%) and Cyanobacteria (14.85%). Pyrrophyta, Cryptophyta, Chrysophyta, and Euglenophyta species were less abundant, accounting for approximately 8%.

Average biomass and density of vertical and surface were higher in phase I than in phase II (Fig. 2). In phase I, vertical biomass ranged from 4.88 to 30.59 mg L–1 (mean 16.09 mg·L–1), reaching the minimum at 6 m during the 3th week and the maximum at 3 m during the 7th week. In phase II, vertical biomass ranged from 2.22 to 21.59 mg L–1 (mean 8.60 mg·L–1), reaching its minimum value at 8 m during the 4th week; the phytoplankton biomass was lower during the first three weeks and achieved its maximum value at 2 m during the 8th week.

Surface biomass ranged from 8.18 to 29.05 mg L–1 (mean 15.40 mg L–1) in phase I, registering the minimum value at the 7th day and increasing from the 28th to the 35th day, when it reached its maximum value, and then decreased. In phase II, surface biomass ranged between 2.46 and 50.61 mg L–1 (mean 17.03 mg L–1), reaching the highest value in the 67th day (May 27) and the lowest in the 50th day (March 27). From the 1th to the 51th day, biomass fluctuated slightly around 10 mg L–1, but after the 52th day (May 12) it increased gradually (Fig. 3).

Phytoplankton density ranged from 41.19x106 to 184.40x106 cells L–1 during phase I; densities were similar before the 28th day, presented a slight increase from the 29th to the 39th day. In phase II, phytoplankton density ranged from 8.26x106 to 473.09x106 cells L–1 density increased slowly until the 25th day, and rapidly from the 50th to 67th day (Fig. 3).

Phytoplankton biomass was positively correlated with DO (r = 0.289, P<0.05, n=55) and negatively correlated with TN (r = -0.315, P<0.05, n=55) and TP (r = -0.307, P<0.05, n=55) in phase I. In phase II, the biomass of phytoplankton was positively correlated with WT (r = 0.480, P<0.01, n=72) and negatively correlated with TN (r = - 0.367, P<0.01, n=72), NO3-N (r = -0.478, P<0.01, n=72), and NH4-N (r = -0.420, P<0.01, n=72).

Short-term succession characteristics

Phytoplankton taxa were sorted into 22 functional groups in phase I and 21 functional groups in phase II (Tab. 1). We found similar changes in the relative biomass of vertical and surface phytoplankton functional groups determined weekly and daily, respectively (Figs. 2 and 3).

In phase I, functional groups S1 (Pseudanabaena limnetica), D (Synedra sp.), Y (Cryptomonas sp.), P (Staurastrum sp.), and Lo (Peridinium sp.) considerably contributed to total biomass. Dominant phytoplankton species of surface phytoplankton in phase I accounted for 70.85% of the total biomass. Among this, S1 biomass ranged from 17.39 to 61.23% of the total biomass, followed by Lo (19%), P (17%), Y (10%), and D (15% before the 29th day and then gradually reduced to around 5%) (Fig. 4a). Early in phase I (1-32th day), relative biomass of S1 was low (<10%), but it gradually increased to 25% from the 33th to the 48th day, and then peaked to about 60% (Fig. 4a).

In phase II, the total biomass of the dominant phytoplankton functional groups S1 (P. limnetica), D (Synedra sp.), Y (Cryptomonas sp.), F (Oocystis sp.), and X2 (Chlamydomonas sp.) at the surface ranged from 29.95 to 89.74% (average 69.82%). The Y and S1 functional groups alternately accounted for 10% to 33% of the total biomass, and groups F and X2 were dominant from March 17 to March 24 with an average relative biomass as high as 30%; Y accounted for 32% of the total biomass, but this was gradually reduced from the 25th to 57th day. Early in phase II (1-32th day), relative biomass of S1 was low (<10%) but it increased gradually to about 25% from the 33th to the 48th day, and then peaked to about 60% (Fig. 4b).

Short-term succession rate (σ)

Rates of change in the phytoplankton community were broader in phase II than in phase I. In summer (phase I) rates ranged from 0 to 0.144 day–1, attaining the highest value from 14th to 15th day (Fig. 5a). Rates of change were higher from 11th to 30th day, and relatively stable from the 1th to 10th day and from 31th to 41th day (Fig. 5a). In phase II, rates of change were between 0.009 and 0.370 day–1 (Fig. 5b), being higher from the 1th to 53th day and relatively stable from the 54th to 67th day (Fig. 5b). In phase I, the rate of change was the largest from the 8th to 15th day (Fig. 5a). In phase II, the rate of change was variable (Fig. 5b), and changes in density and relative biomass were also large (Fig. 3; Fig. 4). The relative density of S1, for example, averaged 31.87% from the 1th to 30th day (Fig. 5b).

Main environmental factors driving at the short-term succession

Average water temperature was higher in phase I than in phase II and there was no thermal stratification of water column in any of the two seasons. In phase I (Fig. 6a), water temperature was usually above 21°C, despite a reduction in 5th and 8th week. In phase II, the temperature ranged between 9.2 and 21.3°C, with an average of 16.5°C (Tab. 2), gradually increasing from the 1th to 5th weeks, and slowly decreasing from the surface to the bottom after the second week; water temperature was always lower at the bottom (6-8 m) than at other depths (Fig. 6a).

In phase I, DO concentrations were lower than in phase II, ranging from 2.14 to 7.75 mg L–1, with an average of 5.17 mg L–1 (Tab. 2). The distributions of DO were similar to that of water temperature, with no obvious vertical stratification in the first two weeks (Fig. 6b). The lowest concentration of DO was registered at the 3th week, and the highest on the 3th and 4th weeks. In phase II, DO ranged from 3.91 to 10.01 mg L–1, with an average of 7.31 mg L–1 (Tab. 2). In phase I, pH varied from 7.39 to 8.64, with an average of 8.16; in phase II, it ranged from 7.29 to 8.5, with a mean value of 8.08. There was a gradual decrease of pH with depth in both phases (Fig. 6c).

No stratification was found for TP and TN (Fig. 6 d,e). The TP concentrations ranged from 10 to 120 μg L–1, with an average of 16 μg L–1 in phase I and 39 μg L–1 in phase II (Tab. 2). The lowest TP concentration occurred on the 3rd week of phase I (Fig. 6d). Because PO4-P concentration was extremely low, it was not a main component of TP. TN was between one and two orders of magnitude higher than TP (Tab. 2). Nitrate (NO3-N) was the main component of TN, its strongest correlations with TN obtained by Pearson’s correlation analysis were r = 0.741 (P<0.01, n=55) and r = 0.891 (P<0.01, n=55). The SD was generally low, with an average of 1.1 m in phase I and an average of 1.36 m in phase II (Tab. 2).

RDA of phytoplankton functional groups and environmental factors

The results of the DCA ordination with phytoplankton functional groups accounted for gradient lengths of 1.267 and 1.542 in phase I and II, respectively. Thus, linear ordination methods such as RDA should be used. In the diagram, the cosine of the angle of environmental factors and species variables represents the correlation between them; an acute angle is a positive correlation and an obtuse angle is a negative correlation. There were 9 environmental factors in both phases, (TP, TN, NO3-N, NO2-N, NH4-N, WT, DO, pH and SD), and 12 functional groups in phase I and 16 functional groups in phase II. In phase I, the Monte Carlo test (P<0.05) demonstrated that the ordination along axis 1 was statistically significant with an eigenvalue of 0.284; 75.8% of the cumulative variance in species distribution was explained by the first two axes (Tab. 3). The most effective explanatory factors were SD (P<0.01, F=8.51, n=499) and TN (P<0.01, F=5.224, n=499), both playing a significant role in the short-term succession. Correlations also showed that WT was the most important variable (P<0.01, F=5.134, n=499). Regarding phytoplankton functional groups, S1 was located towards higher DO and SD values, whereas X2, P, T and LO towards higher SD. Most phytoplankton functional groups (T, Lo, P, X2, LM, Y, F and J) were closely related to four environmental factors, appearing at the bottom left of the diagram. In addition, D and S2 were placed towards higher values of NO2-N, G and B were towards higher NH4-N concentrations, MP was placed towards of water temperature (Fig. 7a).

In phase II, the Monte Carlo test (P<0.05) demonstrated that axis 1 was statistically significant (P<0.01) with an eigenvalue of 0.325; 61.1% of the cumulative variance in species distribution was explained by the first two axes (Tab. 3). The most effective explanatory factors were WT (P<0.01, F=18.27, n=499), NO3-N (P<0.01, F=11.88, n=499), SD (P<0.01, F=11.24, n=499), NH4-N (P<0.01, F=6.66, n=499), and TN (P<0.01, F=5.12, n=499), and played significant roles in the short-term succession of phytoplankton functional groups. While WT (0.684), NO3-N (-0.566) and SD (0.531) were highly correlated with axis 1, TP (0.438), pH (0.317), and SD (- 0.398) were highly correlated with axis 2. Most phytoplankton functional groups (S1, B, X1, Lo, J, SN, F, X2, J and Y) were located towards higher WT, NO2-N and SD. In addition, G, MP, P and D were located towards higher DO, pH, and TP, whereas D, LM, and P towards higher TN and TP concentrations (Fig. 7b).


The short-term succession of phytoplankton response mechanisms to environmental changes

Phytoplankton functional group composition in the Maixi River had particular spatial and temporal distribution characteristics. Our study site is connected to the river and lake ecosystems; thus, its hydrological environment is complex, which influences phytoplankton community dynamics. There was no thermal stratification of water during the study period, phytoplankton community changes were similar between layers, as was observed in shallow water by Nixdorf and Deneke (1997). The dominant functional groups in late spring and early summer were Y and S1, but in summer S1 absolutely dominated. Based on the phytoplankton steady-state conditions, we defined a steady-state in summer and two steady-state periods in late spring and early summer. The steady-state of dominant phytoplankton was controlled by many factors, including nutrients, light, and water stagnation, which changed frequently from deficiency to superfluity (Mischke and Nixdorf, 2003), and by coexistence mechanisms in the relative equilibrium state. There were competition and coexistence mechanism of phytoplankton; both S1 and Y tolerate highly light-deficient conditions, but the first one is sensitive to flushing, contrarily to D that is tolerant to it. In Maixi river S1 is well adapted to the less transparency conditions. Functional group S1 is usually characterized by cyanobacteria, are adapted to high temperature, forming groups, floating, and multiplying fast, which dominate during the summer period (Çelekli and Öztürk, 2014). Some studies have shown that populations of Cyanobacteria, which comprise a few species within the functional group of S1, regularly occur in summer in different water bodies around the world (Rücker et al., 1997; Padisák et al., 2003; Gemelgo et al., 2009; Borics et al., 2012; Li et al., 2013; Stević et al., 2013; Hu et al., 2016). Both Lo and Y have a flagellum that allows their survival in multiple water layers when adverse conditions such as weather changes and water body disturbance are encountered (Padisák et al., 2009; Naselli-Flores and Barone, 2003).

Seasonal variation in phytoplankton is closely related to changes in the water environment, such as sudden changes in temperature, nutritional status, disturbance patterns, zooplankton grazing pressure, and hydrodynamics (Padisák et al., 2003; Wilk-Woźniak and Żurek, 2006; Marija et al., 2007). In fact, changes in water temperature was the main environmental factor affecting the growth of S1 during late spring and early summer (i.e., phase II). Many studies also pointed out that water temperature plays a major role in seasonal changes, growth, and community structure of phytoplankton (Padisák et al., 2003; Salmaso and Zignin, 2010). In summer, with continuously high temperature, this variable did not have a significant effect and DO was the main factor influencing the rate of succession.

Besides, ratios of TN/TP have also been used to examine nutrient limitations on phytoplankton growth (Becker et al., 2010). Given the high TN concentrations observed during the two study periods, our results suggest that TP was probably a limiting factor for phytoplankton growth; the N present in the Maixi River was mainly in the form of NO3-N, which reflected the amount of exogenous pollutants in the water. In addition, this showed that the water self-purification process was ongoing, with the high ammonia utilization ratio leading to increasing nitrate contents. Therefore, these stable and suitable environmental conditions are able to the establishment of a steady state for a longer time, and a less marked rate of phytoplankton change. Our hypothesis was confirmed, the variability in phytoplankton succession to be driven by the physical and chemical factors, which were steady-state phenomena during the short-term duration of the study.

Reaction of phytoplankton functional groups to ecological status

Because the phytoplankton functional groups approach can be used without geographic limitations (Padisák et al., 2009; Crossetti and Bicudo, 2008b; Becker et al., 2009; Pasztaleniec and Poniewozik, 2010), and habitat requirements and phytoplankton are sensitive to changes in water quality (Çelekli and Öztürk, 2014), it is a potential monitoring tool for the assessment of ecological status in the context of the Water Framework Directive (Cellamare et al., 2012). For instance, phytoplankton functional groups estimated reliable water quality states between hypertrophic and oligotrophic conditions in a reservoir of Turkey (Çelik and Sevindik, 2015), indicated eutrophication due to human impacts along the river Loire in France (Abonyi et al., 2012), and allowed an ecological assessment of French Atlantic lakes (Cellamare et al., 2012).

This system is hypereutrophic based on chlorophyll, SD, TN and sometimes due to TP, and also the presence of the dominant functional groups S1, Lo, and P indicated this trophic condition. In addition, we found high values of biomass and density from spring to summer, which is the most likely period for the occurrence of phytoplankton blooms of P. limnetica (in May particularly). Many studies have confirmed that the worse the ecological condition of the water, the higher the thin filamentous cyanobacteria biomass (Hajnal and Padisák, 2008). The water management should take preventive measures against this phenomenon during this month in future. Baihua reservoir is one of the suppliers of drinking water for Guiyang City, and Maixi River has a direct impact on the water quality of this reservoir, threatening the safety of human life. On the other hand, the human activities in this area varied frequently, and a large amount of pollutants and complex hydrological conditions affected phytoplankton composition (Li et al., 2011, 2013). For that situation, we need to pay attention to control pollution in key river valleys and regions, and we should take appropriate measures to reduce nutrient input, especially during the period of the temperature rise and continuous rainfall.


  • This study confirmed the potential utility of the phytoplankton functional groups approach in the Maixi River. Two consecutive short-term sampling phases yielded 23 functional groups; phase I was dominated by the functional group S1, and phase II was dominated by the functional groups Y and S1.

  • Water temperature was a key factor in the selection of phytoplankton species and were likely to be the critical factors affecting phytoplankton communities in the steady-states occurring in later spring to early summer, whereas SD and DO were key factors during summer. Functional groups S1, D, and Y were in a steady-state dynamics in the Maixi River.

  • Based on phytoplankton functional groups and the environmental factors Maixi River can be classified as hypertrophic. Functional group S1 may begin to grow in May, and this period may coincide with phytoplankton blooms. Therefore, water resources management should focus on this period.


We thank two anonymous reviewers for helpful comments. This research was in part supported by the Natural Science Foundation of China (U1612442), the department of science and technology of Guizhou Province ([2015]2001, [2015]10, [2014]7001), Water Resources Department of Guizhou Province KT201401 and the department of science and technology of Baiyun ([2014]69, [2014]70).



A Abonyi, M Leitao, AM Lançon, J Padisák, 2012. Phytoplankton functional groups as indicators of human impacts along the River Loire (France). Hydrobiologia 698:233-249.


DM Anderson, K Rengefors, 2006. Community assembly and seasonal succession of marine dinoflagellates in a temperate estuary: The importance of life cycle events. Limnol. Oceanogr. 51:860-873.


V Becker, V Huszar, L Crossetti, 2009. Responses of phytoplankton functional groups to the mixing regime in a deep subtropical reservoir. Hydrobiologia 628:137-151.


V Becker, L Caputo, J Ordónes, R Marcé, J Armengol, LO Crossetti, VLM Huszar, 2010. Driving factors of the phytoplankton functional groups in a deep Mediterranean reservoir. Water Res. 44:3345-3354.


S Bonilla, D Conde, L Aubriot, MD Pérez, 2005. Influence of hydrology on phytoplankton species composition and life strategies in a subtropical coastal lagoon periodically connected with the Atlantic Ocean. Estuar. Coast. 28:884-895.


G Borics, B Tóthmérész, BA Lukács, G Várbíró, 2012. Functional groups of phytoplankton shaping diversity of shallow lake ecosystems. Hydrobiologia 698:251-262.


LDS Cardoso, DDM Marques, 2003. Rate of change of the phytoplankton community in Itapeva Lake (North Coast of Rio Grande do Sul, Brazil), based on the wind driven hydrodynamic regime. Hydrobiologia 497:1-12.


A Çelekli, B Öztürk, 2014. Determination of ecological status and ecological preferences of phytoplankton using multivariate approach in a Mediterranean reservoir. Hydrobiologia 740:115-135.


K Çelik, TO Sevindik, 2015. The phytoplankton functional group concept provides a reliable basis for ecological status estimation in the Çaygören Reservoir (Turkey). Doga Turk. J. Bot. 39:1226-1230.


M Cellamare, S Morin, M Coste, J Haury, 2012. Ecological assessment of French Atlantic lakes based on phytoplankton, phytobenthos and macrophytes. Environ. Monit. Assess. 184:4685-4708.


LS Costa, VLM Huszar, AR Ovalle, 2009. Phytoplankton functional groups in a tropical estuary: hydrological control and nutrient limitation. Estuar. Coast. 32:508-521.


LO Crossetti, CEM Bicudo, 2005. Structural and functional phytoplankton responses to nutrient impoverishment in mesocosms placed in a shallow eutrophic Reservoir (Garças Pond), Sao Paulo, Brazil. Hydrobiologia 541:71-85.


LO Crossetti, CEM Bicudo, 2008a. Phytoplankton as a monitoring tool in a tropical urban shallow reservoir (Garças Pond): the assemblage index application. Hydrobiologia 610:161-173.


LO Crossetti, CEM Bicudo, 2008b. Adaptations in phytoplankton life strategies to imposed change in a shallow urban tropical eutrophic reservoir, Garças Reservoir, over 8 years. Hydrobiologia 614:91-105.


MT Dokulil, K Teubner, 2003. Steady state phytoplankton assemblages during thermal stratification in deep alpine lakes. Do they occur? Hydrobiologia 502:65-72.


BM Fonseca, CEM Bicudo, 2008. Phytoplankton seasonal variation in a shallow stratified eutrophic reservoir (Garças Pond, Brazil). Hydrobiologia 600:267-282.


MCP Gemelgo, JLN Mucci, D Navaspereira, 2009. Population dynamics: seasonal variation of phytoplankton functional groups in Brazilian reservoirs (Billings and Guarapiranga, Sao Paulo). Braz. J. Biol. 69:1001-1013.


E Hajnal, J Padisák, 2008. Analysis of long-term ecological status of Lake Balaton based on the ALMOBAL phytoplankton database. Hydrobiologia 599:227-237.


RW Howarth, 1998. An assessment of human influences on fluxes of nitrogen from the terrestrial landscape to the estuaries and continental shelves of the North Atlantic Ocean. Nutr. Cycl. Agroecosys. 52:213-223.


R Hu, Q Li, B Han, L Naselli-Flores, J Padisak, N Salmaso, 2016. Tracking management-related water quality alterations by phytoplankton assemblages in a tropical reservoir. Hydrobiologia 763:109-124.


C Kruk, N Mazzeo, G Lacerot, CS Reynolds, 2002. Classification schemes for phytoplankton: a local validation of a functional approach to the analysis of species temporal replacement. J. Plankton Res. 24:901-912.


G Marija, A Plenković-Moraj, K Kralj, I Grigorszky, D PerošPucar, 2007. The relationship between phytoplankton species dominance and environmental variables in a shallow lake (Lake Vrana, Croatia). Hydrobiologia 584:337-346.


J Mieleitner, M Borsuk, HR Bürgi, P Reichert, 2008. Identifying functional groups of phytoplankton using data from three lakes of different trophic state. Aquat. Sci. 70:30-46.


U Mischke, B Nixdorf, 2003. Equilibrium phase conditions in shallow German lakes: how Cyanoprokaryota species establish a steady state phase in late summer. Hydrobiologia 502:123-132.


G Morabito, A Oggioni, P Panzani, 2003. Phytoplankton assemblage at equilibrium in large and deep subalpine lakes: a case study from Lago Maggiore (N. Italy). Hydrobiologia 502: 37-48.


L Naselli-Flores, J Padisák, MT Dokulil, I Chorus, 2003. Equilibrium steady-state in phytoplankton ecology. Hydrobiologia 502:395-403.


B Nixdorf, R Deneke, 1997. Why ‘very shallow’ lakes are more successful opposing reduced nutrient loads. Hydrobiologia 342/343:269-284.


B Nixdorf, U Mischke, J Rücker, 2003. Phytoplankton assemblages and steady state in deep and shallow eutrophic lakes: an approach to differentiate the habitat properties of Oscillatoriales. Hydrobiologia 502:111-121.


Q Li, L Chen, P Xia, S Liu, F Chen, D Yu, C Li, 2011. [Structure of phytoplankton community and its relationship with environmental factors at the estuary of Maixi River in Baihua Reservoir, Guizhou Province].[Article in Chinese]. J. Lake Sci. 23:612-618.


Q Li, L Chen, F Chen, T Gao, X Li, S Liu, C Li, 2013. Maixi River estuary to the Baihua Reservoir in the Maotiao River catchment: phytoplankton community and environmental factors. Chin. J. Oceanol. Limn. 31: 290-299.


J Padisák, G Borics, G Fehér, I Grigorszky., 2003. Dominant species, functional assemblages and frequency of equilibrium phases in late summer phytoplankton assemblages in Hungarian small shallow lakes. Hydrobiologia 502:157-168.


J Padisák, L Crossetti, L Naselli-Flores, 2009. Use and misuse in the application of the phytoplankton functional classification: A critical review with updates. Hydrobiologia 621:1-19.


A Pasztaleniec, M Poniewozik, 2010. Phytoplankton based assessment of the ecological status of four shallow lakes (Eastern Poland) according to Water Framework Directive-a comparison of approaches. Limnologica 40:251-259.


A Pannard, M Bormans, Y Lagadeuc, 2008. phytoplankton species turnover controlled by physical forcing at different time scales. Can. J. Fish. Aquat. Sci. 65:47-60.


CS Reynolds, 2006. Phytoplankton periodicity: the interactions of form, function and environmental variability. Freshwater Biol. 14:111-142.


CS Reynolds, VLM Huszar, C Kruk, L Naselli-Flores, S Melo, 2002. Towards a functional classification of the freshwater phytoplankton. J. Plankton Res. 24:417-428.


J Rücker, C Wiedner, P Zippel, 1997. Factors controlling the dominance of Planktothrix agardhii and Limnothrix redekei in eutrophic shallow lakes. Hydrobiologia 342/343:107-115.


E Wilk-Woźniak, R Żurek, 2006. Phytoplankton and its relationships with chemical parameters and zooplanktonin the meromictic Piaseczno reservoir, Southern Poland. Aquat. Ecol. 40: 165-176.


B Wissel, B Fry, 2005. Tracing Mississippi River influences in estuarine food webs of coastal Louisiana. Oecologia 144:659-672.


N Salmaso, A Zignin, 2010. At the extreme of the physical gradients: phytoplankton in the highly flushed, large rivers. Hydrobiologia 639:21-36.


N Salmaso, L Naselli-Flores, J Padisák, 2015. Functional classifications and their application in phytoplankton ecology. Freshwater Biol. 60:603-619.


TJ Smayda, CS Reynolds, 2003. Strategies of marine dinoflagellate survival and some rules of assembly. J. Sea Res. 49:95-106.


F Stević, M Mihaljević, D Špoljaric, 2013. Changes of phytoplankton functional groups in a floodplain lake associated with hydrological perturbations. Hydrobiologia 709:143-158.


J Sun, DY Liu, SB Qian, 1999. Study on phytoplankton biomass I. Phytoplankton measurement biomass from cell volume or plasma volume. Acta Oceanol. Sin. 21:75-85.


L Xiao, T Wang, R Hu, B Han, S Wang, X Qian, J Padisák, 2011. Succession of phytoplankton functional groups regulated by monsoonal hydrology in a large canyon-shaped reservoir. Water Res. 45:5009-5019.

Fig. 1.

Location of Baihua Reservoir and of the sampling site within Maixi River.

Fig. 2.

Vertical distribution of phytoplankton biomass x-axis: vertical sampling once a week, eight times during each study phase. Phase I, summer, from August to September 2013; phase II, late spring and early summer, from March to May 2014.

Fig. 3.

Changes in surface phytoplankton biomass and density during each sampling period x-axis: phase I and phase II sampling of surface water samples lasted 42 and 63 days.

Fig. 4.

Relative biomass (%) of surface functional groups. a) Others: X2, LM, F, B, C, E, G, H1, J, M, MP, S2, SN, T, TC, W1, W2, X1. b) Others: B, P, Lo, C, E, G, H1, J, LM, M, MP, S2, SN, T, TC, W1, W2, X1.

Fig. 5.

Rate of change in the phytoplankton community (σ day–1) during the two sampling phases in the Maixi River.

Fig. 6.

Spatiotemporal distribution of the values of the main environmental variables. a) Water temperature (°C); b) dissolved oxygen (mg L–1); c) pH; d) total phosphorus (mg L–1); e) total nitrogen (mg L–1).

Fig. 7.

Redundancy analysis of phytoplankton functional groups related to environmental factors in the Maixi River.

Tab. 1.

Phytoplankton functional groups and their typical habitat were used as described by Reynolds et al. (2002) and Padisák et al. (2009).

Functional groups Most important taxa Typical habitat Phase I Phase II
B Cyclotella sp. Vertically mixed, mesothrophic small medium lakes + +
C Asterionella formosa Melosira ambigua Mixed, eutrophic small medium lakes + +
D Synedra acus Shallow, enriched turbid waters, including rivers + +
E Dinobryon divergens Usually small, oligotrophic, base poor lakes or heterothrophic ponds + +
F Oocystis sp., Kirchnecriella sp. Clear epilimnia + +
G Eudorina elegans Short, nutrient-rich water columns + +
H1 Aphanizomenon sp. Dinitrogen-fixing, nostocaleans + +
J Scenedesmus sp., Pediastrum sp., Coelastrum sp. Shallow, enriched lakes ponds and rivers + +
LM Navicula sp., Achnanthes sp. Summer epilimnia in eutrophic lakes + +
LO Dactyloccopsis, Peridinium Summer epilimnia in mesotrophic lakes + +
M Microcystis sp. Daily mixed layers of small eutrophic, low latitude lakes + +
MP Chlorella vulgaris Frequently stirred up, inorganically turbid shallow lakes + +
P Staurastrum sp. Fragilaria sp. Eutrophic epilimnia + +
S1 Pseudanabaena limnetica Turbid mixed layers + +
SN Raphidiopsis, Cylindrospermopsis sp. Warm mixed layers + +
T Quadrigula sp. Deep, well-mixed epilimnia + +
TC Phormidium sp. Eutrophic standing waters, or slow- flowing rivers with emergent macrophytes +  
W1 Euglenophyta sp. Small organic ponds + +
W2 Trachelomononas sp. Shallow mesotrophic lakes + +
X1 Pyrrophyta, Chroococcus sp. Shallow mixed layers in enriched conditions + +
X2 Chlamydomonas sp. Shallow, clear mixed layers in meso-eutrophic lakes + +
Y Cryptophyta Usually, small, enriched lakes + +

[i] Phase I, August 16 - September 26, 2013; Phase II, March 21 - May 22, 2014.

Tab. 2.

Environmental variables in the Maixi River during the two periods.

Phase I Phase II
Unit Mean Range Mean Range
Water temperature °C 23.65 21.1~25.9 16.5 9.2~21.3
DO mg L–1 5.17 4.09~7.75 7.31 3.91~10.01
Conductivity μs m–1 470.11 420~545 532.04 498~689
pH 8.16 7.39~8.64 8.08 7.29~8.5
TP μg L–1 16 10~40 39 10~120
TN mg L–1 1.83 1.55~2.48 2.43 1.14~5.87
TN/TP 126.4 50.9~200.3 72.7 17.0~233.1
NO3-N mg L–1 1.28 1.007~2.17 1.86 0~4.81
NO2-N mg L–1 0.05 0.031~0.08 0.035 0.022~0.056
NH3-N mg L–1 0.011 0~0.05 0.057 0~0.37
PO4-P μg L–1 2 0~4 3 0~7
Chlorophyll a μg L–1 24.7 16.09~32.3 17.1 0.37~33.45
SD m 1.1 0.8~1.4 1.36 0.9~1.75
Tab. 3.

Main results of redundancy analysis and Monte Carlo permutation test for the relationship between biovolume of phytoplankton species and environment variables.

Phase I Phase II
Axes λ1 λ2 λ3 λ4 Total variance λ1 λ2 λ3 λ4 Total variance
Eigenvalues 0.284 0.087 0.038 0.022 1 0.325 0.040 0.028 0.022 1
Species-environment correlations 0.793 0.758 0.580 0.741 0.881 0.611 0.492 0.467
Cumulative percentage variance
of species data 28.4 37.1 40.9 43.1 32.5 36.6 39.3 41.5
of species-environment relation 62.7 81.9 90.2 95.0 74.0 83.1 89.4 94.3
Sum of all eigenvalues 1 1
Sum of all canonical eigenvalues 0.453 0.440
Abstract views:


Article Metrics

Metrics Loading ...

Metrics powered by PLOS ALM

Copyright (c) 2017 Guojia Huang, Qiuhua Li, Xiaoqing Wang, Mengshu Han, Lei Li, Jing Xiao, Yipeng Liu

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
© PAGEPress 2008-2018     -     PAGEPress is a registered trademark property of PAGEPress srl, Italy.     -     VAT: IT02125780185     •     Privacy