Phytoplankton nutrient dynamics and flow cytometry based population study of a eutrophic wetland habitat in eastern India, a Ramsar site

Phytoplankton diversity, their abundance based on flow cytometric (FCM) analysis and seasonal nutrient dynamics were investigated from a waste water fed wetland of Eastern India (88° 24.641′E and 22° 33.115′N). The primary objective of the study was to correlate the seasonal fluctuations in phytoplankton abundance to the environmental variables. Total chlorophyll content and FCM based cell counts were used to characterize and quantify the phytoplankton population. Multivariate statistical methods were employed in predicting the possible relationships between biotic and abiotic variables. Distinct seasonal variations characterized by high abundance during the pre-summer period compared to other seasons were detected. The results indicated that environmental factors like water temperature and nutrients, such as various forms of nitrogen and phosphate, influenced the seasonal phytoplankton accumulation. Cluster analysis and non-metric multidimensional scaling helped analyze the seasonal distribution of phytoplankton based on their composition. The dominant genera among the entire phytoplankton community were Scenedesmus spp. of Chlorophyta, followed by Merismopedia spp. of Cyanoprokaryota. Around 3.7 × 105 phytoplankton mL−1 were recorded during the study period. Due to the very high count of individual species in the community, FCM based counting was applied for determination of Species Diversity Index. The entire population was divided into 13 subpopulations based on the cell sorting method and the seasonal abundance in each sub-population was illustrated.


Introduction
Understanding the various aspects of habitable environments controlling resident community accumulation is a basic research objective in ecology. As phytoplankton are the main organisms responsible for introducing energy into food webs, particular focus has been devoted to understand the factors that control their diversity. [1][2][3] The emerging consensus is that bioticabiotic interactions and ecological dri are the key factors for phytoplankton distribution patterns. 4 Much research has been centered on open ocean ecosystems where it has long been acknowledged that phytoplankton contribute to $50% of global primary production. 1,2 Their vulnerability to alterations in the physico-chemical parameters of their aquatic environment renders them important biological indicators of water quality. 5,6 Phytoplankton are potential mitigators of eutrophication since they assimilate excess nutrients rapidly due to their fast growth processes. 5,[7][8][9][10][11][12][13][14] Elevated nutrient levels, especially nitrogen and phosphorus, created by hydrologic processes including water supply, sewage disposal, sheries, waste water management, recreation, etc., boost phytoplankton populations, leading to extensive blooms. 15 However, freshwater phytoplankton taxa have N and P requirements different from marine ones, and may respond quite differently to altered nutrient composition. 16 Even among freshwater taxa, various phylogenetic groups can respond to nutrient conditions in distinct ways. 17,18 Various methods like microscopic cell counting, chlorophyll estimation, biomass estimation, etc. have been exploited for the quantication and characterization of phytoplankton communities as an index of water quality. During the past few decades, ow cytometry (FCM) has been recognized as a potent tool for the study of phytoplankton ecology, especially for studying spatial and seasonal trends. 10,[19][20][21] Due to the auto-uorescent properties of the phytoplankton, mixed aquatic populations can be discriminated with the help of FCM. 22 Generally, allometric and taxonomic analyses of FCM data contribute to characterization of plankton assemblages. 23,24 Besides open oceans, documentation of standing crop of phytoplankton from different wetlands and their ecological factors have also been carried out throughout the world by various authors in North American Great Plains, 25 southern coastal areas of North America, 26 Eastern England, 14 Eastern Europe, 27 Southern Africa 28 etc. Wetlands are ideal habitats for phytoplankton, which act as nutrient sinks, ood control buffers and breeding grounds for aquatic fauna. 29 Some noteworthy works on the nutrient dynamics study related to phytoplankton productivity from fresh water wetlands are available. 11,13,30 In India wetlands are economically important and are mainly used for sh cultivation; moreover, they have distinct architecture, resulting in extensive purication of waste waters. 31 Pradhan et al., 32 suggested that phytoplankton growth could be an important factor for greater sh production and could also act as a biomonitor for water quality assessment in the wetland ecosystems. The wetlands of eastern India represent one of the world's largest integrated resource recovery practice based on a combination of aquaculture, agriculture and horticulture practices.
The wetland currently under study is a wastewater fed aquaculture pond of East Kolkata Wetlands (EKW) of Eastern Indiaa Ramsar site. Here phytoplankton-nutrient dynamics have direct role in sh production together with natural carbon sequestration. 32 Out of 26 wetlands in India, this wetland is one of the world's largest and oldest integrated resource recovery system based on aquaculture production. 31 There have been some sporadic reports regarding waste water management at the EKW. [31][32][33][34][35] The phytoplankton diversity of the EKW has already been reported by some of the present authors. [36][37][38] In the present investigation an attempt has been taken to explain the seasonal variations in phytoplankton population in response to changes in environmental variables of EKW with special emphasis on Flow Cytometry based cell sorting methods.

