Quantifying photosynthetic performance of phytoplankton based on photosynthesis–irradiance response models

Clarifying the relationship between photosynthesis and irradiance and accurately quantifying photosynthetic performance are of importance to calculate the productivity of phytoplankton, whether in aquatic ecosystems modelling or obtaining more economical production. The photosynthetic performance of seven phytoplankton species was characterized by four typical photosynthesis–irradiance (P–I) response models. However, the differences were found between the returned values to photosynthetic characteristics by different P–I models. The saturation irradiance (Isat) was distinctly underestimated by model 1, and the maximum net photosynthetic rate (Pnmax) was quite distinct from its measured values, due to the asymptotic function of the model. Models 2 and 3 lost some foundation to photosynthetic mechanisms, that the returned Isat showed significant differences with the measured data. Model 4 for higher plants could reproduce the irradiance response trends of photosynthesis well for all phytoplankton species and obtained close values to the measured data, but the fitting curves exhibited some slight deviations under the low intensity of irradiance. Different phytoplankton species showed differences in photosynthetic productivity and characteristics. Platymonas subcordiformis showed larger intrinsic quantum yield (α) and lower Isat and light compensation point (Ic) than Dunaliella salina or Isochrysis galbana. Microcystis sp., especially M. aeruginosa with the largest Pnmax and α among freshwater phytoplankton strains, exhibited more efficient light use efficiency than two species of green algae. The present work will be useful both to describe the behavior of different phytoplankton in a quantitative way as well as to evaluate the flexibility and reusability of P–I models. Meanwhile we believe this research could provide important insight into the structure changes of phytoplankton communities in the aquatic ecosystems.


Background
Phytoplankton are a key functional component of aquatic ecosystems and play a pivotal role in biogeochemical cycles [1]. In particular, marine phytoplankton, as the principal driving force of ocean carbon cycles and energy flows, fix approximately 50 gigatons of inorganic carbon annually, almost half of the total global primary production [2,3]. They show higher CO 2 fixation rates and higher biomass productivity than any other photosynthetic organisms [3]. As the increase of CO 2 concentration in the atmosphere and global warming, an accurate estimate of photosynthetic productivity of phytoplankton becomes ever more important for modelling primary production and structure changes of phytoplankton communities in aquatic ecosystems, especially eutrophic lakes (e.g., Taihu, Erie, Winnipeg lake) and estuaries (e.g., Yangtze River).
Clarifying the relationship between photosynthesis and irradiance is a basis to evaluate the growth performance of phytoplankton. Irradiance acts as a driving force in photosynthesis. The level of irradiance affects the growth, CO 2 fixation efficiency, carbon metabolism, and cell composition of photosynthetic organisms [4][5][6][7][8]. High irradiance causes photoinhibition by the production of reactive oxygen species (ROS) and damages the function of the most light-sensitive complex PSII [5]. While extensive studies have been carried out and many insights have enriched the basis of phytoplankton physiology in recent decades [9][10][11], the relationship remains poorly understood for phytoplankton. Irradiance availability also affects phytoplankton community composition and is one of the key factors causing cyanobacteria blooms [12]. Resource competition theory shows that species with lower "critical light intensity" are often superior, such as Microcystis [13].
On the other hand, phytoplankton cells are rich in proteins, polysaccharides, lipids, vitamins, and polyunsaturated fatty acids, which have stirred up great attention as a promising potential feedstock for biofuel, nutraceuticals, animal and aquaculture feed production [10,14]. Many species have been used for commercial development, such as Dunaliella salina, Isochrysis galbana, Spirulina (or Arthrospira), Haematococcus pluvialis, and Scenedesmus obliquus [2,6,10]. Almost all fishes, bivalve molluscs, and crustaceans primarily graze on phytoplankton to build immunity against diseases during their early larval stages [12]. However, large-scale production of phytoplankton has rarely been successful, with no more than 1 g DW L −1 biomass that is mainly limited by the inefficiency of photosynthesis in high-cell density cultivation [11,14,15]. The photosynthetic parameters can be seen as indicators to achieve sustainable carbon assimilation and TAG accumulation in Isochrysis zhangjiangensis [8]. Therefore, accurately quantifying photosynthetic performance is crucial for more economical integration of production management and operation of industrial-scale phytoplankton culture systems [16].
The response curve of photosynthesis to irradiance (P-I) has been frequently used to characterize photosynthetic performance by fitting experimental data (measured as oxygen evolution or carbon uptake) with P-I models [17]. Obtained photosynthetic parameters, including the maximum net photosynthetic rate (P nmax ), the optimal intensity of irradiance (I sat ), and the dark respiration rate (R d ) can be regarded as indicator to evaluate the response of organism to meet environmental changes. In phytoplankton, some classical models describing the P-I curves were developed by Webb et al. [18], Platt et al. [19][20][21] and Eilers and Peeters [22], respectively. Later, some mechanistic models characterizing P-I curves with considering the underlying biophysical processes of photosynthesis were mainly used to investigate properties of photosynthetic physiology for phytoplankton [23][24][25][26][27][28], for example, photoinhibition, photoacclimation and dynamic down-regulation of photosynthesis. However, most of mechanistic models have many parameters and complex relationship between photosynthetic rate and irradiance. Consequently, the most extensive application is still found in those classical models. For example, an examination of the literature overwhelmingly reveals in excess of 1950 papers on the double exponential model proposed by Platt et al. [20]. This is most probably because these classical models are simpler than those new models with many complex parameters and processes, and thus those new models would certainly take many years to be fully adopted. Higher plant and phytoplankton possess similar photosynthetic systems. Recently, Ye et al. developed a mechanistic model for higher plants that parameterizes the core characteristics of photosynthesis to highlight processes, including solar energy absorption of photosynthetic pigment molecules, energy transfer, and electron transport between photosynthetic apparatuses [29]. The model with four parameters is relatively simple, and has been widely applied in rice, wheat, soybean, sunflower and other plants [30,31].
The objective of this study was to determine the various relationships between the photosynthetic productivity of phytoplankton and irradiance intensity and investigate the reliability of P-I models to estimate the photosynthetic performance for phytoplankton. We selected the rather extensive range of phytoplankton, including three isolated from the ocean and four from lakes, to measure their photosynthetic oxygen evolution under different irradiance intensity. Obtained P-I data were fitted by using P-I model for quantization the photosynthetic performances. The P-I model for higher plants developed by Ye et al. (it was represented as model 4 in this study) was first used to compare against three most widely applied models for phytoplankton (them were represented as models 1, 2 and 3 in this study, respectively).

