Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.

org/content/8/1/26

AQUATIC BIOSYSTEMS
Open Access

RESEARCH

Physicochemical and biological factors controlling water column metabolism in Sundarbans estuary, India
Kaberi Chaudhuri1, Suman Manna1, Kakoli Sen Sarma1, Pankaj Naskar1, Somenath Bhattacharyya1 and Maitree Bhattacharyya2*

Abstract
Background: Sundarbans is the single largest deltaic mangrove forest in the world, formed at estuarine phase of the Ganges - Brahmaputra river system. Primary productivity of marine and coastal phytoplankton contributes to 15% of global oceanic production. But unfortunately estuarine dynamics of tropical and subtropical estuaries have not yet received proper attention in spite of the fact that they experience considerable anthropogenic interventions and a baseline data is required for any future comparison. This study is an endeavor to this end to estimate the primary productivity (gross and net), community respiration and nitrification rates in different rivers and tidal creeks around Jharkhali island, a part of Sundarbans estuary surrounded by the mangrove forest during a period of three years starting from November08 to October11. Results: Various physical and chemical parameters of water column like pH, temperature, conductivity, dissolved oxygen, turbidity, suspended particulate matter, secchi disc index, tidal fluctuation and tidal current velocity, standing crop and nutrients were measured along with water column productivity. Relationship of net water column productivity with algal biomass (standing crop), nutrient loading and turbidity were determined experimentally. Correlations of bacterial abundance with community respiration and nitrification rates were also explored. Annual integrated phytoplankton production rate of this tidal estuary was estimated to be 151.07 gC m-2 y-1. Gross primary productivity showed marked inter annual variation being lowest in monsoon and highest in postmonsoon period. Conclusion: Average primary production was a function of nutrient loading and light penetration in the water column. High aquatic turbidity, conductivity and suspended particulate matter were the limiting factors to attenuate light penetration with negative influence on primary production. Community respiration and nitrification rates of the estuary were influenced by the bacterial abundance. The estuary was phosphorus limited in postmonsoon whereas nitrogen-limited in premonsoon and monsoon period. High algal biomass and primary productivity indicated the estuary to be in eutrophic state in most of the time throughout the year. Our study also indicated a seasonal shifting between autotrophic and heterotrophic conditions in Sundarban estuarine ecosystem and it is a tropical, well mixed (high tidal influx) and marine dominated (no fresh water connection) system. Keywords: Net ecosystem metabolism, Gross primary productivity, Community respiration, Nitrification, Nutrient load, Sundarban estuary

* Correspondence: bmaitree@gmail.com Equal contributors 2 Department of Biochemistry, University of Calcutta, 35, Ballygunge Circular Road, Kolkata 700019, India Full list of author information is available at the end of the article
2012 Chaudhuri et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 2 of 16

Background Primary productivity of marine and coastal phytoplankton contributes to 15% of global oceanic production. Phytoplankton biomass and primary production mainly depend on nutrient dynamics of coastal and estuarine ecosystems [1,2]. Anthropogenic inputs of excess nutrients to coastal waters is a global problem and has increased dramatically over last few decades [1,3,4]. These interventions affect phytoplankton dynamics and primary productions to a great extent [5]. On the other hand enhanced nutrient loading also stimulates primary production [2,3]. Phytoplankton are known for their rapid responses to altered environmental conditions [6], such as anthropogenically introduced eutrophication of coastal waters [7], alteration of coastal configuration, degradation of environmental conditions e.g. sewage outfalls with the consequent water pollution [8,9] etc. The dynamic structure of phytoplankton communities directly reflects the health of aquatic ecosystems [10]. In India, Ganges-Brahmaputra estuary is particularly vulnerable to anthropogenic perturbations due to high nutrient loads from riverine discharge, increasing human population density and rapid economic growth [11,12]. In Indian Sundarbans a huge quantity of leaf litter is loaded to the estuarine water from surrounding mangrove forests. Land mass washes during monsoon and effluents from shrimp culture farms also contribute to this huge nutrient load. Nutrient enrichment or loading has already been identified as a serious problem in recent times in many estuaries [13] in the world namely Baltic sea, Adriatic sea, Gulf of Maxico, Chesapeake Bay and San Francisco Bay especially during rapid growth of population, agriculture and fertilizer production [1,2]. Rates and patterns of nutrient assimilation vary widely among the coastal ecosystems all over the world. While some estuaries exhibits acute eutrophication leading to enhanced algal biomass and primary production, other nutrient rich estuaries maintain low algal biomass and primary production [2]. Estuarine dynamics has been well studied in temperate system such as Chesapeake Bay [14,15], San Francisco Bay [16], and the Baltic sea [17]. Tropical and subtropical estuaries received comparatively less study but are experiencing noticeable anthropogenic alterations [18,19]. Estuaries are highly dynamic systems linking land to the ocean with large seasonal and spatial gradients of biogeochemical compounds and processes [20]. They are often greatly influenced by anthropogenic activities, including enhanced organic matter and nutrient loadings. Net ecosystem metabolism (NEM), the net effect of production and respiration, in coastal and estuarine environments is a means to evaluate whether such environments are sources or sinks of carbon [21,22].

NEM is positive (autotrophic) when production exceeds respiration and negative (heterotrophic) when respiration exceeds production. Most estuaries are net heterotrophic [21,22] and generally a net source of CO2 to the atmosphere [23-25]. The balance between organic matter and nutrient loading is critical in determining the balance between autotrophy and heterotrophy at the ecosystem level [26,27]. Our study area, Jharkhali estuary is a part of HooglyMatla estuary situated in Indian Sundarbans surrounded by the mangrove forests (Figure 1). This macrotidal estuary experiences a subtropical monsoon climate with annual rainfall of about 16001800 mm and is characterized by high tidal energy, moderate to high temperature (22C - 34C) and high nutrient influx. Our study period consists of three seasons namely postmonsoon, premonsoon and monsoon. Postmonsoon period extends from November to February and is characterized by moderate to low temperature (21-27C), moderate salinity (1417 PSU) and intermediate aquatic turbidity (2555 NTU) with little or no rainfall. Premonsoon period extends from March to June and is characterized by moderate to high temperature (2734C), high salinity (2124.8 PSU) and low aquatic turbidity (2025 NTU) with little or no rainfall but occasionally with premonsoon shower. Monsoon period extends from July to October and is characterised by high temperature (29-32C), low salinity (1215 PSU) and high aquatic turbidity (57125 NTU) with moderate to heavy rainfall. Present work reports the productivity of this subtropical estuary and explores its correlation with standing crop, nutrient concentration, turbidity, tidal fluctuation and tidal current velocity. Many other physicochemical parameters like temperature, conductivity, dissolved oxygen concentration (DO), sechhi depth were also monitored throughout the study period. Bacterial abundance was also estimated and its relationship with community respiration and nitrification rates was worked out. In a nutshell, this work explores the estuarine metabolism as a function of both physico-chemical and biological processes.

