Changes in Salinity, Mangrove Community Ecology, and Organic Blue Carbon Stock in Response to Cyclones at Indian Sundarbans

Climate change-induced frequent cyclones are pumping saline seawater into the Sundarbans. Fani, Amphan, Bulbul, and Yaas were the major cyclones that hit the region during 2019–2021. This study represents the changes in the soil parameters, mangrove biodiversity and zonation due to the cyclone surges in the Indian Sundarbans between 2017 and 2021. Increasing tidal water salinity (parts per thousand) trends in both pre-monsoon (21 to 33) and post-monsoon (14 to 19) seasons have been observed between 2017 and 2021. A 46% reduction in the soil organic blue carbon pool is observed due to a 31% increase in soil salinity. Soil organic blue carbon has been calculated by both wet digestion and the elemental analyzer method, which are linearly correlated with each other. A reduction in the available nitrogen (30%) and available phosphorous (33%) in the mangrove soil has also been observed. Salinity-sensitive mangroves, such as Xylocarpus granatum, Xylocarpus moluccensis, Rhizophora mucronata, Bruguiera gymnorrhiza, and Bruguiera cylindrica, have seen local extinction in the sampled population. An increasing trend in relative density of salinity resilient, Avicennia marina, Suaeda maritima, Aegiceras corniculatum and a decreasing trend of true mangrove (Ceriops decandra) has been observed, in response to salinity rise in surface water as well as soil. As is evident from Hierarchical Cluster Analysis (HCA) and the Abundance/Frequency ratio (A/F), the mangrove zonation observed in response to tidal gradient has also changed, becoming more homogeneous with a dominance of A. marina. These findings indicate that cyclone, climate change-induced sea level rise can adversely impact Sustainable Development Goal 13 (climate action), by decreasing organic soil blue carbon sink and Sustainable Development Goal 14 (life below water), by local extinction of salinity sensitive mangroves.


Introduction
Mangroves are a plant group that possesses physiological or anatomical adaptations to survive in a hyper-saline tidal environment. These plants grow between the high tide level (HTL) and low tide level (LTL) along the tropical and subtropical coasts. Mangroves render several ecosystem services, such as protection against natural disasters and storm surges, blue carbon sequestration, a breeding ground for marine and estuarine biota, and consolidating soil against tidal currents [1]. Climate change is destabilizing both terrestrial and marine ecosystems worldwide. In some parts, it is increasing desertification along with the frequency of droughts, whereas on the tropical/subtropical coasts, it is increasing the frequency as well as the intensity of cyclones/typhoons/storm surges [2]. Caught  [15] The hypothesis of this study is that a cyclone-induced salinity rise in surface water will decreases the blue carbon sequestration ecological function of mangroves, change the mangrove community ecology and zonation pattern in the mudflats. Thus, the aim of the present study is to trace the gradual changes in mangrove community ecology due to the changes in surface water salinity because of frequent cyclonic surges in the duration of four years from 2017 to 2021 and to access the alteration in ecosystem services (organic blue carbon sequestration) and ecosystem functions (nitrogen and phosphorus biogeochemical cycles).

Study Site
The study site has been selected in the central part of the Indian Sundarbans (near Amtali and Chotto-mollakhali revenue villages), as detailed in Figure 1. Located on both sides of a major creek, a natural mangrove patch and its location make it a suitable plot for studying the impact of surface water salinity on the mangrove community ecology. The study site covers a length of 1017 m and a maximum breadth of 164.1 m, spreading from 22 • 53 22 E. The location is an extended mudflat, though near the human settlement, and has limited anthropogenic interventions. Hence, the natural growth of mangroves and the changes in mangrove community ecology can be studied. The area is sheltered on all sides by human-inhabited islands and is about 5 km west of the Jhilla River, which is the natural border between India and Bangladesh. This river also connects to the reserve forests (Sajnekhali Sanctuary, Jhilla Reserve Forest, and Sundarban National Park, India). The sites are on both sides of a creek with an approximate width of 126 m.
four years from 2017 to 2021 and to access the alteration in ecosystem services (organic blue carbon sequestration) and ecosystem functions (nitrogen and phosphorus biogeochemical cycles).

Study Site
The study site has been selected in the central part of the Indian Sundarbans (near Amtali and Chotto-mollakhali revenue villages), as detailed in Figure 1. Located on both sides of a major creek, a natural mangrove patch and its location make it a suitable plot for studying the impact of surface water salinity on the mangrove community ecology. The study site covers a length of 1017 m and a maximum breadth of 164.1 m, spreading from 22° 14′ 39" N, 88° 53′ 14" E to 22° 15′ 04" N, 88° 53′ 36" E and 22° 14′ 50" N, 88° 53′ 25" E to 22° 14′ 53" N, 88° 53′ 22" E. The location is an extended mudflat, though near the human settlement, and has limited anthropogenic interventions. Hence, the natural growth of mangroves and the changes in mangrove community ecology can be studied. The area is sheltered on all sides by human-inhabited islands and is about 5 km west of the Jhilla River, which is the natural border between India and Bangladesh. This river also connects to the reserve forests (Sajnekhali Sanctuary, Jhilla Reserve Forest, and Sundarban National Park, India). The sites are on both sides of a creek with an approximate width of 126 m.

