The distribution and diversity of benthic macroinvertebrate fauna in Pondicherry mangroves, India

Background Species distribution, abundance and diversity of mangrove benthic macroinvertebrate fauna and the relationships to environmental conditions are important parts of understanding the structure and function of mangrove ecosystems. In this study seasonal variation in the distribution of macrobenthos and related environmental parameters were explored at four mangrove stations along the Pondicherry coast of India, from September 2008 to July 2010. Multivariate statistical analyses, including cluster analysis, principal component analysis and non-multidimensional scales plot were employed to help define trophic status, water quality and benthic characteristic at the four monitoring stations. Results Among the 528 samples collected over 168 ha of mangrove forest 76 species of benthic macroinvertebrate fauna were identified. Macrofauna were mainly composed of deposit feeders, dominated numerically by molluscs and crustaceans. Statistical analyses yielded the following descriptors of benthic macroinvertebrate fauna species distribution: densities between 140–1113 ind. m-2, dominance 0.17-0.50, diversity 1.80-2.83 bits ind-1, richness 0.47-0.74 and evenness 0.45-0.72, equitability 0.38-0.77, berger parker 0.31-0.77 and fisher alpha 2.46-5.70. Increases of species diversity and abundance were recorded during the post monsoon season at station 1 and the lowest diversity was recorded at station 2 during the monsoon season. The pollution indicator organisms Cassidula nucleus, Melampus ceylonicus, Sphaerassiminea minuta were found only at the two most polluted regions, i.e. stations 3 and 4. Benthic macroinvertebrate fauna abundances were inversely related to salinity at the four stations, Based on Bray-Curtis similarity through hierarchical clustering implemented in PAST, it was possible to define three distinct benthic assemblages at the stations. Conclusions From a different multivariate statistical analysis of the different environmental parameters regarding species diversity and abundance of benthic macroinvertebrate fauna, it was found that benthic communities are highly affected by all the environmental parameters governing the distribution and diversity variation of the macrofaunal community in Pondicherry mangroves. Salinity, dissolved oxygen levels, organic matter content, sulphide concentration were the most significant parameters.


Background
India is a large coastal nation located along the Indian Ocean with 7,517 km of coastline, along which there are many biotopes such as estuaries, lagoons, backwaters, mangroves, salt marshes, coral reefs and creeks. Mangroves in tropical and subtropical intertidal regions of the world support rich faunal resources and play an important role in estuarine and coastal food webs [1]. India has approximately 2.7% of the world's mangroves, covering an estimated area of 4,827 sq km. Almost 80% of the mangrove forests are located along the east coast and the remaining 20% are located on the west coast [2]. The three main benthic faunal components (i.e. microfauna, meiofauna and macrofauna) represent important ecological indicators. Studies on benthic diversity, population dynamics and changes caused by natural or anthropogenic processes are essential for resource management [3]. Understanding the structure of the benthic faunal communities in relation to the impacts of pollution is an important part of monitoring changes in mangrove ecosystems in India [4][5][6][7][8][9][10]. Succession in macrobenthic communities, in relation to organic enrichment and pollution in the marine environment has been reported by Samidurai et al. [11]. "According to [12,13], the distribution of macrobenthic communities is highly correlated with sediment characteristics, which is related to a wider set of environmental condition". Water quality and benthos characteristics have been investigated in coastal ecosystems around the world [10,[14][15][16] and indicate that the health of benthic communities is related to water quality conditions in fringing communities, such as mangroves. Environmental conditions like salinity, oxygen, temperature and nutrients influence the composition, distribution and growth of biota [17]. Total organic carbon of the sediment influences fertility of the soil, thereby enhancing biological activity [18].
Diversity and density of the macrobenthos is dependent on chance settlement of pelagic larval forms of different species, affinity to suitable substratum and also the degree of stress effect caused by strong waves and tide currents. The aim of the present study was to identify the relationships between water quality and benthic macroinvertebrate fauna characteristics in Pondicherry mangroves.