Materials and method

The study area includes Jharkhali Island in the northcentral part. Two major channels namely Matla and Bidya run respectively on eastern and western sides of Jharkhali. These rivers are connected with each other by a prominent creek, locally known as Bara Herobhanga Khal which is about 150 meters in width. Apart from this creek, several smaller creeks and rivulets present within the study area also act as conduits of water from both sides. Herobhanga Reserve Forest is located just south of

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 3 of 16

Figure 1 Map of the Sundarbans showing sampling area.

the Jharkhali Island. Both Bidya and Matla rivers are flanked by inhabited islands and mangrove islands. Altogether 11 (eleven) sample points were selected (Table 1, Figure 1). Out of these eleven points, points 5 and 6 lie in Matla River and Points 1, 2, 3, 4 and 11 are within Bidya River, while points 7, 8, 9 and 10 lie within the Baro Herobhanga Khal (creek) joining these two mighty rivers. These points have been selected keeping in mind the hydrodynamics of the region.
Sample collection

Table 1 Details of the sample points in Sundarban estuary with specific latitude and longitude
Details of Sample Points Point No. 1 Point No. 2 Point No. 3 Point No. 4 Point No. 5 Point No. 6 Point No. 7 Point No. 8 Point No. 9 Point No. 10 Point No. 11 Latitude 22 04' 13.42" 22 04' 23.79" 22 02' 12.82" 21 59' 00.15" 22 02' 49.17" 22 01' 37.96" 22 01' 25.60" 22 01' 07.78" 22 00' 41.28" 22 00' 37.59" 21 56' 19.61" Longitude 88 41' 56.29" 88 45' 03.48" 88 44' 45.11" 88 44' 21.99" 88 39' 39.82" 88 39' 13.55" 88 39' 52.08" 88 40' 55.80" 88 41' 17.43" 88 42' 55.04" 88 40' 56.46"

Study period extended from November08 to October11 postmonsoon (November-February), premonsoon (MarchJune) and monsoon (July-October). Samples were collected from the water surface (0.5m depth) of all the eleven stations along the Matla and Bidya river. Field trips were arranged at each fortnight and approximately sixty trips

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 4 of 16

were conducted to collect the samples. Collected samples were filtered for the study of physicochemical and biological parameters. Samples were preserved in cold condition and transported to the laboratory within three hours of collection to analyze immediately.
Physico-chemical analysis

microscope according to UNESCO PROTOCOL [31]. Total number of phytoplankton present in a litre of water sample was calculated using the following relation: N n x v=V Where, N = Total number of phytoplankton cells per litre of water filtered. n: Average number of phytoplankton cells in 1 ml of plankton sample taken for counting. v: Volume of plankton concentrate (ml) V: Volume of total water filtered (l) Standing crop/ m3 of water = N X 1000 Quantitative estimation of bacterial population was performed using Zeiss confocal fluorescence microscope. Bacterial cells were counted on black polycarbonate nucleopore membrane (pore size: 0.22 m) using acridine orange stain [32,33].
Primary productivity

A suit of environmental variables such as temperature, pH, conductivity, turbidity, light, tidal fluctuation, current speed, nutrients were measured along with primary productivity for three years study period in Sundarban estuary. Water temperature (accuracy: 0.3C), pH (accuracy: 0.01) and salinity (accuracy: 0.5% of reading) were measured in situ with Hach Portable Meters (HQ40d) and turbidity was measured by using portable turbidity meter (Hach 2100P, Turbidity accuracy: 0.1NTU). Dissolved oxygen (DO) concentration was determined according to Winkler method [27], nutrients like total nitrogen (TN), total phosphate (TP), silicate and ammonia were measured according to the standard methodology [28]. The light attenuation coefficient was evaluated by an empirically derived relationship of 1.57/ Secchi disc depth [3]. Subsamples for suspended load (SPM) were filtered through pre-weighed 0.45 pm Millipore filters. They were stored in small petri dishes, dried for 24 h at 60C and reweighed. The difference between the two weights equaled the amount of suspended load in the water [29]. Tide measurement was performed by Valeport MIDAS WTR non directional tide gauge serial no. 34890 (Valeport, U.K). The MIDAS WTR Wave Recorder was used to prove Linear Wave Theory analysis method of measurement. It has high accuracy, piezo-resistive pressure sensors and a fast response PRT temperature sensor as standard. Current speed and direction was measured with Aanderra made Doppler Current Sensor 4420 Serial no. 282 Signal type CANbus.
Biological analysis Phytoplankton biomass (Chlorophyll-a)

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

Phytoplankton standing crop (biomass) was estimated by cell counting method. Phytoplankton were collected from surface water and preserved after treatment with Lugols iodine and buffered formaldehyde. Subsamples (12 liter) were used for quantitative enumeration utilizing a Sedgwick-Rafter counting chamber and Zeiss research

Primary productivity in a water body can be determined by different methods following the guidance of APHA [34]. We adopted the oxygen method where changes in oxygen concentration were measured in light and dark bottles. Planktonic photosynthesis rates were derived in terms of evolved oxygen in the process. These rates were converted into carbon units assuming photosynthetic quotient (PQ) of 1.3 based on C: N: P molar elemental composition of phytoplankton [35], and a respiratory quotient (RQ) of 1 [36]. The advantage of this method is that it estimates gross (GPP) and net (NPP) productivity along with community respiration (CR). One of the sources of uncertainty in estimating gross primary production (GPP) is the metabolic activity of chemoautotrophic bacteria. This bacterial population is very active in estuaries and consumes significant quantity of oxygen in nitrification process. The oxygen consumption should be included to avoid overestimation of community respiration (CR). This problem was resolved by introducing one additional set of incubation containing nitrification (NR) inhibitor (Nitrapyrine and Chlorate, 5 mg/L and 10mM respectively) [37]. Samples were collected at each preselected depth (on the basis of light availability) in light bottles, dark bottles (with and without nitrification inhibitor) and in initial analysis bottles, each in triplicates. Light and dark bottles were suspended at the collection depth and incubated in the middle of the photic zone for 12 hours (from 6 AM to 6 PM). Dissolved oxygen of the initial bottles were fixed with NaI-NaOH and MnCl2 in the beginning of incubation period. At the end of the

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 5 of 16

