Skip to main content

Dynamics of Sundarban estuarine ecosystem: eutrophication induced threat to mangroves



Sundarbans is the largest chunk of mangrove forest and only tiger mangrove land in the world. Compared to the rich species diversity and uniqueness, very few studies have so far been conducted here, mainly due to its inaccessibility. This study explores water quality, density of biomass, species diversity, phytoplankton abundance and bacterial population of a tidal creek in Sunderban estuary during the post and pre monsoon period of 2008-09.


Phytoplankton community was observed to be dominated by diatoms (Biacillariophyceae) followed by Pyrrophyceae (Dinoflagellates) and Chlorophyceae. A total of 46 taxa belonging to 6 groups were recorded. Other algal groups were Cyanophyceae, Euglenophyceae and Chrysophyceae. Species diversity was highest in summer (March) and lowest in winter season (November) in all the sample stations indicating its close correlation with ambient temperature. Species evenness was fairly high in all five stations throughout the study period. Present study indicated that dissolved oxygen, nutrients and turbidity are the limiting factors for the phytoplankton biomass. The estuary was in eutrophic condition (Chlorophyll-a ≥10 μg/L) in winter. During the month of May phytoplankton biomass declined and at high salinity level (21.2PSU) new phytoplankton species take over, which are definitely better resilient to the high saline environment. Bio-indicator species like Polykrikos schwartzil, Dinophysis norvegica and Prorocentrum concavum points to moderately polluted water quality of the estuary.


Eutrophication as well as presence of toxic Dinoflagellates and Cyanophyceae in the tidal creek of Sundarban estuary definitely revealed the deteriorated status of the water quality. The structure and function of the mangrove food web is unique, driven by both marine and terrestrial components. But little attention has been paid so far to the adaptive responses of mangrove biota to the various disturbances, and now our work unfolds the fact that marine status of Sundarban estuary is highly threatened which in turn will affect the ecology of the mangrove. This study indicates that ecosystem dynamics of the world heritage site Sundarban may facilitate bioinvasion putting a question mark on the sustainability of mangroves.


Sundarban is the single largest chunk of mangrove forest in the world. Formed at estuarine phase of the Ganges - Brahmaputra river system, the Sundarban ecosystem is unique in many respects. The area experiences a subtropical monsoon climate with the annual rainfall of about 1600-1800 mm and several cyclonic storms. This mangrove ecosystem of Indian subcontinent is well known not only for the aerial extent, but also for the species diversity. The biodiversity of Sundarban includes numerous species of phytoplankton, zooplankton, micro-organisms, benthic invertebrates, mollusks, amphibians and mammals [1]. It is the only mangrove tiger land on the earth. It has been declared as a world heritage site by International Union for Conservation of Nature (1987). However, the landscape of the Sundarbans have changed remarkably due to neo-tectonic movement compounded with large scale human intervention from the beginning of last centaury, as a result several species have become extinct or are in very much threatened or degraded state [1, 2]. But any systematic approach towards studying the ecosystem dynamics of Sundarban has not been attempted so far [1, 3].

Our work attempts to explore and understand the correlation between different components of Sundarban ecosystem. The study area is located in the reserved mangrove forest of Sundarbans within the 24-Parganas Forest Division. Herobhanga Forest Block, the northern-most block out of seven forest blocks covers over an area of about 200 km2. Out of this, about 16 square kilometers of mangroves were lost due to encroachment and human intervention over a period of about last fifteen years. Total mangrove coverage in this block was 73.05 km2 as per Survey of India topographic sheet (Surveyed in 1988-89) which was dwindled to 57.53 km2 as deciphered through land use land cover study using remote sensing data of the year 2006 (IRS P-6, LISS-IV data). The entire loss of mangrove coverage was recorded on the northern side of the Bara Herobhanga Khal adjacent to the inhabited Jharkhali Island. In fact, out of nine Forest compartments present in Herobhanga Forest Block, compartment numbers 1 (one), 2 (two) and 3(three) have been completely reclaimed and converted into either aquaculture ponds or agricultural land. Effluents from these aquaculture ponds are disposed through another small creek into this Bara Hero Bhanga Khal (creek), thus acting as a point source contaminant into this mangrove ecosystem; while rain washings from the agriculture fields mix up with this Bara Hero Bhanga Khal as a non-point source of impure water. At the same time, this particular creek separates inhabited Jharkhali Island from the Herobhanga Forest, which is a dense mangrove forest having hardly any human intervention.

This creek also joins two mighty tidal rivers of Sundarbans, namely Matla River in the west and Bidya River in the east and plays an important agent for hydrodynamic set up of this area. Thus, this creek along which sampling and analysis of water was carried out in pre and post-monsoon time for one-year period represents both pristine environment along its south bank (the border with dense mangrove forest) and human interfered environment along its north bank (the border with inhabited Jharkhali Island). During different seasons of the year water quality was studied and the concentrations of the nutrients like ammonia, nitrite, silicate and phosphate was quantitated. Primary producers in this estuarine ecosystem was characterized and the dominant species was identified. Phytoplankton is good indicator of trophic states and many species of this community are sensitive to environmental changes. Their presence or absence from the community indicates changes in physio-chemical environment of the estuary [4]. Seasonal distribution patterns of phytoplanktons and primary producers were investigated thoroughly accompanied with the bacterial abundance in the estuarine water. It was envisaged that analysis of different physico-chemical as well as microbial parameters of water samples along different points on both sides of this river might throw some light on the effect of human intervention on the Sundarban eco-system to make the study significant.


Study area

The exact area for this study was in Jharkhali island, a small locality about 130 km from Kolkata (Figure 1). The area is located in Survey of India topographic sheet no. 79 B/12. A small creek (about 150 meter wide) known as 'Bara Herobhanga Khal' separates the Jharkhali Island from Herobhanga Reserve Forest (which is within Project Tiger area). The 'Bara Herobhanga Forest' joins two mighty river of Sundarbans namely the Matla River on the West and the Bidya River in the East. Matla is connected to Bidya and ultimately flows to the Bay of Bengal. The fresh water connection and discharge to this river has been lost in recent times. Salinity of the river water is relatively high owing to fresh water cutoff from upstream region. There is no previous physico-chemical or biological investigation of this tidal creek (River Matla, Herobhangakhal and adjoining places).