Study area
The present study area was within 11°46′03″ to 11°53′40″ North latitude and 79°49′45″ to 79°48′00″ East longitude ( Figure 1). In the study region, mangroves are fringing vegetation covering 168 ha distributed along the Ariankuppam estuary, which is seasonally bar-built and has semi diurnal tides that flow eastwards and empty into the Bay of Bengal at Veerampatinam on south east coast of India, carrying the waste from the adjacent agriculture lands and industries, in addition to domestic municipal and distillery effluents. The details on GPS coordinates, zone, and soil substratum are presented in Table 1. The present investigation was carried out in four Stations: 1 Veerampattinam; 2 Thengaithittu; 3   [20]. The benthos in the sediment samples recovered after sieving through 0.5 mm mesh sieve was brought to the laboratory in polythene bags, transferred to a large, white-bottomed tray, and the animals were hand sorted. After this preliminary examination, the whole sample was treated with 5% buffered formalin and kept for further analysis. Annual rainfall, temperature and relative humidity data was obtained from meteorological department at Chennai. Dissolved oxygen was estimated by Winkler's methods and sulphide by [21] and salinity using a refractometer. pH and temperature were measured using a pH meter, electrical conductivity (EC) was determined by using an EC (Elico) meter. Sediment texture was determined by a pipette analysis method [22]. The organic matter of the sediment was analyzed by a wet oxidation method [23]. Fauna were identified to the lowest practical taxonomic level using standard references; Polychaeta [24]; Crabs [25][26][27], Amphipods [28]; Mollusc [29]. For the sake of interpreting the data, a calendar year was divided into four main seasons, pre monsoon (July-September), monsoon (October-December), post monsoon (January-March), and summer (April-June). The identified samples were expressed as No/m 2 . Biodiversity indices such as species diversity, richness and   evenness were calculated following standard formulae [30][31][32][33].

Statistical analysis
All the statistical analysis methods used were according to Johnson and Wichern [34]. A variety of diversity indices have been used in benthic ecology to assess the environmental quality and the effect of disturbances on benthic communities. In the present study, calculation of two diversity indices were carried out using Shannon Wiener diversity (H'), Margalef diversity and Pielou Evenness (J') reflects the even occurrence of species within a community. Simpson Index species richness (D) was used for univariate measures to assess community structure. Differences in univariate measures between sites were tested using correlation coefficient estimation. Correlations between the benthic macroinvertebrate fauna assemblage and physical-chemical characteristics of water and sediment were evaluated using Pearson correlation analysis. Mean and standard deviations were calculated for each parameter. These statistical analysis programs are part of the SPSS statistical program (Version 13.0 for Windows XP, SPSS, and Chicago, IL, USA). In recent years, multivariate statistical techniques are the preferred tool for a meaningful data reduction and interpretation. Multivariate statistical techniques such as cluster analysis (CA), Non-Multidimensional Scale plot (MDS) and Principle Component analysis (PCA) have widely been used as unbiased methods in analysis between water quality and marine organisms; water quality [35][36][37], phytoplankton characteristics [38,39], benthos characteristics [40][41][42]. Multivariate analysis such as, CA, MDS, and PCA were constructed based on macro faunal abundance, diversity indices and with environmental parameters. Based on the groups obtained from a cluster analysis, species that contributed the most to this distribution were determined using similarity percentage program PAST (Statistical Version 1.93 for Windows XP).

Data treatment
Nearly all the multivariate statistical methods need variables to confirm the normal distribution, thus, the normality of the distribution of each variable was checked by analyzing kurtosis and skewness statistical tests before multivariate statistical analysis was conducted [43]. The original data demonstrated values of kurtosis ranging from −1.82 to 10.56 and skewness values ranging −1.91 to 3.61 indicating that the data was not normally distributed.
Since most of values of kurtosis and skewness were >0, the raw data of all variables were transformed in the form x′=log 10(x). After transformation, the kurtosis and skewness values ranged from −1.081 to 3.01 and 1.01 to 1.66,  respectively, indicating that all the data were normally distributed or close to being normally distributed. In the case of CA, PCA, and MDS, all log-transformed variables were also z-scale standardized to minimize the effects of difference units and variance of variables and to render the data dimensionless [44].

