Mangrove Methane Biogeochemistry in the Indian Sundarbans: A Proposed Budget

Biogeochemical cycling of CH4 was investigated at Lothian Island, one of the relatively pristine islands of Indian Sundarbans and its adjacent Saptamukhi estuary, during June 2010 to December 2012. Intertidal mangrove sediments were highly anoxic and rich in organic carbon. Mean rates of methanogenesis were 3547 and 48.88 µmol m-3 wet sediment d-1, for intertidal (up to 25 cm depth) and sub-tidal sediments (first 5 cm depth), respectively. CH4 in pore-water was 53.4 times more supersaturated than in adjacent estuarine waters. This resulted in significant CH4 efflux from sediments to estuarine waters - via advective and diffusive transport. About 8.2% of the total CH4 produced in intertidal mangrove sediments was transported to the adjacent estuary through advective flux, which was 20 times higher than diffusive CH4 flux. Mean CH4 concentrations in estuarine surface and sub-surface waters were 69.9 and 56.1 nM, respectively, with a dissolved CH4 oxidation rate in estuarine surface waters of 20.5 nmol L-1 d-1. An estimated 0.09 Gg yr-1 of CH4 is released from estuaries of Sundarbans to the regional atmosphere. The mean CH4 mixing ratio over the forest atmosphere was 2 ppmv. On annual basis, only 2.75% of total supplied CH4 to the forest atmosphere was transported to the upper atmosphere via biosphere-atmosphere exchange. Mean CH4 photo-oxidation rate over the forest atmosphere was 3.25 x 10-9 mg cm-3 d-1. Using new and and previously published data we present for the first time, a CH4 budget for Sundarbans mangrove ecosystem which in part, revealed the existence of anaerobic CH4 oxidation in the mangrove sediment column.