Figure 1
figure 1

Map of the Jharkhali tidal creek with marked sampling points.

Description of study sites

Baro Herobhanga Khal is a creek which joins Matla River in the west and the Bidya River in the east. Both these rivers have no upland freshwater discharge now and have been transformed into mere tidal creeks. All these five sample points are distributed along the junctions of these rivers with the Bara Herobhanga Khal and along the stretch of the Bara Herobhanga Khal itself.

Sample Point No. 1: It is on the Matla River itself. A small creek from the Forest Jetty joins Matla at this point.

Sample Point No. 2: It is just at the junction of the Matla River with Bara Herobhanga Khal near the margin of the river. It is about 3.180 km north of 1st Sample point along the river.

Sample Point No. 3: It is near the concrete jetty at Jharkhali village within the Bara Herobhanga River near the northern margin. It is about 2.170 km towards east of sample point no. 2.

Sample Point No. 4: It is located within the small creek going inside mangrove forest on the southern side of the Bara Herobhanga River. It is about 1.2 km south-east of point no.3.

Sample Point no. 5: It is near the junction point of Bara Herobhanga Khal with the Bidya River, another mighty tidal river just like Matla River, running approximately in the North-South Direction.

This Bara Herobhanga River is extremely significant since both Matla and Bidya Rivers are mighty rivers running along Nort-South direction in Sundarban. The tide comes simultaneously through both the rivers. This Bara Herobhanga River joins both the rivers in an East-West direction. It is of interest to investigate which river system dominates over the other through such connection. This study may throw some light on this direction. However, it needs further investigation, especially hydrodynamic studies which may indicate what happens in such cases. It is to be noted here that Sundarban is a network of rivers, where such comparatively smaller creeks join the larger rivers forming a deltaic environment.

Sample collection

Five stations were set up in the tidal creek to capture the overall diversity of the estuarine ecosystem in and around Jharkhali island. The details of the five stations are given in Table 1.

Table 1 Sample points.

Field trips were conducted at fortnight intervals to collect samples at the points that were either measured on board of the launch or were brought back to laboratory for analysis. This study period extended from November to February (post monsoon study) and March to May (pre monsoon study). Samples were collected from the water surface (0.5 m depth) of all the five stations along the Matla and Bidya river. The water was filtered through 25 μm Nitex mesh at 10-15 mm of Hg vacuum to remove larger organisms and debris. Samples were also collected on deeper cast and filtered appropriately (Membrane filters 0.45 μm) using a Millipore suction apparatus for the study of water quality parameters. All samples were preserved in cold condition, transported to the laboratory within three hours of collection to analyze immediately.

Physico-chemical analysis

Water temperature, pH and conductivity were measured in situ with Hach Portable Meters (HQ40d). Turbidity was measured by using portable turbidity meter (Hach 2100P), salinity was determined in practical salinity units by Knudsen method [5], dissolved oxygen concentration was studied according to the method of Winkler [6], nutrients like inorganic nitrogen (ammonia, nitrite), soluble phosphate and reactive silicate were measured according to the same methodology [6].

Biological analysis

Phytoplankton biomass (Chlorophyll-a)

Chlorophyll samples were drawn from all stations, with a maximum vertical spacing of 10 m through the chlorophyll maximum layer, at least one sample was always taken within 5 m of the maximum concentration. Chlorophyll samples were filtered through Whatman GF/F (0.45 μ) filters and extracted in acetone in dark and refrigerated condition. Chlorophyll-a was determined spectrofluorimetrically [7].

Total count of Phytoplankton and Bacteria

Fluorescence microscope was used to estimate the total number of phytoplankton and bacteria. Immediately after sampling, 50 ml of seawater was preserved with 25% gluteraldehyde (0.2-μm-prefiltered) and stored in cold dark environment to prevent reduction of counts. Cells of phytoplankton and bacteria were collected onto a 25-mm black polycarbonate Nucleopore membrane with a 0.45 μm and 0.2 μm pore size respectively and stained with Acridine orange. Slides were preserved at -20°C until they were counted. Twenty random fields were counted in a Zeiss confocal fluorescence microscope coupled with an image analysis system [8, 9]. Direct estimation of phytoplankton cell count was also performed using Sedgwick-Rafter counting chamber [10]. Viable count of Bacterial colonies was also performed using Luria-Bertani medium by serial dilution method [11].

Community structure analysis

Direct estimation of phytoplankton cell abundance and diversity was performed by cell counting method. Surface phytoplankton was collected and the Lugol's preserved subsamples (1-2 liter) were used for quantitative enumeration utilizing a Sedgwick-Rafter counting chamber and Zeiss research microscope according to UNESCO PROTOCOL [10]. Three indices were used to obtain the estimate of the species diversity (H1), species richness (d) and species evenness (J).

Shannon & Wiener [12] diversity index value was obtained using the following equation:

H 1 = I = 1 p i log e p i

Where H1 = Shannon & Wiener diversity index

Pi = Proportion of sample made up by the ith species

S = Total number of species

Species richness (d) was obtained using the equation

d = S 1 / ln N

Where d = Margalef's diversity index [13]

S = Total number of species

N = No of individuals

Species evenness was determined by using the expression of Pielou [14]

J = H 1 / ln S

Where H1 = Shannon and Wiener index

J = Evenness

S = Total number of species

Identification of phytoplankton

Surface water samples were collected using plankton net (20 μm) and immediately fixed with Lugol,s solution and buffered formaldehyde. The preserved samples were kept in refrigerator until analysis. Before identification water samples were allowed to settle for 24 hours and the supernatant was decanted until a concentrate 10 ml was achieved. Few drops of concentrated sample were taken in a common glass slide and observed under Zeiss binocular microscope equipped with phase contrast optics and photographed with Cannon A 1000 camera. In most of the time phase contrast optics was used because it revealed especially well lightly silicified cells of diatoms. For identification of common diatoms examination of raw (without acid cleaned) material in a water mount was done as stated earlier. But for identification of diatoms with specialized structure like striation (Navicula sp) or raphe (Pseudonitzchia sp) acid clearing of samples was adopted which was used to separate diatoms frustules into single valves on which structure diatoms were best seen. Several keys and illustration were consulted to confirm identification [1522].