Biodiversity Assessment
A temporal study on the mangrove biodiversity was conducted in twelve quadrat plots (10 m × 10 m) and transects (5 m × 20 m) at the study site every year in the month of June between 2017 and 2021. The plant species are noted for reference, and the average individual number of each species was calculated based on the surveyed data following the  [17] and Chowdhury et al. [5]. Through relative density, the composition and richness of the mangrove species were determined.

Relative Density (D) =
Number of individuals of a species Number of individuals of all the species × 100 This estimator provides information on a species' status across the board (i.e., across all species present in the region), which aids in the comparison of that species' relative positions in the target environment. Simpson's index for dominance, diversity, and Shannon's diversity index were used to analyze the effect of the species on the health of the ecosystem, their diversity, and dominance [5]. The ratio between abundance and frequency (A/F) has been estimated to shed light into the distribution pattern of the species [7]. An A/F value < 0.025 indicates 'regular distribution', whereas 0.025-0.05 indicates 'random', and values > 0.05 indicate 'contiguous' distribution pattern.

Salinity Assessment of Tidal Water
Salinity assessment has been studied during the pre-monsoon (May) and post-monsoon (October) seasons at the same site every year from 2017 to 2021. Annual salinity is calculated by taking the mean values of pre-and post-monsoon salinity assessments. During each sampling period in a year, five samples (n = 5) have been collected and analyzed. The first sample was collected during the peak of diurnal high tide; the second sample was collected during the peak of diurnal low tide on the same day; the third sample was collected in the night of the first day during the receding tide; the fourth sample was collected during the peak of high tide the next day; and the fifth sample was collected during the second day's peak low tide. This sampling method was followed to minimize the temporal fluctuations in the salinity results. Changes in the salinity profile in the region are assessed using the Mohr-Knudsen method [5,20]. Then, salinity is calculated from the formula Salinity (S) = 1.805Cl − + 0.03 g/kg

Soil Sampling and Analyses
Each year from 2017 to 2021, soil samples were collected from the biodiversity sampling location (10 m × 10 m quadrats). The intact cores were persevered in icebox and bought to the lab and homogenized. The bulk samples were physically sorted from the gravel and coarse roots, and the soils were then lightly pulverized with a mortar and pestle. Three replicate subsamples of the sample (10 cm deep) were taken, and then combined and homogenized to generate a composite soil sample. Before being delivered to the lab, the gathered soil samples were put in plastic bags and marked. To remove large pebbles and plant debris, samples were sieved using a 2 mm sieve after being air-dried at ambient temperature (25-30 • C). The core approach was used to calculate soil bulk density.
According to Fang et al., litter primarily affects the top 10 cm of soil depth, and cyclonic influences can be noticed in the topsoil layers [21]. Studies conducted by Castañeda-Moya et al. to understand the impact of Hurricane Wilma in the mangroves of the Florida Coastal Everglades also considered 0-10 cm depth, as the hurricane-transported sediment was observed to impact 0.5-4.5 cm soil depth [22].
A multiparameter pH probe (HI-2020, Hanna Instruments, India) was used to measure the pH of the soil. After being shaken and allowed to settle for an hour, soil pH was measured in a suspension of soil and deionized water (1:2.5, w/v), and bulk density was calculated according to an established protocol [23]. After digestion with 0.32% KMnO 4 solution and titration with 0.02 N H 2 SO 4 , the semi-automatic nitrogen estimation system calculated the amount of nitrogen (N) that was available to plants. The Olsen method [24] was used to determine how much phosphorus (P) was available to plants. To determine soil organic carbon (SOC EA ) using the elemental analyzer (Euro Vector, EA 3100), the air-dried soils were pretreated with 6 M HCl to eliminate inorganic carbon present in the form of carbonates [25]. The elemental analyzer has been widely used in the carbon estimation of coastal saline or hypersaline soil [26,27].

Wet Digestion Method
Soil organic carbon (SOC w ) was determined through the wet digestion method [23]. An oxidizing agent (K 2 Cr 2 O 7 in an acidic H 2 SO 4 medium) has been used to react with organic matter (reducing agent) present in the soil samples in the presence of a color change indicator [27]. was titrated with ammonium iron (II) sulfate with a visible color change from blue to green, and then, it was stoichiometrically correlated with reduced substrate (SOC). As applied by Chowdhury et al., in the analysis of SOC, 5 gm Ag 2 SO 4 L −1 of H 2 SO 4 was added before the use of H 2 SO 4 to minimize the interference of Cl − in saline mangrove soil [1]. This Walkley and Black chromic acid wet oxidation method has been widely used by researchers estimating the SOC and the blue carbon stock of mangrove, hypersaline edaphic environments [27,28].

Organic Blue Carbon Pool Assessment
The following calculations were made for the carbon pool (0-10 cm) utilizing SOC EA concentration, bulk density, and specific soil depth [27,29]: where BOC S is the blue organic carbon stock in kg C m −2 , SOC EA is the soil organic carbon estimated by the elemental analyzer in g kg −1 , BD = bulk density (g cm −3 ), and T = soil thickness (cm).