Mangrove sediment characteristics
Mangrove sediment substratum was mainly composed of sand with an admixture of silt and clay. The sand  fraction ranged between (39.54-87.31%), followed by silt (9.63-32.37%), clay (3.06-31.20%) and organic matter (0.94 -4.64%) ( Figure 8). Seasonally, Station 1 recorded higher fractions of sand in the summer, higher silt content during the post-monsoon and pre-monsoon period, and higher clay during the post monsoon season at Station 4. Changes in sediment composition were mainly due to transport of sediments by tides and currents.

Dominant taxa
Species dominance in the present study varied from 0.174 to 0.508, the minimum value of dominance was recorded in the pre-monsoon season and maximum value during the monsoon season in 2009 (Table 4). Dominance values showed significant positive correlation between the berger parker (r = 0.843; p < 0.01), and fisher alpha (r = 0.650; p < 0.01). However, there was a negative correlation between evenness (r = − 0.526; p < 0.01), richness (r = −0.586; p < 0.01) and equitability (r = −0.764; p < 0.01 Seven species of amphipods were recorded for the first time from Pondicherry mangroves. The Gammaroidea species Eriopisella sp and Eriopsia chilkensis were collected at stations 1 and 2. The Melitidae species M. dentada and Ampithoidae species Cymadusa. pathyi were observed at stations 3 and 4. The Aoridae species Grandidierella bonnieroides, G. pathyi and Isala Montagui were observed at station 4. A k-dominance plot curve was drawn based on high-and low flow macroinvertebrate community data. In the present investigation the data collected during various seasons and from all four stations was fed into to the dominance plot ( Figure 10). The highest dominance was seen at Station 1, where the macroinvertebrate assemblage had the highest diversity. The curve for Stations 3 and 4 showed the minimum diversity. The highest diversity was recorded in post monsoon season and lowest diversity in the monsoon season. The dominance curve did not show an 'S' shape due to the presence of opportunistic species C. nucleus, M. ceylonicaus, and S. minuta.

Seasonal variation of species diversity
Shannon diversity (H') varied between the stations ranging from 1.80-2.83 and Margalef diversity ranged from 1.67 to 2.70. Diversity was high in the post monsoon season at    There was a negative correlation between Fisher alpha and equitability (r = − 0.643; p < 0.01).

Relationship between benthic fauna and environmental factors
Correlation analysis showed a relationship between benthic diversity indices and abiotic variables ( Table 2)

Multivariate statistical analysis
Cluster analysis revealed three distinct benthic macroinvertebrate fauna groupings, which appeared to reflect differences in sediment/habitat types within Pondicherry mangroves ( Figure 14). Cluster 1 consisted of Stations 1 and 2, outside of the monsoon season, with high diversity, richness and abundance of organisms, particularly C. cingulata; T. telescopium, Anadara rhombea and M. meritrix. Cluster 1 was also characterized by high proportion of coarse sediments and high DO levels. From the resulting dendrogram, it was possible to grade the results according to stations and seasons. They consisted Cluster −1 (PM1, PM2, PM11, PM22, PR1, PR2, SU1-SU3, SU11-SU33), Cluster 2 consisted of all stations during the monsoon season, which is a period characterized by high rainfall and river flow. Cluster 3 consisted of Stations 3 and 4 outside of the monsoon season and was characterized by low species abundance and diversity, along with high organic matter (OM) and sulphide content. Cluster 3 also contained pollution In the PCA analyses, variables associated with principle components 1 and 2 accounted for 70.3% variability among the samples (Figure 16). Stations 1 and 2 grouped on the left side of the plot (correlating with high abundances of macro benthic organisms, high salinity and DO, and low sulphide and OM concentrations). Stations 3 and 4 grouped together on the plot (correlating with low abundances of macro benthos organisms, low DO and salinity, and high sulphide and OM levels). In the multivariate analysis, monsoonal samples were ordinated separately from all other samples.

Water quality parameters
Variations observed in water quality parameters can be attributed to regional patterns in climate and biological activity. Generally, surface water temperature is influenced by the intensity of solar radiation, evaporation, freshwater  influx and cooling and mixing up with ebb and flow from adjoining neritic waters [45]. In the present study, summer peaks and monsoonal troughs in air and water temperatures were similar to those reported in the east and west coast of India [45]. Salinity acts as a limiting factor in the distribution of living organisms, and its variation caused by dilution and evaporation influences faunal distribution in the intertidal zone [46]. In the present study salinity at all the four stations was high in summer and low in the monsoon season indicating that variation in salinity at study sites was affected by freshwater runoff entering the creek systems, as observed in the Gulf of Kachchh [47] and Uppanar estuary [48]. Hydrogen-ion concentration varied little among the four stations and pH was alkaline throughout the study period. Higher pH observed in summer season could be attributed to the removal of CO 2 by the photosynthetic organisms and the lower pH during monsoon season could be due to the dilution of saline water by freshwater inflow [47]. The relatively low DO values observed in the summer are attributable to the prevalence of high saline neritic waters in the mangrove channels, as well as fluctuations in temperature and salinity, which affect dissolution of oxygen [49]. The low DO levels recorded at Stations 3 and 4 may be attributable to the inflow of pollutants and high oxygen demand from elevated levels of hydrogen sulphide. Higher sulphide levels at Station 3 and 4 are in part related l to industries discharging effluents near these two stations. A peak value of H 2 S (0.92 μg g-1) in soil has been reported at Kayamkulam estuary [50]. The presence of high sulphide content in the sediments of polluted regions in this study is likely one of the limiting factors in benthic macroinvertebrate fauna abundance and distribution. By  contrast sediments with a mixture of organic matter, sand and clay, but low sulphide, seemed to support higher abundances of benthic marco fauna. Distribution and ecology of benthic communities in relation to station and season. This study provides a baseline for the distribution, abundance and diversity of benthic macroinvertebrate fauna of Pondicherry mangroves in India. The order of importance of mollusks (gastropods and bivalves), crustaceans, amphipods and polychaetes in the present study is similar to that observed by [51]. Macro benthic faunal densities observed in this study (140-1113 ind. m-2) were higher than that reported by [52] in Zuari estuary (50-1037 ind.m-2) and [53] in Andaman seas (80-998 ind.m-2). The densities were comparable to those observed in the gulf of Arid Zone mangroves of Gulf Kachchh on the west coast of India (i.e. 424-2393 ind.m-2) [8]. The high densities recorded in the postmonsoon season in Pondicherry mangroves could be due to low temperatures and turbidity coupled with stable environmental conditions. Post-monsoon season (Nov-Feb) peaks in density have also been reported for the west coast of India [51].
In the present study, 22 species belonging to 12 genera and 5 families of brachyuran crabs were recorded in Pondicherry mangroves. The distribution of crabs showed relationships to substratum characteristics, salinity, degree of tidal inundation and wave exposure. The distributional patterns are significant from a harvest perspective. Uca annulipes, U. inversa, and U. triangularis are largely caught during the monsoon season. Scylla serrata and Thalamita crenata are mostly caught during the post-monsoon season. Portunus sanguinolentus and Calappa lophos are generally caught mouths of mangrove regions.
In Pondicherry mangroves 16 bivalves and 19 gastropods were recorded, similar to the numbers previously reported for mangroves of Sunderbans [4], Saravanakumar   Schrijvers et al. [54] reported the dominance of polychaetes over mollusks in Kenyan mangrove fauna due to silty clay substratum. Soft mangrove substrates favor tube dwellers over diggers and burrowing animals, such as bivalves. The high diversity of polychaete indicate favorable ecological conditions that exist in mangrove ecosystems [55]. Benthic macro fauna are important components of coastal food webs. Many fish species in estuarine ecosystems are strongly dependent taxa available that reside on the sediment surface, such as amphipods, mysids, and surface deposit feeders (SDF) [56]. Several studies have demonstrated the presence of amphipods in Pitchavaram mangroves of India [5,57]. In the present study seven amphipod species were observed. Due to limited dispersion capabilities and habitat specificity of amphipods, they may be of use in biogeography and environmental monitoring of mangrove ecosystems.
Limited quantitative information has been published on mangrove habitats in India, including the distribution of benthic macroinvertebrate fauna and their relationship to environmental factors [3][4][5][6][7][8][9][10]. The results of cluster analysis in this study show three major groupings, mainly segregated by region and season. Stations 1 and 2 clustered together in all but the monsoon season, largely based on the high density and diversity of organisms. The increased abundance of organisms, species richness and diversity observed in Stations 1 and 2 could be due to the presence of coarse sediments, high DO and salinity [58,59]. Stations 3 and 4 clustered together based on low abundance and diversity, perhaps due to the silty characteristics of the sediments, along with high organic matter content, sulphide concentrations and low DO [16]. From a seasonal perspective, all sites clustered together during the monsoon season, demonstrating the impact of climatic conditions.

Species diversity indices
Pearson and Rosenberg [11] proposed that diversity indices provide important insights into faunal communities at different stages in succession. Snelgrove [60] proposed that the species diversity is mainly controlled by fluctuations in the environment that leads to less diversity. Species diversity tends to be low in physically controlled ecosystems [61]. In the present study, species richness of benthic macrofauna was highest during the post-monsoon season and summer, similar to observations from Cochin backwaters [62]. Low population density recorded in the monsoon season was apparently due to the effect of heavy rainfall, as previously observed by Saravanakumar et al. [8], whoireported a 'severe decline' of macro benthos in the shallow waters during southwest monsoon attributed to lowered salinity [51].
Wilhm and Dorris [63] propose that values less than 1.0 for diversity index (H) in estuarine waters indicate heavy pollution, values between 1.0 and 3.0 indicate moderate pollution, and values exceeding 3.0 indicate non-polluted water. Diversity values in the study area ranged from 1.83-2.83. Thus, these values suggest that the mangroves examined in this study are moderately polluted and the macro benthic community is under stress due to natural and/or anthropogenic factors. Maximum diversity and richness recorded in this study during the post-monsoon season might be due to stable environmental factors, such as high DO and salinity, which play a vital role in faunal distribution. Relatively high species richness, evenness and diversity were observed at Stations 1 and 2, compared to Stations 3 and 4. The low species richness recorded in this study during monsoon might be due to the freshwater runoff containing inadequately treated sewage and low salinity, which in turn affected the distribution of benthos.

Multivariate statistics
CA, MDS and PCA have been widely used in the evaluation of spatial and temporal variations in water quality and benthic characteristics of aquatic ecosystems [35,42,64]. CA was used in this study to examine the differences between the monitoring stations during the four seasons of the year. Based on the cluster analysis, the concentration of sulphide and organic matter in cluster 3 (i.e. Stations 3 and 4) were high compared to clusters 1 (Stations 1 and 2) and 2 (all sites during the monsoon season). In defining environmental factors important in characterizing mangrove condition MDS analyses were used. Low DO, high sulphide, low salinity, high OM and high clay content appear to indicate deteriorating water quality from the standpoint of benthic macroinvertebrate fauna. Human activities have a strong influence on the aquatic environment in the southeast coast of Bay of Bengal [45]. Sources of pollution include agricultural runoff, leaching from solid waste disposal sites, domestic and industrial waste disposal. Similar approaches based on Results of multivariate analysis suggests that for mangrove benthic macroinvertebrate fauna that environmental parameters such as DO, sulphide, salinity, OM and clay have a strong influence on species composition and diversity.

Conclusions
This study provides insights into the effects of a range of environmental parameters on macro benthic communities of Pondicherry mangroves in India. Altogether 76 species of benthic macroinvertebrate fauna, belonging to five major groups, were identified at the four sampling stations. Station 1 was dominated by sandy sediment, high salinity, high DO and relatively low sulphide levels. The region displayed the high species diversity, abundance and species richness. Station 2 was characterized by sandy sediment, low organic matter and relatively high species diversity. Stations 3 and 4 had higher sulphide concentrations, silty sediments and lower DO, with relatively low species diversity. The temporal distribution of benthic macro invertebrate fauna exhibited the highest species density during post-monsoon season. The decrease of benthos during the monsoon may be attributable to low temperatures and salinities. CA, MDS and PCA analyses were useful in helping to define spatial and temporal patterns in mangrove water quality and benthic macro invertebrate fauna characteristics in the Pondicherry mangroves. Temperature, salinity, DO, sulphide, sediment composition and organic matter content all proved to be important descriptive parameters in terms of the abundance and distribution of benthic fauna.