incubation period light and dark bottles were similarly fixed and all the bottles were brought back to laboratory in cold condition for analysis. Then dissolved oxygen concentrations were estimated by Winklers method [28]. In the dark bottles, amount of oxygen was reduced due to its consumption in respiration of plants, animals and bacteria, known as community respiration (CR). But in addition to community respiration, nitrification (NH4 NO3) also consumed oxygen from water. Therefore, two metabolic processes i.e. respiration and nitrification simultaneously consume oxygen from water in the dark bottles. But in the dark bottle with nitrification inhibitor (NI), only respiration was responsible for oxygen consumption and therefore change in oxygen concentration provided the measure of community respiration (CR) only. In the dark bottle without nitrification inhibitor both respiration and nitrification were involved in oxygen consumption. So, the difference in oxygen concentration in the dark bottle with and without NI provided nitrification rate. In the light bottle all the three metabolic processes namely photosynthesis (produced oxygen), respiration (consumed oxygen) and nitrification (consumed oxygen) occurred simultaneously. The difference in oxygen concentration in the light bottle provides net primary productivity (NPP). Gross primary productivity (GPP) is a measure of the sum of net primary productivity (NPP), community respiration (CR) and nitrification (NR). So, NPP, CR, NR and GPP can be estimated using the following relations: NPP Light bottleInitial bottle CR Initial bottleDark bottle with nitrification inhibitor: NR Dark bottle with nitrification inhibitorDark bottle without nitrification inhibitor GPP Net primary productivity community respiration Nitrification rate:
Statistical treatment

Results
Physicochemical parameters

Temperature of the estuary showed predicted range of variability indicating lowest in postmonsoon in the month of January (21.5C) and highest in premonsoon in the month of June (33.5C) (Table 2). The pH of water was weakly alkaline and more or less constant throughout the study period (88.15) (Table 2). The salinity of water increased gradually from postmonsoon (17.3 PSU) to premonsoon period (24.5 PSU) and decreased to a lowest value in monsoon (12.6 PSU). Highest salinity was observed in June (24.5 PSU) and lowest in October (12.6 PSU). Moderate to high DO concentration (6.5-9.8 mg/L) was observed throughout the year, highest in January (9.8 mg/L) and lowest in June (6.5 mg/L) (Table 2). Maximum aquatic turbidity was observed in monsoon in the month of October (125 NTU) followed by postmonsoon in the month of November (55 NTU) and premonsoon in the month of March (25 NTU), also evidenced by SPM concentrations (248.4 mg/L, 172.1 mg/L and 87.8 mg/L) respectively. This observation was also reflected in secchi disc indices (20.0cm, 52.3cm, 137.5cm) and light attenuation coefficients (0.078, 0.030, 0.011). Highest aquatic turbidity and SPM concentration were evidenced in October (125 NTU, 248.4 mg/L) and lowest in June (20 NTU, 57.5 mg/L) (Table 3). Nutrient concentration i.e. total nitrogen (TN), ammonianitrogen, total phosphate (TP) and silicate showed higher concentration in postmonsoon and monsoon compared to premonsoon (Figure 3a). TN and ammonia-nitrogen were estimated to be 34.14 mol/L, 2.06 mol/L in postmonsoon, 20.52 mol/L,1.24 mol/L in premonsoon and 28.22 mol/L, 1.36 mol/L in monsoon (averaged) respectively. Highest TN and ammonia-nitrogen concentration were observed in the month of February (36.25 mol/L) and January (2.3 mol/L) respectively and lowest
Table 2 Monthly variation of physical parameters in Sundarban estuary
Month November December January February March April May June July August September October Temperature (C) pH 27.50.1 24.50.1 21.50.1 25.00.1 27.50.1 32.50.1 33.50.1 34.00.1 33.50.1 32.00.1 30.50.1 29.00.1 7.900.1 8.00.15 8.100.1 8.200.1 8.100.12 8.050.1 7.950.14 7.900.1 7.950.15 7.900.14 7.900.1 7.850.11 Salinity (PSU) DO (mg/L) 14.9 0.2 15.5 0.1 16.6 0.1 17.3 0.1 21.2 0.2 22.3 0.2 23.6 0.1 24.7 0.1 15.0 0.2 14.1 0.1 13.3 0.2 12.7 0.1 8.40.15 8.50.12 9.20.2 9.80.18 7.10.1 7.00.12 6.70.2 6.50.15 6.60.14 6.90.12 7.60.15 7.80.18

During three years study period (November2008 October2011) approximately sixty field trips were conducted (twice in a month). Samplings were made from eleven stations in triplicate and they were treated as subsamples. Reading of each month represented an average of 180 samples. The results were expressed as differences between the groups considered significant at p<0.05. Different statistical analysis and correlation regression analysis were performed using the software STATISTICA and in all cases p<0.01 indicating the relations were highly significant [38].

Each value represents mean of 180 samples.

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 6 of 16

Table 3 Monthly variation of light availability parameters in Sundarban estuary

November December January February March April May June July August September October

555.0 485.5 394.8 255.3 250.6 240.4 230.4 200.8 574.8 1103.2 1142.5 1256.2

Each value represents mean of 180 samples.

in the month of June (14.15 mol/L, 0.97 mol/L). Total phosphate concentration was estimated to be 2.09 mol/L in postmonsoon, 1.44 mol/L in premonsoon, 1.73 mol/L in monsoon (averaged), being highest in February (2.15 mol/L) and lowest in June (1.08 mol/L). Silicate concentration was observed to be 24.23 mol/L in postmonsoon, 14.58 mol/L in premonsoon and 22.47 mol/L in monsoon (averaged) (Figure 3a). Highest silicate concentration was recorded in October (29.58 mol/L) and lowest in June (11.02 mol/L). TN:TP ratio was greater than Redfield ratio (16:1) in postmonsoon (averaged 17.3) and less than redfield ratio in premonsoon (averaged 14.5) and monsoon (averaged 14.7) (Figure 3b). Tidal fluctuation ranged from 5.1 to 6 metre in height depending on spring to neap tide condition (Figure 3a). Tidal variation was observed to be quite high ranging from 56 metre in average and tidal wave travel from sea side which is nearly 50 km towards south to 30 km upstream to Port Canning. Current speed, which is the ultimate driving force for nutrient distribution in an estuary and it was observed to be quite high in this region, ranging from 140180 cm /sec (Figure 2a and 2b). Sundarban estuary is characterized by subtropical monsoon climate with an annual rainfall of about 16001800 mm, maximum recorded in monsoon followed by premonsoon and postmonsoon (Figure 2c).
Biological parameters