Statistical Tests
SPSS 16 (SPSS Inc. Chicago, IL, USA) and Microsoft Excel 2007 have been used for statistical analysis. The trends in temporal changes in salinity regime and biodiversity have been studied using regression analysis. The standard error is used to assess the data's dependability. Analysis of Variance (one-way ANOVA) has been used to analyze the variation between mean values, followed by Duncan's post-hoc test. Hierarchical Cluster Analysis (HCA) was used on the species distribution with 'Ward linkage' and 'Squared Euclidian distance' as the major attributes of determination.

Variations in Tidal Water Salinity Profile
The result showed a drastic change in the salinity profile in the region owing to cyclone impacts. As evident from Figure 2a, the salinity in the pre-monsoon assessment is higher than the post-monsoon salinity depicted in Figure 2b. The pre-monsoon salinity ranges from 21,240 mg L −1 (21.2 ppt) in 2017 to 32,580 mg L −1 (32.6 ppt) in 2021. Post-monsoon salinity was far lower than the pre-monsoon estimates, ranging from 13,320 mg L −1 (13.3 ppt) in 2017 to 19,300 mg L −1 (19.3 ppt) in 2021. The reason behind this is the impact of monsoon rains on the salinity profile of the region. According to previous studies, it is evident that the annual rainfall in the region is high, around 1800 mm y −1 , of which 80% of precipitation occurs during the monsoon months (July-October) [30]. The regression trend line in both pre-monsoon and post-monsoon estimates shows an increasing salinity trend from 2017 to 2021.  (13.3 ppt) in 2017 to 19,300 mg L −1 (19.3 ppt) in 2021. The reason behind this is the impact of monsoon rains on the salinity profile of the region. According to previous studies, it is evident that the annual rainfall in the region is high, around 1800 mm y −1 , of which 80% of precipitation occurs during the monsoon months (July-October) [30]. The regression trend line in both pre-monsoon and post-monsoon estimates shows an increasing salinity trend from 2017 to 2021. An increasing salinity trend has been observed since 2017, as shown in Figure 2c. The one-way ANOVA test revealed a difference in the yearly means from 2017 to 2021, with the salinity profiles of 2017, 2018, and 2019 showing a similar trend, while those of 2019, 2020, and 2021 show similar trends indicating the cyclone-induced sea water intrusion is there from 2019 to 2021, as per the Duncan post-hoc test followed by ANOVA analysis. The F value was found to be significant with a probability level of 0.017. This firmly indicates that the cyclone induced, increase in salinity in the region. The year 2019 witnessed two cyclones, Fani and Bulbul, while the following years recorded one annual cyclonic event, namely Amphan and Yaas (Table 1). Hence, the ANOVA test depicted in Figure 2c showed 2019 having similar salinity data to previous years-2018, 2017, and coming years impacted by frequent cyclones, i.e., 2020 and 2021. This assessment is in concurrence with the existing literature, which has also indicated an increase in salinity in surface water post-cyclone. The Amphan cyclone was associated with an increase in salinity of ∼1.5 PSU An increasing salinity trend has been observed since 2017, as shown in Figure 2c. The one-way ANOVA test revealed a difference in the yearly means from 2017 to 2021, with the salinity profiles of 2017, 2018, and 2019 showing a similar trend, while those of 2019, 2020, and 2021 show similar trends indicating the cyclone-induced sea water intrusion is there from 2019 to 2021, as per the Duncan post-hoc test followed by ANOVA analysis. The F value was found to be significant with a probability level of 0.017. This firmly indicates that the cyclone induced, increase in salinity in the region. The year 2019 witnessed two cyclones, Fani and Bulbul, while the following years recorded one annual cyclonic event, namely Amphan and Yaas (Table 1). Hence, the ANOVA test depicted in Figure 2c showed 2019 having similar salinity data to previous years-2018, 2017, and coming years impacted by frequent cyclones, i.e., 2020 and 2021. This assessment is in concurrence with the existing literature, which has also indicated an increase in salinity in surface water post-cyclone. The Amphan cyclone was associated with an increase in salinity of ∼1.5 PSU in the impacted region [12]. After the impact of the Fani and Bulbul cyclones, an increase in surface water salinity was observed in the area [11,31]. The Yaas cyclone resulted in increased salinity due to flood water influx in the region, which also impacted agriculture [13].

