Seasonal Variations of Carbon Dioxide, Water Vapor and Energy Fluxes in Tropical Indian Mangroves

We present annual estimates of the net ecosystem exchange (NEE) of carbon dioxide (CO 2) accumulated over one annual cycle (April 2012 to March 2013) in the world's largest mangrove ecosystem, Sundarbans (India), using the eddy covariance method. An eddy covariance flux tower was established in April 2012 to study the seasonal variations of carbon dioxide fluxes due to soil and vegetation-atmosphere interactions. The half-hourly maximum of the net ecosystem exchange (NEE) varied from´6 µmol¨m´2¨s´1 during the summer study period, the study area was a carbon dioxide sink with an annual net ecosystem productivity (NEP = ´NEE) of 249 ˘ 20 g¨C m ´2¨year´1. The mean annual evapotranspiration (ET) was estimated to be 1.96 ˘ 0.33 mm¨day´1. The gap-filled NEE was also partitioned into Gross Primary Productivity (GPP) and Ecosystem Respiration (R e). The total GPP and R e over the study area for the annual cycle were estimated to be 1271 g C m ´2¨year´1 and 1022 g C m ´2¨year´1 , respectively. The closure of the surface energy balance accounted for of about 78% of the available energy during the study period. Our findings suggest that the Sundarbans mangroves are currently a substantial carbon sink, indicating that the protection and management of these forests would lead as a strategy towards reduction in carbon dioxide emissions.


Introduction
Carbon dioxide (CO 2 ) has been recognized as the most important greenhouse gas responsible for global warming-related climate change.Intergovernmental panel on climate change (IPCC 2013) reported a significant increase of 40% in the global atmospheric concentration of CO 2 from a pre-industrial value of about 278 ppm in 1750 to 390.5 ppm in 2011.The annual growth rate of 2.0 ˘0.1 ppm per year in the last decade (i.e., 2002-2011) was observed to be higher than any previous measured decadal rates [1].The concern about ever-increasing global warming has pushed global interest to understand the potentials of terrestrial ecosystems such as forests to reduce the rise in CO 2 levels.Further, the quantification of terrestrial carbon stocks and the understanding of temporal and spatial variation in carbon uptake would help in developing mitigation strategies for increasing atmospheric carbon dioxide.
The terrestrial ecosystems, especially forests, are major CO 2 sinks, and they sequester about one-third of the total anthropogenic emissions [2].Forests also play a crucial role in global carbon dynamics by exchanging trace gases between the atmosphere and the biosphere [3].Gross primary production (GPP) of the terrestrial ecosystems is the largest global CO 2 flux which is estimated to be 123 ˘8 Pg C year ´1 [4].
In the tropics, mangroves are the most carbon-rich forests with an average estimate of 1023 Mg C ha ´1 [5] and high rates of primary production estimated at 218 ˘72 Tg C year ´1 globally [6].Thus, due to their tropical distribution [7,8] and rich biogeochemistry [9,10], they are one of the important terrestrial ecosystem at the interface of land and ocean to be understood.
Measurements of CO 2 fluxes between forests and atmosphere are important to understand their response to climate change and predict the vegetation feedback to the climate system [4].Eddy covariance (EC) is the most direct way to measure the exchange rate of CO 2 across the soil and vegetation-atmosphere interface at scales of ~1 km 2 using micrometeorological theory for different time scales, ranging from hours to years [11][12][13][14].The net ecosystem exchange (NEE), defined as the net CO 2 exchange between the ecosystem and the atmosphere, is generally approximated by net ecosystem productivity (NEP) as well in many terrestrial ecosystems (NEE = ´NEP) over short time periods [13].The EC method is most widely used method to measure CO 2 fluxes in over 400 flux research sites in FLUXNET (global flux network) which assimilates information from these sites over a wide range of biomes [15].
Though the importance of mangrove ecosystems in the global carbon budget is significant, limited information exists on the carbon assimilation capacity of mangroves [16].The oft-reported flux tower in mangroves was established in the coastal Florida Everglades in 2003.Results suggested that everglade mangroves are a carbon sink that had an annual NEP of 1170 ˘127 g C m ´2¨year ´1 during 2004 [17].
Sundarbans mangroves are the largest continuous stretch of mangrove forests in the world, covering about 10,000 km 2 (i.e., nearly 6%-7%) of the global mangrove area (15 ˆ10 4 km 2 ) spreading over India and Bangladesh [18].To understand and estimate the source and sink strengths of these mangroves, few studies have been carried out in selected locations of the Sundarbans deltaic system.The studies using the micrometeorological gradient technique has revealed that the ecosystem was a CO 2 sink.Mukhopadhyay et al. had observed the environment of Jambu Island (southwest corner of Indian Sundarbans) to be a sink for CO 2 in the pre-monsoon season at a rate of ´56.5 g CO 2 m ´2¨day ´1 [19].Ganguly et al. have reported a mean net flux of ´48.3 g CO 2 m ´2¨day ´1 at the same Lothian Island and Sajnekhali, and have also estimated a total sink strength of 205 Gg¨day ´1 for entire reserve forest areas (4264 km 2 ) [20].During 2011-2012, a series of experiments using micrometeorological measurements were carried out by Chanda et al. in Sundarbans mangroves at different stations viz.Jharkhali, Henry Island and Bonnie Camp to understand the spatial variability of CO 2 flux [21].They have reported that the ecosystem is a sink for CO 2 during the daytime with varying fluxes between ´33.69 to ´114.91 g CO 2 m ´2¨day ´1.
Though several studies have been attempted to understand the source and sink strengths of the Sundarbans mangroves, there is a huge uncertainty in annual estimates due to variability of CO 2 uptake across days over a year.Since the studies so far are mainly for short time periods, the annual CO 2 sink strength of this complex mangrove system remains unknown.Recent studies by Pan et al. focus on estimates of carbon sources and sink strengths in various global forests [22].A huge data gap for the global carbon sequestration estimate in the Indian subcontinent was evident.The present study will also provide the first continuous and accurate measurements of carbon sequestration from Indian mangroves on an annual basis.
This study is an initiative under the Indian Space Research Organisation-Geosphere Biosphere Programme (ISRO-GBP) under the 'National Carbon Project (NCP)'.One of the main objectives of NCP is to establish a network of EC flux towers in different representative forest ecosystems in the Indian subcontinent to understand their CO 2 source and sink strengths [23][24][25].As part of this, an eddy covariance flux tower was established in mangrove forests in the Sundarbans Biosphere Reserve in April 2012 and has been taking measurements since then.The present EC flux tower stands as one of the few towers established in mangrove ecosystems across the world (FLUXNET).The main objective of the project is to continuously measure and understand vegetation-atmosphere exchanges of CO 2 and energy fluxes over Sundarbans (India).The current study is an extension of measuring CO 2 and H 2 O fluxes over the flux tower site and presents the first annual results (April 2012-March 2013) of NEE measurements over the Sundarbans mangroves, India.The paper also discusses NEE measurements partitioned into GPP and R e (ecosystem respiration, both autotrophic and heterotrophic respiration).Further, we discuss the partitioning of net radiation into sensible and latent heat fluxes during the study period.