Statistical analysis

The results were expressed as differences between the groups considered significant at p < 0.05. Data comparison and influence of the environmental factors on phytoplankton were evaluated by stepwise multiple regression [23] Different statistical analysis and correlation regression analysis were performed using the software STATISTICA.


Physical and Chemical analysis

Temperature of the surface water varied continuously through post monsoon to pre monsoon (Nov'08-May'09) (Figure 2a). The lowest surface water temperature 21°C was recorded in January (2009) while the highest 33°C was recorded in May (2009). pH of the water sample was weakly alkaline and more or less constant in the range of 8.0-8.15 throughout the study period. Salinity level gradually increased from post monsoon to pre monsoon period in the range 10.6-24.6 PSU, lowest value being recorded in November and the highest in May. A steady increase in salinity from station 1 to station 5 was observed in all the sampling months. Average salinity of winter months (Nov-Feb) was 16.17 PSU, lower compared to summer months (March'09-May'89) (23.5PSU) (Figure 2b).

Figure 2
figure 2

Seasonal variation in (a) temperature (°C), (b) salinity (PSU), (c) dissolved oxygen (mg L-1), (d) turbidity (NTU) and (e) nutrients (μg L-1) of estuary. Points represent mean of twenty samples in every month.

Moderate to high dissolved oxygen concentration was observed throughout the sampling months with its maximum in February. Dissolved oxygen concentration was observed to be 8.5 mg L-1 in average in winter and 7.5 mg L-1 in summer (Figure 2c).

Turbidity (NTU) is used to describe water clarity and the values were observed in the range of 6.05 NTU to 21.05 NTU (Figure 2d) during post monsoon to pre monsoon.

Nutrient level i.e. Nitrite-Nitrogen, ammonia-Nitrogen, phosphate and silicate showed higher concentration in winter months compared to summer months. Concentration of Phosphate and Ammonia-Nitrogen were always lower than Silicate and Nitrite. Nitrite-Nitrogen concentration of samples ranged from 7.05 μg/L in summer to 18.25 μg L-1 in winter. Ammonia-Nitrogen concentration was estimated to be 6.17 μg L-1 (average) in winter and 0.94 μg L-1 (average) in summer. Phosphate concentration of water sample ranged from 3.44 μg L-1 (in winter) to 1.41 μg L-1 (in summer). Silicate concentration in surface water sample ranged from 8.17 μg L-1 to 2.86 μg L-1 in winter and summer respectively (Figure 2e).

Biological Analysis


Chlorophyll-a concentration is an index of phytoplankton biomass and the concentration was estimated to be 5.9 μg L-1 to 43.80 μg L-1 through premonsoon to post monsoon months (Nov'08-May'09). Chlorophyll-a concentration was observed to be 33.15 μg L-1 (average) in summer and 19.86 μg L-1 (average) in winter. Chlorophyll-a concentration reached its maximum in February'09 in all the five sampling stations (Figure 3) and similar trend was observed for phytoplankton cell count. Phytoplankton cell count ranged from 1.80 × 104 cells L-1 to 2.05 × 107 cells L-1 with an average of 2.52 × 106 cells L-1 in winter and 9.93 × 105 cells L-1 in summer. Fluorescence count of phytoplankton cell ranged from 2.17 × 106 to 1.42 × 108 in winter and 2.2 × 105 to 4.9 × 106 in summer (Figure 4a, Table 2). The conventional counts for viable bacteria obtained on Luria-Bertani agar medium ranged from 3.68 × 106 CFU L-1 to 2.64 × 107 CFU L-1 during winter season and between 4.48 × 107 CFU L-1and 8.9 × 108 CFU L-1 during summer. Direct counts of bacterial cells ranged from 8.54 × 107 cells L-1 to 9.5 × 108 cells L-1cells in winter and between 1.0 × 109 cells L-1and 4.52 × 1010 cells L-1 in summer (Figure 4b, Table 3).

Table 2 Abundance of phytoplankton in Sundarban estuary.
Table 3 Abundance of bacteria in Sundarban estuary.
Figure 3
figure 3

Seasonal variation of chlorophyll-a concentration (μg L-1) of five sampling points throughout the study period (Nov'08-May'09).

Figure 4
figure 4

Acridine Orange stained cells in microscope fields on black Nucleopore filters (a) Phytoplankton cells and (b) Bacterial cells.

Species composition