Variations in Mangrove Diversity
Previous research has indicated that both Avicennia marina and Aegiceras corniculatum can tolerate higher salinities than similar other mangrove [32][33][34]. Both are sympatric species with similarities in salinity tolerance, and the former has resilience to cyclone disasters [14,35]. In the current study, Figure 3 shows the relationship between the relative densities of A. marina, A. corniculatum and Ceriops decandra with average annual salinity recorded from 2017 to 2021. Figure 3 indicates that there is an observable positive trend between increases in the number of the former two species with rising salinity, while the last species is a true mangrove of the Rhizophoraceae family, which seems to be sensitive to rapid salinity fluctuations. Previous studies on the impact of salinity on these species show that it is a common trend for reduction of propagule density in C. decandra due to rising salinity [14]. Avicennia alba showed a drastic reduction in relative density from 13.3 to 1.6 from 2017 to 2021, respectively. Proteresia corctata is a Poaceae family member that has been observed colonizing the newly formed mudflats across the Indian Sundarbans. The relative density of the mangrove grass increased from 12.6 to 16 from 2017 to 2021. It is also evident from Table 2 that the relative density of the mangrove grass never dropped below 10, indicating its resilience to the salinity rise induced by the saltwater influx due to cyclones [36]. This indicates that the sedimentation in the mudflat has created an opportunistic ground for the ecosystem to bounce back if the salinity regime is ameliorated in the future [1].     Table 2). Relative density reflects a steady reduction in the biodiversity of salt-sensitive species in response increase in salinity. Xylocarpus granatum, Xylocarpus moluccensis, Rhizophora mucronata, Bruguiera gymnorrhiza, and Bruguiera cylindrica used to have less representation in 2017 but went locally extinct from the studied population by 2021. Previous research on the physiological adaptation of B. gymnorrhiza indicated that an increase in salinity results in an augmentation in the evapotranspiration rate of the seedlings [37]. B. cylindrica is a sympatric species like B. gymnorrhiza. X. granatum and X. moluccensis share a similar niche and position in the Sundarbans mangrove ecosystem. X. granatum seedlings show greater growth in freshwater in comparison to higher salinity (around 23 ppt), as per a previous study conducted in Micronesia [38]. Following a similar trend, the species showed extinction in the studied population when the salinity rose to around 26 ppt in 2020 due to the Amphan super cyclone. Even without salinity-increasing impacts, these species showed smaller population proportions (relative density = 0.94 for X. granatum and relative density = 0.56 for X. moluccensis) in 2017, and when salinity stress was impacted due to cyclones, the regeneration of the seedling became difficult due to seedling mortality. Hence, the population would be extinct by 2021. R. mucronata shows tolerance to salinity fluctuations by increasing xylary vessel density [39]. However, if there is a rapid fluctuation in salinity, it may not have enough time to adjust to the changing abiotic factors. The Rhizophoraceae (true mangrove) member has a very small population; beginning in 2017, the population could not adapt significantly to adjust to the changing salinity pattern and became extinct in the regional population. Previous studies have already shed light on the proneness of rare plant species to go extinct in the case of abiotic changes or sudden random events [40]. A similar case was noticed where R. mucronata dominated shorelines, resulting in 100% mortality due to the impact of Cyclone Eline [41]. This indicates that this species is less resilient to cyclone stress. Sonneretia caseolaris showed a complete extinction in the population until 2020, which may be due to the Amphan super-cyclone-induced flushing of saltwater in the inner estuaries. However, as per previous estimates, there are populations of S. caseolaris in the eastern part of the Indian Sundarbans, and the propagules floating through the creek may have settled and started a natural regeneration during 2021. Previous studies also indicated the resilience of this species during instances of cyclones or storm surges [42]. An increase in the relative density of Excoecaria agallocha from 5.6 in 2017 to 10.87 in 2019 and a reduction back to 7.52 in 2021 indicated a variable but salinity-resilient response of this mangrove associate. Studies conducted at other mangrove habitats under cyclone impact show the resilience of this plant [39]. Phoenix padulosa is generally observed in degraded mudflats as an edaphic sub-climax species [5,43]. In the current study, this species shows a decreasing distribution across the years, with a relative density of 3 in 2017  Figure 3 shows the relationship between relative density and salinity rise.
Overall, the study indicates an increasing dominance of salinity-resilient A. marina over the years, which is one of the most resilient mangroves [5,44]. HCA has been used to find associations between the species across the years (Figure 4).  Table 2.
The HCA indicates the same association as is evident from the visual and species diversity analysis. A. marina was found to be distributed without any specific association with other species, which is probably due to its high salinity tolerance. Species such as X. granatum, X. moluccensis, R. mucronata, B. gymnorrhiza, B. cylindrica, R. mucronata, C. decandra, S. caseolaris, and P. padulosa, showed similar associations due to their decreasing diversity because of the salinity rise induced by cyclones. These species are under the maximum threat of local extinction due to the cyclone-induced salinity rise. P. coarctata: as the pioneer species in the mudflat mangrove succession, it is seen to be distantly associated with species that can tolerate and survive the increasing salinity regime, and hence is closely clustered with species that can tolerate the cyclonic disturbances. This is also evident with the increasing value of Simpson's Index of Dominance from 0.08 in 2017 to 0.16 in 2021 (Table 3). S. maritima is a mangrove associate herb with high salinity tolerance and is mostly seen near the low tide level in salt-encrusted mudflats, as evident in previous study [5]. This mangrove associate is showing a weak but positive linear correlation with the rise in annual salinity gradient of tidal water. This indirectly impacted the diversity of the meta-population, as the dominance index is reciprocal to diversity estimation; hence, more dominance indicates lesser diversity [5,45]. Similarly, as Simpson indexes are biased toward dominant species, the information statistic index of the Shannon Index of Diversity reflects the overall health of the ecosystem, considering the rarer species in the equation (Table 3) [45]. The Shannon index generally ranges from 1.5 to 3.5 in a natural ecosystem, and it reflects the overall health of the ecosystem. Table  3 Table 2.
The HCA indicates the same association as is evident from the visual and species diversity analysis. A. marina was found to be distributed without any specific association with other species, which is probably due to its high salinity tolerance. Species such as X. granatum, X. moluccensis, R. mucronata, B. gymnorrhiza, B. cylindrica, R. mucronata, C. decandra, S. caseolaris, and P. padulosa, showed similar associations due to their decreasing diversity because of the salinity rise induced by cyclones. These species are under the maximum threat of local extinction due to the cyclone-induced salinity rise. P. coarctata: as the pioneer species in the mudflat mangrove succession, it is seen to be distantly associated with species that can tolerate and survive the increasing salinity regime, and hence is closely clustered with species that can tolerate the cyclonic disturbances. This is also evident with the increasing value of Simpson's Index of Dominance from 0.08 in 2017 to 0.16 in 2021 (Table 3). S. maritima is a mangrove associate herb with high salinity tolerance and is mostly seen near the low tide level in salt-encrusted mudflats, as evident in previous study [5]. This mangrove associate is showing a weak but positive linear correlation with the rise in annual salinity gradient of tidal water. This indirectly impacted the diversity of the meta-population, as the dominance index is reciprocal to diversity estimation; hence, more dominance indicates lesser diversity [5,45]. Similarly, as Simpson indexes are biased toward dominant species, the information statistic index of the Shannon Index of Diversity reflects the overall health of the ecosystem, considering the rarer species in the equation (Table 3) [45]. The Shannon index generally ranges from 1.5 to 3.5 in a natural ecosystem, and it reflects the overall health of the ecosystem. Table 3 indicated a gradual decline in ecosystem health from 2017 to 2021, with Shannon index values of 2.51 to 2.06. This shows a declining diversity and increasing dominance of salinity-tolerant mangrove species. A phytosociological study also revealed a change in the zonation pattern at Indian Sundarbans, which is depicted in Figure 5.  A phytosociological study also revealed a change in the zonation pattern at Indian Sundarbans, which is depicted in Figure 5.  (2017), the study area was divided into five zones, where Zone 1 has a predominance of Suaeda maritima and Proteresia coactata, followed by Zone 2 with Avcennia marina, Acanthus ilicifolius and Aegiceras corniculatum, followed by Avicennia alba, Avicennia officinalis and Rhizophora mucronata Zone 3. Zone 4 is mostly dominated by Aegialitis rotundifolia and Ceriops decandra, while the last zone (Zone 5) has tree species such as Sonneretia caseolaris, Xylocarpus granatum, Xylocarpus moluccensis and Excoecaria agallocha. (b) Phytosociological assessment of the study area after consecutive cyclone instances and a rise in salinity (2021) indicated a shift in the community composition in the zones. Zone 1 is again predominated by Suaeda maritima and Proteresia coactata followed by Avicennia marina and Acanthus ilicifolius in Zone 2. Zone 3 has an association of Avicennia marina and Aegiceras corniculatum followed by Avicennia officinalis and Avicennia alba, while Zone 5 was not found to be distinctly identifiable. Hence, it is seen that mangrove zonations in response to tidal amplitude (from LTL to HTL) are becoming homogeneous with a dominance of A. marina in most of the associations.
Previous research has highlighted the presence of zonation at the Indian Sundarbans from low tide level (LTL) to high tide level (HTL) [46,47]. In the present study, a similar zonation was observed in the study site. However, over the period of 5 years, A. marina started to dominate both Zones 2 and 3, while the true mangrove tree species was replaced  (2017), the study area was divided into five zones, where Zone 1 has a predominance of Suaeda maritima and Proteresia coactata, followed by Zone 2 with Avcennia marina, Acanthus ilicifolius and Aegiceras corniculatum, followed by Avicennia alba, Avicennia officinalis and Rhizophora mucronata Zone 3. Zone 4 is mostly dominated by Aegialitis rotundifolia and Ceriops decandra, while the last zone (Zone 5) has tree species such as Sonneretia caseolaris, Xylocarpus granatum, Xylocarpus moluccensis and Excoecaria agallocha. (b) Phytosociological assessment of the study area after consecutive cyclone instances and a rise in salinity (2021) indicated a shift in the community composition in the zones. Zone 1 is again predominated by Suaeda maritima and Proteresia coactata followed by Avicennia marina and Acanthus ilicifolius in Zone 2. Zone 3 has an association of Avicennia marina and Aegiceras corniculatum followed by Avicennia officinalis and Avicennia alba, while Zone 5 was not found to be distinctly identifiable. Hence, it is seen that mangrove zonations in response to tidal amplitude (from LTL to HTL) are becoming homogeneous with a dominance of A. marina in most of the associations.
Previous research has highlighted the presence of zonation at the Indian Sundarbans from low tide level (LTL) to high tide level (HTL) [46,47]. In the present study, a similar zonation was observed in the study site. However, over the period of 5 years, A. marina started to dominate both Zones 2 and 3, while the true mangrove tree species was replaced by A. alba, A. officinalis, and C. decandra in Zone 4. The zones are becoming more and more homogeneous with an overall dominance of A. marina. During the 2017 assessment, A. marina was seen to be dominating in Zone 2, whereas in 2021, it was found to be dominating in both Zones 2 and 3, whereas Zone 5 was not observed, which was unlike in the 2017 assessment. The A/F value noted in Table 2 indicates a 'contiguous' (>0.05) distribution pattern for S. maritima, P. corctata, A. marina, E. agallocha, A. ilicifolius, and A. corniculatum between 2017 and 2021, but for C. decandra, the distribution pattern changed from 'contiguous' in 2017 to 'random' in 2021 (i.e., 0.025-0.05). Species such as X. granatum, X. moluccensis, R. mucronata, B. gymnorrhiza, and B. cylindrica showed a 'regular' (< 0.025) distribution in 2017 but became locally extinct by 2021. This indicates that these species are sensitive to random environmental fluctuations due to the cyclones and need further research on identifying the best restoration techniques.

