Study area description
The study was conducted in the Bangladesh Sundarbans that lies in the coast of Bay of Bengal (21o 30´– 22o 30´ N, 89o 00´– 89o 55´E; Fig. 1). Of its total area (6017 km2), about 69 % is land and the remaining areas are waterbodies, such as rivers, small channels and creeks (Wahid et al. 2007). It is inundated regularly twice in a day by high tides and the inundation level varies due to varying elevation and tide levels. The regional hydrology and salinity are regulated by the water enter into the Sundarbans for the sea (the Bay of Bengal) and the Gorai river. The salinity varies spatially and temporally within the ecosystem (Siddique et al. 2021). Based on the spatial salinity variation, the ecosystem is divided into three salinity zones (Sarker et al. 2019), such as low saline zone (LSZ), mid saline zone (MSZ) and high saline zone (HSZ). Salinity also varies substantially between the seasons. For example, salinity drops due to higher monsoonal (June – September) precipitation and freshwater flows from the upstream rivers (Chowdhury et al. 2016a). Precipitation in the study areas decreases during winter (December – February) followed by pre-monsoon (March – May), and consequently increases the salinity, particularly in pre-monsoon (Siddique et al. 2021). During the post-monsoon (October – November), salinity is also low, despite lower precipitation and freshwater flows.
The Sundarbans is rich in biodiversity and home to a diverse range of tree species, including the IUCN listed globally endangered tree species H. fomes (Sarker et al. 2019). The Sundarbans is a part of the Indo-Malaysian mangrove system (Reef and Lovelock 2015). H. fomes is the most dominant and flagship tree species, while E. agallocha is the second most dominant tree species in the Sundarbans (Chowdhury et al. 2016a, 2023a). Other important tree species in this mangrove are S. apetala, X. moluccensis, Ceriops decandra (Griff.) Ding Hou, Amoora cucullate Roxb., A. officinalis, Bruguiera sp., Cynometra ramiflora L., Aegiceras corniculatum (L.) Blanco, Lumnitzera racemose Willd., Rhizophora mucronata Lamk. and X. granatum J. Koenig. Besides a wider range of climbers, creepers, grasses and palms are available. Due to different ecological characters, tree species forms different associations in the Sundarbans which vary across the salinity gradients (Iftekhar and Snegar 2008; Sarker et al. 2016).
Tree measurements and wood sample collection
In this study, five dominant tree species with varying diameters from two distinct shade tolerance categories (Siddiqi 2001), light-demanding (S. apetala, A. officinalis and E. agallocha), and shade tolerant (X. moluccensis and H. fomes), were selected for sampling. In the Sundarbans, each species was sampled in three distinct salinity zones (e.g., LSZ, MSZ, and HSZ). In each salinity zone, 10 to 15 trees were selected from each species for wood core collection. Diseased, leaning trees as well as trees near the edges of river, creeks, and trails were excluded from sampling. Tree height (m) was measured using a clinometer (Sunnto PM5-360), and diameter at breast height (DBH) over bark (cm) of each tree was measured by using a diameter tape. Slenderness ratio of tree is defined as the ratio of total height to diameter at 1.3 m above ground (Zhang et al. 2020) was calculated for each sample tree. A total of 150 (10 x 3 x 5) trees were cored at breast height (1.3 m above from the ground level) with an increment borer (4.15 mm diameter, Haglöf, Sweden), and the extracted radial cores were examined to confirm the presence of the pith. To ensure accurate measurement of green volume, wood cores were immediately placed in an airtight container to prevent them from drying out.
(Basic) wood density measurement
The wood cores were cut into 1 cm long pieces from pith to bark. Inner wood was defined as an arithmetic mean of the innermost wood of two consecutive pieces (1cm/piece) from pith and outer wood as an arithmetic mean of the outermost wood of 2 consecutive pieces (1 cm/piece) from bark (Bastin et al. 2015). The average wood density was defined as an arithmetic mean of all wood segments from pith to bark. The volume of each sample was calculated by measuring the fresh sample with a slide caliper, i.e., volume = πr2h, where r is radius of the increment core and h is the sample length (Chave et al. 2009). Following that, the oven-dry mass of each wood segment was measured using an oven for 24 hours at fixed temperature of 105° C (Chowdhury et al. 2009). Subsequently, basic wood density (g cm-3) was calculated as the ratio of oven-dry mass to green volume.
Above ground tree biomass estimation
Above ground biomass (AGB) was estimated using previously developed multi-species allometric model for the Sundarbans, as follows (Mahmood et al. 2019):
Ln (TAGB) = a + b Ln (DBH) + c Ln (H) + d Ln (WD) ………………………………. (Eq. 1)
where, a = - 6.7189,
b = 2.1634, c = 0.3752,
d = 0.6895,
AGB = Total Above-ground biomass (kg),
DBH = DBH (cm),
H = Tree height (m),
WD = Wood density (kg m-3), and
a, b, c, and d are the constants.2019) for accurate AGB estimations, with values of 1.0119, 1.0222, 1.0119, 1.02 and 1.0103 for S. apetala, A. officinalis, E. agallocha, X. moluccensis and H. fomes, respectively.
Deviations in AGB estimations
AGB calculations differed because different values of wood density were used. The first AGB, for example, was calculated using basic wood density from a repository (Kattge et al. 2011). The second AGB measured in this study was for outer wood basic density, and the third was for inner wood basic density (Bastin et al. 2015). Finally, the reference AGB was calculated using the average basic wood density, which was calculated as mean basic wood density of all segments in each tree from pith to bark. The AGB differences were calculated as follows:
AGB difference repository (%) = 100 * (AGB repository - AGB reference)/ AGB reference …..…. (Eq. 2) AGB difference outer (%) = 100 * (AGB outer - AGB reference)/ AGB reference ……………… (Eq. 3)
AGB difference inner (%) = 100 * (AGB inner - AGB reference)/ AGB reference ………..……….(Eq. 4)
Statistical analysis
The Shapiro-Wilk test was used to check normality of the data. Regression analyses were used to examine the radial wood density trend. The variation in wood density among the species was investigated using the Kruskal-Wallis test, followed by a post-hoc (Dunn's) test. The Wilcoxson test was used to examine differences in average wood density between light and shade tolerant categories. Regression analyses were used to examine the variation in wood density in relation to the slenderness ratio, as well as the relationship between inner and outer wood density. The R platform, version 4.0.5 (R Core Team, 2022) was used for all analyses.