Six major classes were recorded in this study. Bacillariophyceae (Diatoms), Chlorophyceae (Green algae), Cyanophyceae (Blue green algae), Pyrrophyceae (Dinoflagellates) and Chrysophyceae (Figures 5 & 6). The phytoplankton community was dominated by Diatoms. Out of 46 phytoplankton taxa identified, 27 genera (46 species) were Diatoms. Centrales were represented by 13 taxa containing 27 species and Pennales by 14 taxa containing 19 speceis. Centric Diatoms predominated in winter months while Pennates in summer (Figure 7). Most abundant pennate Diatoms were Pinnularia, Navicula, Gyrosigma, Thalassionema and Climacosphenia and most abundant centric Diatoms were Chaetoceros, Coscinodiscus, Cyclotella, Bacteriastrum, Actinocyclus and Planktoniella. Pyrrophyceae (Dinoflagellates) made up of 6 taxa containing 9 speceis. Among Pyrrophyceae ((Dinoflagellates) 6 genera were Prorocentrum, Protoperidinium, Peridinium, Ceratium, Dinophysis and Polykriskos. Protoperidinium, Dinophysis and Ceratium were the most abundant genera in Pyrrophyceae (Table 4). The abundance of Dinoflagellates (Table 5) were slightly higher in pre monsoon (99 × 103) than post monsoon (59 × 103). The green algal forms (Chlorophyceae) comprised of 6 genera namely Cosmarium, Closterium, Netrium, Chorella, Dunaliella, and Drapernaldia. Blue-green algae (cyanophyceae) comprised of 5 genera namely Anabaena, Stigonema, Oscillatoria, Gleocapsa and Trichodesmium. Euglenoids were represented by only one taxa Euglena and Crysophyceae by two species of single genus Dictyocha. (Table 4). In the sampling stations 1 and 2 species composition was similar. Interestingly, in stations 3 and 4 species composition were also similar. This might be due to proximity of stations 1 & 2 and 3 & 4. The stations 3 and 4 received nutrient rich discharges from aquaculture ponds. This area was also subjected to anthropogenic influence due to tourism. As a result the pollution level in water in station 3 and 4 were higher. This was reflected in the presence of higher frequency of Dinoflagellates there. The station 5 showed different species compositions with higher abundance of salinity tolerance species, Chlorophyceae and Cyanophyceae. Station 5 was located towards the sea and had higher salinity level, which influenced the species composition here (Table 6). In general higher planktonic biomass was recorded in winter months than summer. Occurrence of Prorocentrum concavum from the estuary deserves special mention here because its occurrence was reported in sub-tropical mangrove habitats [24] outside India so far our knowledge concern. This is the first reporting when Prorocentrum concavum has been identified in Sundarban, the main mangrove in India.

Table 4 Phytoplankton composition
Table 5 Abundance of Dinoflagellates
Table 6 Abundance and composition of phytoplankton taxa.
Figure 5
figure 5

Percentage of algal division in Sundarban estuary in the total study period (Nov'08-May'09).

Figure 6
figure 6

Relative abundance (cells L-1) of planktonic classes over the study period (Nov'08-May'09).

Figure 7
figure 7

Representative phytoplankton taxa identified in Sundarban estuary. a) Bacteriastrum sp b) Coscinodiscus sp c) Chaetoceros curvesetus d) Bacteriastrum hyalinum e) Closterium f) Asteronella g) Grammatophora marina h) Nitzschia sp i) Navicula penata j) Gyrosigma baticum k) Anabaena l) Asterionellopsis gracilis m) Dinophysis norvegica n) Polykrikos schwartzil.

Community structure

Phytoplankton cell count (cell density) was higher in winter months than summer. Lowest count was observed in May, 2009 and highest in Feb, 2009 in all the five sampling locations. Species richness (d) and species diversity (H1) showed more or less higher value (>0.8) throughout the study period. Species richness was the highest (4.36) in March (Station 3). Diversity index showed higher value in Station 3 and 5 in the month of March. Species evenness (J) i.e., Species equitability recorded highest value (3.25) in Station 1 in the month of April. Species evenness was found to be higher (>1) throughout the study period (Figure 8).

Figure 8
figure 8

Seasonal variation of diversity index (H), evenness (J) and richness (d) of the estuary. Points represent mean of twenty samples in every month.


Chlorophyll-a as biomass index

Phytoplankton is a good indicator of trophic states; each and every change in environment affects this community. Many species of this community are very sensitive to changes and also they respond very quickly. We attempted to assess and predict the trophic state of the tidal creek on the basis of phytoplankton data. The phytoplankton biomass in water of Sundarban estuary was measured by chlorophyll-a concentration when typical seasonal growth pattern was observed in all stations. The phytoplankton biomass increased steadily from November 2008 to February 2009 and thereafter dropped sharply and ultimately reached at a minimum value in May 2009 (Figure 3). The phytoplankton cell density evidenced exactly similar pattern as phytoplankton biomass.

Similar community structure was also reported in Adriatic Sea [25] and in Boka Kotorska Bay [26] in Europe. Phytoplankton cycle reached its maximum in winter and slowed down gradually in summer season, which is perfectly similar to the universal pattern.

The factors that regulate biomass of planktons include nutrients like Nitrogen [27] Phosphorus [28] and Silica [29] but simultaneously high nutrient concentration also enhances the risk of Eutrophication.

Dynamic relationship was noticed between level of these nutrients and Chlorophyll a concentration (Figure 9), correlation coefficients between these nutrient levels and chlorophyll a was determined. (Table 7). Coefficients with p-value < 0.2 had a significant relationship with chlorophyll-a [23]. Correlation coefficient determination yielded value of p much less than 0.2 for all the parameters and showed positive correlations signifying that these nutrients levels of nitrite + ammonia-nitrogen (N), phosphate (P) and silicate (Si) regulated phytoplankton biomass production in the estuary (Figure 9).

Table 7 Correlation coefficients
Figure 9
figure 9

Seasonal variation in chlorophyll-a concentration (μg L-1) with nutrients throughout the study period (Nov'08-May'09). Points represent mean of twenty samples in every month.

A regression equation was explored to describe the interrelations amongst the parameters like nutrients, dissolved oxygen and turbidity and their combined effect on chlorophyll-a level. The pH and salinity level of water were not included in this equation as these variables showed negative correlation with chlorophyll level and p values were greater than 0.2. There were 150 observations in the data sets that include chlorophyll-a for possible use in the equation.

A sensitivity test was performed to identify the most dominant parameters. According to their importance in the equation the parameters were organized in descending order: phosphate, DO, turbidity, nitrite + ammonia-nitrogen, silicate.

[Chl a] = 10.34 + 5.31 [P] + 2.28 [DO]

-1.10 [Tur] - 0.18 [N] - 0.02 [Si]

R 2 = 0.81

Where, Chl a = concentration of Chlorophyll-a

DO = concentration of dissolved oxygen (mg L-1)

N = concentration of nitrite + ammonia-nitrogen (μg L-1)

P = concentration of total phosphorous (μg L-1)

Si = concentration of silicate (μg L-1)

Tur = turbidity (NTU)

The coefficient of determination (R2 = 0.81) was relatively high and the relationship described by this equation was highly significant (p < 0.0001). Actual versus predicted values for Chl a concentration using the data set in conjunction with the equation developed in this study provided a good prediction (Figure 10).