Changes in Carbon Pool and Soil Parameters
The concentrations of soil organic carbon, soil salinity, available nitrogen, and available phosphorus and the nutrient cycling in sediments can be critical decisive factors controlling the density and development of mangrove forests. The growth pattern of mangrove forest is highly dependent on the above-mentioned abiotic factors. However, the change in the mangrove biodiversity also plays a crucial role in the change in these soil parameters as evident from earlier studies [48,49]. Thus, the changes in biodiversity have an impact on ecological services.
Research has shown a decrease in blue carbon pools with degrading or deteriorating ecological functions and biodiversity [34,50]. Table 4 elucidates the changes in soil chemical constituents at the sampling site during the study period. The pH and bulk density do not show significant changes with the year, but one-way ANOVA indicates a significant temporal variation in SOC EA , SOC W , available nitrogen, and available phosphorous. Soil salinity significantly increased from 2019 onwards after the Fani, Bulbul, Amphan and Yaas cyclones hit the region. However, a slight increase in salinity was there already in 2018. The post-hoc test is also indicating the similarity between the salinity in year 2017 and 2018 and between 2019 and 2021. Soil salinity is also dependent on tidal water salinity, which is evident from Figure 6. The letters in the superscripts indicate the DMRT post-hoc test, and the significant 'F' results (with significance bellow 0.5) have been put in bold.
The concentrations of available N and available P were recorded as being approximately 1.5 times lower in 2021 in comparison to 2021 ( Table 4). The gradual decrease in the concentration of the available N and the P with increasing soil salinity indicates the impact of salinity on the micro-organisms as well as on the reducing mangrove density. The cyclone Fani hit west Bengal in 2019, and a significant decrease in SOC, available N and available P can be seen in Table 4. Previous studies also suggested that that the concentration of organic C, available N and P was lower in the disturbed sites, which indicates the lower microbial activities [34].
do not show significant changes with the year, but one-way ANOVA indicates a significant temporal variation in SOCEA, SOCW, available nitrogen, and available phosphorous. Soil salinity significantly increased from 2019 onwards after the Fani, Bulbul, Amphan and Yaas cyclones hit the region. However, a slight increase in salinity was there already in 2018. The post-hoc test is also indicating the similarity between the salinity in year 2017 and 2018 and between 2019 and 2021. Soil salinity is also dependent on tidal water salinity, which is evident from Figure 6. The concentrations of available N and available P were recorded as being approximately 1.5 times lower in 2021 in comparison to 2021 ( Table 4). The gradual decrease in SOC W is prone to interferences due to manual error, while titration or changing strength of K 2 Cr 2 O 7 during analysis can influence the outcome. An elemental analyzer has been widely used in the carbon estimation of coastal saline or hypersaline soil [25][26][27]. Although a positive relationship was observed between the concentration of SOC W and SOC EA with a 'R 2 ' value of 0.9 (Figure 7), a lower SOC W value has been observed in comparison with the SOC EA value, which may be due to interferences or manual errors while conducting the wet digestion analysis. As evident from work of Passos et al. and Nóbrega et al., C estimation with an elemental analyzer has better accuracy in case of coastal and mangrove soil [26,27].
SOCW is prone to interferences due to manual error, while titration or changing strength of K2Cr2O7 during analysis can influence the outcome. An elemental analyzer has been widely used in the carbon estimation of coastal saline or hypersaline soil [25][26][27] Although a positive relationship was observed between the concentration of SOCW and SOCEA with a 'R 2′ value of 0.9 (Figure 7), a lower SOCW value has been observed in com parison with the SOCEA value, which may be due to interferences or manual errors while conducting the wet digestion analysis. As evident from work of Passos et al. and Nóbrega et al., C estimation with an elemental analyzer has better accuracy in case of coastal and mangrove soil [26,27]. Blue carbon sequestration is one of the key services of mangrove forests [1]. The changes in biodiversity patterns do impact the carbon pool, as evident from Figure 8. I indicates a steady decline in BOCS from 2017 to 2021 with a negative correlation with years and an 'R 2′ value of 0.88. A rise in soil salinity or electrical conductivity is coinciding with cyclonic disturbances in the area. Higher salinity can change the nitrogen dynamics of soil, which may be the reason for lower N over the years in the present study [51]. cyclonic disturbances in the area. Higher salinity can change the nitrogen dynamics of soil, which may be the reason for lower N over the years in the present study [51]. The work of Chambers et al. also highlights an increase in dissolved organic carbon (DOC) in the mangrove soil pore water due to increased inundation and salinity rise [52]. An increase in DOC means more allochthonous transport to a distant location and less accretion of carbon in the mudflat. A. marina residing along the fringe of the forest has a leaf carbon content of 39.1%, while species such as B. gymnorrhiza (45.6% C), E. agallocha (44.8% C), R. mucronata (41.2% C), and X. granutum (40.7% C) have a higher carbon content in their leaves [53,54]. The mangrove ecosystem is evergreen. Most of the energetics of the ecosystem is managed through the detritus food chain and the decomposition of nutrients in the leaves. With low microbial diversity, it takes more time for the waxy mangrove leaves to be decomposed in the soil and release carbon and nutrients. The salinity-tolerant A. marina has lower carbon content in the leaves, resulting in low carbon sequestration in the soil in comparison with a diverse ecosystem with other species. Salinity also impacts the biodiversity of microbes and hence can alter the time required for releasing nutrients back to the soil, including C, in the case of an increased salinity increment.
This study indicated a lower concentration of the major macromolecules in soil, proving a degradation of ecological functions and services that can be linked with changes in mangrove community ecology in response to natural disaster-induced salinity intrusions. The work of Mace et al. also indicates a relationship between ecosystem structure, function, and services [55].