in June (4.85 g L-1, 1.07 x108 cells m-3 ) (Figure 4, Figure 5a and 5b).Chlorophyll-a concentration was observed to be 33.15 g L-1 in postmonsoon , 16.02 g L-1 in premonsoon and 19.05 g L-1 in monsoon (averaged) (Figure 4). Phytoplankton standing crop ranged from 3.9X1010 cells m-3 in postmonsoon, 1.34X109 cells m-3 in premonsoon and 1.39X109 cells m-3 in monsoon (averaged) (Figure 5a and 5b). Bacterial abundance ranged from 5.14X105 cells ml-1 in postmonsoon, 1.77X107 cells ml-1 in monsoon and 2.68X107 cells ml-1 in premonsoon (averaged), being the highest in June (4.9X107 cells ml-1) and lowest in January (1.18 X 104 cells ml-1) (Figure 5c). Phytoplankton primary productivity also followed a seasonal cycle being highest (597.3 mgC m-2 d-1) in the month of February (postmonsoon) and lowest (311.0 mgC m -2 d-1) in the month of August (monsoon) (Figure 6). Primary productivity showed high degree of positive correlations with temperature (r=0.79, p=0.002) (Figure 7a), DO (r=0.90, p= 0.00008) (Figure 7b), TN (r=0.80, p=0.0016) (Figure 7d), TP (r=0.80, p=0.0017) (Figure 7e), silicate (r=0.81, p=0.0011) (Figure 7f ) and phytoplankton standing crop (r=0.71, p=0.0095) (Figure 7g) and negative correlation with turbidity (r=0.66, p=0.013) (Figure 7c). Photosynthesis depends on large number of external conditions, a regression equation was derived to describe the influence of various parameters like dissolved oxygen, turbidity, total nitrogen, total phosphate silicate and ammonia on NPP. NPP 239:76 39:04DO 1:07TUR 11:21TN 149:96 TP 25:41AMMONIA 14:19SILICATE Where R2 = 0.85 and NPP = Net primary productivity (mgC m-2 d-1) DO = concentration of dissolved oxygen (mg L-1) Tur = turbidity (NTU) TN = concentration of total nitrogen (M) TP = concentration of total phosphate (M) The coefficient of determination (R2 = 0.85) was relatively high and the relationship described by this equation was highly significant (p < 0.0001). The pH and salinity level of water were not included in this equation as no significant correlations were detected between them and NPP. There were more than 500 observations in the data sets that include NPP for possible use in the equation. A sensitivity test was performed to identify the most dominant parameters. According to their importance the parameters can be

Annual cycle of phytoplankton biomass (chlorophyll-a concentration) and abundance (phytoplankton standing crop) evidenced seasonal variation being highest in February (43.80 g L-1, 1.37 x1011 cells m-3 ) and lowest

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 7 of 16

Figure 2 (See legend on next page.)

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 8 of 16

(See figure on previous page.) Figure 2 a A representative profile of Tidal fluctuation at Jharkhali over the period of one month, showing mixed semidiurnal tidal pattern, with two consecutives spring tide nearly 15 days interval with a neap phase in between. Maximum tidal range on Day 07 (6.1 metre) coincides with Blue moon event on Day 06. Subsequent highest peak on Day 22, having average tidal range 5.1metre corresponds to next new moon event. b Tidal current at Jharkhali for the period of one month, showing highest speed on Day 07 with velocity 140 cm/ sec, coinciding with maximum tidal range of Blue moon period and next highest speed attains on Day 22, matches on New moon event on the same day. c Distribution of monthly precipitation (mm) in Sundarbans. The precipitation data was averaged over 9 years (20022011) showing highest in July (21%) and lowest in December (0.4%).

organized in descending order: total phosphate> DO> silicate> turbidity> total nitrogen> ammonia. Community respiration and nitrification rates were found to be the highest in premonsoon in the month of June (126.3 mgC m2 d-1 and 126.6 mgC m2 d-1) and lowest in postmonsoon in the month of January (80.3 mgC m2 d-1 and 60.5 mgC m2 d-1) for Sundarban estuary (Figure 6). Community respiration and nitrification rates indicated positive correlations with temperature (r= 0.94,

p=0.000007 and r=0.98, p=0.000000008) (Figure 8a and Figure 9a). Both community respiration and nitrification evidenced positive correlations with bacterial abundance (r=0.72, p=0.009 and r=0.82, p=0.0012) (Figure 8b and Figure 9d). Nitrification showed negative correlations with DO and ammonia concentration (r=0.86, p=0.0003 and r= 0.93, p =0.00001) (Figure 9b and 9c). There was a seasonal shifting between autotrophic and heterotrophic conditions in the estuary, net autotrophic in

Figure 3 a Monthly variation of nutrient concentrations such as total nitrogen (mol/L), silicate (mol/L), total phosphate (mol/L) and ammonia (mol/L) in Sundarbans estuary. Each value represents mean of 180 samples. b Monthly variation of redfield ratio (Total Nitrogen/ Total Phosphorus) ina Monthly variation of nutrient concentrations such as total nitrogen (mol/L), silicate (mol/L), total phosphate (mol/L) and ammonia (mol/L) in Sundarbans estuary Sundarbans estuary.

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 9 of 16

Table 4 Gross primary production for major world estuaries and bays: a comparative assessment
Estuary/Bay Dellware Estuary Chesapeake Bay Hudson river estuary Sundarban estuary Ems Dollard (middle) estuary San Francisco Bay Fraser river estuary Columbia river estuary Wassaw Sound estuary Primary Production Referrences (g C m-2 y-1) 307 200 180 151.07 100 140 130 120 90 90 Coljn 1983 [49] Peterson 1979 [50] Parson et al. 1970 [51] Anderson 1975 [52] Turner et al. 1979 [53] Pennock and Sharp 1986 [47] Cloern 1979 [48] Coljn 1983 [49]