Study Area
Sundarbans is the largest single block of tidal halophytic mangrove forest in the world (Figure 1).It lies on the Ganga-Brahmaputra delta at the point where it merges with the Bay of Bengal.The Sundarbans mangrove ecosystem is spread across several islands with 34% in Indian territory and 66% in Bangladesh.The total area of the Indian Sundarban region is about 9630 km 2 , out of which reserve forest area occupies 4264 km 2 while the rest of the area is covered by water bodies [20].The average altitude of the Sundarbans mangroves is 7 m and it has a wide biodiversity of tree and shrub species.The name Sundarbans is derived from the dominant plant species Heritiera fomes (locally "sundari").The mangrove system with an unstratified dense canopy is largely undisturbed with no major changes in last few decades (FAO 2005) and it has been a natural habitat for many endangered species including the Bengal tiger.The annual temperature fluctuations in Sundarbans are between 30 ˝C and 21 ˝C, with high temperatures during summer (mid-March to mid-June) and low temperatures during winter (December and January).The Sundarbans region is generally very humid and receives precipitation ranging from1500 mm to 1800 mm.discusses NEE measurements partitioned into GPP and Re (ecosystem respiration, both autotrophic and heterotrophic respiration).Further, we discuss the partitioning of net radiation into sensible and latent heat fluxes during the study period.

Study Area
Sundarbans is the largest single block of tidal halophytic mangrove forest in the world (Figure 1).It lies on the Ganga-Brahmaputra delta at the point where it merges with the Bay of Bengal.The Sundarbans mangrove ecosystem is spread across several islands with 34% in Indian territory and 66% in Bangladesh.The total area of the Indian Sundarban region is about 9630 km 2 , out of which reserve forest area occupies 4264 km 2 while the rest of the area is covered by water bodies [20].The average altitude of the Sundarbans mangroves is 7 m and it has a wide biodiversity of tree and shrub species.The name Sundarbans is derived from the dominant plant species Heritiera fomes (locally "sundari").The mangrove system with an unstratified dense canopy is largely undisturbed with no major changes in last few decades (FAO 2005) and it has been a natural habitat for many endangered species including the Bengal tiger.The annual temperature fluctuations in Sundarbans are between 30 °C and 21 °C, with high temperatures during summer (mid-March to mid-June) and low temperatures during winter (December and January).The Sundarbans region is generally very humid and receives precipitation ranging from1500 mm to 1800 mm.

Eddy Covariance Flux Tower Setup
On an uninhabited island of about 11 km 2 in Sundarbans, the EC flux tower of 15 m was established at 21°49′ N latitude and 88°37′E longitude on 1 April 2012 (Figure 1) to measure CO2 and energy fluxes.The tower is surrounded by a continuous mangrove patch with dominant species viz.Avicennia alba, Bruguiera gymnorrhiza, Rhizophora, with a mean canopy height of about 5 m [23].
A canopy walkway was constructed with cement blocks of 1 m height across the path to the tower in order to reach the EC tower through the marshy land and water-inundated areas (Figure 2).The tower base is established on a 2 m × 2 m cement base to prevent it from sinking.The tower of 15 m is composed of triangular aluminum sections.The well-established guy wires hold the tower straight to perform measurements.All the major tower electronics are housed in a water-proof box at a height of 2 m above the ground.The tower is powered by three solar panels (120 W each) with rechargeable batteries.

Eddy Covariance Flux Tower Setup
On an uninhabited island of about 11 km 2 in Sundarbans, the EC flux tower of 15 m was established at 21 ˝49 1 N latitude and 88 ˝37 1 E longitude on 1 April 2012 (Figure 1) to measure CO 2 and energy fluxes.The tower is surrounded by a continuous mangrove patch with dominant species viz.Avicennia alba, Bruguiera gymnorrhiza, Rhizophora, with a mean canopy height of about 5 m [23].
A canopy walkway was constructed with cement blocks of 1 m height across the path to the tower in order to reach the EC tower through the marshy land and water-inundated areas (Figure 2).The tower base is established on a 2 m ˆ2 m cement base to prevent it from sinking.The tower of 15 m is composed of triangular aluminum sections.The well-established guy wires hold the tower straight to perform measurements.All the major tower electronics are housed in a water-proof box at a height of 2 m above the ground.The tower is powered by three solar panels (120 W each) with rechargeable batteries.Micrometeorological and trace gas (CO2) measurements were made at a height of 10 m.A three-dimensional (3D) sonic anemometer, CSAT-3 (Campbell Scientific Inc., Logan, UT, USA), was used to measure high frequency wind velocity components (i.e., 3D orthogonal wind components ux, uy, uz) and sonic temperature (T).Considering the low power and low maintenance requirement at a remote site location such as Sundarbans (uninhabited and marshy land), an open path infrared gas analyzer system (IRGA), LI-7500A (LI-COR, Lincoln, NE, USA) was installed to measure the densities of carbon dioxide and water vapor.The IRGA was slightly angled to allow surface moisture drainage and prevent direct solar radiation from possibly contaminating readings [26].The raw measurements from the sonic anemometer and IRGA were recorded using an analyzer interface unit, LI-7550 (LI-COR, Lincoln, USA), at a frequency of 10 Hz.The data is collected from 1 April 2012 to 31 March 2013 on a half-hourly basis with short time gaps (about 6%) due to low solar power (mostly occurring during rainy season due to continuous rains).
The additional environmental variables were measured above the canopy at 10 s interval, averaged over 10 min and logged on to CR3000 (Campbell Scientific Inc. Logan, UT, USA) data logger.These measurements include net radiation (CNR 4; Kipp and Zonen, Delft, Netherlands) and incoming photosynthetic active radiation (PAR, LI-190SB; LICOR, Lincoln, USA) at height of 10 m.Measurements also include air temperature and relative humidity (HMP45C-L; Campbell Scientific Logan, UT, USA), wind speed and wind direction (Wind Monitor 5103; RM Young, Traverse, MI, USA) at heights of 2 m, 4 m and 8 m.In addition, soil heat flux (HFT-3; Campbell Scientific, Logan, USA) is measured at two levels (2.5 and 5 cm) and soil temperature (107-B; Campbell Scientific, Logan, UT, USA) at four levels (5, 10, 20 and 40 cm).Precipitation at the study site has been recorded using a tipping bucket TB4-L (Campbell Scientific Inc. Logan, UT, USA).