Cyclones, Salinity Intrusions and Impact on Global Blue Carbon Sink
Global greenhouse gas emissions are responsible for trapping heat in the atmosphere. In a positive feedback loop, this entrapment of heat results in a rise in the overall oceanic temperature and accelerates the surface evaporation process. The accelerated evaporation cycle in the tropical oceans has increased the frequency and intensity of cyclones in the last decade. In the Indian Ocean and nearby oceanic regimes, the temperature is warming the oceans, disrupting the general oceanic circulations, and increasing the risk of the development of tropical disturbances or cyclones [56][57][58]. Similarly, the Bay of Bengal has also been impacted by the global problems of climate change and global warming. El Niño The work of Chambers et al. also highlights an increase in dissolved organic carbon (DOC) in the mangrove soil pore water due to increased inundation and salinity rise [52]. An increase in DOC means more allochthonous transport to a distant location and less accretion of carbon in the mudflat.
A. marina residing along the fringe of the forest has a leaf carbon content of 39.1%, while species such as B. gymnorrhiza (45.6% C), E. agallocha (44.8% C), R. mucronata (41.2% C), and X. granutum (40.7% C) have a higher carbon content in their leaves [53,54]. The mangrove ecosystem is evergreen. Most of the energetics of the ecosystem is managed through the detritus food chain and the decomposition of nutrients in the leaves. With low microbial diversity, it takes more time for the waxy mangrove leaves to be decomposed in the soil and release carbon and nutrients. The salinity-tolerant A. marina has lower carbon content in the leaves, resulting in low carbon sequestration in the soil in comparison with a diverse ecosystem with other species. Salinity also impacts the biodiversity of microbes and hence can alter the time required for releasing nutrients back to the soil, including C, in the case of an increased salinity increment.
This study indicated a lower concentration of the major macromolecules in soil, proving a degradation of ecological functions and services that can be linked with changes in mangrove community ecology in response to natural disaster-induced salinity intrusions. The work of Mace et al. also indicates a relationship between ecosystem structure, function, and services [55].