postmonsoon and net heterotrophic in premonsoon and monsoon (Figure 10).

Discussion Estuaries are dynamic ecosystems, forming a transition zone between river and ocean environments. They are subject to both marine influence, such as tides, waves and the influx of saline water and other riverine influences. The primary productivity of an ecosystem also depends on the favorable status of the physical characteristics and hydrodynamic conditions of the estuary [39]. Phytoplankton species composition, richness, population density, and primary productivity will vary from coast to coast and sea to sea depending upon the varying hydro biological features [40]. Our work expresses the metabolism of Sundarbans estuary as a function of physico-chemical and biological processes. Gross primary production (GPP) predictability evidenced high degree of positive correlation with phytoplankton abundance and standing crop i.e. biomass. This

Figure 4 Monthly variation of chlorophyll-a concentration (g/L) in Sundarbans estuary. Each value represents mean of 180 samples.

is expected as higher level of plankton biomass enhances the sites for photosynthesis. A moderate degree of correlation was observed between GPP and temperature as indicated by a negative slope. GPP was found to be declined above an optimal temperature of 25.5C and it was also positively correlated with DO level of water. Fairly high concentration of DO throughout the year indicated consistently high wave action (Figure 2a and 2b) and steady state primary productivity of the estuary. Turbidity of water influenced the primary productivity through its effect on light attenuation indicating negative exponential correlation with gross primary production [41]. The pH value of water did not bear any significant correlation with gross primary production and remained almost constant throughout the year. This is contrary to general notion that removal of CO2 from water during photosynthesis enhances the pH value [42]. Primary productivity increased steadily from November to February, despite the fact that salinity of water was increased during the study period. Thus, primary production was independent of salinity as long as the salinity regime did not exceed a certain tolerance level. But this pattern broke down at high salinity values in the month of April. At salinity level 24.7 PSU primary productions as well as planktonic biomass was reduced indicating that the threshold tolerance level of majority of planktonic species was crossed [38] in the estuary. The tidal fluctuation is the main driving force behind distribution of nutrients along the entire estuary. The distribution pattern and dispersion of the nutrients and pollutants within the estuary is dependent on the different motions of the water body. In estuaries, nutrient availability is generally adequate to support primary productivity [43]. Gross primary production bore high degree of positive correlation with total nitrogen concentration. Similarly, primary productivity was positively correlated with total phosphate content. Highest phytoplankton standing crop (biomass) and productivity was observed during post monsoon and magnitude of these maxima was closely related to nitrogen and phosphorus availability in the aquatic system. These findings were in good agreement with Boynton et al. [14]. This estuary was characterized by high level of total nitrogen and moderate to low level of total phosphate (Figure 2). High leaf litter loading from surrounding mangrove forest was the possible reason for nitrogen enrichment of the estuary. Primary growth limiting nutrient for phytoplankton appeared to be predictable from total nitrogen: total phosphorous (TN:TP) ratios [44]. During postmonsoon TN:TP ratio was greater than Redfield ratio (16:1), so the estuary was phosphorus- limited in this time whereas it was nitrogen-limited during monsoon and pre monsoon (Figure 3b) since TN:TP value was less than Redfield ratio [44-46]. Silicate concentration was

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 10 of 16

Figure 5 a Monthly variation of phytoplankton standing crop (cells/m3) in Sundarbans estuary. Each value represents mean of 180 samples. b Monthly variation of phytoplankton abundance (x 105cells/L) in Sundarban estuary. Each value represents mean of 180 samples. c Monthly variation of bacterial abundance (cells/ ml) in Sundarbans estuary. Each value represent mean of 180 samples.

observed to be higher throughout the year (Figure 3a), the highest silicate concentration being recorded in monsoon (October) and late monsoon (November). Silicate is a primary growth limiting nutrient for diatoms. Silicate concentration was found to be highest during monsoon and late

monsoon due to run-off and seasonal opening of upstream channels of Sundarbans. Silicate concentration also evidenced a positive correlation with primary productivity. Phytoplankton production in the estuary displayed marked seasonal variability, lowest value being recorded

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 11 of 16

Figure 6 Monthly variation of net primary production (NPP, mg C m-2d-1), community respiration (CR, mgC m-2d-1) and nitrification (NR, mgC m-2d-1) in Sundarbans estuary. Each value represents mean of 180 samples.

in August (monsoon) (314.0 mg m-2d-1) and highest in February (postmonsoon) (597.3 mg m-2d-1) with an average of 151.07 gC m-2 Y-1 compared to other nutrient rich estuarine estuaries of the world (Table 4) [47-53]. The NPP of the estuary reached its peak in the post monsoon months. An examination of physical and chemical parameters revealed that most of them stay at a comfort-zone favouring high level of photosynthesis during this time. Total nitrogen, phosphate and silicate levels of estuarine water reached their maximum level, making the estuarine water nutrient rich for metabolic activity. This nutrient enrichment of water was achieved by land mass wash off from adjoining agricultural lands, drainage of waste from shrimp culture farms and upwelling of sediments during monsoon. Among physical factors, light penetration and temperature remained conducive for photosynthesis. The cyclonic activity died down during this time allowing estuarine water to become calm and clean, thereby allowing greater light penetration. The temperature also remained at the optimal level. At the premonsoon season, this favourable condition for high productivity gets disturbed due to drastic change in a few parameters. The dominant parameters adversely affecting NPP were high temperature and high salinity. As a result, the estuary showed marked decline in NPP. During monsoon months, the turbidity of water became severely limiting due to high suspended load to block light penetration towards sub-surface levels. High nutrient concentration and favourable conductivity level or high DO of water could hardly help in elevating the NPP value. The community respiration (CR) of the estuary ranged from 91.35 mgC m-2d-1 to 132.3 mgC m-2d-1. The CR value was found to be highest in the month of June and