Phytoplankton cultivation
The three strains of marine phytoplankton (Isochrysis galbana, Dunaliella salina and Platymonas subcordiformis) isolated from East China Seas were grown aseptically in f/2 medium. The four strains of freshwater phytoplankton (Microcystis aeruginosa FACHB-905, Microcystis wesenbergii FACHB-1112, Scenedesmus obliquus FACHB-116 and Chlorococcum sp. FACHB-1556) were purchased from the Freshwater Algae Culture Collection (FACHBcollection) of the Institute of Hydrobiology, Chinese Academy of Sciences (Wuhan, China) and cultivated in BG11 medium. The cultures were illuminated by cool white fluorescent bulbs (60 µmol photons m −2 s −1 ) with a photoperiod of 12 h per day at 26 ± 1 °C.

Measurement of photosynthetic oxygen evolution
After 7 to 10 days of incubation, the photosynthetic oxygen-evolving rate of microalgal cells reaching the exponential growth phase was determined using a bio-oxygen meter (Yaxin-1151, Beijing Yaxinliyi Science and Technology Co., Ltd., China). Eight mL cell suspensions of each strain were exposed to increasing orders of irradiance intensity (0, 25, 50, 100, 150, 200, 300, 400, 500, 600, 800, 1000, and 1200 µmol photons m −2 s −1 ), given by a digital LED light source (YX-11LA, Beijing Yaxinliyi Science and Technology Co., Ltd., China), at 25 ± 1 °C. The meter took reads once every 3 s for 5 min in each irradiance measurement point, during which a linear relationship varying with time in oxygen concentration was obtained. Triplicate samples were prepared and measured for each test. The response of the photosynthetic oxygen-evolving rate to irradiance (P n -I) was fitted with four P-I models [19][20][21][22]29].