INTRODUCTION
Methane (CH 4 ) is the key gas produced in anaerobic environments and represents the second most abundant greenhouse gas associated with climate change (Forster et al., 2007). Moreover, about 1% of the annually-fixed CO 2 via photosynthesis is converted back to CO 2 through microbial methanogenesis-with an estimated 1 billion tones of CH 4 cycled this way each year (Rudolf and Seigo, 2006). The global atmospheric CH 4 mixing ratio has increased from 722 ppbv in the year 1,750 to 1,840 ppbv in 2016 (http://cdiac.ornl.gov/pns/current_ghg.html; http://www.esrl. noaa.gov/gmd/ccgg/trends_ch4/). The cause of this large increase in atmospheric CH 4 is not fully understood, but is probably related to a surge in CH 4 emissions from wetlands that contribute approximately 20-39% of the annual global atmospheric CH 4 budget (Hoehler and Alperin, 2014). Since CH 4 is 26 times more effective than CO 2 in its radiative forcing as a greenhouse gas (Lelieveld et al., 1993), we need to remain vigilent about the continually-evolving changes in the cycling of CH 4 during the Anthropocene (Bridgham et al., 2013).
Mangroves are one of the most productive coastal environments, characterized by high turnover rates of organic matter and as well as exchange with the adjacent marine system (e.g., Donato et al., 2012;Alongi, 2014). Organic matter mineralization in sediments is a multi-step process that depends on the availability of oxygen and presence of terminal electron acceptors (TEAs). At the terminal step of organic matter decomposition, when all the TEAs have been consumed and electron donors are in surplus, CH 4 is produced (requisite redox potential: <150 mV; Wang et al., 1993) by a fermentative disproportionation reaction of low-molecular weight compounds (e.g., acetate), or by reduction of CO 2 via hydrogen or simple alcohols (e.g., Canfield et al., 2005). In wetlands, sedimentary-derived CH 4 can escape to the adjacent water/atmosphere via diffusive evasion, ebullition, and plantmediated transport (Chanton and Dacey, 1991). Past work has shown that diffusion is the least effective pathway (Chanton and Dacey, 1991;Laanbroek, 2010;Bridgham et al., 2013), with plant-mediation being the most globally-important source to the atmopshere from shallow-water ecosystems-which are largely freshwater and with emergent rooted plants (Bastviken et al., 2011). In coastal wetlands, the role of plant-mediated transport is significantly reduced because methanogenesis is largely suppressed by sulfate reduction. However, tidal effects can introduce another pathway of transport whereby during low-tide conditions, CH 4 -rich pore water is transported to adjacent creeks and estuaries through hypsometric gradients (Deborde et al., 2010;Santos et al., 2012;Stieglitz et al., 2013).
In estuarine water columns, CH 4 can be partly oxidized to CO 2 by methanotrophs (e.g., Hanson and Hanson, 1996), which reduces CH 4 fluxes across water-atmosphere interface. In fact, for stratified lakes, methanotrophs can consume up to 90% of the dissolved CH 4 (Utsumi et al., 1998;Kankaala et al., 2006). However, in well-mixed estuaries, the activity of methanotrophs has been shown to be significantly less, allowing for the possible escape of CH 4 to the atmospshere (Abril et al., 2007). The amount of dissolved CH 4 that escapes microbial oxidation in estuaries depends in large part on the CH 4 concentration gradient at the air-water interface and the gas transfer velocity. The remaining CH 4 can be exported to adjacent continental shelf regions where it only plays a very minor role in the larger global coastal carbon budget (Bauer et al., 2013). Once CH 4 enters the atmosphere, across the sediment-atmosphere and water-atmosphere interfaces, there is enrichment of the atmospheric CH 4 mixing ratio at a regional level. At this point, an entirely different set of complex atmospheric chemical transformations will determine its fate. Any emitted CH 4 from mangroves will exchange across biosphere-atmosphere interface largely as a function of atmospheric turbulence. Another significant fraction of this emitted CH 4 will undergo photo-oxidation-depending upon ambient NO x and OH radical concentrations (Wayne, 1991). Thus, the role of photo-oxidation of CH 4 across the steep light gradients that exist in a mangrove forest canopy, which remain largely unknown, need to be considered when examining fluxes across the sediment-water-atmosphere continum. Finally, another interesting factor controlling the fate emitted CH 4 in freshwater forested systems such as Cypress Swamps, is the role of plant-mediated CH 4 efflux through cypress knees (composed of lenticels and aerenchyma tissue), where significant efflux of CH 4 has been obeserved across the sediment-air interface (Bianchi et al., 1996).
The main objective of this study was to measure methanogenesis (in intertidal and sub-tidal sediments) and methanotrophy (in intertidal sediment) in the Indian Sundarbans, the largest tidal mangrove forest in the world. This builds on our previous work (Dutta et al., , 2015a only with more focus on the potential influences of sediment physicochemical parameters on pore-water CH 4 distribution, the effects of pneumatophores and bioturbation density on the variability of sediment-atmosphere CH 4 fluxes, and estuarine physicochemical parameters on CH 4 oxidation. Perhaps more importantly, we present for the first time to our knowledge, a comprehensive quantitative CH 4 budget for this globally-important wetland region.

SAMPLING LOCATION
The Sundarbans-located across India and Bangladesh at the land-ocean boundary of Ganges-Brahmaputra delta and the Bay of Bengal, is the largest single tidal mangrove forest in the world. This extensive natural mangrove forest has been established as a UNESCO world heritage site and covers an area of 10,200 km 2 , of which 4,200 km 2 of reserved forest resides in India, with the remainder in Bangladesh. More specifically, the Indian Sundarbans Biosphere Reserve (SBR), which extends over an area of 9,600 km 2 , is comprised of 1,800 km 2 of estuarine waterways and 3,600 km 2 of reclaimed areas, in addition to the aforementioned mangrove reserve forest. The forest is about 140 km in length from east to west, and extends approximately 50-70 km from the southern margin of the Bay of Bengal toward the north.
The Indian Sundarban mangrove delta is fed by the sediment-rich waters of seven rivers, namely the Mooriganga, Saptamukhi, Thakuran, Matla, Bidya, Gosaba, and Haribhanga. These rivers form a sprawling archipelago of 102 islands, of which 54 have been reclaimed for human settlement with the others essentially uninhabited and relatively prisitine. Lothian Island, is one of these pristine islands situated at the buffer zone of the Sundarbans Biosphere Reserve-covering an area of 38 km 2 (Figure 1). The island is completely intertidal and is occupied by thick, robust, and resilient mangrove trees, with a mean height of ca. 10 m. The dominant species of mangroves are Avicennia alba, Avicennia marina, and Avicennia officinalis, with Excoecaria agallocha marginally distributed, and Ceriops decandra sparsley scattered across the island. The mangrove sediments are a silty clay largely composed of quartzo-feldspathic minerals such as quartz, albite, and microline (Ray et al., 2013). Lothian Island also borders the Saptamukhi River estuarine system (Figure 1), which has a mean depth ∼6 m (Dutta et al., 2015a) with no perennial source of freshwater, but receives significant amounts of agricultural and anthropogenic runoff-especially during monsoon season. Like entire Sundarbans, the regional climate of the study area is characterized by pre-monsoon (February-May), southwest monsoon (June-September), and northeast monsoon (or postmonsoon) (October-January) conditions. Although, Lothian Island, and its associated Saptamukhi estuary, represent only a small fraction of the larger Sundarbans mangrove ecosystem, we belive that the hydrologic, ecological, geological and climatic settings of the smaller sampling location is an ideal location for studying a subcomponent of CH 4 cycling in the vast mangrove region.

Field Sampling
The present study was conducted during June 2010 to December 2012. Sediment and air samples were collected from intertidal mangrove sediments and a watch tower, located at the center of the Lothian Island (21 • 42.58 ′ N: 88 • 18 ′ E), respectively (Figure 1). Water samples were collected in duplicate, on monthly basis from the Saptamukhi estuary along an estuarine stretch of 5 km (with 1 km intervals) near the Lothian Island jetty. Water temperature and pH were recorded in situ using a field thermometer and a portable pH meter (Orion Star A211), with a Ross combination electrode calibrated on the NBS (US National Bureau of Standards) scale (Frankignoulle and Borges, 2011). Reproducibility was ±0.005 pH units. Salinity and dissolved oxygen concentrations in surface and sub-surface waters were measured onboard, following the Mohr-Knudsen and Winkler titration methods, respectively (Grasshoff et al., 1983). Transparency of the water column was measured with a 15 cm diameter secchi disc. Water samples for dissolved inorganic nitrogen (DIN) analyses (e.g., nitrite, nitrate, and ammonia) were filtered in the field through Whatman GF/F filter paper, stored in HDPE bottles, and transported on ice to the laboratory.

Analysis of Sediment Samples and Flux Calculations across Sediment-Atmosphere and Sediment-Water Interfaces
Intertidal sediment cores (0-25 cm depth) were collected seasonally from different littoral zones (upper, mid, and lower) on the island using hand-held stainless-steel corers (diameter: 10 cm). Estuarine bottom sediments (sub-tidal) were collected using grab samplers. All sediment samples were sealed in polythene bags and stored on ice during transport to the laboratory.
Cores were sliced in the laboratory (using a N 2 -filled glovebag) at the following 5 cm intervals: 0-5; 5-10; 10-15; 15-20, and 20-25 cm. Sediment redox potential (E h ) was measured immediately using a platinum electrode and Ag/AgCl reference electrode (Fiedler et al., 2003). Wet-sediment samples (both intertidal and sub-tidal) were processed for measurement of CH 4 concentrations. Nitrite concentrations as well as sulfate, acidvolatile sulfide (AVS), and organic carbon percentages were also measured in sediments samples at all depth intervals .
CH 4 production rates of collected sediment samples (except 0-5 cm surface layer) were measured using anaerobic sediment incubations. Approximately 10 g of wet sediment were placed in an incubation bottle (1.2 cm i.d. and 10 cm long) fitted with a teflon septum. The bottles were then flushed with high purity N 2 for 1 min. to create completely anaerobic conditions, and incubated in dark at a fixed ambient temperature (i.e., field temperature at that time) for 24 h. At the end of the incubation, a 1 mL gas sample was withdrawn from the headspace through the rubber stopper using a gas-tight glass syringe (Lu et al., 1999). The collected gas sample was then analyzed for CH 4 concentration, based on the purity of nitrogen as blank. CH 4 production (in µmol CH 4 m −3 wet sediment d −1 ) was calculated based on CH 4 enrichment in the headspace and volume of the sediment.
CH 4 oxidation was measured for intertidal surface sediment following incubation of wet sediment (∼6 cm 3 ) in a 60 mL rubber septum-fitted flask, with a headspace that was spiked with a CH 4 standard (10.9 ppmv; procured from Chemtron Science Laboratories Pvt. Ltd.) (Saari et al., 1997). These flasks were then incubated at ambient temperature (i.e., field temperature) for 4 days in the dark. Gas samples from the headspace were drawn at the onset of incubation and over 24 h intervals until the termination of the experiment, and then analyzed for CH 4 concentrations. CH 4 oxidation was calculated by converting the loss of CH 4 in the headspace with respect to volume of sediment.
During low tide conditions, CH 4 emission rates from the intertidal sediment surface to the atmosphere, were measured using a static perspex chamber method (Purvaja et al., 2004;Dutta et al., 2015b). Basically, the chamber was placed in the sediment for a particular duration and CH 4 emission rate was calculated based on the enrichment of the CH 4 mixing ratio inside the chamber-relative to ambient air. Based on the precision of the air sample measurement method (discussed later), the minimum detectable flux using this chamber method would be 23.30 µmol m −2 d −1 . Advective CH 4 fluxes from intertidal forest sediment to the estuarine water column (F ISW ) were computed according to the method used by Reay et al. (1995): F ISW = × ν × C; where, = porosity of sediment = 0.58 , ν = mean linear velocity = d −1 (d = specific discharge), C = pore water CH 4 concentration in intertidal sediment. Porewater specific discharge was measured using a traditional, albeit rather "crude" method, which was based on the accumulation of pore water in an excavated pit of known surface area over time (Dutta et al., 2015a). This was performed during extreme low tide conditions on an intertidal flat at 100 m intervals during receding ebb flow. Diffusive CH 4 fluxes were calculated using Fick's law of diffusion (e.g., Sansone and Graham, 2004).

Analysis of Estuarine Water Samples and Water-Atmosphere CH 4 Flux Calculation
Collection and analysis of dissolved CH 4 concentration and CH 4 oxidation rates, in the Saptamukhi River estuarine waters, were performed according to the methods described by Dutta et al. (2013) and Dutta et al. (2015a), respectively; replicates were found to be within 2.30-3.14%. CH 4 fluxes across the air-water interface were calculated based on the wind speed parametization of the gas transfer velocity (k), using the following expression: (Liss and Merlivat, 1986). Although, water current speed can play a significant role in generating water-side turbulence (and therefore "k"), as recently reported by Ho et al. (2016), we were not able to include this parameter. A positive value denotes flux from water to the atmosphere and vice versa.

Analysis of Atmospheric Samples, Micrometeorology, and Biosphere-Atmosphere Flux Calculation
Samples for measurement of atmospheric CH 4 mixing ratios were collected on monthly basis using an air sampling bulb, at heights of 10 and 20 m above the forest floor of Lothian Island, and then transported to laboratory for analysis. Two reference gas standards (10.9 and 5 ppmv, procured from Chemtron Science Laboratories Pvt. Ltd.) were run before and after every measurement to check for instrument. Duplicate samples were analyzed periodically and the replicate measurements were found to be within 2-3.2%. Throughout the observation period the measured mixing ratios were corrected for water vapor. Meteorological parameters like air temperature and wind velocity were simultaneously recorded at 10 and 20 m heights above the forest floor of the island using a portable weather station (Model: Davis 7440). Mangrove biosphere-atmosphere CH 4 exchange fluxes (F BA ) were calculated using the following micrometeorological relationship (Barrett, 1998;Ganguly et al., 2008): where, χ = difference of CH 4 mixing ratio between 10 and 20 m heights, and V C = exchange velocity = 1/(r a + r s ) (r a = aerodynamic resistance and r s = surface layer resistance). Negative flux values indicate influx from the atmosphere to the biosphere, while positive flux indicates an efflux emission from the biophere.
The aerodynamic resistance (r a ) was calculated as follows (Wesely and Hicks, 1977): where, Z o is the roughness height, c is a correction function for atmospheric stability. c was computed as follows (Wesely and Hicks, 1977): Stable condition: c = −5 (Z/L) when 0 < Z/L < 1 (Z = height and L = Obukhov scale length, Z/L = atmospheric stability parameter).
The friction velocity (u * ) was calculated as follows: u * = k (u 1 − u 2 )/ln(Z 2 / Z 1 ) [k is the Von Karman constant, u 2 and u 1 are wind speeds at two heights, Z 2 and Z 1 ] Z o was determined by plotting the wind profile as lnZ versus u. The slope of the resulting straight line is k/u * and the intercept is lnZ o . For forest cover, a displacement length (d) equal to 80% of the mean height of the roughness element (mean height of mangrove plant = 10 m) was considered (Panofsky and Dutton, 1984). The scale length (L) was evaluated by the use of Pasquill stability classes A-F (Pruppacher and Klett, 1978) and is related to Z o (Golder, 1972) as: 1/L = a + b lnZ o , where "a" ranges between 0.035 and −0.096 and "b" ranges between 0.029 and −0.036.
at 10 m height was computed as: C d = r s / (ρ × V 102 ); ρ = air density and V 10 = wind velocity at 10 m height.
Sensible heat flux (H) was calculated as follows (Ganguly et al., 2008): H = ρ t C p (T 10m -T 20m )/(r a + r s ); ρ t and C p are density and specific heat of air, respectively.
Height (h) of planetary boundary layer (PBL) was computed as follows (Pal Arya, 2001): h = (0.25 u * )/| F |. F is the Coriolis parameters related to the rotational speed of the earth ( ) and latitude ( ) as F=2 sin .
CH 4 photo-oxidation rates (P) in the lower forest atmosphere were calculated based on the reaction (CH 4 + OH → CH 3 + H 2 O) and equation of P = k [CH 4 ] [OH] (Dutta et al., 2015b), where, k = rate constant of the reaction between CH 4 , and OH = 1.59 × 10 −20 T 2.84 exp (−978/T) cm 3 molecule −1 s −1 (Vaghjiani and Ravishankara, 1991), [CH 4 ] = mean of all CH 4 mixing ratio measurements during the day time at 10 m height in the diurnal cycle, and [OH] = mean of all OH radical concentrations during the day time at 10 m height in the diurnal cycle in molecules cm −3 . OH radical concentrations were computed using photolysis frequency of O 3 (O 3 mixing ratio in the mangrove forest atmosphere ranged between 14.66 ± 1.88 and 37.90 ± 0.91 ppbv; Dutta et al., 2015b), based on the empirical relation given by Ehhalt and Rohrer (2009).

Analytical Measurements
All CH 4 concentrations were measured according to method described Knab et al. (2009)-whereby, measurement of headspace CH 4 was performed using gas chromatograph (Varian CP3800 GC), fitted with chrompack capillary column (12.5 m × 0.53 mm), and a flame ionization detector (FID); the relative uncertainty for measurments were ± 2.9%.

Statisical Methods
All the statistical analyses were performed using Sigma Plot Statistical Software version 11. A simple t-test was used to test for significance acorss different seasons, while multiple regression analyses were used to identify key controlling factors for variability in pore water CH 4 concentrations as well as sediment and water column CH 4 oxidation.

Sediment Geochemistry in the Mangroves of Lothian Island
Spatial patterns in sedimentary CH 4 production in the mangrove links significantly with redox potential (E h ). E h values of mangrove surface sediment varied between −119.8 and −103.2 mV having a mean of −111.4 ± 6.78 mV (Figure 2). Mangrove sediment redox potential values decreased significantly with sediment depth-reaching a minimum at 25 cm (varied between −167.8 and −221.4 mV, mean = −186.7 ± 24.6 mV). These negative redox potential values indicated strong reducing conditions in mangrove sediments; even the oxic surface was not easily detectable. These observed redox condidions, which are highly conducive for anearobic microbial metabolism favors significant production of biogenic trace gases such as N 2 O, CH 4 , and H 2 S in this mangrove system . In contrast to the Sundarbans mangrove ecosystem, mangrove sediments of New Caledonia were found to have positive redox potential values, to a sediment depth of 30 cm, with negative E h values only found in considerably deeper sediments (Deborde et al., 2015). Such extreme differences sediment redox potential across different mangrove systems can typically be attributed to the geological setting (sediment porosity, organic matter supply), ecological factors (floral and faunal activities), hydrologic regime (tidal inundation and flushing). Past work has shown that methanogensis typically occurs at a redox of < −150 mV (Wang et al., 1993) similar to values observed in Lothian Island sediments (at depth of 10-25 cm), consistent with our previous work Dutta et al. (2013). Surface sediment SO 2− 4 -S concentrations ranged from 15.36 to 25.07 mM, with a mean of 20.54 ± 4.89 mM (Figure 3). The range of SO 2− 4 -S concentrations between sediment depths of 20-25 cm was 27.16-27.70 mM, were 1.34 times higher than surface sediment values. Once again, if we compare these SO 2− 4 -S concentrations to the mangrove system in New Caladonia, they are considerably lower (Deborde et al., 2015). One other difference was that in the Lothian Island Sundarbans sediment, SO 2− 4 -S concentrations had a significant decrease only within upper 15 cm depth of sub-surface mangrove sediment consistent  with patterns of peak sulfate reduction, with methanogenesis occurring deeper in the sediments (more discussion on this later).
AVS concentrations in the Lothian Island Sundarbans sediment varied between 13.32 and 13.95 mM with a mean of 13.63 ± 0.32 mM (Figure 3). In contrast to sulfate profiles, downcore AVS concentrations showed distinctly higher values within the 10-15 cm region of the sediment profile. Mangrove surface sediment NO − 2 -N concentrations varied between 1.61 and 4.14 µM, with a mean of 2.51 ± 1.41 µM. However, downcore NO − 2 -N concentrations reached a maximum between 10 and 15 cm depth in sediments (varied between 1.99 and 18.53 µM), possibly indicating the dominance of denitrifying bacterial activity at that depth (Figure 3). Percent OC in mangrove surface sediments on Lothian Island varied between 2.1 and 2.51, with a mean of 2.33 ± 0.21 (Figure 2), and decreased in concentration by ∼34% at sediments depths of 20-25 cm, where concentrations varied between 1.41 and 1.63%. The %OC measured here was lower than values previously reported by Jennerjahn and Ittekkot (1997) for the mangrove area in the Paraiba do Sul river mouth (4.82%). However, they were within the range of those reported by Bouillon et al. (2004) for mangroves located near Godavari estuary (0.6-31.7%), near the Indian and southwest coast of Srilanka. The decreasing trend of %OC coupled with decreasing redox potential across deep mangrove sediment indicates significant OC mineralization by anaerobic microbial respiration.

CH 4 Production Rates in Sediment
While only marginal CH 4 production (215-1,176 µmol CH 4 m −3 wet sediment d −1 ) was observed in sediments between 5 to 10 cm in depth, a significant [p (= 0.02) < 0.05] increase (2,089 to 4,033 µmol CH 4 m −3 wet sediment d −1 ) was observed at sediment depths of 10 to 15 cm (Figure 4). This trend of increasing CH 4 production rates with sediment depth continued even deeper in mangrove sediments with increases between 2,919 to 5,903 and 4,290 to 7,353 µmol CH 4 m -3 wet sediment d The production rates measured for this study site varied between 215 and 7,353 µmol CH 4 m −3 wet sediment d −1 , mean = 3,547 ± 2,205 µmol CH 4 m −3 wet sediment d −1 , and were within the range reported for pristine mangrove forests at Balandra, Mexico (Strangmann et al., 2008). While we did not measure CH 4 in the upper 0 to 5 cm in this study, due to the low methanogenic and intense methanotrophic activities in that layer, these rates are likely an underestimate-but were still quite high. However, in the bordering Saptamukhi estuarine system near our study site, CH 4 production rate in the top 5 cm of sub-tidal sediment varied between 18.72 and 85.74 µmol CH 4 m −3 wet sediment d −1 , with a mean of 48.88 ± 26.04 µmol CH 4 m −3 wet sediment d −1 .
Seasonally, mangrove sediment CH 4 production rates were higher [p (= 0.21) > 0.05] in post-monsoon (4,616 ± 2,666 µmol CH 4 m −3 wet sediment d −1 ) compared to pre-monsoon (2,378 ± 1,799 µmol CH 4 m −3 wet sediment d −1 ) seasons. The higher post-monsoon CH 4 production rate may be attributed to the presence of high mangrove leaf litter-fall at this time (58.79 g dry wt C m −2 month −1 ; Ray et al., 2011), which likely mediated by the supply of labile organic matter in these sediments-supporting higher CH 4 production. In contrast, higher salinity conditions during pre-monsoon may have partially inhibited sediment CH 4 production rate. For example, Marton et al. (2012) also reported lower CH 4 production rates during high salinity conditions. Other work in tidal mangrove forests have shown that intrusion of high salinity waters results in higher pore-water salinities, along with availibility of terminnal electron acceptors (e.g., NO − 3 , SO 2− 4 )-which constrain CH 4 production (Craft et al., 2009;Larsen et al., 2010).

Utilization of Organic Carbon and Linkages to Methanogenesis
Methanogens can only utilize a limited number of substrates and with the major pathways being (I) fermentation of acetate, commonly referred to as acetoclastic methanogenesis (AM) (Equation 1), and (II) acetate oxidation (Equation 2), where CO 2 is reduced with H 2 (hydrogenotrophic methanogenesis [HM]) (Equation 3) (Zinder, 1993;Mayumi et al., 2013). The aforementioned equations are shown below: Mayumi et al. (2013) calculated change in Gibbs free energy ( G • ) values for the aforementioned reactions in well-controlled CO 2 -injected microcosms. This work showed that acetate oxidation, under these microcosm conditions, was endergonicindicating a predominance of AM over the HM in a CO 2rich environment. Moreover, Lessner (2009) reported that ∼70% of biologically-produced CH 4 originated from the conversion of the methyl group of acetate to CH 4 . Thus, considering the CO 2 -rich character of the Sundarban mangrove environment (Biswas et al., 2004), we speculate that AM likely predominated in these sediments. If we make this assumption, the stoichiometric equation of acetoclastic methanogenesis (Equation 1) would require the utilization of 2 moles of OC per mole of CH 4 production. Along with production of CH 4, the reaction also produces one mole of CO 2 as by-product. Using the above equations, we calculated that in the intertidal sediment (from 0 to 25 cm depth), 4,966 µmol m −3 d −1 of OC was transformed through an AM pathway producing both CH 4 and CO 2 as 2,483 µmol m −3 d −1 . Similarly in the estuarine bottom sediment, adjacent to the island, 68.33 µmol m −3 d −1 of OC was transformed through the AM pathway producing 34.22 µmol m −3 d −1 of CH 4 and 34.22 µmol m −3 d −1 of CO 2 . The release of CH 4, and CO 2 as by-products, possibly enhancing potential regional effects on climate change. However, mangroves are also well-known to be blue carbon habitats, which sequester and store CO 2 (Bouillon et al., 2008;Alongi, 2012). Thus, further work is needed to establish the overall role of these mangroves as net sinks or sources of greenhouse gases. Finally, although we did not examine the role of sulfate reduction in this study, based on the relatively low %OC in these mangrove sediments (Prasad et al., 2017), it is likely to be another very important anaerobic process for sedimentary OC processing in the Sundarbans mangrove ecoystem.

Pore-Water CH 4 Levels and its Fluxes to the Adjacent Saptamukhi River Estuary
Pore-water CH 4 concentrations in intertidal mangrove sediments of Lothian Island were relatively constant within the 10 to 15 cm sediment layers (varied between 1,881 and 2,429 nM), and then abruptly [p (= 0.001) < 0.05] increased to concentrations of 3,804 and 4,542 nM between 15 and 20 cmwith further significant increases [p (= 0.02) < 0.05] to 5,274 and 6,670 nM between 20 and 25 cm (Figure 4) PW and E h (p = 0.017)-but not with others. This suggests the potential impact of sediment redox potential on the variability of mangrove sediment pore-water CH 4 concentrations. Mean sub-tidal sediment pore-water CH 4 concentrations were significantly higher in pre-monsoon (3,980 ± 1,227 nM) than in post-monsoon conditions (2,770 ± 1,039 nM) (Dutta et al., 2015a), which followed a similar seasonal trend for intertidal sediments. Overall, sediment pore waters in the intertidal mangrove sediments of Lothian Island were ca. 53.4 times more supersaturated than the adjacent estuarine sediments. This clearly indicates a significant source CH 4 to the estuarine CH 4 sediment pool, via advective (ranged between 115.81 ± 31.02 and 199.15 ± 47.89 µmol m −2 d −1 ; Figure 5) and diffusive fluxes (ranged between 7.06 ± 1.95 and 10.26 ± 2.43 µmol m −2 d 1 ; Figure 6) (Dutta et al., 2015a). Unfortunately, to the best of our knowledge no estimates have been made on advective CH 4 fluxes from intertidal mangrove to Saptamukhi River estuarine system. It should be noted that these are the the first reported measurments of advective CH 4 fluxes from Indian Sundarbans, albeit using a traditional rather "crude" method. Although, our estimated pore water specific discharge rates (3.60-5.76 cm d −1 ) were in the range of that measured for other mangrove regions located at temperate and tropical regions (2.10-35.50 cm d −1 ; Tait et al., 2016) by combined natural tracer Radon ( 222 Rn) with hydrodynamic models but as such, more accurate measurements are needed in near future to better constrain such measurements in the this globally-important ecosystem.
When compared to a temperate system, diffusive CH 4 fluxes from sub-tidal sediment of the Saptamukhi estuary were comparable to those found in the Arcachon tidal lagoon (France) (11.97 µmol m −2 d −1 ; Deborde et al., 2010). Similarly, our diffusive CH 4 fluxes were also higher than in the Yangtze temperate estuary, China (1.7-2.2 µmol m −2 FIGURE 5 | Seasonal variations of mangrove sediment pore water CH 4 concentration ([CH 4 ] IS ) and its advective flux to estuary (F ISW ) FIGURE 6 | Seasonal variations of sub-tidal sediment pore water CH 4 concentration ([CH 4 ] SS ) and its diffusive flux to estuary (F SSW ). d −1 ; Zhang et al., 2008), yet lower than White Oak river estuary, in northern California (17.1 µmol m −2 d −1 ; Kelley et al., 1990). While such cross-system comparisons are frought with the complexity of regional differences in ecological and physicochemical drivers of fluxes, we provide them for a simple geographic perspective here. Nevertheless, we can conclude that in the Indian Sundarbans fluxes measured during post-monsoon period were higher compared to premonsoon conditions, but were only significant for advective fluxes [p (= 0.0006) < 0.05] -not for diffusive fluxes [p (= 0.21) > 0.05]. The peak post-monsoon advective CH 4 flux may be due to maximal intertidal sediment pore water CH 4 concentrations, as well as pore water specific discharge (5.76 cm d −1 ). Peak diffusive CH 4 flux may be attributed to maximal subtidal sediment pore water CH 4 concentrations, which resulted in a greater CH 4 concentration gradient at sediment-water interface.

CH 4 Oxidation and Emission at/from Mangrove Surface Sediment
Mean surface sediment CH 4 oxidation rates (1.76 ± 0.34 mg m −2 d −1 ) measured in the Lothian Island mangrove forest were within the range of that reported by Bradford et al. (2001) and Jang et al. (2006). The mean peak pre-monsoon CH 4 oxidation rate (2.149 ± 0.16 mg m −2 d −1 ; Figure 7) are likely linked with maximum soil surface temperatures. Seasonal CH 4 exchange across mangrove sediment-atmosphere interface in Lothian Island mangroves revealed an annual mean flux of 452 ± 126 µmol m −2 d −1 (Figure 7), indicative of a net source of CH 4 from mangrove sediments to the atmosphere. These fluxes were within the range of those reported by Purvaja et al. (2004) in the Pichavaram mangrove, located at the northern end of the Cauvery delta, India, and along the coast of Puerto Rico (Sotomayor et al., 1994). Soil temperatures and salitinties have been show to have a significant impact on CH 4 emissions (F SA ) from mangroves (Bartlett et al., 1987;Lekphet et al., 2005). Lothian Island mangroves soil temperature and pore water salinity ranged between 18.25 ± 0.22 and 28.36 ± 1.02 • C (t) and 22.55 ± 0.31 and 28.88 ± 0.13 (s), respectively. Both soil temperature and salinity were maximal during pre-monsoon and minimal during post-monsoon seasons, respectively. "F SA " was best fitted linearly with "t" (R 2 = 0.35, F = 5.33, p = 0.041, n = 12), however, a second order polynomial relationship with "s" (R 2 = 0.77, F = 7.71, p = 0.029, n = 12), indicated a cumulative effect of temperature and salinity on variability of sediment CH 4 emission fluxes. A similar effect was previously reported in Ranong Province mangrove area, Thailand (Lekphet et al., 2005) and in a salt marsh at Queen's creek (Bartlett et al., 1987).
In Lothian Island mangroves, we observed higher spatial variability on CH 4 emissions (432-761 µmol m −2 d −1 ) in the upper littoral zone compared to mid and lower littoral zones. This was likely due to the higher pneumatophore density in that region (42 number m −2 ) and the associated plant-mediated diffusion of CH 4 through them . This further supports the role of plant-mediated transport of CH 4 across the sediment atmosphere interface in aquatic system (Chanton and Dacey, 1991;Laanbroek, 2010;Bridgham et al., 2013). In fact, a significant correlation was found between "F SA " with both pneumatophore (P no ) and bioturbation (B o ) densities, with positive correlation with pneumatophore (F SA = 400.13 + 0.56 P no ; R 2 = 81.9%, F = 6.94, p = 0.032, n = 20) and negative with bioturbation (F SA = 444.38-0.028 B no ; R 2 = 61.9%, F = 5.94, p = 0.041, n = 20). The positive correlation between sediment CH 4 emission flux with pneumatophore density via plant-mediation, and negative correlation via bioturbation e.g., burrowing crabs (Kristensen and Alongi, 2006), reflect the complex and changing dynamics of plant-and-animal populations on intertidal sediment-atmosphere CH 4 emission fluxes.

Physicochemical Properties of the Saptamukhi River Estuary
Physicochemical parameters of the Saptamukhi River estuarine waters are presented on seasonal basis in Figures 8, 9. For both temperature and salinity, the values were highest during the pre-monsoon and lowest during the post-monsoon (for temperature) and monsoon (for salinity) months. The lack of variability between estuarine surface and sub-surface water temperature and salinity clearly indicated the presence of a well-mixed water column (Dutta et al., 2015a). Seasonal differences of pH were not significant and ranged between 8.10 ± 0.03 and 8.17 ± 0.16. Dissolved oxygen (DO) concentrations in estuarine surface and sub-surface waters were high (0.19 ± 0.023 to 0.23 ± 0.036 and 0.17 ± 0.001 to 0.19 ± 0.025 mM, respectively; Figure 8), and DO % of saturation ranged from 94.8 to 99.3. These high oxygen values reflect a well-mixed oxygenated water column, that should restrain most anaerobic microbial metabolism of organic matter within estuarine water column. However, it should also be noted that CH 4 production in the isolated anoxic microhabitats of sinking particulate organic matter (POM), in well-oxygenated water column, have been observed in the open ocean (see Reeburgh, 2007).

Dissolved CH 4 Oxidation in Saptamukhi River Estuary
The Saptamukhi River estuarine water column is a welloxygenated system that favors the activity of methanotrophs.
Other studies have shown the mean dissolved CH 4 oxidation rate in these estuarine surface waters to be 20.59 nmol L −1 d −1 , and varies between 12.96 ± 2.88 and 30.24 ± 6.48 nmol L −1 d −1 ( Table 1) (Dutta et al., 2015a). Seasonally, CH 4 oxidation rates during the pre-monsoon period were significantly [p (= 0.002) < 0.05] higher than post-monsoon periods (Figure 10). Aquatic CH 4 oxidation has been shown to be significantly impacted by other physicochemical parameters such as temperature (thermal), salinity (tonicity), oxygen (oxidative), DIN (nutrient), and turbidity (surface). In the case of turbidity, estuarine methanotrophs, associated with suspended particulate matter, may be more effeicient in oxidizing CH 4 in the water column (Abril et al., 2007). but not with others. This further supports the need to examine broader spectrum and variability of drivers on CH 4 oxidation in these complex estuarine waters.

Saturation Percentage of CH 4 and Air-Water CH 4 Exchange
Percent saturation of dissolved CH 4 in the Saptamukhi River estuary ranged from 2,483 ± 950 to 3,525 ± 1,054, indicating that estuarine waters were supersaturated in CH 4 , highly conducive for CH 4 efflux across the water-atmosphere interface. Seasonal variability in CH 4 exchange across air-water interface of this estuary are shown in Figure 10 (Dutta et al., 2015a). Flux values calculated during the monsoon period (10.67 ± 6.92 µmol m −2 d −1 ) were significantly higher [p (= 0.007) < 0.05] compared to pre-monsoon conditions (6.27 ± 1.61 µmol m −2 d −1 ) ( Table 1). Air-water CH 4 fluxes were in the range of that calculated for other tropical estuaries, like the Hooghly estuary (0.88-148.63 µmol m −2 d −1 ) and Yangtze River estuary (6-25 µmol m −2 d −1 ) (Biswas et al., 2007;Zhang et al., 2008). The lowest CH 4 fluxes during pre-monsoon season were likely attributed to the coupled impact of low wind speeds, resulting low gas transfer velocities, as well as low dissolved CH 4 concentrations in surface waters. Estimated CH 4 flux values estimated for the Saptamukhi River estuary showed asignificant correlations with water temperature and salinity (water temperature: R 2 = 61%, F = 6.71, p = 0.029, n = 24; salinity: R 2 = 54%, F = 5.31, p = 0.037, n = 24), indicating cumulative effects of temperature and salinity on CH 4 emission fluxes (Dutta et al., 2015a). On an annual basis, the mean CH 4 emission flux across the wateratmosphere interface was 8.88 µmol m −2 d −1 , indicating that the estuary acted as a source of CH 4 to the regional atmosphere during this study period.

Micrometeorology and Atmospheric CH 4 Mixing Ratio
Previous work has shown that the lowest air temperatures and wind velocities, on Lothian Island, occurred during the post-monsoon season and were at their maximum during premonsoon conditions (Dutta et al., 2015b). Friction velocity (u * ), which controls stability of the atmosphere (e.g., atmospheric turbulence), varied between 0.01 and 1.2 m s −1 . PBL heights over the mangrove forest atmosphere (702.45 to 936.59 m), and were highest and lowest during pre-monsoon and monsoon periods, respectively. The values of drag coefficient (0.16-0.39) and roughness height (varied between 1.63 ± 1.02 and 3.77 ± 3.01 m) also followed the same seasonal trends. Dutta et al. (2015b) measured atmospheric CH 4 mixing ratios at 10 and 20 m heights in the mangrove forest on Lothian Island (Figure 11). The lowest atmospheric CH 4 mixing ratios were observed during pre-monsoon (at both 10 and 20 m heights) season, while maximum ratios were found during monsoon and post-monsoon periods for 10 m and 20 m heights, respectively ( Table 1). The significantly higher CH 4 concentrations in air [p (= 0.001) < 0.05] at 10 m height during monsoon season, compared to pre-monsoon, may be attributed to peak monsoon CH 4 emissions from sediments and aquatic surfaces, which primarily contribute to lower atmospheric CH 4 pool. We have also shown peak CH 4 concentrations in air occurring during early morning, which may be attributed to CH 4 accumulation within a stable boundary layer during that period (Dutta et al., 2015b). Concentrations of CH 4 in the lower atmosphere have been shown to decrease with advancing daylight due to increases in atmospheric turbulence-which tends to breakup a stable boundary layer (Mukhopadhyay et al., 2002). The atmospheric stability parameter (Z/L) and CH 4 mixing ratio (at 10 m height) were significantly correlated (R 2 = 74.8%, p < 0.001, F = 29.73, n = 30), this further supports that micrometeorology plays an important role in the variability of CH 4 mixing ratio in the lower atmosphere of mangrove forests in the Indian Sundarbans (Dutta et al., 2015b).

Exchanges of CH 4 across Biosphere-Atmosphere Interface
On an annual basis, atmosphere CH 4 mixing ratios measured at 10 m height over Lothian Island were 2% higher compared to 20 m, indicating CH 4 exchanges across the biosphere-atmosphere interface that were modulated by atmospheric turbulence. Our previous work (Dutta et al., 2015b) has shown monthly CH 4 exchange fluxes of across mangrove biosphere-atmosphere interface (Figure 11), that have maximal fluxes during monsoon and minimal during post-monsoon periods ( Table 1). These values also confirmed that mangroves were a CH 4 source to the upper atmosphere during the monsoon period (flux values positive), when χ was a significantly positive sink, in contrast to periods during pre-and post-monsoon seasons (flux values negative), when χ was negative. During the observation period, the mean mangrove biosphere-atmosphere CH 4 exchange flux was estimated to be 5.5 µmol m −2 d −1 , which indicated that on an annual basis these mangroves were dominant sources of CH 4 to the upper atmosphere. The mean compensation point, where the net biosphere-atmosphere CH 4 flux is zero, was found to be 1.997 ppmv. Sensible heat flux (H), which moderately controls atmospheric transport of energy and mass, was significantly correlated with biosphere-atmosphere CH 4 flux (F BA ). This was defined by a second order polynomial relationship (R 2 = 0.53, F = 7.72, p = 0.002, n = 12), further supporting the significant influence of sensible heat flux on the variability of biosphere-atmosphere CH 4 exchange fluxes (Dutta et al., 2015b).

Atmospheric CH 4 Photo-Oxidation
On annual basis, the mean daytime CH 4 mixing ratio was estimated to be 1.03 times lower than at nighttime; this variability is expected to be cumulatively governed by photo-oxidation and diurnal changes of PBL height. Dutta et al. (2015b) reported an insignificant correlation between daytime and nighttime CH 4 mixing ratios with PBL height ( CH 4 = −0.0118 PBL + 0.3045; R 2 = 16.2%, F = 0.321, p = 0.987, n = 12), and for the first time, showed the presence of CH 4 photo-oxidation within this tropical mangrove forest atmosphere. More specifically, this work reported CH 4 photo-oxidation rates in the forest atmosphere that varied between 6.05 × 10 10 and 1.67 × 10 11 molecules cm −3 d −1 , and were maximal and minimal during monsoon and postmonsoon periods, respectively ( Table 1). The significantly higher monsoon CH 4 photo-oxidation rate [p (= 3.2 × 10 −7 ) < 0.05] compared to post-monsoon, may be attributed to the combined effect of maximum CH 4 influx to atmosphere through emissions across sediment-atmosphere and water-atmosphere interfaces-as well as high UV indices and UV erythermal doses of irradiance that occur during this period in subtropical latitudes (Panicker et al., 2014). When considering a mean day light period of 12 h and that 6.023 × 10 23 molecules equals to 1 mole or 16,000 mg CH 4 , the mean CH 4 photo-oxidation rate to the atmosphere in this tropical mangrove-dominated island was calculated to be 3.25 × 10 −9 mg cm −3 d −1 .

Quantitative CH 4 Budget for Sundarbans Mangrove Ecosystem
A box model was developed on the biogeochemical cycling of CH 4 -as well as a quantitative CH 4 budget for the Sundarbans mangrove ecosystem (Figure 12). In the model different subecosystems are designated as separate reservoirs. The CH 4 pool in each reservoir and exchange fluxes between different reservoirs are presented on an annual mean basis. The major outputs from the box model are described below.

Mangrove/Intertidal Sediment Methane Budget
Total annual CH 4 production in Sundarbans mangrove sediment, within an intertidal sediment depth of 25 cm, was estimated to be 21.7 Gg year −1 , with a daily rate of 3,547 µmol m −3 d −1 . The CH 4 sediment reservoir pool was 0.03 Gg, with a mean pore-water CH 4 concentration 3,451 nM. Pore water was ∼55 times more supersaturated than adjacent estuarine waters (63.0 nM). This suggests that during low tide, depending upon the hypsometric gradient, there is a significant outflux of CH 4 -rich pore water from the intertidal mangrove sediments to estuary, via advective transport. Extrapolating a mean advective CH 4 influx from mangrove sediment to the adjacent estuary (159.5 µmol m −2 d −1 ) for entire intertidal area of Indian Sundarbans (45% of total forest area; http:// www.sundarbanbiosphere.org/html_files/sunderban_biosphere_ reserve.htm), we estimate that about 8.2% of the total mangrove sediment produced CH 4 is advectively transported to the adjacent estuarine system, with ∼12.4% undergoing oxidation at sediment surface by methanotrophs. Total annual CH 4 emission, from intertidal mangrove sediment to the forest atmosphere, was estimated to be 10.8 Gg year −1 (∼49.6% of total produced CH 4 in sediment), with a mean daily rate of 7.1 mg m −2 d −1 . When examining the ratio of [sediment emission]: [aerobic oxidation]: [advective flux], which computes to 6.05:1.51:1, it reveals that CH 4 emission is major CH 4 removal mechanism in these intertidal mangrove sediments. Excess CH 4 remaining ([ CH 4 ] S ) in these mangrove sediments, beyond the aforementioned "sink" mechanisms, was computed as follows: [ CH 4 ] S = Total CH 4 production − Total outflux (sediment emission + aerobic oxidation + advective flux). When using this equation, we find sink value [ CH 4 ] S of 6.4 Gg year −1 . This estimated value reflects the removal of CH 4 , via anaerobic methane oxidation (AMO), within sub-surface mangrove sediments, that were not considered in this study.

Sub-Tidal Sediment CH 4 Budget
Total CH 4 production at the 0-5 cm depth of sub-tidal sediment was 0.03 Gg year −1 , with a mean of 48.88 µmol m −3 d −1 . Mean sub-tidal sediment pore-water CH 4 concentration was 3,286 nM; which was ∼ 52.1 times supersaturated, compared to overlying estuarine water CH 4 levels. This supports the notion of diffusive CH 4 transport mechanism from sub-tidal sediments to the overlying estuary. On annual basis total diffusive CH 4 influx from sub-tidal sediment to the overlying estuary was 0.09 Gg, with a mean of 8.4 µmol m −2 d −1 .

Estuarine CH 4 Budget
The total CH 4 input to the Saptamukhi River estuary (by both advectice and diffusive transport) was 1.8 Gg year −1 . The ratio between advective and diffusive CH 4 flux was 20:1, indicating that advective CH 4 flux from intertidal mangrove sediments was the dominant driver in the buildup of the estuarine CH 4 pool. The estuary stands as a reservoir pool of 0.01 Gg CH 4, having a mean concentration of 63.0 nM. An estimated 74.5% of total CH 4 supplied to the estuary was removed from the estuarine systemvia microbial oxidation, with only 5.0% lost by an air-water CH 4 exchange flux. The ratio between microbial oxidation and airwater CH 4 exchange flux was 14:1, indicating microbial oxidation was the principal CH 4 removal pathway in this estuary. Mean turnover time of CH 4 in the water column relative to oxidation and emission was 3.7 days. After oxidation and emission, excess CH 4 left ([ CH 4 ] E ) in the estuarine system was computed as: [ CH 4 ] E = Total influx (Advective + Diffusive) − Total outflux (microbial oxidation + estuarine emission) Using the equation, [ CH 4 ] E was computed as 0.5 Gg year −1 , which is estimated CH 4 exported to the adjacent continental shelf region, thereby enriching the CH 4 pool in northern Bay of Bengal.

Atmospheric CH 4 Budget
The total annual CH 4 influx for Sundarbans mangrove ecosystem, via mangrove sediment and estuarine emissions, was 10.9 Gg year −1 , of which 99.1% was from sediments. Atmospheric CH 4 mixing ratios at 10 and 20 m heights in the forest atmosphere were 2.03 and 1.98 ppmv, respectively having mean of 2 ppmv. Extrapolating over entire Sundarbans, up to the mean PBL height (811.7 m), the atmosphere stands as a reservoir pool of 11.2 Gg CH 4 . The annual mean mangrove biosphere-atmosphere CH 4 exchange flux was estimated to be 5.38 µmol m −2 d −1 . However, when extrapolated for entire forest area it was found that only 2.75% of total annual CH 4 influx to the forest atmosphere was transported to the upper forest atmosphere. Moreover, ∼85% of total annual CH 4 influx to the forest atmosphere experienced photo-oxidation within the atmospheric boundary layer of Sundarbans, with a mean rate of 3.25 × 10 −9 mg cm −3 d −1 . Excess CH 4 that was from the forest atmosphere ([ CH 4 ] A ) beyond the photooxidation and biosphere-atmosphere fluxes was calculated as follows: [ CH 4 ] A = Total influx (estuarine emission + sediment emission) − Total outflux(photo-oxidation + biosphere-atmosphere flux) Using this mass balance equation [ CH 4 ] A , we estimated that 1.3 Gg year −1 enriches the regional atmospheric CH 4 mixing ratio and further contributes to regional climate change scenarios. The impact of these gases need to be considered as a delineated reservoir, within the context of a regional climate and atmospheric boundary layer height (ABL) (mean height = 811.7 m; Dutta et al., 2015b). Considering a radiative forcing efficiency of CH 4 in the atmosphere as 3.7 × 10 −4 W m −2 ppb −1 (https://www.ipcc.ch/publications_and_data/ar4/wg1/en/ ch2s2-10-2.html), [ CH 4 ] A in this mangrove forest atmosphere resulted 0.11 Wm −2 year −1 of radiative forcing to the regional atmosphere.

CONCLUSIONS
In intertidal mangrove sediment column CH 4 production rate, down to depth of 25 cm in the sediments, was estimated to be 21.75 Gg year −1 with a mean pore water CH 4 concentration was 3,541 nM. CH 4 emission across intertidal sediment-atmosphere interface acted as major sink for the CH 4 produced in intertidal sediments over surface layer CH 4 oxidation and advective CH 4 transport to Saptamukhi River estuary. The estuary, which was well-oxygenated, which constrained methanogenesis within estuarine water column, where dissolved CH 4 was considered to be entirely exogenous in nature. Advective CH 4 flux, which was 20 times higher than diffusive flux, was the major source for CH 4 to the estuary. CH 4 oxidation which is 14 times higher than the water-atmosphere exchange, was considered the principal CH 4 removal mechanism in this estuary. On an annual basis, total CH 4 emissions from sediments and waters of the Sundarbans mangrove biosphere was 10.9 Gg of which sediment contributed principally 99.1%. Compared to total CH 4 supply to the forest atmosphere, about 85% photooxidized within atmospheric boundary layer of Sundarbans while 2.75% is transported to the upper atmosphere through mangrove biosphere-atmosphere CH 4 exchange flux. Based on our proposed CH 4 budget the mangrove forest atmosphere resulted in a radiative forcing 0.11 Wm −2 year −1 to the regional atmosphere.

AUTHOR CONTRIBUTIONS
MD and SM: Have desinged the investigations, data genaration as well as manuscript preparation. The author TB has corrected the manuscript and made significant improvement of the manuscript.