lowest in January. The CR value is observed to be correlated positively with temperature and bacterial abundance. These findings are in agreement with Iriarte et. al. [54] for a temperate estuary in Northern Spain. These high rates of CR in summer months in the estuary must be associated to enhanced organic carbon input from sources other than the primary production since the later was low during summer. Fresh water flow cannot explain this unknown carbon source to the estuary as there was hardly any rainfall in summer and connection to fresh water river is cut off due to siltation. Thus high CR value during summer may only be explained by an increase of lateral organic carbon input into the estuary possibly from mangrove vegetation in the catchment. This organic carbon consumed and stored into the sediments probably balances the carbon budget of the estuary [55]. Thus planktonic CR in the Sundarbans estuary was primarily driven by allochtonous organic matter rather than local production. Similar feature was also observed in other estuaries of the world e.g. Scheldt estuary [56]. The nitrification rate (NR) of the estuary ranged from 60.5 mg C m-2 d-1 to 123.6 mg C m-2 d-1, attaining the highest in June and lowest in January indicating a positive correlation with ambient temperature [54]. This can be explained on the basis of bacterial abundance in the estuary showing a positive correlation with NR. During the January-June period ammonia-nitrogen content also declined steadily, revealing a negative correlation with NR. Since ammonia is the substrate on which nitrifying bacteria metabolize, the observed pattern is self explanatory. Reduced nitrification level should ideally improve the DO level of water. Elevated DO level in response to decreased nitrification activity in postmonsoon and vice versa in premonsoon establishes this hypothesis i.e. DO level is negatively correlated with nitrification activity.

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 12 of 16

Figure 7 Correlation of gross primary production (GPP, mg C m-2d-1) with ( a) temperature (C) and (b) dissolved oxygen concentration (DO, mg/L), (c) turbidity (NTU), (d) total nitrogen (TN, mol/L) (e) total phosphate (TP, mol/L ), (f) silicate (mol/L) and (g) phytoplankton abundance (cells/m3).

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 13 of 16

Figure 8 Correlation of community respiration (CR, mgC m-2d-1) with (a) temperature (C) and( b) bacterial abundance (cells/ ml).

Figure 9 Correlation of nitrification (NR) with (a) temperature (C) (b) dissolved oxygen concentration (DO, mg/L), (c) ammonia concentration (mol/L) and (d) bacterial abundance (cells/ml).

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 14 of 16

Figure 10 Net ecosystem metabolism in Sundarban estuary.

Estuaries function as a biogeochemical reactor since many autochthonous processes occur here [57,58]. The metabolic state of an ecosystem maintains a dynamic balance between primary production and community respiration [26]. This estuary remained autotrophic for five months (November to March) of the year. During this time the primary production is greater than community respiration resulting in export or burial of organic matter through conversion of inorganic matter and carbon dioxide. The estuary remained heterotrophic during remaining seven months (April to October). Community respiration was greater than primary production during that period and allochthonous materials are re-mineralized leading to production of inorganic nutrients and carbon dioxide. Thus Sundarbans estuary was a net sink of CO2 for five months of the year and net source of CO2 for remaining seven months. If net ecosystem metabolism of the entire year was taken as a whole, the estuary can be designated as a net source of CO2 [22-24,59]. Trophic state of the estuary was worked out using chlorophyll-a concentration (Jones and Fredly, 1982). This estuary showed chlorophyll-a concentration greater than 10 g/L throughout the year , except for the month of May and June (9.07 and 4.85 g/L respectively). Thus, the estuary was identified as mesotrophic only in the months of May and June but eutrophic in the remaining ten months [38,60].

was a function of nutrient loading and light penetration in the water column. This result was the outcome of a nutrient rich ecosystem which also potentiated the growth of phytoplankton. High aquatic turbidity, conductivity and suspended particulate matter were also the other limiting factors to attenuate light penetration with negative influence on primary production. The growth rate of phytoplankton populations was a function of light availability and the variability in biomass can be justified by variation of light exposure. Community respiration and nitrification rates were also influenced by the bacterial abundance in the estuary. Sundarban estuary was identified as a tropical, well mixed (due to high tidal influx) and marine dominated (as there is no fresh water connection) system.
Competing interests The authors declare that they have no competing interests. Authors contributions KC and SM 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. KSS prepared the map and performed the tide and current experiments. PN participated in the field experiments and assisted in preparing manuscript. SB planned the project, was involved in acquisition of funds, selected the site and field stations including the geo-referencing of the stations, prepared the map, guided the field study, sample collection and the tide and current experiments. All authors read and approved the final manuscript. Acknowledgement We are grateful to Ministry of Earth Science, Government of India for financial support. Thanks are due to Dr. B.R. Subramanian, Director, ICMAMPD, Chennai for his constant support. We express our sincere gratitude to Mr. Arijit Banerjee, Director, IESWM for his continuous support and interest in our work. We acknowledge the instrumental facility provided by UGC, DST and World Bank in the Department of Biochemistry. We acknowledge Mr. Ankan Dey, IESWM, Kolkata for working out the rainfall profile of Sundarbans. It would not be possible to carry out this work without the help and support of the local people of Sundarbans. We express our inability to acknowledge them individually. Author details 1 Institute of Environmental Studies and Wetland Management, DD-24, Sector-I, Salt lake, Kolkata 700064, India. 2Department of Biochemistry, University of Calcutta, 35, Ballygunge Circular Road, Kolkata 700019, India.

Conclusion Phytoplankton production in the tidal creeks in Sundarbans estuary averaged 151.07 g C m-2 Y-1 designating this ecosystem to be moderately productive compared to other nutrient rich estuarine ecosystems of the world (Table 4). Different indices and limiting relations achieved from the experiments support the generalized concepts regarding light availability, nutrient load, conductivity and bacterial abundance in an ecosystem. Average primary production

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 15 of 16

Received: 4 April 2012 Accepted: 15 October 2012 Published: 19 October 2012 References 1. Nixon SW: Coastal marine eutrophication: a definition. social causes and future concerns. Ophelia 1995, 31:199219. 2. Cloern JL: The relative importance of light and nutrient limitation of phytoplankton growth: a simple index of coastal ecosystem sensitivity to nutrient enrichment. Aquatic Ecology 1999, 33:316. 3. Murrell MC, Hagy JD III, Lores EM, Greene RM: Phytoplankton production and nutrient distributions in a subtropical estuary: importance of freshwater flow. Estuaries and Coasts 2007, 30(3):390402. 4. Howarth RW, Paerl HW: Coastal marine eutrophication: control of both nitrogen and phosphorus is necessary. Proc Nat Acad Sci USA 2008, 105: E103E103. 5. Cloern JE, Cole BE, Oremland RS: Autotrophic processes in meromictic big soda lake. Nevada. Limnol Oceanogr 1983, 28(6):10491061. 6. Stolte W, Mccllin T, Noordeloos A, Reigman R: Effect of nitrogen-source on the size distribution within marine-phytoplankton populations. J. Exp. Mar. Biol. Ecol. 1994, 184:8397. 7. Richardson K: Harmful or exceptional phytoplankton blooms in the marine ecosystem. Adv Mar Biol 1997, 31:301385. 8. Pereira LCC, Jimnez JA, Medeiros C: Environmental Degradation of the Littoral of Casa Caiada and Rio Doce, Olinda-PE (Brazil). J Coastal Res 2003, 35:502508. SI. 9. Pereira LCC, Jimnez JA, Koening ML, Neto FP, Medeiros C, Da Costa RM: Effect of coastline properties and wastewater on plankton composition and distribution in a stressed environment on the North Coast of Olinda-PE (Brazil). Braz Arch Biol Technol 2005, 48:10131026. 10. Valdes-Weaver LM, Piehler MF, Pinckney JL JL, Howe KE, Rossignol K, Paerl HW: Long-term temporal and spatial trends in phytoplankton biomass and class-level taxonomic composition in the hydrologically variable Neuse-Pamlico estuarine continuum, North Carolina U. S. A. Limnol Oceanogr 2006, 51(3):14101420. 11. Seitzinger SP, Harrison JA, Dumont E, Bensen AHW, Bouwman AF: Sources and delivery of carbon, nitrogen and phosphorous to the coastal zone: an overview of global nutrient export from watersheds (NEWS) models and their application. Global Biogeochem Cycles 2005, 19:GB4S01. doi:10.1029/ 2005GB002606. 12. Biswas H, Dey M, Ganguly D, De TK, Ghosh S, Jana TK: Comparative analysis of phytoplankton composition and abundance over a twodecade period at the land-ocean boundary of a tropical mangrove ecosystem. Estuaries and Coasts 2010, 33:384394. 13. Smith VH: Eutrophication of fresh water and coastal marine ecosystem: a global problem. Environ Sci Pollut Res 2003, 10:114. 14. Boynton WR, Kemp WM, Keefe CW: A cornparative analysis of nutrients and other factors influencing estuarine phytoplankton production. In Estuarine comparisons. Edited by Kennedy V. New York: Academic Press; 1982:6990. 15. Kemp WM, Boynton WR, Adolf JE, Boesch DF, Boicourt WC, Brush G, Cornwell JC, Fisher TR, Glibert PA, Hugy JD, Harding LW, Houde ED Jr, Kimmel DG, Millar WD, Newell RIE, Roman MR, Smith EM, Stevenson JC: Eutrophication of Chesapeake Bay: historical trends and ecological interactions. Mar eco Pro Ser 2005, 303:129. 16. Cloern JE: Phytoplankton bloom dynamics in coastal ecosystems: a review with some general lessons from sustained investigation of San Francisco Bay. California. Rev Geophysics 1996, 33:127168. 17. Conley DJ, Kaas F, Mohlenberg F, Rasmussen B, Windolf J: Characteristics of Danish estuaries. Estuaries 2000, 23:820837. 18. Eyre BD: Water quality changes in an episodically flushed sub-tropical Australian estuary: a 50 year perspective. Mar chem 1997, 59:177187. 19. Turner RE, Rabalais NN: Linking landscape and water quality in the Mississippi River basin for 200 years. Bioscience 2003, 53:563572. 20. Heip CHR, Goosen NK, Herman PMJ, Kromkamp J, Middelburg JJ, Soetaert K, Ansell AD, Gibson RN, Barnes M: Production and consumption of biological particles in temperate tidal estuaries. Oceanography and Marine Biology. An Annual Review 1995, 33:1149. 21. Smith SV, Hollibaugh JT: Coastal metabolism and the oceanic organic carbon balance. Rev Geophys 1993, 31:7589. 22. Gattuso JP, Frankignoulle M, Wollast R: Carbon and carbonate metabolism in coastal aquatic ecosystems. Annu Rev Ecol Syst 1998, 29:405434.

23. Borges AV: Do we have enough pieces of the jigsaw to intergrate CO2 fluxes in the coastal ocean? Estuaries 2005, 28:125. 24. Cai WJ: Estuarine and coastal ocean carbon paradox; CO2 sinks or sites of terrestrial carbon incineration. Ann Rev Mar Sci 2011, 3:123145. 25. Laruelle GG, Drr HH, Slomp CP, Borges AV: Evaluation of sinks and sources of CO2 in the global ocean using spatially-explicit typology of estuaries and shelves. Geophys Res Lett 2010, 37:L15607. 26. Kemp WM, Smith EM, Marvin-Dipasquale M, Boynton WR: Organic carbon balance and net ecosystem metabolism in Chesapeake Bay. Mar Ecol Prog Ser 1997, 150:229248. 27. Eyre BD, Mckee LJ: Carbon, nitrogen, and phosphorus budgets for a shallow subtropical coastal embayment (Moreton Bay, Australia). Limnol Oceanogra 2002, 47:10431055. 28. JGOFS Protocol: Protocol for the Joint Global Flux Study (JGOFS): Core Management. 1994. 29. Harrison PJ, Khan N, Yin K, Saleem M, Bano N, Nisa M, Ahmed SI, Rizvi N, Azam F: Nutrient and phytoplankton dynamics in two mangrove tidal creeks of the Indus River delta, Pakistan. Mar Ecol Prog Ser 1997, 157:1319. 30. Ventrick EL, Hayward TL: Determining chlorophyll on the 1984 CalCOFI surveys. CalCOFI Rep 1984, 25:7479. 31. UNESCO: Monograph and Oceanographic methodology. In Phytoplankton Manual, UNESCO Technicalpaper. 28th edition. Edited by Sournia A. PARIS: UNESCO; 1978:337. 32. Hobbie JE, Daleyr RJ, Jasper S: Use of nuclepore filters for counting bacteria by fluorescence microscopy. Appl Environ Microbiol 1976, 33:12251228. 33. Bianchi A, Giuliano L: Enumeration of viable bacteria in the marine pelagic environment. Appl Enviro Microbiol 1995, 62:174177. 34. APHA: Standard Methods for the Examination of Water and Wastewater. Washington, DC: American Public Health Association; 1998. 35. Redfield AC, Ketchum BH, Richards FA: The influence of organisms on the composition of seawater. In The sea. Interscience. 2nd edition. Edited by Hill MN. 1963:2663. 36. Hopkinson CSK, Smith EM: Estuarine respiration: an overview of benthic, pelagic and whole system respiration. In Respiration in aquatic ecosystem. Edited by Williams PJLB, del Giorgio PA. Oxford: Oxford University Press; 2005:123147. 37. Berounsky VM, Nixon SW: Rates of nitrification along an estuarine gradient in Narraganset Bay. Estuaries 1993, 16:718730. 38. Manna S, Chaudhuri K, Bhattacharyya S, Bhattacharyya M: Dynamics of Sundarban estuarine ecosystem: eutrophication induced threat to mangroves. Saline Systems 2010, 6:8. 39. Webster IT, Ford PW, Robson B, Margvelashvili N, Parslow JS: Conceptual models of the hydrodynamics, fine-sediment dynamics, biogeochemistry,and primary production in the fitzroy estuary. Coastal CRC Project CM-2. Technical Report 2003, 8:141. 40. Prabhahar C, Saleshrani K, Enbarasan R: Studies on the ecology and distribution of phytoplankton biomass in Kadalur coastal zone Tamil nadu. India. Curr. Bot. 2011, 2(3):2630. 41. Cloern JE: Turbidity as a control on phytoplankton biomass and productivity in estuaries. Continental shelf Research 1987, 7:13671381. 11/12. 42. Sharma SK: Adsorptive iron removal from groundwater. PhD dissertation. The Netherlands: Wageningen University / IHE Delft; 2001a:187198. Sharma SK: Comparison of physicochemical iron removal mechanisms in filters. In Journal of water supply: Research and Technology aqua. 2001b. 43. Cole BE, Cloern JE: An empirical model for estimating phytoplankton productivity in estuaries. Mar Ecol Prog Ser 1987, 36:295305. 44. Smith VH: Responses of estuarine and coastal marine phytoplankton to nitrogen and phosphorus enrichment. Limnol Oceanogr 2006, 51(1, part 2):377384. 45. Downing JA: Marine nitrogen: phosphorus stoichiometry and the global N: P cycle. Biogeochemistry 1997, 37:237252. 46. Sarma VVSS, Prasad VR, Kumar BSK, Rajeev K, Devi BMM, Reddy NPC, Sarma VV, Kumar MD: Intra-annual variability in nutrients in the Godavari Estuary, India. Cont Shelf Res 2010, 30(19):20052014. 47. Pennock JR, Sharp JH: Phytoplankton production in Delaware estuary: temporal and spatial variability. Mar Ecol Prog Ser 1986, 34:143155. 48. Cloern JE: Phytoplankton ecology of the San Francisco Bay system: the statues of our current understanding. In San Francisco Bay: the urbanized estuary. Edited by Conomos JT. San Francisco: Allen Press; 1979:247264.

Chaudhuri et al. Aquatic Biosystems 2012, 8:26 http://www.aquaticbiosystems.org/content/8/1/26

Page 16 of 16

49. Colijn F, Ludden E: Primary production of phytoplankton in the EmsDollard estuary. In Primary production in the Ems-Dollard estuary. Edited by Colijn F. Groningen: Drukkerij Van Denderen; 1983:3899. 50. Peterson DH: Sources and sinks of biologically reactive oxygen, carbon, nitrogen, and silica in Northern San Francisco Bay. In San Francisco Bay, The urbanized estuary. Edited by Conomos TJ. San Francisco: Allen Press; 1979:175193. 51. Parson TR, Le Brasseur RJ, Barracloug WE: Levels of production in the pelagic environment of the Strait of Georgia, British Columbia: A review. J Fish Res Board of Can 1970, 27:12511264. 52. Anderson GC: Aspects of marine phytoplankton studies near the Columbia River, with special reference to a subsurface chlorophyll maximum. In The Columbia River estuary and adjacent oceanwaters. Edited by Prljter AT, Alverson DL. Seattle: University of Washington; 1972. 53. Turner RE, Woo SW, Jitts HR: Phytoplankton production in a turbid, temperate salt-marsh estuary. Est Coast Mar Sci 1979, 9:603613. 54. Iriarte A, de Madariaga I, Diez-Garagarza F, Revilla M, Orive E: Primary plankton production, respiration and nitrification in a shallow temperate estuary during summer. J Exp Mar Biol Ecol 1996, 208:127151. 55. Gocke K, Cortes J, Murillo MM: The annual cycle of primary productivity in a tropical estuary. the inner regions of the Golfo de Nicoya, Costa Rica. Rev Biol Trop 2001, 49:289306. 56. Gazeau F, Gattuso J, Middelburg JJ, Brion N, Schiettecatte L, Frankignoulle M, Borges AV: Planktonic and whole system metabolism in a nutrient-rich estuary (the Scheldt Estuary). Estuaries 2005, 28(6):868883. 57. Savoye N, David V, Morisseau F, Etcheber H, Billy I, Charlier K, Abril G, Oggian G, Derriennic H, Sautour B: Origin and composition of particulate organic matter in a macrotidal turbid estuary: The Gironde Estuary, France. Est, Coast Shelf Sci 2012, 30:13. 58. Ram ASP, Nair S, Chandramohan D: Seasonal shift in net ecosystem production in a tropical estuary. Limnol Oceanogr 2003, 48(4):16011607. 59. Gupta GVM, Shoji D, Thottathil K, Balachandran K, Madhu NV, Madeswaran P, Nair S: CO2 supersaturation and net heterotrophy in a tropical estuary (Cochin, India): Influence of Anthropogenic effect. Ecosystems 2009, 12(12):11451157. 60. Jones RA, Fred Lee G: Recent advances in assessing impact of phosphorus load on eutrophication-related water quality. Water Res 1982, 16:503515.
doi:10.1186/2046-9063-8-26 Cite this article as: Chaudhuri et al.: Physicochemical and biological factors controlling water column metabolism in Sundarbans estuary, India. Aquatic Biosystems 2012 8:26.

Submit your next manuscript to BioMed Central and take full advantage of:
Convenient online submission Thorough peer review No space constraints or color gure charges Immediate publication on acceptance Inclusion in PubMed, CAS, Scopus and Google Scholar Research which is freely available for redistribution
Submit your manuscript at www.biomedcentral.com/submit