Figure 10
figure 10

Plot of observed versus predicted values of chlorophyll a from predicted equation.

Biomass production gradually increased starting from November 2008 to February 2009 despite the fact that salinity of water also increased during the period; indicating that the salinity regime was within the tolerance level of planktonic biomass. During the month of March 2009 this pattern was discontinued and at salinity level 21.2 PSU the phytoplankton biomass declined in all stations indicating the tolerance level for majority of phytoplankton species. At this higher salinity level, a new set of euryhaline phytoplankton species (e.g. Dunalielle salina, Trichodesmium, Chlorella salina etc.) was observed which were resilient to higher salinity level [3033]. The bacterial population (Table 3) was the lowest in January'09 (3.68 × 106 CFU L-1) and the highest in May'09 (8.9 × 108 CFU L-1) showing exponential relation with temperature (p = 0.006, r = 0.896) (Figure 11). The bacterial cell count obtained by fluorescence microscopy was higher than that of conventional method (plate count), as both culturable as well as nonculturable bacterial population were observed in fluorescence microscopy.

Figure 11
figure 11

Variation of bacterial abundance (CFU ml-1) with temperature (°C). Points represent mean of twenty samples in each month.

Eutrophication of Estuary

The balance of water ecosystem is disturbed by eutrophication i.e. excessive fertilization, which, in turn, leads to increases in phytoplankton quantity and primary production. Eutrophication can have significant deleterious effects on the beneficial uses of estuarine and marine waters. Eutrophication also initiates changes in phytoplankton community structure, decrease in diversity and frequency of harmful algal blooms. Metrices based on phytoplankton quantity and productivity are widely used indicators of eutrophication in the status assessment of surface waters [34, 35]. There are a number of ways in which eutrophication of estuary manifest itself: increase in phytoplankton biomass [36] and macroalgae [37] anoxia and hypoxia [38, 39] even fish and benthos kill [40]. The most commonly used indicator of eutrophication in waterbody, however, is chlorophyll a [41].

During winter months (November 2008 - February 2009), the chlorophyll a concentration remained very high (>10 μg L-1) in all the five stations. This indicates that the estuary was in eutrophic condition during that time [42]. High nutrient input in coastal estuary water may be responsible for successful establishment and enhancement of nuisance algal species like Cyanophyceae and Dinophyceae [43]. The chlorophyll a level dropped rapidly at the onset of summer and reached a comparatively lower value (8.66 μg L-1) in May. Thus the estuary was mesotrophic-eutrophic (<10 μg L-1) in summer.

The poor water quality of the estuary could be ascertained from the presence of toxic Dinoflagellates like Dinophysis [44], Polykrikos [45], Prorocentrum [46], Cyanophyceae like Anabaena, Oscillatoria [47] and diatom like Nitzschia [47]. Many studies on mangrove sediments of Sundarban also indicated the occurrence of pollutant chemicals in Sundarban region [4851].

Eutrophication seems to be a global problem. Nutrient-enrichment of the coastal zone increases the mortality of mangroves by enhancing shoot growth relative to root which makes them vulnerable to environmental stresses like salinity, drought that adversely affect plant water relationships [52]. The Eutrophication of this tidal creek may have detrimental effect on the mangrove vegetation.

Community Structure and biodiversity

The Sundarbans estuarine phytoplankton community is rich in species diversity and species richness [53]. In all the five stations, throughout the study period, species diversity and species richness showed a value >0.8 indicating healthy phytoplankton community in the estuary. However, the diversity index increased steadily from November 2008 to March 2009 and ultimately reached its peak (2.5) in March and then declined slightly. The species richness index also followed the same pattern but attained the maximum in the month of April, 2009.

A high value of diversity index generally implies healthy ecosystem while a low value indicates degraded state. Investigation was made to find out a relationship between the trophic state of ecosystem and species diversity. It was observed that during post monsoon (November 2008 - February 2009) when species diversity increased steadily, the phytoplankton biomass (chlorophyll-a > 10 μg L-1) also increased in all the stations and the estuary was in a eutrophic state [54]. During premonsoon (March-May, 2009), the species diversity index declined and so also the phytoplankton biomass.

However, phytoplankton diversity depends on a number of factors other than nutrient supply. For instance, at a higher trophic level, the impact of predation of fish constitutes a strong top down control on phytoplankton assemblage [55]. The pollution status of water also strongly influences the phytoplankton species diversity. Thus, phytoplankton diversity index should not be accounted as a good indicator of trophic state of an estuary. The diversity indices very well characterize the differences between assemblages and associations, though the factors which influence diversity are seldom governed by trophic state [23].


The dyamics of Sundarban is mainly maintained by sedimentations from the three major rivers Ganga, Bramhaputra and Meghna. Both the spatial and temporal influences have been demonstrated on the salinity in this region. While describing mangrove forest of Sundarbans, general tendency is to describe this forest as 'Pristine' and totally devoid of any human intervention [56]. However, detail regional study using remote sensing technique on a temporal scale coupled with extensive ground truth survey points to the fact that in the fringe areas of these forested islands with inhabited islands, human intervention is not very rare; rather in some cases, the change in land cover condition is so pronounced that it may be termed as to cross the threshold limit. The most glaring examples of such human interventions have been identified in the Herobhanga Forest Block along the creek. That is precisely the reason behind identification of this area as study area so as to estimate the effect of human intervention on the overall mangrove eco-system of Sundarbans.

The estuary remained eutrophic for most time of the year and mesotrophic-eutrophic during the summer months. Overload of nutrients flowing into the estuary resulted in high levels of dissolved Nitrogen, Phosphorus and Silicate and supporting high algal growth. The Correlation coefficient (p 0.2) for all the parameters signify that these components play crucial role to drive the trophic level of the estuary. The coefficients of the regression equation indicate that phosphate, DO, turbidity, nitrogen and silicate control the production of algal biomass (and hence eutrophic state) in decreasing order of influence.