Cyclones, Salinity Intrusions and Impact on Global Blue Carbon Sink
Global greenhouse gas emissions are responsible for trapping heat in the atmosphere. In a positive feedback loop, this entrapment of heat results in a rise in the overall oceanic temperature and accelerates the surface evaporation process. The accelerated evaporation cycle in the tropical oceans has increased the frequency and intensity of cyclones in the last decade. In the Indian Ocean and nearby oceanic regimes, the temperature is warming the oceans, disrupting the general oceanic circulations, and increasing the risk of the development of tropical disturbances or cyclones [56][57][58]. Similarly, the Bay of Bengal has also been impacted by the global problems of climate change and global warming. El Niño Southern Oscillation (ENSO) has been proven to also influence the Indian Ocean dipole and influence cyclogenesis throughout the southern oceans [58]. Kang and Elsner estimated an average global cyclone intensity of 1.3 m s −2 over the last 30 years due to global warming impacts [2]. Although predictive models are present, it is difficult to comprehend the impact of climate change or global warming at regional levels, particularly in coastal ecosystems. Saltwater flows in through the tides, while the freshwater flow is maintained due to riverine transport in a deltaic or estuarine system. Cyclones destabilize ecosystems. As per previous research, cyclones and natural disasters are one of the main reasons behind declining mangrove cover across the globe [59]. Due to the sediment load of the Ganga, Brahmaputra, and Meghna River systems, along with the earth's movement some 200 years ago, the Bengal Basin is tilting eastward. Recent trends indicate a reduction of freshwater from Himalayan rivers coupled with an increase in sea level impacting the Sundarban mangrove ecosystem [60]. This geological reason can also be responsible for increased salinity in the delta by facilitating seawater intrusion. Around 123 countries located in the tropical and subtropical regions of the globe have mangrove vegetation along their coastlines, estuaries, or deltas [61]. Mangrove ecosystems are regarded as one of the best natural models to study the impact of climate change and global warming-related vagaries such as rises in salinity and sea level and cyclone impacts, and they give a glimpse into understanding the impact of these global problems at regional levels [62]. Hence, research on understanding climate change and natural disaster impacts on mangrove biodiversity also paves the way to understanding the impact of these global problems on different ecosystems across the world. Blue carbon sequestration by mangroves and sea grasses plays a crucial role in the fight against climate change (SDG 13), and a healthy ecosystem conserves the species diversity in the wetlands (SDG 14). It is observed by scientists that mangroves have the potential to withstand the vagaries of climate change and sea level rise according to Menéndez et al., but this resilience may have a limit and will eventually result in irreparable changes in mangrove community structure and composition, as is also evident from this study [63].

