The present study investigates the effect of different habitat attributes on brachyuran crab diversity in two different study sites in the Sundarban mangrove, India. The two sites differ in the level of anthropogenic intrusion and in the age of the mangrove forest. Seasonal changes in the environment and in brachyuran faunal abundance were recorded for three years. Species composition varied between the two habitats irrespective of season. The habitat heterogeneity and the recorded crab community was analysed by several univariate and multivariate statistical techniques. The newly replanted mangrove site showed lesser diversity than the natural one. Ocypodid crabs, mainly Uca rosea, dominated both study sites, whereas Uca triangularis was totally absent from the replanted site. Canonical correspondence analysis showed that the total acidity, total alkalinity, pH content of water, total dissolved solids, inorganic phosphate content of water, soil specific gravity, soil density and the physical constructions of the habitat play a crucial role in moderating the crab community structure. This study reveals that brachyuran crab diversity can be used as a potential indicator of the alterations of mangrove habitats.
RESUMEN
El presente estudio investiga el efecto de las diferentes características de hábitat sobre la diversidad de los cangrejos braquiuros en dos lugares diferentes del manglar Sundarban, India. Los dos sitios difieren en el nivel de intrusión antropogénica, así como en la edad del bosque de manglar. Se registraron cambios estacionales en el medio y en la abundancia de la fauna de braquiuros durante tres años. A pesar de las diferencias estacionales, la composición de especies difiere en los dos hábitats. Se analizó la heterogeneidad del hábitat y la comunidad de cangrejos mediante técnicas estadísticas univariantes y multivariantes. La zona de manglares replantada recientemente mostró menor diversidad que la zona natural. Los cangrejos ocipódidos, principalmente Uca rosea, dominaron en ambas zonas de estudio, mientras que Uca triangularis estuvo totalmente ausente en la zona replantada. El análisis de Correspondencia Canónica (CCA) mostró que la acidez total, la alcalinidad total, el contenido de pH del agua, los sólidos totales disueltos (TDS), el contenido de fosfato inorgánico del agua, el peso específico del suelo, la densidad del suelo, junto con las construcciones físicas del hábitat desempeñan un papel fundamental en la estructura de la comunidad de cangrejos. Este estudio revela que la diversidad de los cangrejos braquiuros puede ser utilizada como un potencial indicador de las alteraciones de los hábitats de manglares.
KEYWORDScrab diversitycanonical correspondence analysisedaphic qualitieshabitat heterogeneityhydrological parametersmangrovesPALABRAS CLAVEdiversidad de cangrejosanálisis de correspondencia canónicacualidades edáficasheterogeneidad del hábitatparámetros hidrológicosmanglaresINTRODUCTION
Mangroves are the most productive and dynamic ecosystems in tropical and subtropical intertidal areas (Robertson and Blaber 1992, Kathiresan and Bingham 2001). The Sundarban mangrove (21°30’ to 22°40’N, 88°05’ to 89°55’E) in India is the single largest mangrove in the world and has the greatest variability in the hydrological features both at a spatial and temporal scale. This mangrove was declared as a world heritage site by IUCN in 1987 for its rich faunal diversity. The topography and quality of the substratum, the salinity variations and the interactions among these factors result in a high habitat heterogeneity in the mangrove, ensuring a similarly diverse faunal community (Gopal and Chauhan 2006).
The habitat parameters (physical structure and edaphic and hydrological parameters) normally have a strong influence on the diversity and abundance of the associated organisms (MacArthur and MacArthur 1961, Lawton 1983, Ellner 2001). The complex system of prop-roots and pneumatophores in mangrove forests, in association with an often dense layer of leaf litter and detritus on the substratum, provides generous cover from predators for many invertebrate species (Acosta and Butler 1997, Primavera 1997, Laegdsgaard and Johnson 2001). In addition, variation in the density and structure of mangrove roots produces an array of shelters (Primavera 1997, Macia et al. 2001). Such variation normally affects the benthic fauna in a species-specific manner (Kon et al. 2009). Although some studies have examined effects of mangrove loss on benthic and epibenthic animals (Levings et al. 1994), few studies have examined how smaller-scale modification of habitat complexity affects the biodiversity and abundance of mangrove-associated fauna (reviewed by Hatcher et al. 1989).
Among all benthic macro fauna in the intertidal zone, brachyurans are the most significant because of their great diversity, comprising about 1271 genera and 6793 species (Ng et al. 2008). Brachyurans are also the most important macro-invertebrates as they make up 80% of the total faunal biomass (Golley et al. 1962) and attain high density (Macintosh 1984). Their role in the ecological functioning of the mangrove ecosystem (Lee 1998) through litter turnover (Robertson 1986) and bioturbation of soil (Smith et al. 1991) provides crucial food sources for coastal fisheries (Macintosh 1984). For all these reasons, Smith et al. (1991) recommended crabs as possible keystone species in the mangrove habitat. Except for the distribution and abundance of some of the dominant species of the Sundarban mangrove, no detailed information is available on the crab community. Therefore, the objectives of the present study were (i) to compare the community structure of brachyuran crabs at two mangrove sites with different habitat attributes because each habitat functions according to different physical, edaphic and hydrological properties, and (ii) to evaluate whether brachyuran crabs can be used as an indicator of habitat status in mangrove ecosystem.
MATERIALS AND METHODSStudy sites
The present study was conducted in the central part of Sundarban mangrove, namely at Jhorkhali Island (22°01’14.55”N, 88°41’30.39”E) located between the two large rivers Matla and Bidya. Two study sites were selected; site-1 along the mudflats of the Matla River and site-2 along the Boro Herobhanga creek (Fig. 1). Both sites experience a subtropical monsoonal climate with an annual rainfall of about 1600-1800 mm (Manna et al. 2010). At site 1, the mangrove forest is truly natural without much anthropogenic intrusion, and with four dominant arboreal species (The dominant mangrove species are Avicennia marina, A. officinalis, Rhizophora apiculata and Excoecaria agallocha). Site 2 has a docking site for local ferry boat services and a mangrove replantation plot dating back to the 1990s, with a mixed vegetation composition (R. apiculata, A. marina, A. officinalis, Bruguiera cylindrica, Haritiera fomes, Ceriops sp., Sonneratia alba, etc.). The replantation plot serves as a source for replantation programmes throughout Sundarban.
Map showing the Study sites as black dots in Jhorkhali Island, Sundarban.
The sampling area of each site was 3 km2, within which three 1-km transects were laid (two along the two sides of the square and one along the middle). About sixty 2×2 m2 quadrats were placed randomly along the transects for the collection of the data. In all the quadrats physical, hydrological and edaphic parameters were thoroughly analysed twice in a month for three seasons (pre-monsoon, monsoon and post-monsoon) in three consecutive years. For each of the variables, 10 replicates were made during each sampling. The number of the pneumatophores, senescent and decaying leaves, shrub and herb density were recorded in each quadrat. Percentage of canopy cover was determined using a locally built densiometer, while canopy height and shrub height were determined using measuring tape.
Soil sampling design
Soil samples were collected from each quadrat with the help of a soil corer and carried to the laboratory in plastic zipper pouches. In the laboratory, the initial weight of the soil samples was taken, the samples were oven dried (60°C for 48 hours) and the final weight of the samples was taken for the measurement of soil moisture. Soil pH and soil specific gravity were determined following the method of Brower et al. (1998). Soil organic carbon and organic matter were assessed using a modified Walkey and Black method (Nelson and Sommers 1996).
Water sampling design
At both sites, subsurface water samples were collected twice every month simultaneously with crab sampling during high tide from a depth of 20-30 cm. Overall 13 hydrological variables were considered. Air temperature and water temperature were measured instantly with a mercury thermometer with 0.1°C graduations. Dissolved oxygen was measured at the site using a Lutton, DO-5509 Dissolved Oxygen Meter. Total Dissolved Solute and pH were also recorded in the field using portable testers (Eutech instrument, cyber scan). Water samples were collected in amber-coloured, labelled glass bottles for analysing the rest of the hydrological parameters. Samples were preserved at 4°C without freezing and taken to the laboratory within 3 h of collection for analysis as per standard methods (APHA 1998). In the laboratory, total alkalinity was measured following the acid titrimetric method, using H2SO4 as titrant and methyl orange as an indicator. Salinity was determined in practical salinity units (PSU) by the Knudsen argentometric method (Strickland and Parsons 1972). Total acidity was measured following the alkaline titrimetric method using NaOH as the titrant and phenolphthalein as the indicator. Total hardness was assayed by the EDTA titrimetric method. Inorganic phosphate, reactive silicate, nitrate-nitrogen and nitrite-nitrogen were measured using a UV spectrophotometer following standard sea water analysis methods (Grasshoff 1983, Grasshoff et al. 1983, Strickland and Parsons 1972).
Crab sampling design
During both low spring tide and low neap tide crabs were collected quantitatively by performing three individual time-based samplings in each of the 2×2 m2 quadrats. The collection was based on the assumption that brachyuran crabs generally do not change significantly over the time scale (Ashton et al. 2003a). Precisely, three independent time-based samplings were employed in each quadrat. One sampling (considered as one replicate) represented one person collecting crabs for 15 min. This method, first devised by Ashton et al. (2003a), was found to offer a good representative sample of the crab population residing at each site. Crabs were collected mostly by hand, sometimes with the help of a trowel. To reduce the bias in favour of slow-moving species, great efforts were given to catching the larger and fast-moving elusive species. After collection, the crabs were kept in 70% alcohol inside labelled plastic containers and carefully transported back to the laboratory. In the laboratory the crabs were identified and preserved in 4% formalin for future identification.
Data analysis
As the study focuses on the temporal factors (seasons and years) at two different study sites, one-way analysis of variance (ANOVA) at 5% level of significance was used to test for significance for the environmental parameters followed by Duncan’s post hoc tests to check for differences in the environmental factors between every pair of seasons and year. These parts of the statistical analyses were carried out by means of SPSS (Statistical Package for Social Sciences, version 17.0 Norušis 2000).
A total of 2055 individuals were collected during the period of three years. Total crab abundance was calculated at two different sites for three seasons in each year: pre-monsoon, monsoon and post-monsoon. From the sampled data, Pielou evenness (J) (Pielou 1969) and species dominance (D) index (Berger and Parker 1970) were calculated for two sites in three seasons each year. Univariate measures were transformed using square-root to normality, if necessary. The similarity in seasonal species composition at each site was studied by calculating the Bray-Curtis coefficient based on the square-root-transformed species abundance data. The result was displayed by non-metric multidimensional scaling (nMDS) plot (Clarke 1993). For differences between the study sites, formal significance tests for similarity were performed using the one-way ANOSIM permutation test on the crab similarity matrix. The crab species which are contributing to the dissimilarities between the two study sites were investigated using the similarity percentage (SIMPER) (Clarke and Warwick 1994). All mentioned statistical treatments were performed using PRIMER-E Software (v. 6).
In order to assess the interaction between different hydrological, edaphic and physical parameters and the crab species assemblage, unimodal distribution of species was used to explain the abundance of species with environmental variables (water temperature, pH, salinity, total acidity, total alkalinity, total hardness, inorganic phosphate, nitrate-nitrogen and nitrite-nitrogen, soil organic carbon, soil organic matter, soil pH, soil density, soil specific gravity, % of canopy cover, canopy height, pneumatophore density, % of leaf litter, herb density, etc.). Canonical correspondence analysis (CCA) was employed with biplot type of scaling on inter-species distance and down weighting of rare species. The species abundance data were square-root transformed. A test of significance of the first canonical axis and all canonical axes was made using the Monte Carlo test with 499 permutations under the reduced model (Ter Braak 1986, 1988, Jongman et al. 1995). The relationships of crab species abundance with environmental variables were presented in ordination diagrams (biplots) for the first two axes using CANOCO 4.5 software package.
RESULTS
Seasonal and annual trends in hydrological and edaphic parameters for the two sites are given in Table 1 and Table 2, respectively. Most of the parameters, and specially inorganic phosphate, nitrite-nitrogen, nitrate-nitrogen, salinity of water and organic carbon along with organic matter and moisture content of soil, differed significantly (p<0.05) between the two sites. At site 1, all the hydrological parameters and most of the edaphic factors varied significantly (p<0.05) between three seasons of a year and also between years of a particular season. At site 2, hydrological and edaphic parameters of the habitat differed (p<0.05) mostly in all the seasons and in all three years of a particular season with the exception of soil density in the pre-monsoon and monsoon seasons. Most of the environmental parameters were correlated with each other according to Spearman’s correlation coefficient. Soil moisture was negatively correlated with number of senescent or decaying leaves (r=–0.61, p<0.01), whereas the number of fallen leaves was positively correlated with pneumatophore density, canopy height and shrub height (r= 0.764, p<0.01, r=0.639, p<0.01 and r=0.625, p<0.01 respectively). Soil pH was positively correlated with soil moisture (r=0.522, p<0.01) but negatively correlated with soil organic carbon and soil organic matter (r=–0.335, p<0.05 for both the cases).
Hydrological parameters of study sites in the Indian Sundarbans. SAL, salinity (PSU); ACI, total acidity (mg L–1); ALK, total alkalinity (mg L–1); TH, total hardness (mg L–1); WT, water temperature (°C); AT, air temperature (°C); pH; TDS, total dissolved solids (ppm); inorganic phosphate (µM L–1); NO3-N, nitrate-nitrogen (µM L–1); NO2-N, nitrite-nitrogen (µM L–1). Different letters indicate significant differences at 0.05% level of significance.
SAL
ACI
ALK
TH
WT
AT
pH
TDS
Ip
NO3-N
NO2-N
Year 1
PMo
St1
27.42a
±0.32
18.15a1
±0.10
139.18a2
±0.11
120.24a3
±0.04
30.38a4
±0.06
33.33a5
±0.08
8.20a6
±0.04
1448.48a7
±0.27
0.26a8
±0.00
0.93a9
±0.03
0.48a10
±0.00
St2
25.91a’ ±0.23
31.87a1’
±0.35
165.00a2’
±1.86
181.14a3’
±3.08
32.58a4’
±0.09
36.35a5’
±0.06
7.63a6
±0.06
1437.22a7’
±0.90
0.07a8’
±0.00
0.87a9’
±0.01
0.18a10’
±0.01
Mo
St1
15.49b ±0.88
44.88b1
±0.23
117.99b2
±0.34
56.43b3
±0.14
27.52b4
±0.19
31.52b5
±0.12
7.95b6
±0.03
1675.24b7
±3.00
2.30a8
±0.07
2.55a9
±0.07
0.57a10
±0.01
St2
21.27b’ ±0.61
63.06b1
±0.53
160.52b2’
±0.25
76.16b3’
±0.85
31.60b4’
±0.09
35.28b5
±0.05
7.45b6’
±0.03
1653.94b7
±1.57
1.05b8’
±0.01
1.62b9’
±0.12
0.18b10’
±0.01
PoMo
St1
23.81c ±1.01
36.85c1
±0.23
125.68c2
±0.47
62.31c3
±0.08
22.65c4
±0.09
26.45c5
±0.13
7.17c6
±0.01
2334.63c7
±1.59
1.35a8
±0.01
3.37a9
±0.07
0.86a10
±0.01
St2
27.99c’ ±0.21
57.97c1
±0.33
195.77c2
±1.26
133.50c3’
±16.3
26.45c4’
±0.06
30.20c5’
±0.07
6.78c6’
±0.02
2434.52c7
±10.24
0.94c8’
±0.01
1.37c9’
±0.03
0.28c10’
±0.01
Year 2
PMo
St1
31.35a ±0.42
27.99d1
±0.26
112.58d2
±0.24
43.50d3
±1.19
33.08d4
±0.81
36.48a5
±0.03
7.30d6
±0.07
1161.62d7
±3.25
1.45a8
±0.01
3.20d9
±0.00
0.84d10
±0.01
St2
26.93a’
±0.33
64.01d1’
±0.60
156.81d2’
±1.15
91.27d3’
±0.98
36.63a4’
±0.10
38.45a5’
±0.06
6.78d6’
±0.09
1133.23d7’
±0.71
0.48a8’
±0.00
2.17d9’
±0.01
0.45b10’
±0.01
Mo
St1
13.99e ±0.50
57.92b1
±0.21
100.71b2
±0.43
38.48b3
±0.41
28.35e4
±0.03
34.23e5
±0.10
7.20b6
±0.04
1345.91e7
±1.31
3.75a8
±0.03
3.08e9a9
±0.03
1.34d10a10
±0.01
St2
21.99e’ ±0.26
95.96b1’b1
±0.38
118.42e2’
±0.48
64.50b3’
±0.47
32.23e4’
±0.06
36.28e5
±0.09
6.40e6’
±0.04
1338.13e7
±0.24
0.88a8’b8’
±0.00
1.64e9’b9’
±0.01
0.66b10’
±0.02
PoMo
St1
15.03f
±0.47
56.96c1
±0.37
129.75f2
±1.60
27.12f3
±0.47
26.58f4
±0.08
29.23f5
±0.06
7.28f6
±0.05
2489.97f7
±4.50
3.45a8
±0.04
3.53f9a9
±0.03
0.97d10
±0.01
St2
28.28f’ ±0.22
86.21c1
±0.29
135.59f2’
±1.12
88.52f3’
±0.35
29.58f4
±0.14
34.85f5’
±0.24
6.75f6’
±0.03
2440.80f7
±0.25
1.19a8’c8’
±0.00
2.45f9’c9’
±0.02
0.53b10’c10’
±0.00
Year 3
PMo
St1
27.80a
±0.47
63.16g1
±0.43
88.27g2
±0.06
22.71g3
±0.28
29.50a4’
±0.04
34.38a5
±0.10
7.21g6
±0.00
1557.02g7
±0.47
2.84a8
±0.02
3.91g9
±0.02
1.70g10
±0.15
St2
32.00a’ ±0.27
77.99g1
±0.18
144.62g2
±1.14
107.97g3’
±1.44
33.28a4’
±0.07
32.28a5’
±0.09
7.04g6
±0.07
1180.24g7
±7.35
0.85a8’
±0.00
1.69g9’
±0.09
0.63g10’
±0.00
Mo
St1
9.96h
±0.43
86.66b1
±0.24
90.26b2
±0.04
32.29b3
±0.13
27.65h4
±0.12
32.33h5
±0.11
7.53b6
±0.11
1473.10h7
±5.32
4.07h8a8
±0.01
8.46h9a9
±0.07
2.78h10a10
±0.07
St2
12.30h’ ±0.11
94.43g1
±0.32
122.06b2
±1.48
110.83b3’
±3.49
31.50h4’
±0.10
34.83h5
±0.36
6.73h6’
±0.09
1348.62h7’
±3.58
0.98g8’b8’
±0.00
3.25h9’b9’
±0.00
1.45b10’
±0.18
PoMo
St1
15.59i
±0.40
81.16c1
±0.41
101.57i2
±0.26
18.26i3
±0.02
23.33i4
±0.11
28.05i5
±0.29
7.13i6
±0.03
1520.93i7
±2.89
2.84i8
±0.02
8.50i9a9
±0.05
3.45i10
±0.09
St2
21.36i’ ±0.04
84.10g1
±0.30
118.49i2’
±0.09
114.80i3’
±1.29
29.08i4
±0.30
33.58i5’
±0.09
6.45i6’
±0.13
1144.96i7
±1.60
1.27g8’c8’
±0.01
5.35i9’c9’
±0.11
1.39c10’
±0.05
Edaphic parameters of study sites in the Indian Sundarbans; soil organic carbon (%); soil organic matter (%); soil moisture (%); soil pH; soil density and soil specific gravity. Different letters indicate significant differences at 0.05% level of significance.
Soil organic carbon
Soil organic matter
Soil moisture
Soil pH
Soil density
Soil specific gravity
Year 1
PMo
St1
1.00a1±0.00
1.73a2±0.00
42.47a3±0.17
7.58a4±0.11
2.73a5±0.07
1.39a6±0.04
St2
0.44a1’±0.00
0.76a2’±0.00
46.91a3’±0.45
7.75a4±0.01
2.58a5’±0.05
1.48a6’±0.08
Mo
St1
1.03b1±0.00
1.77a2±0.00
53.10b3±3.80
7.57b4±0.10
2.90b5±0.02
1.46b6±0.02
St2
0.76a1’±0.01
1.30a2’±0.02
59.54b3±0.42
7.47b4a4±0.02
1.69b5’±0.08
1.23b6±0.00
PoMo
St1
1.04c1±0.00
1.79a2±0.00
43.73c3±0.58
7.54c4±0.03
2.98c5±0.05
1.36c6±0.02
St2
0.94a1’±0.01
1.62a2’±0.01
34.87c3’±0.29
6.71c4’a4±0.08
1.73c5’±0.01
1.09c6’±0.03
Year 2
PMo
St1
0.95d1±0.02
1.63d2±0.02
40.44a3±3.80
7.67d4±0.04
1.80d5±0.13
1.24d6±0.00
St2
0.54d1’±0.10
0.92d2’±0.17
48.41d3’±0.40
7.88d4±0.04
2.67a5’±0.15
1.24a6’±0.01
Mo
St1
1.31e1±0.08
2.26e2±0.15
53.65e3b3±0.67
7.67b4±0.01
2.35b5±0.01
1.05e6b6±0.00
St2
0.88a1’±0.01
1.52a2’e2±0.01
54.89e3b3±0.24
7.63b4±0.02
1.49b5±0.09
1.17e6’b6±0.01
PoMo
St1
1.36a1’±0.01
2.35f2±0.02
48.30f3±0.16
7.73f4±0.05
3.03f5±0.08
1.58c6±0.02
St2
1.19a1’±0.00
2.04f2a2’±0.00
32.56f3’±0.09
6.61f4c4’±0.05
1.68f5’±0.01
1.07f6’±0.02
Year 3
PMo
St1
0.95g1±0.02
1.64g2±0.04
45.77a3±7.78
7.10g4±0.20
1.78g5±0.08
1.10g6±0.03
St2
0.76g1’±0.01
1.32g2’±0.02
48.44g3’±0.64
7.77g4’±0.05
2.37a5’±0.08
1.24c6a6’±0.00
Mo
St1
1.24h1±0.01
2.14h2±0.01
57.12h3b3±1.05
7.67h4b4±0.01
2.38b5±0.02
1.05h6b6±0.00
St2
0.98a1’±0.00
1.69a2’±0.01
36.15h3b3±0.41
6.77h4’b4±0.07
1.54h5’±0.01
1.02h6’c6±0.00
PoMo
St1
1.33i1±0.07
2.29i2±0.12
43.19i3±0.20
7.90i4±0.03
2.65i5±0.12
1.60i6c6±0.03
St2
1.27a1’±0.01
2.19a2’±0.02
38.39i3’±0.08
7.22i4c4’±0.09
1.97i5’±0.01
1.34i6’c6±0.01
A total of six species of crabs were recorded from the two study sites. These included three ocypodid species [Uca rosea (Tweedie, 1937), Uca triangularis (A. Milne-Edwards, 1873) and Uca dussumieri (H. Milne-Edwards, 1852)] and one species from each sesarmid [Episesarma mederi (A. Milne-Edwards, 1854)], portunid [Scylla serrata (Forskål, 1775)], menippid [Myomenippe hardwickii (Gray, 1831)]. Out of the three ocypodid species, U. triangularis was present only at site 1, whereas U. rosea was found to be the most dominant species at both the sites. When an independent sample t-test was performed between the two sites, they differed significantly in terms of species evenness (J) and dominance (D) indices accordingly (p<0.05).
Species density was found to be much lower at site 2 (27-32 ind. m–2) than site 1 (58-68 ind. m–2). Site 1 showed a seasonal pattern in diversity indices. The pre-monsoon season represented significantly higher species evenness (J) (p=0.000) than the post-monsoon and monsoon seasons (Fig. 2), whereas at site 2, J and D were found to be significantly different (p=0.000) among three respective seasons. Throughout the study period post-monsoon was found to be the most favourable season for abundance and even distribution of crabs.
Seasonal and inter-annual variation of Pielou’s evenness (J) and species dominance (D) indices for two study sites.
One-way ANOSIM represented significant differences in brachyuran crab community structure between the two sites (R=0.05, p<0.001). Similarity percentage analyses (SIMPER) of square-root-transformed brachyuran crab abundance data revealed that site 2 had the highest similarity (62.6%) constituted mainly due to the contribution of only one species: U. rosea (cumulative contribution to abundance 62% compared with 56.6% at site 1). The discriminating species that contributed most to the dissimilarity between the two were U. rosea (40.9%) and U. triangularis (21.4%). The nMDS plot (Fig. 3) based on the species abundance in three seasons in three years at both study sites represented a clear division of the crab community into five different groups (derived from 60% of Bray-Curtis similarities; stress value, 0.02). The first group was composed of the abundance of U. rosea at both sites; the second was formed by the abundance of Myomenippe hardwicki (Gray, 1831) at site 1; the third group was formed by the absence of U. triangularis at site 2; and the fourth group was mainly formed by the cumulative abundance of Uca dussumieri (H. Milne-Edwards, 1852) and M. hardwicki at site 2. The fifth group was formed by all the other remaining species at both sites.
Two-dimensional nMDS plot of the crab assemblages (based on crab abundances) according to Bray-Curtis similarity. Stress value (2D): 0.02.
With regard to the percentage abundance of individuals, one-way ANOVA showed significant differences (p<0.05) between seasons in each year and between years for each season at two study sites, as shown in Figure 4. When seasons were compared, U. rosea was found to be the highest in percentage of individuals representing maximum abundance during all seasons at both sites. A significantly higher percentage abundance of U. triangularis was recorded in the pre-monsoon season than in the other two seasons in every year at site 1; E. mederi was also comparatively higher in percentage during the monsoon than in the other two seasons. At site 2, U. rosea was the most dominant species throughout the year, followed by S. serrata and E. mederi.
Seasonal and interannual species composition at site 1 (A) and Site 2 (B).
Figure 5 shows the CCA analysis for the brachyuran crab assemblage based on 36 samples at each of the stations, including all the seasons in the study period. It indicates that the environmental variables explain significantly about 96% and 97.5% of the crab assemblages at site 1 and site 2, respectively, considering a cumulative percentage variance of four axes. The brachyuran crab species at each site are plotted as the points in relation to the environmental parameters plotted as arrows. Although a number of axes were determined within the study, only axes 1 and 2 were taken into consideration as they accounted for 77.5% and 75.2% of the total variability of the four axes.
The CCA Graph showing species vs. environmental correlation at site 1 (A) and site 2 (B). Abbreviations: U. rosea, Uca rosea; U. triang, Uca triangularis; U. dussum, Uca dussumieri; E. meder, Episesarma mederi; M. hardwi, Myomenippi hardwicki; S. serrat, Scylla serrata; C.Cvr, canopy cover; C. Ht, canopy height; S. Den, shrub density; S. Ht, shrub height; Pneu, pneumatophore density; H. Den, herb density; FL, fallen leaves; Soil OC, soil organic carbon; Soil OM, soil organic matter; Soil M, soil moisture; Soil den, soil density; Soil S G, soil specific gravity; WT, water temperature; AT, air temperature; ALK, alkalinity; ACI, Acidity; TH, total hardness; SAL, Salinity; Ip, inorganic phosphate; NO3-N, nitrate-nitrogen; NO2-N, nitrite-nitrogen.
Considering their vector length, acidity (r=0.58), pH of water (r=0.64), alkalinity (r=0.64), TDS (r=0.55) and soil specific gravity (r=–0.59) might be the most important environmental parameters influencing the crab assemblage in site 1(Fig. 5A). On the other hand, inorganic phosphate content of water (r=0.51), soil moisture (r=0.42), soil density (r=–0.46), soil specific gravity (r=0.42) and canopy cover (r=–0.42) were the important environmental variables that influenced the crab assemblage at site 2 (Fig. 5B). At site 1, axis 1 separated the pH-sensitive, acidity-tolerant species (on the right) from the species preferring a nutrient-rich environment. Axis 1 also explained S. serrata to be sensitive to changes in the moisture content of soil. Again, both axes separated the species with highest acidity tolerance (the upper half) and lower acidity tolerance (the lower half). At site 2, axis 1 separated the species which preferred heavier soil (contributing to the higher values of soil density and soil specific gravity) than the others. Axis 2 separated mostly the species having affinity towards inorganic phosphate (the upper half) from the species sensitive to it (the lower half).
DISCUSSION
Among all the benthic macro fauna which dwells in the intertidal zone, brachyurans are the most prominent (Jones 1984, Macintosh 1988) as they play a fundamental role in the structuring and functioning of the mangrove ecosystem (Lee 1999). Although extensive literature exists on crab diversity in southeast Asia and even in the coastal region of India (Tan and Ng 1994, Ashton et al. 2003b, Sakthivel and Fernando 2012), very few (Chakraborty and Choudhury 1992a,b) have dealt recently with crab diversity in Sundarban.
The number of crab fauna recorded in different climatic zones of the world varies considerably. Ashton et al. (2003a) have recorded 31 species of crab from Sarawak mangrove in Indonesia. The number of species found during the present study is low compared with other previous studies from the country (Table 3). However, the crab diversity in Sundarban has not been studied extensively. Moreover, the effect of habitat on the crab diversity has not yet been explored in this mangrove. The present study establishes that the brachyuran crab fauna reflects habitat status when natural mangrove forests and replanted forests are compared. Within the limitations of the sampling, two interesting facts appear to affect the brachyuran crab community structure most strongly: the age of the forest stand and the level of anthropogenic turbulence.
Previous reported works on mangrove crabs in India.
Geographical Location
Authors
No. of brachyuran species reported
Sagar island (21°31’N to 21°53’N and 88°02’E to 88°15’E)
Chakraborty et al. 1992
14
Prentice Island (21°40’N to 21°48’N and 88°16’E to 88°21’E)
Chakraborty and Chaudhury 1992a
18
Sagar island (21°31’N to 21°53’N and 88°02’E to 88°15’E)
Chakraborty and Chaudhury 1992b
4
Pichavaram Mangrove
Ravichandran et al. 2000
23
Pichavaram Mangrove (11°27’N and 79°47’E)
Vellar Mangrove (11°29’N to 79°46’E)
Ajmal Khan et al. 2005
38
8
Pichavaram Mangrove (11°27’N and 79°46’E)
Ravichandran et al. 2007
46
Karawar mangrove, west coast of India (74°18’N and 14°88’E)
Bandekar et al. 2011
13
Gulf of Kutch
Trivedi et al. 2012
19
Mudasal Odai (11°29’N and 79°47’E) and Nagapattim (10°46’N and79°59’E)
Sakthivel and Fernando 2012
38
Six species of brachyuran crabs were recorded from two different study sites. Among the two, all the species were abundant at site 1 in all seasons. The higher values of D at site 2 (Fig. 2) can be directly attributed to the maximum predominance of a particular species, U. rosea. It seems from the present study that most of the crab species preferred forested mudflats over replanted site. The assemblage groups revealed from the nMDS analysis (Fig. 3) showed clear separation according to species abundance. The high dominance of U. rosea, widely distributed at the replanted site and in the forested mangrove, contributed the overall lower diversity values. Absence of U. triangularis at the replanted site and low abundance of M. hardwicki placed them in the separate plots from the others. At site 2, microhabitat differentiation was not prominent due to low canopy cover (up to 20%) and root ramification. Hence, absence of U. triangularis, a spatially inferior competitor, is a foregone conclusion due to presence of all other competitively superior species.
Tan and Ng (1994) suggested that the maintenance of high crab species diversity is integral to the health of the mangroves. Natural and human-induced disturbances pose serious threats to the functioning of mangrove ecosystems (Osborn and Polsenberg 1996). The reaction of natural ecosystem to disturbances is often discontinuous, determined inadvertently and hence unpredictable (Holling et al. 1995). The present study shows that brachyuran crab community structure changes with the age of the stand. High ocypodid crab abundance dominates disturbed and young plantations (site 2), whereas a large abundance of grapsoid crabs, exclusively sesarmids, dominates mature forests where the canopy closes. Similar observations have been established elsewhere in different mangroves of the world (Macintosh 1984, Aksornkoae 1993, Sasekumar and Chong 1998, Ashton et al. 2003b). The result of independent sample t-test of the two study sites varied significantly (p<0.05) in terms of pneumatophore density, herb density, canopy cover and number of decaying leaves. Site 2 has a lesser number of Avicennia sp. facing the river, resulting in a lesser number of breathing roots. Reduction in the pneumatophore density and the attached algal mass would lead to higher near-bottom flow rates, with a concomitant reduction in deposition of nutrients (Abelson and Denny 1997). The present study establishes U. rosea to be an opportunistic species which has a negative correlation with the pneumatophore density and therefore a lesser need for deposited nutrients, a finding which is corroborated by Macintosh (1984). Open mudflats provide a unique habitat for many species of brachyuran crabs; 10 species of brachyuran crabs have been recorded from the open mudflat of Mahi river estuary (Pandya and Vachhrajani 2013). However, in the current study, out of the six species of crab recorded, U. rosea dominates the replanted open mudflat habitat, whereas U. triangularis mostly prefers forested mangrove.
Elucidating the relationship between species occurrence and environmental variables allows us to identify the ecological processes that regulate different crab populations. Regardless of the changes in environmental conditions, many species change their habitats according to their ontogeny or seasonal rhythms, which means that relations between different species and their environment or habitat are spatially and seasonally dynamic (Morrison et al. 2002, Kanou et al. 2005). Mangrove crabs are generally affected by the abiotic and biotic factors of the habitat (Macintosh 1988, Lee 1999, Ashton et al. 2003b). The CCA in the present investigation recognized that hydrological and edaphic factors and the physical architecture of the habitat have a great impact on the brachyuran assemblage. The findings were similar to earlier observations (Soundarapandian et al. 2008, Chakraborty and Choudhury 1992b). The crab assemblage of the study area of Sundarban mangrove is influenced seasonally by a variety of factors, viz. acidity, pH content of water, alkalinity, soil moisture content, soil pH, nutrient content of water and habitat attributes such as canopy cover and canopy height.
High dominance of a single species might indicate a stressful environment (Macintosh et al. 2002 and Ashton et al. 2003b). As the ecosystem has not been restored to the normal mature condition by replantation, only one species (U. rosea) dominates the crab fauna at site 2. Further time may change the crab community structure at that location or the anthropogenic effects may irreversibly prevent the ecosystem from going back to its original state (Ashton et al. 2003b). The crab community composition may give an indication of the habitat health. Ghost crabs have already been used as a tool for rapid assessment of human impacts on the exposed sandy beaches (Barros 2001). In mangrove, given the environmental complexity, it is uncertain whether one or two species could be sufficiently sensitive to all major factors affecting a complex ecosystem to serve as indicator species. Rather, monitoring the changes in community structure may provide a more prominent and effective tool for evaluating ecosystem health (Cairns et al. 1993 and Ashton et al. 2003b). As the brachyuran crab communities are omnipresent throughout the Sundarban mangrove and are an important part of the mangrove food web, they can be used as a good indicator of the local condition. The fact that they are primarily sedentary and have a reasonably long life cycle (Ashton et al. 2003b) also supports their candidature. Changes in the crab community thus reflect the mangrove habitat status, which can be used as a tool for management purposes.
In conclusion, the determinant role of the habitat characteristics in moderating crab abundance and diversity is indicated. More meticulous approaches to the assessment of past and current restoration efforts must be developed (Stevenson et al. 1999) in accordance with the study of the health of the mangrove at undisturbed natural sites. An understanding of the ecological effects of disturbance to mangrove forests and the changes in the structure and dynamics of the crab meta-community will be of great importance for conservation of this unique ecosystem.
ACKNOWLEDGEMENTS
The authors are thankful to the local villagers of Sundarbans, and to the Head of the Department of Zoology, University of Calcutta, India for the facilities provided. Financial support from the University Grant Commission-National Eligibility Test project is gratefully acknowledged.
REFERENCES Abelson A. Denny M. Settlement of marine organisms in flow 199728317339http://dx.doi.org/10.1146/annurev.ecolsys.28.1.317 Acosta C.A. Butler M.J. Role of mangrove habitat as a nursery for juvenile spiny lobster, Panulirus argus, in Belize 199748721728http://dx.doi.org/10.1071/MF96105 Ajmal Khan S. Rafii S. Lyla P.S. Brachyuran crab diversity in natural (Pitchavaram) and artificially developed mangroves (Vellar estuary) 20058813161324 Aksornkoae S.1993 Bangkok, Thailand. IUCN, AIT APHA 1998 Washington, DC Ashton E.C. Macintosh D.J. Hogarth P.J. A baseline study of the diversity and community ecology of crab and molluscan macrofauna in the Sematan mangrove forest, Sarawak, Malaysia 200319127142http://dx.doi.org/10.1017/S0266467403003158 Ashton E.C. Hogarth P.J. Macintosh D.J. A Comparison of Brachyuran Crab Community Structure at Four Mangrove Locations under Different Management Systems along the Melaka Straits-Andaman Sea Coast of Malaysia and Thailand 20032614611471http://dx.doi.org/10.1007/BF02803654 Bandekar P.D. Neelkantan K. Kakati V.S. Biodiversity of crabs in Karawar mangrove environment West Coast of India. 2011315 Barros F. Ghost crabs as a tool for rapid assessment of human impacts on exposed sandy beaches 200197399404http://dx.doi.org/10.1016/S0006-3207(00)00116-6 Berger W.H. Parker F.L. Diversity of planktonic foraminifera in deep sea sediments 197016813451347http://dx.doi.org/10.1126/science.168.3937.1345 Brower J.E. Zar J.H. von Ende C.N.1998 New York McGraw-Hill, Cairns J. Mccormick P.V. Niederlehner B.R. A proposed framework for developing indicators of ecosystem health 1993263144http://dx.doi.org/10.1007/BF00006084 Chakraborty S.K. Choudhury A. Ecological studies on the zonation of brachyuran crabs in a virgin mangrove island of Sundarbans, India 199234189194 Chakraborty S.K. Choudhury A. Population Ecology of Fiddler Crabs (Uca spp.) of the mangrove estuarine complex of Sundarbans 1992337888 Chakraborty S.K. Poddar T.K. Choudhury A. Species diversity of macrobenthos of Sagar Island, Sundarbans, India 199245Suppl. A435444 Clarke K.R. Non-parametric multivariate analyses of changes in community structure. 199318117143 http://dx.doi.org/10.1111/j.1442-9993.1993.tb00438.x Clarke K.R. Warwick R.M.1994 Plymouth, U.K Plymouth Marine Laboratory Ellner S.P. Pair approximation for lattice models with multiple interaction scales 200121043547http://dx.doi.org/10.1006/jtbi.2001.2322 Golley F. Odum H.T. Wilson R.L. The structure and metabolism of Puertorican and mangrove forest in May196243918http://dx.doi.org/10.2307/1932034 Gopal B. Chauhan M. Biodiversity and its conservation in the Sundarban Mangrove Ecosystem. 200668338354http://dx.doi.org/10.1007/s00027-006-0868-8 Grasshoff K. Determination of nitrite, nitrate, oxygen, thiosulphate. Grasshoff K. Ehrhardt M. Kremling K.1983 New York Verlag Chemie Weinheim 61150 Grasshoff K. Ehrhardt M. Kremling K.1983 New York Verlag Chemie Weinheim Hatcher B.G. Johannes R.E. Robertson A.I. Review of research relevant to the conservation of shallow tropical marine ecosystems 198927337414 Holling C.S. Schindler D. Walker W. Biodiversity in the functioning of ecosystems: An ecological synthesis. Perrings C. Maler K.G. Folke C.1995 Cambridge, U.K., Cambridge University Press 4483http://dx.doi.org/10.1017/CBO9781139174329.005 Jones D.A. Crabs of the mangal ecosystem Por F.D. Dor I.1984 The Hague, Netherlands W. Junk Publishers 89109 Jongman R.H.G. Ter Braak C.J.F. Van Tongeren O.F.M.1995 Cambridge Cambridge University Press http://dx.doi.org/10.1017/CBO978051152557 Kanou K. Sano M. Kohno H. Larval and juvenile fishes occurring with flood tides on an intertidal mudflat in the Tama River estuary, central Japan 200552158164http://dx.doi.org/10.1007/s10228-005-0267-5 Kathiresan K. Bingham B.I. Biology of mangroves and mangrove ecosystems 20014081251http://dx.doi.org/10.1016/S0065-2881(01)40003-4 Kon K. Kurokura H. Tongnunui P. Do mangrove root structures function to shelter benthic macrofauna from predators? 200937018http://dx.doi.org/10.1016/j.jembe.2008.11.001 Laegdsgaard P. Johnson C. Why do juvenile fish utilise mangrove habitats?2001257229253http://dx.doi.org/10.1016/S0022-0981(00)00331-2 Lawton J.H. Plant architecture and the diversity of phytophagous insects 1983282339http://dx.doi.org/10.1146/annurev.en.28.010183.000323 Lee S.Y. Ecological role of grapsid crabs in mangrove ecosystems: A review. 199849335343http://dx.doi.org/10.1071/MF97179 Lee S.Y. Tropical mangrove ecology: Physical and biotic factors influencing ecosystem structure and function 199924355366http://dx.doi.org/10.1046/j.1442-9993.1999.00984.x Levings S.C. Garrity S.D. Burns K.A. The Galeta oil spill. III. Chronic reoiling, long-term toxicity of hydrocarbon residues and effects on epibiota in the mangrove fringe 199438365395http://dx.doi.org/10.1006/ecss.1994.1026 Macarthur R.H. Macarthur J.W. On bird species diversity 196142594598http://dx.doi.org/10.2307/1932254 Macia A. Quincardete I. Paula J. A comparison of alternative methods for estimating population density of the fiddler crab Uca annulipes at Saco Mangrove, Inhaca Island (Mozambique)2001449213219http://dx.doi.org/10.1023/A:1017502817087 Macintosh D.J. Ecology and productivity of Malaysian mangrove crab populations (Decapoda: Brachyura). Soepadmo E. Rao A.N. Macintosh D.J.1984 KualaLumpur, Malaysia, and UNESCO, Paris University of Malaysia 354377 Macintosh D.J. The ecology and physiology of decapods of mangrove swamps 198859315341 Macintosh D.J. Ashton E.C. Havanon S. Mangrove rehabilitation and intertidal biodiversity: A study in the Ranong mangrove ecosystem, Thailand 200255331345http://dx.doi.org/10.1006/ecss.2001.0896 Manna S. Chaudhuri K. Bhattacharyya S. Dynamics of Sundarban estuarine ecosystem: eutrophication induced threat to mangroves 201068 http://dx.doi.org/10.1186/1746–1448-6–8 Morrison M.A. Francis M.P. Hartill B.W. Diurnal and tidal variation in the abundance of the fish fauna of temperate tidal mudflat 200254793807http://dx.doi.org/10.1006/ecss.2001.0857 Nelson D.W. SommersL.E. Page A.L.Total carbon, organic carbon, and organic matter19969611010Madison, WISoil Science Society of America, American Society of Agronomy Ng P.K.L. Guinot D. Davie P.J.F. Systema Brachyurorum: Part I. An Annotated Checklist of Extant Brachyuran Crabs of the World 2008171286 Norušis M.2000 Upper Saddle River Prentice-Hall Osborn J.G. Polsenberg J.F. Meeting of the mangrovellers: The interface of biodiversity and ecosystem function 199611354356http://dx.doi.org/10.1016/0169-5347(96)30035-9 Pandya P.J. Vachhrajani K.D. Brachyuran crab diversity of lower estuarine mud flats of Mahi River with new record of two species from Gujarat, India. 20132242250 Pielou E.C 1969 New York. Wiley-Interscience Primavera J.H. Fish predation on mangrove-associated penaeids: The role of structures and substrate 1997215205216http://dx.doi.org/10.1016/S0022-0981(97)00046-4 Ravichandran S. Soundarapandian P. Kannupandi T. Zonation and distribution of crabs in Pichavaram mangrove swamp, southeast Coast of India 2000477378 Ravichandran S. Anthonisamy S. Kannupandi T. Habitat Preference of crabs in Pichavaram Mangrove Environment, Southeast Coast of India 200724755http://dx.doi.org/10.3923/jfas.2007.47.55 Robertson A.I. Leaf-burying crabs: Their influence on energy flow and export from mixed mangrove forests (Rhizophora spp.) in northeastern Australia 1986102237248http://dx.doi.org/10.1016/0022-0981(86)90179-6 Robertson A.I. Blaber S.J.M. Plankton epibenthos and fish communities Robertson A.I. Alongi D.1992 New York. Springer-Verlag173224http://www.agu.org/books/ce/v041/CE041p0173/CE041p0173.pdf Sakthivel K. Fernando A. Brachyuran crabs diversity in Mudasal Odai and Nagapattinam coast of south east India. 20121136143 Sasekumar A. Chong V.C. Faunal diversity in Malaysian mangroves 199875760http://dx.doi.org/10.2307/2997697 Smith III T.J. Boto K.G. Frusher S.D. Keystone species and mangrove forest dynamics: The influence of burrowing crabs on soil nutrient status and forest productivity. 199133419432http://dx.doi.org/10.1016/0272-7714(91)90081-L Soundarapandian P. John Samuel N. Ravichandran S. Biodiversity of Crabs in Pichavaram Mangrove Environment, South East Coast of India 20084113118http://dx.doi.org/10.3923/ijzr.2008.113.118 Stevenson N.J. Lewis R.R. Burbridge P.R. Disused shrimp ponds and mangrove rehabilitation Streever W.1999 Dordrecht, Netherlands Kluwer Academic Publishers 277297http://dx.doi.org/10.1007/978-94-011-4683-8_28 Strickland J.D.H. Parsons T.R. A practical handbook of sea-water analysis 19721671311 Tan C.G.S. Ng P.K.L. An annotated checklist of mangrove brachyuran crabs from Malaysia and Singapore 19942857584http://dx.doi.org/10.1007/BF00005655Ter BraakC.J.F.Canonical correspondence analysis: A new eigenvector technique for multivariate direct gradient analysis19866711671179 Ter Braak C.J.F. CANOCO – A Fortran program for canonical community ordination by (partial) (detrended) (canonical) correspondence analysis, principal component analysis and redundancy analysis (version 2.1)6700 AC Wageningen, The Netherlands Institute of Applied Computer Science, Statistical Department Wageningen 1988195 Trivedi J.N. Gadhavi M.B. Vachhrajani K.D. Diversity and habitat preference of brachyuran crabs in Gulf of Kutch, Gujarat, India 201211323