The trend of biomass production of phytoplanktons represents that specific threshold was reached at salinity level around 21 PSU. Not only the biomass production decreased beyond this salinity level, the species assemblage also changed in favour of more salinity tolerant species. Although the bacterioplankton count did not exhibit similar profile. Phytoplankton abundance was in the range of 1.80 × 104 cells L-1 - 2.05 × 107 cells L-1 and bacterial population was in the range of 8.54 × 107cells L-1- 4.52 × 1010 cells L-1 in the estuary throughout the study period. Phytoplankton community was observed to be dominated by diatoms (Bacillariophyceae) followed by Pyrrophyceae (Dinoflagellates) and Chlorophyceae A total of 46 taxa belonging to six algal groups were identified from the estuarine water of Sundarban.

The high species diversity and species richness of phytoplanktons in the estuary throughout the study period indicated overall good health of the producers in the ecosystem. However, during pre-monsoon period, the species diversity and species richness showed declining trend that possibly resulted from stress in the abiotic environment.

Mangroves are the only woody halophytes dominated ecosystem situated at the confluence of land and sea, they occupy a harsh environment, being daily subject to tidal changes in temperature, water and salt-exposure and varying degree of anoxia [57]. Eutrophication as well as presence of toxic Dinoflagellates and Cyanophyceae in the tidal creek of Sundarban estuary definitely revealed the deteriorated status of the water quality.