Conservation Consideration and Recommendation
The Indian part of the Sundarbans comprises 102 islands, out of which 54 are inhabited. The rest of the islands are under mangrove cover. The river Ganges carries down 262million tons of sediment per year, but research indicates that only 7% is diverted through its southern distributaries [64]. The low discharge restricts this sediment in the mangrove forested islands. It creates a unique mud island diurnally washed by tides and soil consolidated by mangrove roots. It is recognized as the 27th Ramsar wetland site of India in 2019, making it the largest wetland in the nation. Ecotourism is booming in the region. There has been a rise in ecotourism income (INR) from 1,508,000 in 2003-2004 to 11,772,000 between 2012 and 2013 [65].
Mangroves hold sediment, and their long-extended rhizosphere helps in consolidating sediment during storm surges and cyclonic disturbances. As per previous studies, mangroves have acted as a bio-shield against natural disasters along the coastlines by acting as a wind break, consolidating soil, reducing wave force, and protecting the coastlines [47,59,63]. Although recent research has focused on the importance of blue carbon and disaster resilience ecosystem service of mangroves for mitigating climate change threats, the ecological investigations on the climate change impacts on mangrove communities are less explored [66,67]. Friess et al. opined that a rise in temperature can destabilize the ecological function and services of global mangrove ecosystem, substantiating a research focus on the climate change-mangrove community ecology nexus [67]. Segaran et al. have indicated that the publication trends in this domain mostly focus on carbon sequestration, nature-based solutions, and disaster resilience, but few works focus on ecological changes in the mangrove community due to climate change vagaries [66]. This necessitates the need for focused ecological research to formulate effective solutions to tackle the issue.
The solution lies in the protection of the areas that are outside the reserve forest and the restoration of mudflats with multispecies mangroves [1,68]. These meta-populations also act as a natural repository of propagules and are fundamental for a sustainable ecosystem by facilitating the colonization of newly formed or degraded mudflats [69,70]. Thus, the changes in community structure can also reduce the availability of propagules of salinity-sensitive mangrove species; hence, they will impact the overall biodiversity of the ecosystem in the unforeseen future.
This research sheds light on the impact of salinity rises in surface water on mangrove ecological composition and services, which can hinder achieving the goals of SDG 13 and 14. Hence, ecological engineers need to comprehend the role of mangrove biodiversity in its carbon sequestration potential. According to this case study, a 45% reduction in the blue carbon pool is observed due to a 47.6% increase in surface water salinity from 2017 to 2021. This study highlights the negative impact of salinity intrusion on blue carbon sequestration as well as mangrove biodiversity. A similar situation may exist at mangrove patches scattered across the globe and requires a similar investigation.

Conclusions
The research proved the increase in surface salinity trend due to cyclones in the innercentral part of the Indian Sundarbans. This is resulting in the dominance of salinity-tolerant mangrove species while causing localized extinctions of salinity-sensitive members. It is also negatively impacting the health of the ecosystem and can disrupt the ecological services of the ecosystem. The changes in community composition are associated with a reduction in organic blue carbon sequestration by mangrove soil along with reduction in other major nutrients such as available nitrogen and phosphorus. This study indicates a deterioration of ecological health, diversity, changes in mangrove species zonation in response to tidal gradient and functions due to a sudden change in the surface water salinity profile due to natural disturbances. Other causes of salinity rise, such as sea level rise, climate change, and a reduction in freshwater flow in the delta, can also contribute to increased salinity, which needs further investigation. Future research focus needs to be directed toward establishing the connection between climate change, global warming, cyclogenesis, and changes in estuarine salinity profiles to formulate any future conservation policies for coastal habitats. Conservation plans also need to be put in place to manage this localized extinction, as it can destabilize the entire Sundarban mangrove ecosystem.