Primary Raw Data Processing
The initial flux processing was done from 1 April 2012 to 31 March 2013 using Eddypro post-processing software (version 5.0, LI-COR, USA).CO2 and energy fluxes are calculated at a standard averaging interval of 30 min using the fast response measurements from the sonic anemometer (3D wind speed components and sonic temperature) and IRGA (concentrations of CO2 and H2O).Every 30 min, raw data is corrected for standard flux corrections for the installed open Micrometeorological and trace gas (CO 2 ) measurements were made at a height of 10 m.A three-dimensional (3D) sonic anemometer, CSAT-3 (Campbell Scientific Inc., Logan, UT, USA), was used to measure high frequency wind velocity components (i.e., 3D orthogonal wind components u x , u y , u z ) and sonic temperature (T).Considering the low power and low maintenance requirement at a remote site location such as Sundarbans (uninhabited and marshy land), an open path infrared gas analyzer system (IRGA), LI-7500A (LI-COR, Lincoln, NE, USA) was installed to measure the densities of carbon dioxide and water vapor.The IRGA was slightly angled to allow surface moisture drainage and prevent direct solar radiation from possibly contaminating readings [26].The raw measurements from the sonic anemometer and IRGA were recorded using an analyzer interface unit, LI-7550 (LI-COR, Lincoln, USA), at a frequency of 10 Hz.The data is collected from 1 April 2012 to 31 March 2013 on a half-hourly basis with short time gaps (about 6%) due to low solar power (mostly occurring during rainy season due to continuous rains).
The additional environmental variables were measured above the canopy at 10 s interval, averaged over 10 min and logged on to CR3000 (Campbell Scientific Inc. Logan, UT, USA) data logger.These measurements include net radiation (CNR 4; Kipp and Zonen, Delft, Netherlands) and incoming photosynthetic active radiation (PAR, LI-190SB; LICOR, Lincoln, USA) at height of 10 m.Measurements also include air temperature and relative humidity (HMP45C-L; Campbell Scientific Logan, UT, USA), wind speed and wind direction (Wind Monitor 5103; RM Young, Traverse, MI, USA) at heights of 2 m, 4 m and 8 m.In addition, soil heat flux (HFT-3; Campbell Scientific, Logan, USA) is measured at two levels (2.5 and 5 cm) and soil temperature (107-B; Campbell Scientific, Logan, UT, USA) at four levels (5, 10, 20 and 40 cm).Precipitation at the study site has been recorded using a tipping bucket TB4-L (Campbell Scientific Inc. Logan, UT, USA).

Primary Raw Data Processing
The initial flux processing was done from 1 April 2012 to 31 March 2013 using Eddypro post-processing software (version 5.0, LI-COR, USA).CO 2 and energy fluxes are calculated at a standard averaging interval of 30 min using the fast response measurements from the sonic anemometer (3D wind speed components and sonic temperature) and IRGA (concentrations of CO 2 and H 2 O).
Every 30 min, raw data is corrected for standard flux corrections for the installed open path system to calculate CO 2 flux (Fc), termed as the mean covariance between deviations in instantaneous vertical wind speed (w') and density of CO 2 in air (c') [11].
The standard flux corrections include the removal of bad data in fast response measurements by de-spiking, de-trending and correcting for time delays [27].The tilt corrections for sonic measurements are made by axis rotation such that mean vertical wind (w) is null [28].Frequency response corrections are applied to compensate for flux losses at different frequencies (eddy sizes) [29,30].The density fluctuations of CO 2 and H 2 O due to temperature (thermal expansion) and water vapor (dilution) were compensated by Webb-Pearman-Leuing (WPL) Corrections [31].The EC system measures the fastest and smallest eddies accurately when the terrain is flat and homogenous.

Secondary Data Processing
The half-hourly computed flux data after standard corrections have spikes due to unrealistic meteorological measurements.These erroneous values affect the annual estimates and thus it is important to further correct the computed fluxes by removal of spurious data and to fill the gaps with more reasonable estimates.Further, the data has gaps due to shut-down of the tower because of the limitations of the solar power supply.Over annual measurements from April 2012 to March 2013, we had an actual data gap of about 11% of the total data due to low power and instrumental errors.
The secondary data processing includes the removal of CO 2 spikes, Latent Heat (LE) and Sensible Heat (H) spikes and corrections for nighttime flux based on statistical analysis.The corrected data is then gap-filled and partitioned into gross primary productivity (GPP) and ecosystem respiration (Re).The secondary data processing methodology used is similar to the proposed corrections by Thomas et al. [32].

CO 2 Flux Spike Removal
The computed flux data may have unrealistic spikes though the raw data are filtered using standard corrections.It is assumed that the extreme spikes are due to non-localized advection events or other non-stationary phenomena and not because of the natural phenomena [32].
The spikes were removed using statistical bounds on the data using a two-stage approach.Firstly, the data was grouped into positive and negative fluxes.Then half-hourly positive flux values which are greater than the respective month average plus two times the standard deviation (x > x + 2 ˆσ) are removed.The opposite was used for removal of absurd data in the negative fluxes (x < x ´2 ˆσ).In the second stage, positive (or negative) flux values greater (or lesser) than the respective monthly diurnal average plus (or minus) two times the standard deviation were removed.After the first stage about 13% of data was removed and an additional 2% was removed after the second stage.

LE and H Flux Spike Removal
The spikes in the energy fluxes are removed simply by using thresholds.Erroneous values of latent heat flux and sensible heat fluxes are removed from the dataset using similar thresholding as that described in Section 3.2.1.These may have resulted from strong advection or sensor interference which affects the reliability of associated CO 2 and energy fluxes, which were subsequently removed.Additional data of about 1% was removed after thresholding.

Removal of Negative Nighttime CO 2 Fluxes
This includes the removal of negative nighttime fluxes in the primary processed flux data post-spike removal.Identification of nighttime fluxes was made using PAR (Photosynthetic Active Radiation) less than 20 µmol¨m ´2¨s ´1.During the night, there is no chance of net CO 2 assimilation by the forest, thus these fluxes are related to the advection events which are likely at this site under stable nighttime conditions.Removal of such data caused a data gap of around 23%.There is no universal method to account for advection [33] and our method is in agreement with the several authors recording such errors using open path systems [26,32,34].

Friction Velocity Correction
Though the eddy covariance method is perhaps the best and most direct method for measuring mass and energy exchanges for terrestrial ecosystems, during the night it becomes complicated and uncertain because of the low turbulence and relatively stable air stratification.The standard method to process nighttime CO 2 flux data is to primarily assess the quality of the nighttime flux data using the friction velocity threshold, then to establish statistical models using reliable nighttime CO 2 flux data with meteorological data.Finally, we correct the incorrect data using the established model [35].
Since low friction velocity (u*) corresponds to calm nights, the u* threshold is used to identify the low turbulent flux data and the reliable flux data as described by Reichstein et al. [36].This method is the statistically based, site-independent determination of the lower u* threshold.The CO 2 flux data is divided into six temperature classes with equal sample size.Each temperature class is subdivided into 20 u* classes.The u* threshold for each temperature class was defined as the u* class where the nighttime flux reached more than 95% of the average flux within the higher u* classes.In order to account for seasonal variation in vegetation, the final u* threshold was calculated in three-month subsets as the median of the six temperature classes.The obtained u* thresholds for three-month subset data were 0.1, 0.13, 0.12 and 0.08, respectively.If the u* threshold was unable to be determined, a default value of 0.2 m¨s ´1 was used.The low turbulent flux data were subsequently removed, which amounts to 12% of the total data.

Gap-Filling and Uncertainty Analysis
In order to form a continuous dataset for annual carbon budgeting, the deleted data is replaced using gap-filled values.Several gap-filling methods have been developed based on interpolation, probabilistic filling, look-up tables, non-linear regression and artificial neural networks [37].In the present study, a total of 35% of data is deleted after spike removal of energy fluxes and CO 2 fluxes (including the deletion of negative nighttime values and low turbulence fluxes).To form a continuous dataset, we used standard procedures such as the marginal data sampling method (MDS) and mean diurnal variation (MDV) to substitute for the deleted values [36,38].Both MDS and MDV methods have proven to be unbiased approaches for gap-filling in comparison with other techniques [37].
The error and bias in fluxes caused by gap-filling methods (MDS and MDV) were estimated using the random creation of artificial gaps and then refilling these gaps using the valid data periods [37].During the study period, around 150 half-hourly fluxes were randomly removed to create an artificial gap and then filled with values to compare with the original values.The root mean square error (RMSE) and the bias error (BE) were determined using the flux estimations from 100 simulations.The computed RMSE of CO 2 flux 1.81 ˘0.24 µmol¨m ´2¨s ´1 was in the reported range of gap-filling errors in six forested types in Europe [37].The low BE (0.13 ˘0.1 µmol¨m ´2¨s ´1) value suggests the minimal bias in the annual NEP estimates.Similarly, the computed RMSE of H 2 O flux was also found to be 0.73 ˘0.12 µmol¨m ´2¨s ´1.

Flux Partitioning
The CO 2 flux data after u* threshold correction were partitioned into gross primary productivity (GPP) and ecosystem respiration.The reliable nighttime data (i.e., after removal of low turbulent flux data) is regressed with air temperature using the Lloyd-Taylor regression model [39].The air temperature is used as a proxy for ecosystem respiration to estimate daily respiration rates over a 14-day period.GPP is subsequently calculated as the difference between NEE and respiration.During nights, GPP cannot be negative or non-zero; such values are treated as errors and are forced to zero.During such cases, the ecosystem respiration was estimated to be the same as NEE due to the uncertainty in the regression model.For the secondary corrections, gap-filling and flux partitioning procedures, we used the REddyproc package [40] in an open source R environment along with MATLAB scripts to form the continuous datasets.REddyproc implements the same standardized methods as those described by Reichstein et al. [36] and is currently employed by CarboEurope.

Meteorological Data
Figure 3 shows the annual daily means of temperature, relative humidity, soil temperature, and solar radiation.The climate of the study area is hot and humid during most of the year.The maximum temperature of 32 ˝C occurred during June 2012 and the minimum temperature of 13 ˝C occurred during December 2012.The relative humidity was recorded to be greater than 70% for most of the year (Figure 3).Lower humidity values (nearly 50%) were recorded during winters (December 2012 to February 2013).Annual variations of temperature imply that May and June were hotter months with an average daily temperature of 30 ˝C whereas December and January were the colder months with an average daily temperature of 18 ˝C.The total precipitation at the study site was recorded to be 1267 mm during the study period, with nearly 80% of the total during monsoons (July to September). Figure 4 shows the monthly precipitation at the tower site.procedures, we used the REddyproc package [40] in an open source R environment along with MATLAB scripts to form the continuous datasets.REddyproc implements the same standardized methods as those described by Reichstein et al. [36] and is currently employed by CarboEurope.

Meteorological Data
Figure 3 shows the annual daily means of temperature, relative humidity, soil temperature, and solar radiation.The climate of the study area is hot and humid during most of the year.The maximum temperature of 32 °C occurred during June 2012 and the minimum temperature of 13 °C occurred during December 2012.The relative humidity was recorded to be greater than 70% for most of the year (Figure 3).Lower humidity values (nearly 50%) were recorded during winters (December 2012 to February 2013).Annual variations of temperature imply that May and June were hotter months with an average daily temperature of 30 °C whereas December and January were the colder months with an average daily temperature of 18 °C.The total precipitation at the study site was recorded to be 1267 mm during the study period, with nearly 80% of the total during monsoons (July to September). Figure 4 shows the monthly precipitation at the tower site.

Variations in CO2 Flux and H2O Flux
The footprint (or fetch) which is attributed to 90% of the total turbulent flux computations was estimated using the Kljun et al., model [41].The average footprint during the day was estimated to be about 160 m with a standard deviation of 22 m during the study period (Figure 5).

Variations in CO 2 Flux and H 2 O Flux
The footprint (or fetch) which is attributed to 90% of the total turbulent flux computations was estimated using the Kljun et al., model [41].The average footprint during the day was estimated to be about 160 m with a standard deviation of 22 m during the study period (Figure 5). Figure 6

Variations in CO2 Flux and H2O Flux
The footprint (or fetch) which is attributed to 90% of the total turbulent flux computations was estimated using the Kljun et al., model [41].The average footprint during the day was estimated to be about 160 m with a standard deviation of 22 m during the study period (Figure 5).Water vapor fluxes (or evapotranspiration (ET) rates), defined as the sum of evaporation and transpiration over an ecosystem, varied between a maximum of 5.5 mmol•m −2 •s −1 in October 2012 to a minimum of 2.5 mmol•m −2 •s −1 in July 2013 (Figure 8).The majority of H2O fluxes were found to be positive throughout the year with a steady increase during the daytime (photosynthesis and high ET rates), peaking during mid-day (1200 to 1300 h) and decreasing to the lowest during the night (no photosynthesis, lowest ET rates).The mean annual ET during the study period was estimated to be 1.96 ± 0.33 mm•day −1 .
The maximum daytime NEE obtained from the Florida Everglades mangroves is about −20 to −25 μmol•m −2 •s −1 [17].The NEE values obtained are only one-third of the NEE observed at the Florida Everglades mangroves.The tree heights in the mangroves of the Florida Everglades are 15-20 m [17] which, compared to the mean height of 5 m in the mangroves of Sundarbans, is the possible explanation of the huge variation in NEE.Water vapor fluxes (or evapotranspiration (ET) rates), defined as the sum of evaporation and transpiration over an ecosystem, varied between a maximum of 5.5 mmol•m −2 •s −1 in October 2012 to a minimum of 2.5 mmol•m −2 •s −1 in July 2013 (Figure 8).The majority of H2O fluxes were found to be positive throughout the year with a steady increase during the daytime (photosynthesis and high ET rates), peaking during mid-day (1200 to 1300 h) and decreasing to the lowest during the night (no photosynthesis, lowest ET rates).The mean annual ET during the study period was estimated to be 1.96 ± 0.33 mm•day −1 .
The maximum daytime NEE obtained from the Florida Everglades mangroves is about −20 to −25 μmol•m −2 •s −1 [17].The NEE values obtained are only one-third of the NEE observed at the Florida Everglades mangroves.The tree heights in the mangroves of the Florida Everglades are 15-20 m [17] which, compared to the mean height of 5 m in the mangroves of Sundarbans, is the possible explanation of the huge variation in NEE.Water vapor fluxes (or evapotranspiration (ET) rates), defined as the sum of evaporation and transpiration over an ecosystem, varied between a maximum of 5.5 mmol¨m ´2¨s ´1 in October 2012 to a minimum of 2.5 mmol¨m ´2¨s ´1 in July 2013 (Figure 8).The majority of H 2 O fluxes were found to be positive throughout the year with a steady increase during the daytime (photosynthesis and high ET rates), peaking during mid-day (1200 to 1300 h) and decreasing to the lowest during the night (no photosynthesis, lowest ET rates).The mean annual ET during the study period was estimated to be 1.96 ˘0.33 mm¨day ´1.
The maximum daytime NEE obtained from the Florida Everglades mangroves is about ´20 to ´25 µmol¨m ´2¨s ´1 [17].The NEE values obtained are only one-third of the NEE observed at the Florida Everglades mangroves.The tree heights in the mangroves of the Florida Everglades are 15-20 m [17] which, compared to the mean height of 5 m in the mangroves of Sundarbans, is the possible explanation of the huge variation in NEE.

Heat Fluxes and Energy Balance
The maximum half-hourly measured net solar radiation (Rn) was nearly 700 W•m  9).During July 2012 to November 2012, energy partitioning into latent heat flux was slightly higher than the sensible heat flux.This is largely due to the availability of water due to rain events, which triggered evapotranspiration at the study site.Closure of the surface energy balance provides insight into the performance of the eddy covariance measurement system [11,42].For each half-hour period, the total net radiation (Rn) should be approximately equal to the sum of the energy distributed between the LE, H, and soil heat flux (G) and canopy heat storage (S).Since the canopy heat storage (S) is expected to be small in short canopies with minimal biomass [42], the surface energy balance at the surface can be expressed as:

Heat Fluxes and Energy Balance
The maximum half-hourly measured net solar radiation (R n ) was nearly 700 W¨m ´2 during the summer season (April 2012 and March 2013) and 500 W¨m ´2 during the winters (November 2012 to January 2013).The partitioning of the net available energy into sensible heat flux (H) and latent heat flux (LE) was almost the same throughout the year.During the study period from April 2012 to March 2013, both H and LE showed a linear relationship with net radiation with values peaking during midday.
The monthly diurnal average values of sensible heat flux (maximum of 240 W¨m ´2 during April 2012) and latent heat flux (maximum of 240 W¨m ´2 during October 2012) depicts that the net radiation was equally divided into both energy fluxes for most months of the year (Figure 9).During July 2012 to November 2012, energy partitioning into latent heat flux was slightly higher than the sensible heat flux.This is largely due to the availability of water due to rain events, which triggered evapotranspiration at the study site.

Heat Fluxes and Energy Balance
The maximum half-hourly measured net solar radiation (Rn) was nearly 700 W•m −2 during the summer season (April 2012 and March 2013) and 500 W•m −2 during the winters (November 2012 to January 2013).The partitioning of the net available energy into sensible heat flux (H) and latent heat flux (LE) was almost the same throughout the year.During the study period from April 2012 to March 2013, both H and LE showed a linear relationship with net radiation with values peaking during midday.
The monthly diurnal average values of sensible heat flux (maximum of 240 W•m −2 during April 2012) and latent heat flux (maximum of 240 W•m −2 during October 2012) depicts that the net radiation was equally divided into both energy fluxes for most months of the year (Figure 9).During July 2012 to November 2012, energy partitioning into latent heat flux was slightly higher than the sensible heat flux.This is largely due to the availability of water due to rain events, which triggered evapotranspiration at the study site.Closure of the surface energy balance provides insight into the performance of the eddy covariance measurement system [11,42].For each half-hour period, the total net radiation (Rn) should be approximately equal to the sum of the energy distributed between the LE, H, and soil heat flux (G) and canopy heat storage (S).Since the canopy heat storage (S) is expected to be small in short canopies with minimal biomass [42], the surface energy balance at the surface can be expressed as: Closure of the surface energy balance provides insight into the performance of the eddy covariance measurement system [11,42].For each half-hour period, the total net radiation (R n ) should be approximately equal to the sum of the energy distributed between the LE, H, and soil heat flux (G) and canopy heat storage (S).Since the canopy heat storage (S) is expected to be small in short canopies with minimal biomass [42], the surface energy balance at the surface can be expressed as: R n ´G " H `LE Figure 10 shows the energy closure of the Sundarbans study area during the day from April 2012 to March 2013 with measured half-hourly flux values.Figure 10 also depicts that the slope of the regression (m = 0.78) was less than one and the intercept (c = 18.7) was slightly greater than zero with an R 2 = 0.74.The surface energy balance of about 78% of the available energy is accounted for by the 30-min fluxes in the present mangrove ecosystem.Although the energy balance closure is not perfect, it is the typical range reported at other flux sites.Wilson et al., 2002, used 22 FLUXNET sites of different vegetation heights and compositions and reported the ranges of slope, intercept and R 2 to be in between 0.53 and 0.99, ´32.9 and 36.9W¨m ´2, and 0.64 and 0.36, respectively.Other long-term eddy covariance studies also indicate that the measurements do not show a closure of surface energy balance, although this discrepancy is not well understood [42][43][44].It is also considered that the imbalance may be due to the sensor accuracy, data processing and under-estimations of the advections due to homogenous terrain.Figure 10 also depicts that the slope of the regression (m = 0.78) was less than one and the intercept (c = 18.7) was slightly greater than zero with an R 2 = 0.74.The surface energy balance of about 78% of the available energy is accounted for by the 30-min fluxes in the present mangrove ecosystem.Although the energy balance closure is not perfect, it is the typical range reported at other flux sites.Wilson et al., 2002, used 22 FLUXNET sites of different vegetation heights and compositions and reported the ranges of slope, intercept and R 2 to be in between 0.53 and 0.99, −32.9 and 36.9W•m −2 , and 0.64 and 0.36, respectively.Other long-term eddy covariance studies also indicate that the measurements do not show a closure of surface energy balance, although this discrepancy is not well understood [42][43][44].It is also considered that the imbalance may be due to the sensor accuracy, data processing and under-estimations of the advections due to homogenous terrain.Over one year (April 2012 to March 2013), the estimated NEP was 249 ± 20 g C m −2 •year −1 .The estimated GPP and Re were 1271 g C m −2 •year −1 and 1022 g C m −2 •year −1 , respectively (Table 1).The estimated NEP values (249 g C m −2 •year −1 ) during April 2012 to March 2013 fall in the reported range for terrestrial ecosystems [45][46][47].In comparison with the mangroves of the Florida Everglades (NEP = 1170 ± 145 g C m −2 •year −1 ), the estimated NEP of Sundarbans is quite low and is attributed to the fact that the Everglades mangroves are dense with trees of greater heights, have a high leaf litter rate and also have a greater wood production [17].Also, the species composition which contributes to biomass and, in turn, NEP is different in the wet evergreen system of Sundarbans (predominant species viz.Avicennia alba, Bruguiera gymnorrhiza and Rhizophora) in comparison with the mangroves of the Everglades (predominant species viz.Rhizophora mangle, Avicennia germinans and Laguncularia racemosa).The presence of significantly less biomass (10-50 t•ha −1 ) in the Sundarbans mangrove forest [48] in comparison with the high biomass of around 150-200 t•ha −1 in the Everglades mangroves [49] explains the huge difference in the NEP.Over one year (April 2012 to March 2013), the estimated NEP was 249 ˘20 g C m ´2¨year ´1.The estimated GPP and R e were 1271 g C m ´2¨year ´1 and 1022 g C m ´2¨year ´1, respectively (Table 1).The estimated NEP values (249 g C m ´2¨year ´1) during April 2012 to March 2013 fall in the reported range for terrestrial ecosystems [45][46][47].In comparison with the mangroves of the Florida Everglades (NEP = 1170 ˘145 g C m ´2¨year ´1), the estimated NEP of Sundarbans is quite low and is attributed to the fact that the Everglades mangroves are dense with trees of greater heights, have a high leaf litter rate and also have a greater wood production [17].Also, the species composition which contributes to biomass and, in turn, NEP is different in the wet evergreen system of Sundarbans (predominant species viz.Avicennia alba, Bruguiera gymnorrhiza and Rhizophora) in comparison with the mangroves of the Everglades (predominant species viz.Rhizophora mangle, Avicennia germinans and Laguncularia racemosa).The presence of significantly less biomass (10-50 t¨ha ´1) in the Sundarbans mangrove forest [48] in comparison with the high biomass of around 150-200 t¨ha ´1 in the Everglades mangroves [49] explains the huge difference in the NEP.The measured respiration rates (2 to 4 µmol¨m ´2¨s ´1) were low throughout the year in comparison with other forest types.It was reported by the measurements across terrestrial AmeriFlux and EuroFlux sites that the respiration rates range from 3.72 ˘2.20 to 5.92 ˘4.40 µmol¨m ´2¨s ´1 at different soil temperatures of 15 ˝C and 20 ˝C ,respectively [38].In the study area, though high soil temperatures (20 ˝C to 30 ˝C) were recorded, the observed respiration rates were less than half of the globally reported values.Due to the consistent high water level around the tower site, the soils are generally anoxic and would have contributed to low respiration rates.

NEE and Photosynthetic Active Radiation (PAR)
Since the photosynthetic activity is limited to the availability of photons rather than other limiting factors, daily variations in GPP were also correlated with daily averaged photosynthetic active radiation (PAR).GPP shows a positive non-linear relation with daily averaged PAR, indicating that GPP increases with PAR up to a certain value and saturates further (Figure 12).
Though Sundarbans mangroves exhibit less phenological changes, a relative increase in NEE sequestration was observed during September 2012 to December 2012 compared to other months (Figure 6).Accordingly, relationships of monthly diurnal averaged NEE with monthly diurnal averaged PAR were analyzed in two phases.During September 2012 to December 2012, NEE increased incrementally with PAR until reaching 900 µmol¨m ´2¨s ´1, thereafter leveling off at ´8 µmol¨m ´2¨s ´1 incrementally with PAR (Figure 13).During the remaining period, NEE showed a similar trend in response with PAR, but saturated at lower levels of sequestration around 5.5 µmol¨m ´2¨s ´1 (Figure 13).The measured respiration rates (2 to 4 μmol•m −2 •s −1 ) were low throughout the year in comparison with other forest types.It was reported by the measurements across terrestrial AmeriFlux and EuroFlux sites that the respiration rates range from 3.72 ± 2.20 to 5.92 ± 4.40 μmol•m −2 •s −1 at different soil temperatures of 15 °C and 20 °C ,respectively [38].In the study area, though high soil temperatures (20 °C to 30 °C) were recorded, the observed respiration rates were less than half of the globally reported values.Due to the consistent high water level around the tower site, the soils are generally anoxic and would have contributed to low respiration rates.

NEE and Photosynthetic Active Radiation (PAR)
Since the photosynthetic activity is limited to the availability of photons rather than other limiting factors, daily variations in GPP were also correlated with daily averaged photosynthetic active radiation (PAR).GPP shows a positive non-linear relation with daily averaged PAR, indicating that GPP increases with PAR up to a certain value and saturates further (Figure 12).
Though Sundarbans mangroves exhibit less phenological changes, a relative increase in NEE sequestration was observed during September 2012 to December 2012 compared to other months (Figure 6).Accordingly, relationships of monthly diurnal averaged NEE with monthly diurnal averaged PAR were analyzed in two phases.During September 2012 to December 2012, NEE increased incrementally with PAR until reaching 900 μmol•m −2 •s −1 , thereafter leveling off at −8 μmol•m −2 •s −1 incrementally with PAR (Figure 13).During the remaining period, NEE showed a similar trend in response with PAR, but saturated at lower levels of sequestration around 5.5 μmol•m −2 •s −1 (Figure 13).

Conclusions
This paper presents the first annual assessment of the CO 2 flux of the Sundarbans mangroves ecosystems, India, using the eddy covariance method, and it aids in providing more direct defensible estimates of the carbon sequestration over the Indian sub-continent.The present eddy flux tower is one of the few towers measuring CO 2 fluxes in mangroves across the globe.Eddy covariance-based NEE estimates in the present study area are more reliable in comparison with previously reported estimates due to their continuous measurements throughout the year.The Sundarbans mangroves had proven to be a consistent carbon sink over different months with an annual NEP of 249 ˘20 g C m ´2¨year ´1 during April 2012 to March 2013.A seasonal pattern of CO 2 flux has been observed over the study period, indicating higher NEE during October to December 2012 and lower NEE during April to June 2012.Surface energy balance closure, which is a key indicator to the performance of the eddy covariance measurement system [11], is estimated to be 78% in the present study and is in agreement with the reported values [42].The small bias error estimate obtained after gap-filling the missing NEE values suggests the reliable annual estimates.Our eddy covariance estimates of NEP, GPP and Re were in the range of the globally reported values from different terrestrial ecosystems.PAR was found to be a strong driver for CO 2 sequestration with a maximum quantum efficiency of 0.0398 ˘0.0038 (95% confidence) mols CO 2 per mol of absorbed photons in September 2012.It was observed that the partitioning of net radiation into latent and sensible hear flux was uniform throughout the majority of the year, whereas during July to November 2014 the latent heat flux was higher than the sensible heat flux due to the availability of water during the monsoon and post-monsoon months.The total evapotranspiration was estimated to be 715 ˘120 mm¨year ´1 over the study period.
The present study reports a conservative yet reliable NEP estimate of 249 g C m ´2¨year ´1 purely based on half-hourly computed CO 2 fluxes within the fetch of EC instruments and extrapolation of the same over the entire study area with the assumption of uniform homogeneous vegetation.The study does not, however, account for other important elements of aquatic pathways within which factors such as tidal and groundwater CO 2 export and dissolved inorganic and organic carbon (DIC) exports can influence the NEP to a considerable degree ( [17,53]).Information on such holistic factors which drive the carbon fluxes can add value to our current estimates, which the authors envisage in future studies.Long-term research and measurements would help in building continuous and consistent eddy covariance data sets, which will help to understand the ecosystem dynamics (such as seasonal patterns), environment and climate change.Also, the understanding of inter-relationships with bio-physical variables in future studies helps in up-scaling the results the regional scales.

Figure 1 .
Figure 1.Sundarbans flux tower location on the false color composite (FCC) of IRS AWIFS.

Figure 1 .
Figure 1.Sundarbans flux tower location on the false color composite (FCC) of IRS AWIFS.

Figure 3 .
Figure 3. Daily averages of meteorological variables over Sundarbans from April 2012 to March 2013.Bold lines represent the averages and the shaded area represents the min-max range.For net radiation, only average value is depicted.

Figure 3 .
Figure 3. Daily averages of meteorological variables over Sundarbans from April 2012 to March 2013.Bold lines represent the averages and the shaded area represents the min-max range.For net radiation, only average value is depicted.
Figure 6 depicts the variations of monthly diurnal averaged carbon dioxide flux over a year (April 2012 to March 2013).Even though mangroves are less prone to phenological changes, there is significant seasonal variation in CO2 flux when observed annually.Variations in ecosystem respiration during nights (positive CO2 flux) were minimal with values ranging between 2 to 4 μmol•m −2 •s −1 throughout the year, whereas CO2 flux due to photosynthesis during the day (negative CO2 flux) peaked at −6 μmol•m −2 •s −1 during the summer season (April 2012 to June 2012) and decreased to −10 μmol•m −2 •s −1 during the winter period (October 2012 to December 2012).The half-hourly maximum daytime NEE of −11.8 μmol•m −2 •s −1 was observed during October 2012.The diurnal and seasonal variation of CO2 fluxes over a year are depicted in Figure 7, which shows the gap-filled NEE (in μmol•m −2 •s −1 ) with respect to time of day and month.

Figure 5 .Figure 4 .
Figure 5. Sundarbans flux tower foot print of 90% of the flux contribution (represented by X 90%) over the study period ranged from 90 m to 300 m.
depicts the variations of monthly diurnal averaged carbon dioxide flux over a year (April 2012 to March 2013).Even though mangroves are less prone to phenological changes, there is significant seasonal variation in CO 2 flux when observed annually.Variations in ecosystem respiration during nights (positive CO 2 flux) were minimal with values ranging between 2 to 4 µmol¨m ´2¨s ´1 throughout the year, whereas CO 2 flux due to photosynthesis during the day (negative CO 2 flux) peaked at ´6 µmol¨m ´2¨s ´1 during the summer season (April 2012 to June 2012) and decreased to ´10 µmol¨m ´2¨s ´1 during the winter period (October 2012 to December 2012).The half-hourly maximum daytime NEE of ´11.8 µmol¨m ´2¨s ´1 was observed during October 2012.The diurnal and seasonal variation of CO 2 fluxes over a year are depicted in Figure 7, which shows the gap-filled NEE (in µmol¨m ´2¨s ´1) with respect to time of day and month.

Figure 4 .
Figure 4. Total monthly precipitation during April 2012 to March 2013.
Figure 6 depicts the variations of monthly diurnal averaged carbon dioxide flux over a year (April 2012 to March 2013).Even though mangroves are less prone to phenological changes, there is significant seasonal variation in CO2 flux when observed annually.Variations in ecosystem respiration during nights (positive CO2 flux) were minimal with values ranging between 2 to 4 μmol•m −2 •s −1 throughout the year, whereas CO2 flux due to photosynthesis during the day (negative CO2 flux) peaked at −6 μmol•m −2 •s −1 during the summer season (April 2012 to June 2012) and decreased to −10 μmol•m −2 •s −1 during the winter period (October 2012 to December 2012).The half-hourly maximum daytime NEE of −11.8 μmol•m −2 •s −1 was observed during October 2012.The diurnal and seasonal variation of CO2 fluxes over a year are depicted in Figure 7, which shows the gap-filled NEE (in μmol•m −2 •s −1 ) with respect to time of day and month.

Figure 5 .Figure 5 .Figure 6 .
Figure 5. Sundarbans flux tower foot print of 90% of the flux contribution (represented by X 90%) over the study period ranged from 90 m to 300 m.

Figure 7 .
Figure 7. Gap-filled net ecosystem exchange (NEE) fingerprint of Sundarbans mangroves showing time of the day (x-axis) and month-wise (y-axis) CO2 fluxes in μmol•m −2 •s −1 during the study period.

Figure 6 .Forests 2016, 7 , 35 9 of 18 Figure 6 .
Figure 6.Monthly mean diurnal variation of half-hourly CO 2 flux during April 2012 to March 2013.The error bars show one standard deviation of the diurnal average fluxes during respective months.

Figure 7 .
Figure 7. Gap-filled net ecosystem exchange (NEE) fingerprint of Sundarbans mangroves showing time of the day (x-axis) and month-wise (y-axis) CO2 fluxes in μmol•m −2 •s −1 during the study period.

Figure 7 .
Figure 7. Gap-filled net ecosystem exchange (NEE) fingerprint of Sundarbans mangroves showing time of the day (x-axis) and month-wise (y-axis) CO 2 fluxes in µmol¨m ´2¨s ´1 during the study period.

Figure 8 .
Figure 8. Monthly mean diurnal variations of half-hourly H2O flux during April 2012 to March 2013.The error bars show one standard deviation of the diurnal average fluxes during the respective months.
−2 during the summer season (April 2012 and March 2013) and 500 W•m −2 during the winters (November 2012 to January 2013).The partitioning of the net available energy into sensible heat flux (H) and latent heat flux (LE) was almost the same throughout the year.During the study period from April 2012 to March 2013, both H and LE showed a linear relationship with net radiation with values peaking during midday.The monthly diurnal average values of sensible heat flux (maximum of 240 W•m −2 during April 2012) and latent heat flux (maximum of 240 W•m −2 during October 2012) depicts that the net radiation was equally divided into both energy fluxes for most months of the year (Figure

Figure 9 .
Figure 9. Annual variation of monthly mean diurnal variation of sensible heat and latent heat fluxes.

Figure 8 .
Figure 8. Monthly mean diurnal variations of half-hourly H 2 O flux during April 2012 to March 2013.The error bars show one standard deviation of the diurnal average fluxes during the respective months.

Figure 8 .
Figure 8. Monthly mean diurnal variations of half-hourly H2O flux during April 2012 to March 2013.The error bars show one standard deviation of the diurnal average fluxes during the respective months.

Figure 9 .
Figure 9. Annual variation of monthly mean diurnal variation of sensible heat and latent heat fluxes.

Figure 9 .
Figure 9. Annual variation of monthly mean diurnal variation of sensible heat and latent heat fluxes.

ForestsFigure 10
Figure10shows the energy closure of the Sundarbans study area during the day from April 2012 to March 2013 with measured half-hourly flux values.Figure10also depicts that the slope of the regression (m = 0.78) was less than one and the intercept (c = 18.7) was slightly greater than zero with an R 2 = 0.74.The surface energy balance of about 78% of the available energy is accounted for by the 30-min fluxes in the present mangrove ecosystem.Although the energy balance closure is not perfect, it is the typical range reported at other flux sites.Wilson et al., 2002, used 22  FLUXNET sites of different vegetation heights and compositions and reported the ranges of slope, intercept and R 2 to be in between 0.53 and 0.99, −32.9 and 36.9W•m −2 , and 0.64 and 0.36, respectively.Other long-term eddy covariance studies also indicate that the measurements do not show a closure of surface energy balance, although this discrepancy is not well understood[42][43][44].It is also considered that the imbalance may be due to the sensor accuracy, data processing and under-estimations of the advections due to homogenous terrain.

Figure 10 .
Figure 10.Surface energy closure in daytime in the Sundarbans mangroves during April 2012 to March 2013 with half-hourly fluxes.Soil heat flux includes energy storage between the surface and the 5 cm soil.

Figure 11 shows
Figure 11 shows the daily 30-min average of the NEP partitioned into GPP and respiration (Re).The bold lines represent the weekly average for the partitioned fluxes.Positive fluxes of NEP and GPP represent net sequestration (sink) by the forest ecosystem whereas negative values indicate a net loss (source) of CO2 to the atmosphere.Respiration represents the net release of CO2 from the ecosystem to the atmosphere.During the study period, daily GPP estimates ranged from 3 to 5.5 g C m −2 •day −1 during October and November 2012, and between 2.5 to 4 g C m −2 • day −1 during January and February 2013.Greater fluctuations (1 to 4.5 g C m −2 • day −1 ) were recorded during monsoon periods (July and August 2012), which may be attributed to less availability of photosynthetic radiation during cloudy/rainy days.The annual variation of the day-wise ecosystem respiration was also similar to GPP, with maximum values of around 5.5 g C m −2 • day −1 during October and November 2012 and minimum values of 2 g C m −2 •d −1 during January and February 2013.Figure 11 also indicates that the daily NEPs are positive throughout the majority of days during the study period, suggesting that the ecosystem is a sink for CO2.During September to December 2012, the values of NEP were measured to be nearly 2.5 g C m −2 • day −1 and they decreased to around 0.5 g C m −2 • day −1 during January to March 2013.
Figure 11 shows the daily 30-min average of the NEP partitioned into GPP and respiration (Re).The bold lines represent the weekly average for the partitioned fluxes.Positive fluxes of NEP and GPP represent net sequestration (sink) by the forest ecosystem whereas negative values indicate a net loss (source) of CO2 to the atmosphere.Respiration represents the net release of CO2 from the ecosystem to the atmosphere.During the study period, daily GPP estimates ranged from 3 to 5.5 g C m −2 •day −1 during October and November 2012, and between 2.5 to 4 g C m −2 • day −1 during January and February 2013.Greater fluctuations (1 to 4.5 g C m −2 • day −1 ) were recorded during monsoon periods (July and August 2012), which may be attributed to less availability of photosynthetic radiation during cloudy/rainy days.The annual variation of the day-wise ecosystem respiration was also similar to GPP, with maximum values of around 5.5 g C m −2 • day −1 during October and November 2012 and minimum values of 2 g C m −2 •d −1 during January and February 2013.Figure 11 also indicates that the daily NEPs are positive throughout the majority of days during the study period, suggesting that the ecosystem is a sink for CO2.During September to December 2012, the values of NEP were measured to be nearly 2.5 g C m −2 • day −1 and they decreased to around 0.5 g C m −2 • day −1 during January to March 2013.

Figure 10 .
Figure 10.Surface energy closure in daytime in the Sundarbans mangroves during April 2012 to March 2013 with half-hourly fluxes.Soil heat flux includes energy storage between the surface and the 5 cm soil.

18 Figure 11 .
Figure11shows the daily 30-min average of the NEP partitioned into GPP and respiration (Re).The bold lines represent the weekly average for the partitioned fluxes.Positive fluxes of NEP and GPP represent net sequestration (sink) by the forest ecosystem whereas negative values indicate a net loss (source) of CO 2 to the atmosphere.Respiration represents the net release of CO 2 from the ecosystem to the atmosphere.During the study period, daily GPP estimates ranged from 3 to 5.5 g C m ´2¨day ´1 during October and November 2012, and between 2.5 to 4 g C m ´2¨day ´1 during January and February 2013.Greater fluctuations (1 to 4.5 g C m ´2¨day ´1) were recorded during monsoon periods (July and August 2012), which may be attributed to less availability of photosynthetic radiation during cloudy/rainy days.The annual variation of the day-wise ecosystem respiration was also similar to GPP, with maximum values of around 5.5 g C m ´2¨day ´1 during October and November 2012 and minimum values of 2 g C m ´2¨d ´1 during January and February 2013.Figure11also indicates that the daily NEPs are positive throughout the majority of days during the study period, suggesting that the ecosystem

Figure 11 .
Figure 11.Daily average of partitioned flux into gross primary productivity (GPP), respiration (Re) and net ecosystem productivity (NEP) over the study period.Bold lines represent seven-day running mean.

Table 1 .
Summary table of annual estimates of CO 2 fluxes and evapotranspiration rates (ET) over the study period.

Table 1 .
Summary table of annual estimates of CO2 fluxes and evapotranspiration rates (ET) over the study period.