Mangrove communities are recognized as highly productive ecosystems that provide large quantities of organic matter to adjacent coastal waters in the form of detritus and live animals [58]. There is a close microbe-nutrient-plant relationship that functions as a mechanism to recycle and conserve nutrients in the mangrove ecosystem. Nutrient enrichment is one of the most serious threats to near shore coastal ecosystems. Increasing nutrients availability introduces an instability into mangrove forests that lowers their resilience to environmental variability. The instability arises because nutrients, particularly nitrogen, stimulate growth of shoots relative to roots [59], thereby enhancing productivity during favorable periods but increasing vulnerability to water stress during drought. Enhanced instability with coastal eutrophication has far reaching consequences for many aspects of mangrove ecosystem function under contemporary and future climate conditions [52]. Apart from other aspects this study gains its significance for throwing light on future of one of the most precious natural resource of this biogeographic region.


  1. Gopal B, Chauhan M: Biodiversity and its conservation in the Sundarban Mangrove Ecosystem. Aquatic Sciences. 2006, 69: 338-354. 10.1007/s00027-006-0868-8.

    Article  Google Scholar 

  2. Sodhi NS, Brok BW, Bradshaw CJA: Diminishing Habitats in Regions of High Diversity. Trophial conservation Biology. Edited by: Sodhi NS, Brok BW, Bradshaw CJA. 2007, Blackwell Publishing Ltd, 1-32. ISBN 978-1-4051-5073-6

    Google Scholar 

  3. Alongi DM: Life in Tidal Water. The Energetics of Mangrove Forests. Springer Science. Edited by: Alongi DM. 2009, 65-87.

    Google Scholar 

  4. Rissik D: Plankton. A Guide to their Ecology and Monitoring for Water Quality. Edited by: Suthers IM, Rissik D. 2009, CSIRO Publishing, Collingwood, Australia, ISBN: 9780643090583

    Google Scholar 

  5. Knudson M: Hydrographical tables. G.E.C. Gad Copenhagen. 1901, 63-

    Google Scholar 

  6. JGFOS Protocol: Protocol for the Joint Global Flux Study (JGFOS) Core Management. 1994, 178-

    Google Scholar 

  7. Ventrick EL, Hayward TL: Determining chlorophyll on the 1984 CalCOFI surveys. CalCOFI Rep. 1984, 25: 74-79.

    Google Scholar 

  8. Hobbie JE, Daleyr RJ, Jasper S: Use of Nuclepore Filters for Counting Bacteria by Fluorescence Microscopy. Appl Environ Microbiol. 1976, 33: 1225-1228.

    Google Scholar 

  9. Bianchi A, Giuliano L: Enumeration of Viable Bacteria in the Marine Pelagic Environment. Appl Enviro Microbiol. 1995, 62: 174-177.

    Google Scholar 

  10. UNESCO: Monograph and Oceanographic methodology. Phytoplankton Manual. UNESCO Technicalpaper No.28. Edited by: Sournia A. 1978, 337-

    Google Scholar 

  11. Cappuccino JG, Sherman N: Cultivation of Microorganisms. Microbiology: A Laboratory Manual. Edited by: Cappuccino JG, Sherman N. 2007, Dorling Kindersley Publishers, India, ISBN 81-317-1437-3 129-134

    Google Scholar 

  12. Shannon CE, Wiener W: The mathematical theory of communication. 1949, University of Illinois Press. Urbana, Illinois, 125-

    Google Scholar 

  13. Margalef DR: Information theory in ecology. Gen Syst. 1958, 3: 36-71.

    Google Scholar 

  14. Pielou EC: Mathematical ecology. 1977, John Wiley & Sons, New York, 385-

    Google Scholar 

  15. Subramanyan R: A systemic account of the marine plankton diatoms of the Madras coast. Proc Indian Acad Sci. 1946, 24: 85-197.

    Google Scholar 

  16. Hendy NI: An introductory account of smaller algae of British Coastal Water, Part; V. Bacillariophyceae(Diatoms). Min Ag Fish Food Fishery Invest Ser. 1964, 4: 317-

    Google Scholar 

  17. Patrick R, Reimer CW: The diatoms of the United states exclusive of Alaska and Hawaii. Vol 1, Monographs Academy of Natural Sciences, Philadelphia. 1966, 13: 688-

    Google Scholar 

  18. Patrick R, Reimer CW: The diatoms of the United states exclusive of Alaska and Hawaii. Vol. 2, Part 1- Monographs Academy of Natural Sciences, Philadelphia. 1975, 13: 213-

    Google Scholar 

  19. Nwanko DI: Contribution to the diatom flora of Nigeria. I Diatoms of lagos lagoon and the adjacent sea. Nig J Bot. 1990, 3: 53-70.

    Google Scholar 

  20. Tomas CR: Identifying marine phytoplankton. 1997, Academic Press, Harcourt Brace and Company, Torpnto, 858-

    Google Scholar 

  21. Wehr JD, Sheath RG: Freshwater Algae of North America: Ecology and Classification. Edited by: JD Wehr RG. 2003, Academic Press. USA and UK

    Google Scholar 

  22. Perry R: A Guide to Marineplankton of southern California. Mar Sci Center, Ocean Globe. 2003, 1-17. 3

    Google Scholar 

  23. Rakocevic-Nedovic J, Hollert H: Phytoplankton Community and Chlorophyll a as Trophic State Indices of Lake Skadar (Montenegro, Balkan). Environ Sci & Pollut Res. 2005, 12 (3): 146-152.

    Article  CAS  Google Scholar 

  24. Faust MA: Morphologic details of six benthic species of Prorocentrum (Pyrrophyta) from a mangrove island, Twin Cays, Belize, including two new species. J Phycol. 1990, 26 (3): 548-558. 10.1111/j.0022-3646.1990.00548.x.

    Article  Google Scholar 

  25. Nincevic Z, Marasovic I: Chlorophyll a and primary production of size fractionated phytoplankton in the middle Adriatic Sea. Rapp Comm Int mer Médit (CIESM), Dubrovnik. 1998, 35: 472-473.

    Google Scholar 

  26. Krivokapic S: Chlorophyll a As Biomass Indicator In: The Boka Kotorska Bay. Ohrid, Republic of Macedonia. 2008, 27: 1-4.

    Google Scholar 

  27. Kivi K, Kaitala S, Kuosa H, Kuparinen J, Leskinen E, Lignell R, Marcussen B, Tamminen T: Nutrient limitation and grazing control of the Baltic planktoncommunity during annual succession. Limnol Oceanogr. 1993, 38: 893-905. 10.4319/lo.1993.38.5.0893.

    Article  Google Scholar 

  28. Andersson A, Hajdu S, Haecky P, Kuparinen J, Wikner J: Succession and growth limitation of phytoplankton in the Gulf of Bothnia (Baltic Sea). Mar Biol. 1966, 126: 791-801. 10.1007/BF00351346.

    Article  Google Scholar 

  29. Zimba PV: The use of nutrient enrichment bioassays to test for limiting factors affecting epiphytic growth in Lake Okeechobee, Florida: Confirmation of nitrogen and silica limitation. Arch Hydrobiol. 1998, 141: 459-468.

    CAS  Google Scholar 

  30. Simon VA, Geoffrey AC, Geoffrey MG: Salt-stimulation of caesium accumulation in the euryhaline green microalga Chlorella salina: potential relevance to the development of a biological Cs-removal process. J Gen Microbiol. 1993, 139: 2239-2244.

    Article  Google Scholar 

  31. Fei-Xue Fu, Bell PRF: Effect of salinity on growth, pigmentation, N2 fixation and alkaline phosphatase activity of cultured Trichodesmium sp. Mar Ecol Prog Ser. 2003, 257: 69-76. 10.3354/meps257069.

    Article  Google Scholar 

  32. Qinghua He, Dairong Q, Zhang Q, Shunji He, Yin Li, Linhan B, Zhirong Y, Cao Yi: Cloning and expression studies of the Dunaliella salina UDP-glucose dehydrogenase cDNA. DNA Sequence. 2005, 16 (3): 202-206.

    Article  Google Scholar 

  33. Oren A: Diversity of Organic Osmotic Compounds and Osmotic Adaptation in cyanobacteria and algae. Algae and Cyanobacteria in Extreme Environments. Edited by: Seckbach J. 2007, Springer, 639-655. full_text.

    Chapter  Google Scholar 

  34. HELCOM: Environment of the Baltic Sea area 1994-1998. Balt Sea Environ Proc No. 82B. 2002, 214-

    Google Scholar 

  35. EEA: Europe's environment - The fourth assessment. State of the environment report No. 1/2007. 2007, 452-

    Google Scholar 

  36. Harding LW, Perry ES: Longterm increase of phytoplankton biomass in Chesapeake Bay, 1950-1994. Mar Ecol Prog Ser. 1997, 157: 39-52. 10.3354/meps157039.

    Article  Google Scholar 

  37. Valiela I, McClelland J, Hauxwell J, Behr PJ, Hersh D, Foreman K: Macroalgal blooms in shallow estuaries: Controls and ecophysiological and ecosystem consequences. Limnol Oceanogr. 1997, 42: 1105-1118. 10.4319/lo.1997.42.5_part_2.1105.

    Article  Google Scholar 

  38. Rosenberg R, Elmgren R, Fleischer S, Johnsson P, Persson G, Dahlin H: Marine eutrophication studies in Sweden. Ambio. 1990, 19: 102-108.

    Google Scholar 

  39. Kiirikki M, Lehtoranta J, Inkala A, Pitkänen H, Hietanen S, Hall POJ, Tengberg A, Koponen J, Sarkkula J: A simple sediment process description suitable for 3D-ecosystem modelling - Development and testing in the Gulf of Finland. J Mar Syst. 2006, 61: 55-66. 10.1016/j.jmarsys.2006.02.008.

    Article  Google Scholar 

  40. Norkko A, Bonsdorff E: Population responses of coastal zoobenthos to stress induced by drifting algal mats. Mar Ecol Prog Ser. 1996, 140: 141-151. 10.3354/meps140141.

    Article  Google Scholar 

  41. Kauppila P: Phytoplankton quantity as an indicator of eutrophication in Finnish coastal waters. Monogr Boreal Environ Res. 2007, 31: 1-57.

    Google Scholar 

  42. Jones RA, Fred Lee G: Recent Advances in Assessing Impact of Phosphorus load on eutrophication-related water quality. Water Research. 1982, 16: 503-515. 10.1016/0043-1354(82)90069-0.

    Article  CAS  Google Scholar 

  43. Paerl HW: Nuisance phytoplankton blooms in coastal, estuarine and inland waters. Limnol Oceanogra. 1988, 33: 823-847. 10.4319/lo.1988.33.4_part_2.0823.

    Article  CAS  Google Scholar 

  44. France J, Patricija M: Ecological characterization of toxic phytoplankton species (Dinophysis spp., Dinophyceae) in Slovenian mariculture areas (Gulf of Trieste, Adriatic Sea) and the implications for monitoring. Mar Pollut Bull. 2006, 52: 1504-1516. 10.1016/j.marpolbul.2006.05.012.

    Article  CAS  Google Scholar 

  45. Jackson JBC: Ecological extinction and evolution in the Brave New Ocean. In the Light of Evolution Vol II: The National Biodiversity and Extinction. Edited by: Avise JC, Hubbleand SP, Ayala FJ. 2008, The National Academics Press. Washington. D.C

    Google Scholar 

  46. Dickey RW: The extraction, purification, and characterization of toxins from the marine dinoflagellates Gambierdiscus toxicus and Prorocentrum concavum. PhD thesis. 1984, Southern Illinois University, Carbondale, Illinois, USA, 125-

    Google Scholar 

  47. Klisch M, Häder DP: Mycosporine-Like Amino Acids and Marine Toxins - The Common and the Different. Mar Drugs. 2008, 6: 147-163.

    Article  CAS  Google Scholar 

  48. Binelli A, Sarkar SK, Chatterjee M, Riva C, Parolini M, Bhattacharya BD, Bhattacharya AK, Satpathy KK: Concentration of polybrominated diphenyl ethers (PBDEs) in sediment cores of Sundarban mangrove wetland, northeastern part of Bay of Bengal (India). Mar Pollut Bull. 2007, 54: 1220-1229. 10.1016/j.marpolbul.2007.03.021.

    Article  CAS  Google Scholar 

  49. Binelli A, Sarkar SK, Chatterjee M, Riva C, Parolini M, Bhattacharya B, Bhattacharya AK, Satpathy KK: A comparison of sediment quality guidelines for toxicity assessment in the Sunderban wetlands (Bay of Bengal, India). Chemosphere. 2008, 73: 1129-1137. 10.1016/j.chemosphere.2008.07.019.

    Article  CAS  Google Scholar 

  50. Saha M, Cobelo-Garcia A, Sarkar SK, Prego R, Bhattacharya B: Distribution of metals in representative biota of sundarban mangrove wetland, northeast India. Bull Environ Contam Toxicol. 2006, 76: 656-662. 10.1007/s00128-006-0970-6.

    Article  CAS  Google Scholar 

  51. Ghosh A, Dey N, Bera A, Tiwari A, Sathyaniranjan KB, Chakrabarti K, Chattopadhyay D: Culture independent molecular molecular analysis of bacterial communities in the mangrove sediment of Sundarban, India. Saline Systems. 2010, 6: 1-10.1186/1746-1448-6-1.

    Article  Google Scholar 

  52. Catherine EL, Marilyn CB, Katherine CM, Feller IC: Nutrient Enrichment Increases Mortality of Mangroves. PLoS ONE. 2009, 4 (5): e5600.-

    Article  Google Scholar 

  53. Biswas H, Dey M, Ganguly D, De TK, Ghosh S, Jana TK: Comparative Analysis of Phytoplankton Composition and Abundance over a Two- Decade Period at the Land-Ocean Boundary of a Tropical Mangrove Ecosystem. Estuaries and Coasts. 2010, 33: 384-394. 10.1007/s12237-009-9193-5.

    Article  CAS  Google Scholar 

  54. Nixon S, Trent Z, Marcuello C, Lallana C: Europe'swater: An indicator-based assessment. European Environment Agency Topic Report 1/2003. 2003

    Google Scholar 

  55. Komarkova J: Fish stock as a variable modifying trophic pattern of phytoplankton. Hydrobiol. 1998, 369/370: 139-152. 10.1023/A:1017091020475.

    Article  Google Scholar 

  56. Ray S: Comparative study of virgin and reclaimed islands of Sundarban mangrove ecosystem through network analysis. Eecol Model. 2008, 215: 207-216. 10.1016/j.ecolmodel.2008.02.021.

    Article  Google Scholar 

  57. Alongi DM: Present state and future of the world's mangrove forests. Environ Cons. 2002, 29: 331-349.

    Article  Google Scholar 

  58. Holguin G, Vazquez P, Bashan Y: The role of sediment microorganisms in the productivity, conservation, and rehabilitation of mangrove ecosystems: an overview. Biol Fert Soils. 2001, 33: 265-278. 10.1007/s003740000319.

    Article  CAS  Google Scholar 

  59. Naidoo G: Differential effects of nitrogen and phosphorus enrichment on growth of dwarf Avicennia marina mangroves. Aqua Bot. 2009, 90: 184-190. 10.1016/j.aquabot.2008.10.001.

    Article  CAS  Google Scholar 