Determination of chlorophyll a concentration and cell counts
The cells for photosynthetic oxygen-evolving measurement were collected by centrifugation (5600×g) for 10 min at 4 °C. Chlorophyll a (Chl a) was extracted from microalgal cells in 90% (v/v) acetone and left overnight at 4 °C in darkness. The extracts were then centrifuged at 3600×g for 10 min. The Chl a concentration was determined spectrophotometrically in the supernatant with a SP752 UV-Vis spectrophotometer (Spectrum Instruments, Shanghai, China) according to the method of Jeffrey and Humphrey [32]. One-mL cultures of each strain were taken and preserved in Lugol's iodine solution for counting algal cells by a haemocytometer. Each test was conducted in triplicate.

Model 1
The light dependence of the net photosynthetic rate (P n ) is expressed as [19]: where P n (μmol O 2 mg −1 Chl a h −1 ) is the chlorophyll a-normalized net photosynthetic rate at irradiance I, P nmax (μmol O 2 mg −1 Chl a h −1 ) is the light-saturated maximum rate of photosynthesis, α (μmol O 2 mg −1 Chl a h −1 /μmol photons m −2 s −1 ) is the light-limited initial (1) P n = P nmax tanh αI P nmax − R d , slope of P n -I curve, and R d (μmol O 2 mg −1 Chl a h −1 ) is the dark respiration rate. As Eq. (1) is an asymptotic function, the saturation irradiance cannot be directly calculated. Therefore, the saturation irradiance (I sat , μmol photons m −2 s −1 ) is obtained by drawing a line from tangent of the initial slope with the plateau of the P n -I curve onto the x-coordinate [1]. I sat is calculated by the following calculation formula: But the analytic solution of the light compensation point (I c , μmol photons m −2 s −1 ) cannot be directly obtained by Eq. (1). In order to obtain I c , Kok effect [33] must be ignored here, and then I c can be calculated as [21]: The photosynthetic quantum efficiency (P n ′, μmol O 2 μmol photons −1 ) is calculated as:

Model 2
The light dependence of P n is expressed as [20,21]: where P n is the chlorophyll a-normalized net photosynthetic rate at irradiance I; α is the light-limited initial slope of P n -I curve; β is the dimensionless parameter reflecting the photoinhibition process; Without photoinhibition, P s is the maximum photosynthetic output; P s is the parameter reflecting the maximum, potential, lightsaturated, rate of photosynthesis at β > 0; and R d is the dark respiration rate.
The I sat is calculated as: The P nmax can be calculated as: However, the analytic solution of I c cannot be directly obtained by Eq. (5). To obtain I c , the Kok effect must be ignored here, and then I c can be calculated as: (3) The photosynthetic quantum efficiency is calculated as:

Model 3
The light dependence of P n is expressed as [22]: Here, P n is the chlorophyll a-normalized net photosynthetic rate at irradiance I; α and β are the fundamental parameters, nondimensional; and R d is the dark respiration rate. The reciprocal of γ is the light-limited initial slope of P n -I curve.
I sat is calculated as: P nmax is given by: When P n = 0, I c is given as follows: The photosynthetic quantum efficiency is calculated as:

Model 4
The light dependence of P n is expressed as [29]: Here P n is the chlorophyll a-normalized net photosynthetic rate at irradiance I, α is the initial slope of the P n -I response curve, β and γ are the nondimensional parameters reflecting photoinhibition and light saturation, respectively, and R d is the dark respiration rate.
I sat is calculated as: P nmax is obtained by: When P n = 0, I c is given as follows, The photosynthetic quantum efficiency is calculated as:

Statistical analysis
P n -I data were fitted using SPSS version 24.0 using nonlinear, least-squares fitting based on the Levenberg-Marquardt algorithm. Duncan's post hoc tests (p < 0.05) were performed to establish differences among fitted results from model 1, model 2, model 3 and model 4. Data were reported as the means and standard errors in the calculations. Goodness of fit of the mathematical models to experimental data was assessed using the adjusted coefficient of determination (R 2 ). Akaike information criterion (AIC) is a standard to measure the best-fit of statistical models. When sample size (n) is small compared to the number of parameters (i.e., n/k < 40), the use of a second order, AIC c (= AIC + 2k(k + 1)/(n − k − 1)) is recommended [34]. In this paper, AIC c of each model was calculated because of n/k = 1 for model 1, yet n/k = 0.75 for models 2, 3 and 4.

Comparison of different P-I models of production curves
Applying different values of the fundamental parameters to the model, the differences in the characteristics of production curves among models 2, 3 and 4 were compared, save for model 1, without consideration of lightinhibition at high irradiant intensity. Assuming that the initial slope α was 0.5 (the initial slope of the curve equals the reciprocal of γ in model 3), increasing values of the light-saturated or photoinhibition parameters decreased P nmax of the curve and increased the magnitude of inhibition in three types ( Fig. 1b-f ), which indicated that they could closely reproduce the trend of the P n -I curve.
However, although P s is defined as being associated with P nmax in model 2, the given value of P s was over 30 ~ 125% of P nmax , for which the biological implication is difficult to understand (Fig. 1a). In Fig. 1c, I sat was kept constant value versus the change of β because I sat was barely related to α or γ, according to Eq. 16. And in fact, greater β values were associated with greater bends of the curve, indicating saturation occurred more easily. Thus, Fig. 1c is clearly contradictory to the basis of photosynthetic physiology.    photo-inhibition, by using models 2, 3 and 4, but in addition to those estimated by model 1.

The morphological and growth characteristics of phytoplankter
M. aeruginosa and M. wesenbergii are two different species of Microcystis sp., and despite having nearly identical P n -I curves, there were some differences in the photosynthetic parameters obtained by the different models ( m −2 s −1 , respectively. No significant differences were found between the I sat calculated by model 4 and the measured data (p > 0.05). The photosynthetic parameters obtained by model 1 were far from the measured data for all strains; above all, the I sat were seriously underestimated. The initial slope of the P n -I curve α, namely, the intrinsic quantum yield, the estimated by model 4 was the highest for all strains among the other three models. There were no significant differences in the estimation of I c or R d among each model. Figure 2B shows that the quantum yield calculated by models 2, 3 and 4 for M. aeruginosa, M. wesenbergii, S. obliquus and Chlorococcum sp. decreased as I increased, until it was equal to zero at the I sat point. Subsequently, it became negative as I increased, which also reveals why P n decreased as I increased above I sat . However, the values of P´n obtained by model 1 were always greater than zero with increasing I due to the asymptotic function in this model.

P n -I curve and P´n-I curve of marine phytoplankton
The P n -I curves of I. galbana, D. salina and P. subcordiformis are shown in Fig. 3A, and obvious differences were observed among strains. P n increased gradually with I towards saturation, which was at 800 μmol photons m −2 s −1 for I. galbana. However, for D. salina and P. subcordiformis, P n increased steeply, almost linearly, within low irradiance intensity (below 200 μmol photons m −2 s −1 ), and it decreased rapidly when it reached the maximum value. As was observed for freshwater phytoplankton, all curves stopped above the I sat , excluding those produced by model 1, which indicates the presence of photoinhibition.
Differences were also observed in photosynthetic characteristic parameters calculated by the four types of models (Table 3). Model 1 either overestimated P nmax or underestimated I sat , and these values showed significant differences with their measured values (p < 0.05). The P nmax obtained by models 2, 3 and 4 for I. galbana were 97.45 ± 3.02, 97.55 ± 3.37 and 98.33 ± 3.20 μmol O 2 mg −1 Chl a h −1 , respectively. The I sat corresponding  Chl a h −1 , respectively; however, the calculated I sat were significantly higher than the measured values (p < 0.05), likely because of the rapid increase of P n during low-intensity irradiance. For α, estimated by model 4 was higher for all strains than those estimated by models 1, 2, and 3, with significant differences (p < 0.05) for D. salina and P. subcordiformis. The photosynthetic quantum yield represents the efficiency of carbon dioxide fixation or oxygen evolution by a photosynthetic apparatus driven by absorbed photon energy, that is, the conversion efficiency of absorbed solar energy into chemical energy. Figure 3B indicates that the nonlinear change of P´n as I in three species of marine phytoplankton was similar to that in freshwater phytoplankton.

Discussion
Photosynthesis is not only a biochemical process achieved by photosynthetic apparatuses, it also contains a biophysical process [5,9,35]. As shown in Fig. 4, photosynthetic pigment molecules (Chl), such as chlorophyll a and b and carotenoids, absorb solar energy, which induces them into an excited state (Chl*). The largest amount of exciton binding energy is transferred to   the photochemical reaction centres (e.g., P 680 and P 700 ), where charge separation occurs and produces electrons (e − ) and accompanied by the splitting of water into P 680 *. Other energy is transformed into fluorescence and heat [5,17,23,24,27]. Chl* conducts de-excitation by photochemistry, non-radiation heat dissipation, and chlorophyll fluorescence is then able to accept new photons, yet the process depends on the lifetime of Chl in the excited state [35,36]. The released electrons pass through pheophytin to the primary electron acceptor Q A and are ultimately transferred via a series of electron carriers to photosystem I, thereby generating ATP and reducing power NADPH to driving photosynthetic carbon fixation and respiratory carbon oxidation [5,26]. Although classical P-I models have been widely used to fit the P n -I curve for estimating photosynthetic performance and responses to environment changes for phytoplankton [18][19][20][21][22], many of them were not built based on the photosynthetic mechanism. The exponential model established by Webb et al. [18] and model 1 are still applied extensively for phytoplankton [37][38][39][40][41] even though they lack photoinhibition function. For example, Ma et al. [40] indicated that the P nmax calculated by model 1 for M. aeruginosa FACHB-905 and M. aeruginosa FACHB-469 were 253.92 ± 6.79 and 231.32 ± 6.40 μmol O 2 mg −1 Chl a h −1 , respectively, at 25 °C, yet the corresponding I sat were only 92.71 ± 7.86 and 88.61 ± 3.22 μmol photons m −2 s −1 , respectively. Furthermore, the shape of their P-I curves did not appear to decline above I sat . In present study, the results are as we expected that model 1 overestimated P nmax and could not directly estimate the saturation irradiance (I sat ) because of it being an asymptotic function. The values of α, P nmax and I sat fitted by model 1 showed significant differences with those obtained by other models (p < 0.05); either P nmax or I sat were distinct from their measured data for seven species of phytoplankton (including M. aeruginosa FACHB-905), which suggests that an insufficient irradiance would be supplied to the cultivation if the I sat obtained by model 1 was used as the optimal intensity of irradiance.
To describe the entire range of light levels of phytoplankton, Platt et al. [20,21] proposed another empirical model with a photoinhibition function (model 2 in this study). Superficially, the P n -I curves fitted by model 2 seem to be perfect as other studies [42,43], but the value of P s among the fitted results was notably higher than the value of P nmax in seven phytoplankton species, specially, for I. galbana, M. aeruginosa and S. obliquus when β > 0 (Table 4). To calculate P nmax , Eq. 7 must satisfy the condition of β > 0, otherwise it is unable to estimate directly P nmax because model 2 will degenerate Fig. 4 Schematic representation of the mechanism of photosynthesis consisting of biophysical and biochemical processes into an exponential function without an extreme value when β = 0, where there was no photoinhibition, and the fitting curves were similar to model 1 (Fig. 5). This reveals a clear disagreement with the definition of P s that characterizes the output of dark reactions of photosynthesis in model 2. The analytic solution of I c cannot be directly obtained by Eq. 5. To obtain I c , the Kok effect must be ignored. However, the Kok effect was widespread observed in phytoplankton [33]. Currently, most researchers ignore these problems when they use model 2 to investigate and fit the P n -I curves of phytoplankton [1,[43][44][45]. Therefore, model 2 may be treated carefully in explaining biological implication of P s and estimating photosynthetic parameters.
Compared with previous models, model 3 is no longer just a mathematical equation describing the dependence of the photosynthetic rate on irradiance intensity. Its foundation is an assumption of "photosynthetic factories" (PSF) on physiological mechanisms proposed by Crill [46]. A PSF that is regarded as a combination of photosystem I (PSI) and PSII conducts one unit of light to generate one unit of photosynthetic product. And Eilers and Peeters assumed that the process of photosynthesis is modeled by changes of the states of PSF from the resting state to the activated and inhibited state [17,22]. Model 3 yielded a good-fitting curve for the P n -I data of all phytoplankton species in this study, and the returned values for P nmax , I c , and R d were close to their measured values, except for I sat , which showed a large deviation (p < 0.05). Meanwhile, Fig. 1c shows that the I sat of curve did not change with the value of β. This may be because the model neglects the detail process of the capture of solar energy, energy transfer process, and electron transport from PSII to Cytb6f and then to PSI.
Although differences between higher plants and phytoplankton are observed in photosynthetic antenna system and photosynthetic components [10,16], in present study the P n -I curves of all phytoplankton species fitted by model 4 were good and the returned values were also close to the measured data. This reveals  5 The P n -I curves produced by the model 2 at β = 0 and β > 0 that model 4 for higher plants is applicable for phytoplankton. Acquiring an accurate and optimal parameter for irradiance intensity is essential to achieve high biomass of phytoplankton in production. Irradiance is rapidly attenuated during high-cell density cultivation of phytoplankton [14,27]. Variation in the pigment composition of light harvesting complexes with irradiance intensity has been observed in most species of phytoplankton [4,5]. Irradiance intensity also regulates the accumulation of triacylglycerols and carbohydrates [6,7]. Note that obtained I sat by model 4 was closer to the measured value than other three models. The differences between the returned values for P nmax , I c , and R d by model 4 and their measured values were slightly larger than those by model 3, without significant differences (p > 0.05). For model 3 and model 4 with the same four parameters, they can be regarded as a semi-mechanistic model and mechanistic model in describing P n -I curve, respectively. The former mainly characterizes and simulates P n -I curves of phytoplankton, the latter mainly focuses on photosynthetic performance of photosynthetic organisms. In addition, the biological significance of parameters in model 3 is implicit except for γ. On the contrary, the biological significance of four parameters in model 4 is distinct. Consequently, model 4 may provide an alternative method in study photosynthetic characteristics of phytoplankton. The fitting curves by model 4 for P. subcordiformis, M. aeruginosa, M. wesenbergii, and Chlorococcum sp. exhibited some deviations under low intensity of irradiance, likely because the model targeted higher plants, which showed higher light dependence than phytoplankton.
In meso-and eutrophic water bodies, irradiance or temperature is a key factor affecting changes of phytoplankton community composition, especially for those that become the dominant population between cyanobacteria and green algae [47]. The results of this study explicitly demonstrate that M. aeruginosa and M. wesenbergii had high intrinsic quantum efficiency (α), while their Chl a content per cell was lower than that of both S. obliquus and Chlorococcum sp., indicating the efficient light harvesting and use for M. aeruginosa and M. wesenbergii because the intrinsic quantum yield represents the numbers of photosynthetic electrons required to assimilate one CO 2 molecule [8]. In addition, almost two times less α than both S. obliquus and Chlorococcum sp., and the largest P nmax were found in M. aeruginosa. However, M. aeruginosa is the main contributor of notorious bloom-forming cyanobacteria in global freshwater bodies, such as Dianchi Lake in China [48]. These results reveal the underlying physiological basis of photosynthesis of Microcystis with lower "critical light intensity", and provide important insights into the management and control of cyanobacteria in changing lakes and estuarine waters.
I. galbana and D. salina are applied world-wide to generate biofuels due to their rich lipids (lipid levels between 23 and 55% by weight of dry biomass), and they are also commonly cultivated with P. subcordiformis (lipid levels between 20% and 30% by weight of dry biomass) for aquaculture in China, Japan, Australia, and southeast Asia [14,49]. To meet nutritional requirements, mixed cultures of two or more species of phytoplankton are often fed to larvae in seed farming of aquatic products [50]. It is critical that the photosynthetic productivity of each strain reach as high as possible during production. The comparison revealed that, although the P nmax lay between I. galbana and P. subcordiformis, other photosynthetic characteristic parameters showed great differences. The smallest α and highest I c were found in I. galbana, which meant a low efficiency of light capture and use for I. galbana. In contrast, the largest α and lowest I sat and I c were in P. subcordiformis, although it possesses lower Chl a content per cell compared to D. salina. Consequently, the ranking of light-dependence in descending order was P. subcordiformis, D. salina, and I. galbana. Under co-culture conditions, a gradient of irradiance from low to mid to high can be supplied in one photoperiod.

Conclusions
Our study showed that significant differences were found between the returned values to photosynthetic characteristics by models 1, 2 and 3, some parameters (e.g., I sat ) were distinctly different to the measured data. Model 4 for higher plants reproduced the irradiance response trends of photosynthesis well, was applicable for phytoplankton, but more studies are required to investigate its flexibility and reusability. Differences in photosynthetic performance were observed among phytoplankton species. P. subcordiformis showed higher light-dependence than D. salina and I. galbana, while M. aeruginosa and M. wesenbergii exhibited more efficient light use than S. obliquus and Chlorococcum sp. These findings could contribute to a better understanding of structure changes of phytoplankton communities in the aquatic ecosystem, especially in those eutrophic lakes and estuaries.