Study area
The East Kolkata Wetlands (EKW) situated along the eastern fringes of Kolkata metropolitan city (West Bengal), make up one of the largest assembly of sewage fed shponds in the world. It collects sewage wastes from municipalities, agricultural practices and industries of urban and semi-urban areas. It was declared as "Ramsar site number 1203" and based on its ecological and socio-cultural importance, the government of India declared EKW as "Wetland of international importance" under Ramsar convention in the year 2003. 39 The present study site, called Captain Bhery, was selected within EKW (88 24.641 0 E and 22 33.115 0 N as determined by Garmin GPS map 76 CSx device) (Fig. 1). It spreads over an area of 125 km 2 and is1.2-1.5 m deep. 35

Water analysis
Samplings for phytoplankton population and water analyses were carried out at bi-weekly intervals for two years from Oct' 13 to Sep' 15, covering different seasons, namely a short postmonsoon (October-November), winter (December-January-February), a pre-summer period (March), a prolonged summer (April-May-June) and monsoon (July-August-September). Sample water was collected in a 1 L PVC (polyvinyl chloride) bottle immersed beneath the surface water from four different transects of the water body. The sample bottles were stored in a cool place and brought back to the laboratory for determination of nutrient concentrations within 30 min of sampling.
The water temperature (temp.) and transparency (transp.) were recorded in situ using a Celsius thermometer and a secchi disc respectively. In the laboratory, pH was measured using an electronic pH meter. Different chemical parameters including nitrate (NO 3 À ) (phenol disulphonic acid method), nitrite (NO 2 À ) (diazotization method), dissolved inorganic phosphate (DIP) (ammonium molybdate method), dissolved inorganic silicate (DSi) (molybdosilicate method), ammonium nitrogen (NH 4 + ) (phenate method) and hardness (EDTA titration method) were measured spectrophotometrically following the standard protocols of APHA. 40 Dissolved oxygen (DO) was measured using Winkler iodometric titration method 41 using the formula: where V 1 ¼ total volume of sample taken (125 mL), V 2 ¼ volume taken for titration (100 mL) v ¼ 2 mL (1 mL MnSO 4 + 1 mL alkaline KI), x ¼ volume of sodium thiosulphate consumed in the titration (APHA 1975(APHA , 1998. The gross primary productivity (GPP), net primary productivity (NPP) and community respiration rate (CRR) were determined following light and dark bottle method aer 3 hours incubation. Productivity rates were determined by converting DO to carbon equivalence using photosynthetic quotient of 1.2 and respiration quotient of 1.0. Productivity values were determined from the following formulae:

Flow cytometric cell counting and diversity index study
The sample water was ltered and collected using phytoplankton net (100 mm pore size) to remove unwanted large soil particles, zooplankton or other debris which could block the machine during sample run. The cell count of the natural heterogenous phytoplankton community was carried out with the help of Fluorescence Activated Cell Sorting (FACS) on a BD FACSAria™III cell sorter based on pigment uorescence, using a method standardized for EKW samples by Roy et al. 21 This FACS machine was equipped with a 100 mm nozzle, with the ow rate of the sample being maintained at 6 mL min À1 (lowest lters indicated cell granularity, size and refractive index by measuring scattering of incident radiation at 488 nm. Two FSC lters were used to channel and record events corresponding to individual cells within a size limit. The recorded unicelled population for each sample was produced as a cytogram represented as population 1 (P1) (Fig. 2). 21 The optical lters other than FSC and SSC transmitted autouorescence wavelengths emitted by different major photopigments excited by lasers. For instance, chlorophyll-a on excitation at 488 nm, emits red uorescence at around 685 nm which can be transmitted through PE-Texas Red (Filter-616/23) lter. Again, phycoerythrin on excitation at 496 nm, emits orange uorescence at 560-585 nm transmitted through PE (Filter-585/42) lter. 43,44 Similarly, another important photopigment, phycocyanin, emits blue-green uorescence at around 670 nm on being excited at 650 nm, transmitted through APC (Filter-660/20) or PE-Texas Red (Filter-616/23) lters. 45 Red autouorescence (at 695 nm) of peridinine-chlorophyll-protein complex within photosynthetic apparatus excited by 482 nm radiation was transmitted through Per CP-CY5-5-A (695/40) lters. 46 Each P1 was gated and sorted into smaller entities (P4, P5, P6 and P7) based on two-colour pigment uorescence (Fig. 2). 21 Sorted entities (P4-P7) were subdivided on the basis of two-colour pigment uorescence using a different set of lters, so that 13 distinct and consistent sub-populations (P8, P9, P11-P18, P25-P27) were obtained. These sub-populations were studied for a 24 month period (October 2013-September 2015) to nd out seasonal variation patterns among the phytoplankton assemblages.
The gating of the entire population into subpopulations was followed by cell sorting of those gated populations in 4 way sorting precision. Aer cell sorting each sub-population (P8, P9, P11-P18, P25-P27) was collected on a microplate. The microplate with the sorted samples was identied under light microscope (BD Pathway 855), showing the different phytoplankton taxa obtained. Flow cytometric determination of abundance of each phytoplankton sub-population in each sample was estimated using the formula: , where q is the ow rate (mL min À1 ), t is the duration of analysis, n is the number of events counted by the FACS, and N is the number of cells per milliliter.
The diversity index of the phytoplankton community study was determined using Shannon Wiener's Index (H 0 ), species richness and species evenness (e). 47

Statistical analysis
The statistical analysis was performed in the free and opensource statistical soware R version 3.2 (R Core Team 2014). Datasets included environmental variables (pH, temp., transp., DO, BOD, GPP, NPP, CRR, NO 3 À , NO 2 À , NH 4 + , DIN, DSi, DIP, hardness) against a biotic variable (total chl content). Multivariate analyses such as principal component analysis (PCA) and cluster analysis (CA) were performed to envisage the possible relationships between biotic and abiotic variables and to group these variables based on similarity. For the multivariate analysis, excessively rare species (below 1% abundance) were removed from the original dataset as the inclusion of very rare species weakens correlation analysis. 48 The correlation structure within the set of environmental variables was visualized by plotting the correlations of the variables against PCA axes. Since each principal component (PC) is a linear function of the given variables, the loadings represent the correlations of the variables to PC1 and PC2. For CA, unweighted pair group average (UPGA) method was followed for better understanding of species succession where seasonal assemblages were separated according to the similarity of species composition. However, the level of abundance among the phytoplankton taxa uctuated according to seasonal variation. Thus in order to get rid of biasness, we scaled the data before clustering by subtracting with the mean followed by dividing with standard deviation. The new scaled data was then used to draw the dendogram and other plots. In general, the goal of the analysis is to detect meaningful underlying dimensions to explain observed similarities or dissimilarities between the investigated objects.

Physicochemical parameter analysis
The two-year long study investigated the seasonal variations in physicochemical parameters and phytoplankton productivity in Captain Bhery ( Fig. 3 and 4). The average values of each parameter recorded as two datasets at een days intervals in each month were taken into consideration. The water was found to be alkaline with pH values ranging from 7.6 to 8.6 units (Fig. 3a). Minimum pH values were recorded during postmonsoon seasons, while maximum pH was obtained during summer for both years. The water temperature followed a seasonal cycle with a maximum mid-summer value at 36 C and a minimum in mid-winter value at 15.7 C (Fig. 3b). The water transparency remained almost same with minute uctuations between 0.20 m and 0.447 m throughout the year (Fig. 3c). The dissolved inorganic nitrogen (DIN) concentration of the habitat water (Fig. 4a) (Fig. 4g). Phytoplankton biomass was estimated in terms of chl (chlorophyll) content. The chl content of sampled water ranged from 0.095 to 0.501 mg L À1 (Fig. 5a) with values higher during the pre-summer and lower during the post-monsoon for both years. In addition to the estimation of the phytoplankton abundance, their photosynthetic activity was also determined in terms of DO. The DO level of the sample water varied from 3.24 mg L À1 to 8.14 mg L À1 following similar trend in seasonal variation as that of chl content (Fig. 5b), and showing maximum values in winter. The BOD value ranged from 1.44 to 6.0 mg L À1 (Fig. 5b). As evident from Fig. 5c, the values of GPP, NPP and CRR ranged from 283.61 to 2147.71 mgC m À3 h À1 , from 114.11 to 1471.86 mgC m À3 h À1 and from 26.00 to 976.67 mgC m À3 h À1 respectively. The GPP values were higher than NPP, thereby establishing a positively productive ecosystem (Fig. 5c). Maximum productivity (GPP) was recorded in the month of March 2014 and minimum productivity was recorded in August 2015.
Correlation matrix (Table 1) (Table 1). A signicant positive correlation of DO with GPP (r ¼ 0.514) suggested an increased oxygen concentration with higher photosynthetic activity. However, increased BOD levels were observed with increasing values of chl, DO and GPP. The DO , DIN and DIP were signicantly correlated with temp., suggesting a probable higher nutrient concentration during the warm season.

Phytoplankton community study
The phytoplankton community study from the investigated wetland revealed the presence of diverse groups of unicellular or colonial, free-oating autotrophic microplanktonic taxa belonging to mainly four different algal phyla, namely Cyanoprokaryota, Chlorophyta, Bacillariophyta and Euglenophyta. Among the total 165 taxa recorded (Table 2), the chlorophytic members were found to account for the bulk phytoplankton population with 84 species followed by Cyanoprokayotes (30 species), Euglenophytes (28 species) and Bacillariophytes (23 species) ( Table 2). Taxonomic identication of the above recorded taxa had been previously carried out. [36][37][38] In the EKW a total of 3.7 Â 10 5 phytoplankton per mL was recorded during the entire study period, where pre-summer season was the most productive (9.2 Â 10 4 cells per mL) and post monsoon was the least (6.5 Â 10 4 cells per mL). From the total cell count of individual phytoplankton genera, it was also evident that Scenedesmus sp. from phylum Chlorophyta followed by Merismopedia sp. of Cyanoprokaryotes showed the maximum profusion (Fig. 6).
Seasonal distribution of the phytoplankton community composition was revealed from the FACS study. Populations of cells (P4-P7) sorted on the basis of two-color pigment uorescence intensity at two FACS lters were further sorted based on uorescence intensity at two other lters (Fig. 2). The phytoplankton abundance in terms of cell count was mapped based on these sorted sub-populations (P8, P9, P11-P18, P25-P27) representing mixed assemblages of taxa with similar pigment proles (Table 3), and plotted on cytograms for different seasons (pre-summer, summer, monsoon, post-monsoon and winter) (Fig. 7). Each of these subpopulations was tagged by a specic colour ( Table 3). The change in colour intensity of the subpopulations during different seasons indicated their variations in abundance at different seasons in terms of abundance. Colour intensity of any one sub-population was directly proportional to number of cells in the corresponding assemblage (Fig. 7). In general, comparatively higher cell counts were obtained for all assemblages during the pre-summer, while lower counts occurred in post-monsoon (for P11, P12, P8 and P9), monsoon (for P15-P18) and winter (for P13, P14, P25-P27). Along with the phytoplankton abundance, types of phytoplankton recorded from each population were also studied through microscopic identication ( Table 3). Most of the subpopulations obtained by FACS consisted of mixtures of microplanktonic phyla, although a few (P14, P17, P25 and P26) contained members of single phyla only.
The comparison between the seasonal variations of total cell count and chl content of the recorded planktonic algal phyla evidenced that there were almost similar seasonal uctuations in total chl content and cell count of individual groups (Fig. 8a). A positive relation between total cell count and total chl content was established (Fig. 8b), indicating the phytoplankton's contribution to the chlorophyll concentration of the present habitat water.

Species diversity index
Diversity of phytoplankton population at the present site showed distinct variations on a seasonal basis in terms of different biotic indices (Fig. 9). Aer FACS based cell counting method, the diversity was measured by Shannon-Wiener's Index (H 0 ) 47 which varied from 3.23 to 3.37, suggesting an intermediate diversity range as compared to other aquatic habitats like streams and lakes. Highest values were recorded during the pre-summer (H 0 ¼ 3.37) and least during the postmonsoon (H 0 ¼ 3.23) seasons. Likewise, seasonal variations were evident for species richness as well. Seasonally, species richness was lowest in post-monsoon (25.45) whereas it was highest in the winters (28.29). Species evenness (e), a measure of the contribution of individual taxa to the phytoplankton population showed insignicant seasonal variation (Fig. 9).

Statistical analysis
In an attempt to aid data interpretation, inferential statistics like Principal Component Analysis (PCA) were used to envisage the possible relationships between biotic and abiotic variables and group these variables on the basis of similarity (Fig. 10). Principal components with eigen values greater than 1 (Kaiser Guttman criterion) were taken into consideration. PCA among environment and biotic variables extracted four signicant factors (eigen values greater than 1). PC1 and 2 jointly contributed to 60.4% of the variation within the data. The distribution of the observations along the axes reects their respective correlations with the variables. Other ordination approaches yielded similar relations between variables (ESI Fig. 1 and 2 †). Among the variables, temperature has the largest but most negative loading for PC1 and thus its variability is explained almost totally by PC1, which accounted for 33% of the total variance (Fig. 10a). The variable chl has equivalent but positive loadings for PC1. This conrmed that chlorophyll content was strongly correlated (negatively) to water temperature. Almost  ), total dissolved inorganic nitrogen (DIN) and dissolved inorganic phosphate (DIP) occurred along PC2 that represented 27.4% of the variance. Thus, these components in the loading plot largely corresponded to high nutrient condition, a possible indication towards eutrophic nature of the habitat. DO and BOD showed high factor loadings along PC1 and intermediate loadings along PC2, which not only indicated similar patterns of variance but established the interdependence between them. The hardness of the water appeared to negatively relate with chl as well. Again, length and direction of the arrows of NO 3 À , DIN, DIP and pH suggest their similar pattern of seasonal variation. This was already conrmed from correlation study as well. The scores for the 24 sample months pointed towards the variables driving phytoplankton abundance in each season (Fig. 10b). While the winter and presummer months (Nov-Mar) showed high scores along PC1, summer (Apr-Jun) showed high scores along PC2 and monsoon (Jul-Oct) showed low scores along both axes. Comparing with the loadings of different variables, it was inferred that high phytoplankton abundance in winter and presummer was encouraged by low water temperature and led to high GPP, low nutrient (DIN and DIP) status, high BOD, high DO and low transparency. In summer, high temperature led to waning of phytoplankton bloom and accumulation of nutrients (eutrophication) as indicated by  lower GPP, chl, DO and BOD but higher DIN and DIP. Arrival of monsoon diluted the nutrient concentrations, raised transparency but hindered phytoplankton bloom.
The cluster analysis (CA) (Fig. 11) and NMDS (Fig. 12) ordinations provided a better insight into the seasonal pattern of species based upon their abundance data. The NMDS algorithm ranks distances between objects, and uses these ranks to map    Table 4 List of recorded genera along with their abbreviated names clustered on the basis of cluster analysis

Cluster
Cyanoprokaryota Chlorophyta Bacillariophyta Euglenophyta pre-summer followed by winter, summer, monsoon and the least in post-monsoon. Cluster III included the winter predominating phytoplankton groups. The Table 4 includes the genera grouped into Cluster I, II and III. Thus from the above two plots it was evident that Euglenophytes were abundant mostly in winters followed by presummer while the Cyanoprokaryotes dominated in the winter season. The Chlorophytes' abundance was maximum during the pre-summer. The maximum chlorophyll content was recorded previously during this season. This could in turn indicate that Chlorophytes accounted for maximum chlorophyll content and phytoplankton productivity.

Discussion
Wetlands being dynamic ecosystems continually undergo natural changes due to inltration of sediments and nutrients that signicantly inuence the phytoplankton assemblages. The EKW is renowned for carrying out sewage treatment in conjunction with aquaculture. The water body under investigation, Captain Bhery, had been found to be associated with sewage treatment facility experiencing high ushing rates. 31,34 It was noted that the environmental variables recorded at Captain Bhery nearly covered the value ranges of corresponding parameters obtained at other different EKW sites. [50][51][52] The Captain Bhery was thus considered to be sufficiently stable and dynamic to be a representative water body for the EKW, and the complex interplay of the different hydrological factors and their impact on the phytoplankton population observed at this site could be extended for majority of the EKW. Here the phytoplankton biomass was found to uctuate seasonally and appeared to be highly inuenced by the nutrient regime as well as environmental conditions. The average phytoplankton abundance during the entire study period was found to be 3.7 Â 10 5 cells per mL with maximum abundance during pre-summer and minimum during monsoon. The dominant taxa recorded amongst the entire phytoplankton community were Scenedesmus spp. followed by Merismopedia spp. of Cyanoprokaryota. The FACS study was found to quite helpful in calculating phytoplankton cell count efficiently. Eminent researchers like Crosbie et al., 53 Tijdens et al., 54 Toepel et al., 55 applied FACS to study the phytoplankton community especially in freshwater environment. Cellamare et al., 56 also sorted 175 algal taxa from different freshwater ecosystems including seasonal dynamics of Synechococcus spp using FACS method. The division Cyanoprokaryota and Euglenophyta mainly predominated during the winter seasons. The CA demonstrated progressive change of dominance through a warmer to cooler temperature gradient. The chlorophyll content and the phytoplankton count including FACS observations suggested a seasonal trend in phytoplankton assemblages with monsoon being least abundant due to seasonal precipitation, and presummer being the maximum. Similarities between variation patterns of total chlorophyll content and phytoplankton cell count suggested that the autotrophic productivity of the present ecosystem was primarily regulated by the phytoplankton biomass as indicated from GPP value also. It has already been observed that the chlorophyll values accorded with that of eutrophic ecosystem. 57 The pH recorded was slightly alkaline which, from the above results, was evidently contributed by the different nutrients present therein. Besides, the higher pH values obtained validated the occurrence of eutrophication. 58,59 An important ecological factor regulating phytoplankton growth is water temperature. 32 Temperature seemed to be primarily responsible for the shis in phytoplankton assemblages as signicant negative correlation of temperature with chlorophyll, DO, GPP and NPP was obtained from the study. This inference was supported by PCA results, which indicated that water temperature and annual precipitation were major determinants of phytoplankton abundance, which was highest in winter and lowest in monsoon ( Fig. 8 and 10). Although chlorophyll shows strong negative correlation with temperature, it is not winter (when temperature is minimum that productivity peaks. It is in presummer when the mix of high nutrient levels and optimum temperature causes maximium productivity. Both abundance and diversity of phytoplankton taxa diminished from presummer to summer despite eutrophic conditions in summer ( Fig. 3-5, 8 and 9). Mesotrophic status was achieved in winter and presummer by assimilation of excess nutrients by phytoplankton, which also raised GPP and DO by oxygenic photosynthesis, and in turn encouraged high microbial abundance as indicated by high BOD and CRR.
The DO is essential to all forms of aquatic life, including those organisms responsible for the self-purication processes in natural waters. 59 A regulatory network of DO along with photosynthetic activity and primary productivity therein equilibrates the ecological balance of the ecosystem. The present ndings accorded with the previous reports of Hardy 60 where DO shows a positive correlation with the phytoplankton biomass (Table 1). According to Ganf and Horne, 61 if in a productive aquatic ecosystem respiration accounted for large proportion of GPP, it would be a measure of eutrophic nature. In general, decomposition of the sewage, dead plankton etc. along with respiration of the inhabitants are responsible for creating avenues for high CRR and ultimately BOD of the water column. The present investigation recorded higher BOD levels, which probably emphasized on the eutrophic status of the selected site. A positive correlation obtained between DO and BOD is indicative of higher heterotrophic microbial community along with planktonic autotrophs supported by higher CRR values and lower percentages of oxygen saturation. However, the obtained DO values ranged up to of 8.14 mg L À1 which according to WHO, 60 supported the survival of biological communities including sh production. Higher values of GPP and its signicant positive correlation with DO were useful in supporting the fact that phytoplankton contributes a natural method of biological purication for the sewage treatment in EKW. This was further supported by other studies from other EKW sites such as Dasgupta et al. 62 During sewage treatment the microbial degradation of the sewage releases the nutrients stored in it, creating eutrophic conditions, which in turn support high rates of primary productivity. 32,35 The excretion of nitrogenous compounds by sh is also a source of NO 3 À , NO 2 À , NH 4

+ and other inorganic
This journal is © The Royal Society of Chemistry 2018 substances. 63 Major nutrient like nitrogen occurs in natural waters in various forms, including NO 3 À , NO 2 À and NH 4 + . The NO 3 À is the essential nutrient for many photosynthetic autotrophs and has been identied as the growth limiting nutrient. 59 However, in municipal and industrial waste-waters or effluents including biological treatment plants, NO 3 À concentrations are enhanced resulting in eutrophication. 57 Similarly NO 2 À concentrations higher than 21.74 mM and NH 4 + greater than 11.11 mM could be an indication of organic pollution. 59 High availability of NO 3 À in EKW not only encourages phytoplankton abundance but also raises anaerobic metabolism, particularly under oxic conditions, while NO 2 À plays only a transient role in N cycling. 64 Another major nutrient of aquatic systems for phytoplankton development is phosphate mostly in form of DIP which in general ranges from 0.053 to 0.21 mM in most natural surface waters. 57 The present records also showed the positive role of NO 3 À , NH 4 + and DIP on phytoplankton biomass growth. 65 However, loadings of water parameters, particularly temperature, transparency, chlorophyll content, BOD and DO, obtained by PCA showed that biomass-induced opacity of the water column was mainly caused by photosynthetic microorganisms, indicating a productive ecosystem occurring in winter and presummer ( Fig. 7-10). Since this productivity was seasonal, an effective way of sewage treatment would be periodic harvesting of the nutrient-containing photosynthetic biomass from the littoral zone throughout winter and presummer, leading to an oligotrophic and oxygenated habitat for pisciculture.

Conclusion
In summary, it may be said that the EKW station is suitable for phytoplankton population growth. Elevated nutrient levels, particularly species of nitrogen and to a lesser extent DIP, contribute to the growth of phytoplankton resulting in high levels of DO, GPP and NPP. These in turn favor the purication process, provided that the photosynthetic biomass is removed at the end of the productive season to prevent decay and anoxia. The levels of phytoplankton biomass, pH, BOD, major nutrient concentrations of the water body and the heterotrophic microbial respiration were in support for the wetland under study to be eutrophic. Even though BOD levels are high, signicantly higher DO along with elevated levels of GPP and NPP support the growth of sh making it suitable for aquaculture. FACS based phytoplankton community study was found to be useful to account for their diversity and abundance. Seasonal variation in the phytoplankton standing stock was primarily regulated by environmental variables like nutrient availability, temperature variations and annual precipitation. From the total cell count of individual phytoplankton genus, it could be inferred that Scenedesmus sp. (Chlorophyta) followed by Merismopedia sp.

Conflicts of interest
There are no conicts to declare.