Download references


We acknowledge the instrumental facility provided by UGC-CAS (phase I) and DST-FIST in Department of Biochemistry, University of Calcutta. We are grateful to Ministry of Earth Science, Government of India for financial support and providing fellowship to S Manna and K Chaudhuri. We thank Dr. B.R. Subramanian, Director, ICMAM-PD, Chennai for his constant support. We express our sincere gratitude to the Director, IESWM, Kolkata for his continuous encouragement and interest in our work. It would not be possible to carry out this work without the help and support of the local people of Sundarban. We express our inability to acknowledge them individually.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Maitree Bhattacharyya.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SM and KC performed all experiments, calculated results, prepared the tables, graphs and diagrams, and composed the draft manuscript in consultation with MB. MB designed the experiments, analyzed and interpreted data and results, modified the manuscript in the final form. SB planned the project, was involved in acquisition of funds, selected the site and field stations according to GPS measurement, prepared the map, supplied remote sensing information, and guided the field study and sample collection. All the authors read and approved the final manuscript.

Suman Manna, Kaberi Chaudhuri contributed equally to this work.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Manna, S., Chaudhuri, K., Bhattacharyya, S. et al. Dynamics of Sundarban estuarine ecosystem: eutrophication induced threat to mangroves. Aquat. Biosyst. 6, 8 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: