Abstract

A quantitative and qualitative study of interstitial fauna and environmental variables was carried out on five selected sandy beaches of the west coast of India. Species of nine interstitial taxa abound the beaches. Nematodes, harpacticoid copepods, turbellarians, and polychaetes constituted the bulk of the population. The available energy in the beaches ranged from 0.2245 to 16.08 joules/mg and the grain size varied from 0.93 to 2.88. Organic matter correlated significantly with coarse sand (Pearson correlation ; ). Organic carbon, particle size, and dissolved oxygen determined the abundance and distribution of interstitial fauna as per multivariate BIOENV analysis. Shannon Wiener diversity index was maximum at Cherai (2.027) and minimum at Sakthikulangara (1.144) beach. The value of nematode/copepod ratio (/) indicated at Sakthikulangara beach validates the increased sensitivity of harpacticoids to environmental stress.

1. Introduction

The coastline of Kerala, which extents up to 590 kms is intercepted by fishing harbours, ports, tourist resorts, mangrove beds, sandy beaches, and so forth [1]. Sandy beaches harbour diverse and abundant assemblage of interstitial organisms and the different taxonomic groups have complex distribution patterns [2]. Cyclic sedimentary processes make the beaches of Kerala fragile. A beach, maintained for a particular period, may not remain there forever. The wave climate and the beach processes along the Kerala coast are found to vary considerably in both spatial and temporal frame [3]. Erosion adversely affects the density of intertidal interstitial populations, which live in the interstices of beach sediments. Biodiversity assessment is the key to understand the relationship between biodiversity and ecosystem functioning [4]. Hydrodynamics changes the physicochemical sedimentary environment and food availability which are directly related to meiofauna spatial distribution patterns [5]. Meiofauna plays a major role in pollution monitoring studies [6, 7]. The meiofaunal species are vulnerable to abiotic and hydrodynamic disturbance [8]. Though arguments were traditionally put forth against the use of meiofauna as biological indicator, underlining difficulties in identification, the high rate of sampling frequency and the microscopic size of the organisms, and new technologies and tools, such as standardized methodologies, electronic identification keys, molecular approaches, and the creation of new indices, currently allow for and promote the use of meiofauna in ecological studies [9]. Though several meiofaunal studies have been conducted on the west coast of India [1016], there is scanty information on the ecology and taxonomy of interstitial fauna from the sandy beaches of Kerala. Further, ecological studies of meiofauna are necessary to understand the trophodynamics of sandy beaches. Although a great deal of uncertainty remains in our understanding of meiofaunal biodiversity, ecology, and evolutionary biology, this area of research remains among the most challenging, the most neglected, and potentially the most enlightening frontiers of discovery in biology [17].

2. Materials and Methods

2.1. Study Area

Study Area. Five sandy beaches (Figure 1) such as Cherai (lat.10°09′N; long.76°02′E; St.1), Fort Kochi (lat.9°18′N; long.76°05′E; St.2), Arthungal (lat.9°10′N; long.76°23′E; St.3), Sakthikulangara (lat.8°45′N; long.76°38′E; St.4), and Veli (lat. 8°29′ N; long.76°59′ E; St.5) located along the coast of Kerala, India, were identified by site survey to undertake a study of the psammophilous fauna. Cherai and Arthungal beaches are relatively stable and less disturbed. Fort Kochi and Veli beaches are affected by industrialization and Sakthikulangara is characterized by natural radioactive pollution of thorium and presence of large quantities of faecal waste. Seasonal collections were made for one year (2005-06). Samples were collected at two points (100 m apart) along a transect from the subtidal area of each beach by using a galvanized iron core of 8 cm inner diameter. Three undisturbed sediment samples each were immediately transferred into bottles. 7% MgCl2 prepared in filtered sea water was added to the sediment sample to anaesthetize the animals. The samples were preserved in buffered 4% formalin prepared with filtered seawater. 0.1% Rose Bengal [18] was added to the sample for efficient extraction of the fauna. Meiofauna was separated by suspension-decantation method [19] with a few modifications [20].

2.2. Granulometry

Sediment samples of approximately 150 gms were used from each station for carbon analysis and textural analysis. Sediment was subjected to repetitive rinsing and decanting employing distilled water to remove the salinity and then dried in an oven under a temperature of 40–50°C. The percentage of organic carbon was estimated with a CHN analyzer (Elementar Vario III) and the organic carbon values were converted to organic matter using a conversion factor of 1.72 [21]. Energy content was obtained using the method adopted by Barnes [22].

Dry sieving method [23] was followed for the textural analysis of the sediment sample, since the beach sand contained less than 5% silt-clay fraction. The equipment for the analysis consisted of a stacked set of Wentworth grade sieves with. 5 phi () intervals within the range of 62–2000  . 80 gm of oven dried salt free sand sample was introduced into the 2000  sieve at the top of the stack and the stack transferred to the mechanical shaker and agitated for 15 minutes. After this, the material in each sieve was taken and weighed. Material (<62 m) that had fallen into the pan at the bottom of the sieve stack was also weighed. Fraction of the sediment (calculated in percentage weight of the total sample) was transferred into cumulative frequency series and then plotted as cumulative frequency curve. From the resulting sigmoid curve, the particle diameter was estimated and interpreted. From the cumulative frequency curve, graphic mean and standard deviation were calculated [24]. The nature of sediment was identified based on the verbal classification measure of Folk [24].

2.3. Hydrography

Interstitial water was siphoned out from the pit that was made by the core sampler employed to collect the sediment. These samples were analysed for temperature (with a high precision mercury thermometer), salinity (using a calibrated salinometer), pH (with Pocket pH Tester, Eutech) and dissolved oxygen by Winkler’s method[25].

2.4. Ecological Indices

Diversity Indices. Species diversity was estimated according to Shannon-Wiener (log2) [26]), Margalef’s richness (), [27] and Pielou’s evenness () [28] indices. The nematode-copepod ()/ ratio was calculated by dividing the number of nematodes ( trophic group or epigrowth feeders) in a sample by the number of copepods [29].

2.5. Statistical Analysis

The software programmes namely, SPSS (Version 11.0) and PRIMER V6 (Plymouth Routines In Multivariate Ecological Research, Version 6) were used for univariate and multivariate analysis of data. Pearson product-moment correlations were performed to test the relationship between meiofaunal distributions and environmental variables. Correlations were considered significant at .

BEST Analysis. The BEST routine available in PRIMER v6 combines the BIO ENV and BV STEP procedures of PRIMER v5. This routine uses all the available environmental variables to find out the combination that “best explains” the patterns in the biological data. PRIMER v6 for windows was used for the analysis of community structure.

Similarity Indices. Multivariate analysis consisted of estimating Bray-Curtis similarity after suitable transformation of sample abundance data. The similarity matrix was subjected to both clustering (hierarchical agglomerative method using group average linking) and ordination (nonmetric multidimensional scaling, MDS) using PRIMER v6. The significance of the clusters formed was tested by the similarity profile (SIMPROF) test.

3. Results and Discussion

In marine benthic ecology, sediment granulometry, along with environmental parameters is considered essential to determine the composition and characterization of benthic organisms [30]. The sediment characteristics and the interstitial faunal density are presented in Tables 1 and 2, respectively. The nature of sediments varied from very well sorted to moderately sorted. Moderately well sorted medium sand harbours the most diverse meiofauna [9]. Coarser grains have small retention capacity than fine ones, but, on the other hand, they have larger capacity of circulation [31]. Therefore, any major change in the sedimentary structure results in a change in the community structure [32]. Absence of meiofaunal taxa in post-monsoon at Fort Kochi beach (Table 2) accompanied by a drastic reduction in the medium sand fraction percentage corroborates the statement.

But a single factor alone cannot be held responsible for meiofaunal reduction. The pore volume in poorly sorted sediments is often reduced to 20%. The poor sorting of sediments at Sakthikulangara beach (Table 1) implies that the transportation agencies were incapable of separating the sediments into different class sizes and differential settling rates may be prevalent in the area. The sandy beaches of Kerala are subjected to seasonal erosion and accretion. Due to various oceanographic phenomena like upwelling and storm surges accompanied by the south west and north east monsoon, the morphology, beach profile, and the physical characteristics of the beaches alter drastically [3]. The variations noticed in the textural characteristics of the beach sediment are due to the physical changes that happen to these sandy beaches periodically.

Besides texture, organic content of sediment is an important factor that determines the infaunal distribution, since it is a dominant source of food, directly for deposit feeders and indirectly for suspension feeders. It is generally believed that sediments with coarse particles are deficient in organic matter whereas fine grain sediments show organic enrichment [30]. However, organic matter correlated significantly with coarse sand (Pearson correlation ; ) at the studied locations. Unlike the earlier findings that coarse sediments contain lesser organic carbon than fine sediments [33], the present study revealed that the angularity of sediments could play a significant role in the amount of organic carbon adsorbed to the sediments. More angular sand particles can harbor more carbon than smooth grains [9, 34]. Consumption of organic matter bound sediments and subsequent generation of faecal pellets will alter the mechanical composition of sediments. The amount of energy available in different beaches during the period of study ranged from 0.2245 to 16.08 joules/mg (Figure 2). A significant relation was not evident between organic carbon and the numerical abundance of interstitial organisms (Table 4). Ansari and Parulekar could not observe any consistent pattern between organic carbon and benthic biomass [15]. It is therefore clear that the quality of organic matter matters.

The hydrographic parameters such as temperature, pH, salinity, and dissolved oxygen are presented in Table 3. The fluctuations in temperature of the interstitial water medium are not very significant in tropical habitats and a temperature fluctuation of 5°C is easily endured by the interstitial animals [35]. Of the several hydrographic parameters of the medium that influence the occurrence and distribution of the interstitial forms, salinity is by far the most important in a tropical habitat. An important feature of the study area is the influence of south west monsoon, which affects the hydrographic condition in a remarkable manner. Wide fluctuations in salinity exert great stress upon interstitial organisms. Low salinity characterized Fort Kochi beach (Table 3).

The beach here is narrow and the lower salinity is due to the proximity of the channel leading to the backwaters and the consequent mixing with brackish water [36]. The inflow of fresh water greatly affects the interstitial fauna inhabiting the area. Increased salinity is favourable to harpacticoids [37]. Low salinity at Fort Kochi beach is marked by the absence of harpacticoids in the monsoon and post-monsoon seasons of the year (Tables 2 and 3). Great reduction of meiofauna abundance has been recorded during periods with reduced salinity, after heavy rain [10, 38]. In tropical beaches, although the seasons are not markedly delineated, periods of “good weather” and “storms” influence vertical migration of meiofauna [39]. In the marine biome, pH plays only a minor role for interstitial fauna, since the slightly alkaline sea water is well buffered against pH fluctuation. Dissolved oxygen of the interstitial water ranged from 3.68 to 8.05 mg/L. Dissolved oxygen in a sandy beach is dependent on the drainage of capillary water and the evaporation of interstitial water in the surface layers. Several reports are available on the correlation between the meiofaunal abundance and dissolved oxygen [40, 41].

Beaches are highly dynamic ecosystems and the study of meiofauna in beach ecology is important for the understanding of trophodynamic processes [42]. The relative abundance of interstitial organisms (annual mean) in different beaches is illustrated in Figure 3. While nematodes were the most abundant taxon in Cherai and Sakthikulangara beaches, the rest of the beaches under study were dominated by harpacticoid copepods. Other groups included turbellarians, polychaetes, oligochaetes, gastrotrichs, isopods, Kinorhyncha, and cnidarians. It could be seen that the distribution of total interstitial fauna and the different taxa was by no means uniform and the order of domination of taxa varied with beaches.

The Pearson correlation matrix between the various meiofauna taxa and the environmental variables is presented in Table 4. Total meiofaunal density is positively correlated with medium sand at 5% level of significance. Harpacticoids and turbellarians correlated at 0.05 levels with medium sand.

Rodríguez et al. [8] found that exposure time, desiccation, availability of food, sediment granulometry, tidal zonation, and interstitial water quality are the physical parameters that regulate the abundance of intertidal meiofauna. A multitude of environmental factors shape the meiofaunal abundance and distribution of a habitat. BEST analysis was done to find the best possible combination of variables that determined the faunal abundance and distribution. Biota-Environmental matching (BEST), as per the primer analysis (Figures 4 and 5), showed that organic carbon, particle size, and dissolved oxygen were the best combination of variables that determined the abundance and distribution of interstitial fauna in the present study.

Species diversity was estimated according to the Shannon-Wiener (log2) diversity, Margaleff richness (), and Pielou’s evenness () indices (Table 5).

Shannon diversity index was maximum (2.027) at Cherai and lowest at Sakthikulangara (1.144). Richness was at its peak in Cherai (0.988) and lowest in Sakthikulangara beach (0.599). The evenness component () varied in conformity with . Multivariate analysis consisted of estimating Bray-Curtis similarity after suitable transformation of sample abundance data. The similarity matrix was subjected to both clustering (hierarchical agglomerative method using group average linking) and ordination (nonmetric multidimensional scaling, MDS). MDS plot gave a good ordination with a stress value of 0.08 (Figure 6). Significant differences are observed in the distribution of taxa between beaches.

The nematode-copepod ratio was first used by Raffaelli and Mason [43] as a fast and reliable tool in monitoring the level of organic pollution. Despite several criticisms on this simplified relationship [4446], several researchers [4750] still use nematode-copepod ratio in assessing and identifying eutrophically enriched or polluted sediments. Benthic harpacticoid copepods are known to be sensitive to sediment metal concentration [51]. The nematode/copepod ratio [43] with modifications [29, 52] was used as an index to monitor pollution of selected beaches. Warwick [29], based on his studies, has suggested that an indication of pollution might be given by ratios around 40 ( ()/) for fine sediments and 10 for sands. The ()/ ratio value of around 12 indicated at Sakthikulangara beach (Figure 7) is indicative of some disturbance of this beach.

Sakthikulangara beach is characterized by natural radioactive pollution by thorium (based on secondary data) and faecal waste disposal. Smol et al. [54] have reported higher percentage of nematodes and lower percentage of harpacticoids within the sewage (TiO2) dumping area of Dutch coast. Ansari et al. [55] also found that copepods are more sensitive to environmental stress than nematodes in their study of interstitial fauna of the shelf region of south east coast of India. Though pollution tests have not been conducted and N/C value cannot be taken as a sole means of indicator of pollution, various features observed and analysed during the course of this study depict a different scenario of Sakthikulangara beach.

The assessment of ecological role and the spatial and temporal changes of meiobenthos can be used in environmental monitoring programs whose main goal is often to identify patterns in community structure and to relate them to measured environmental variables including pollutants [42]. Cherai beach supported maximum diversity and abundance of interstitial taxa. Further analysis with the help of available primary data would be carried out to explain the trophodynamics and ecology of the nematodes and harpacticoids of Cherai and Sakthikulangara beaches.

4. Conclusion

Nine interstitial taxa dominated by nematodes and harpacticoids abound the five ecologically vibrant beaches of Kerala. The organic matter in the beaches ranged from 0.01 to 0.745 and the grain size varied from 0.93 to 2.88. Organic matter correlated significantly with coarse sand (Pearson correlation ; ). The angularity of sediments plays a significant role in the amount of organic carbon adsorbed to the sediments. Total meiofaunal density correlated positively with medium sand (; ). Organic carbon, particle size, and dissolved oxygen determined the abundance and distribution of interstitial fauna as per multivariate BIOENV analysis. Diversity index (Shannon Wiener ) was maximum at Cherai (2.027) indicating the stability of environmental parameters and minimum at Sakthikulangara (1.44) characterized by regular physical changes. The value of nematode/copepod ratio ( ()/) obtained for Sakthikulangara beach indicate disturbed nature of this beach. Nonmetric multidimensional scaling, MDS (stress = 0.08), and cluster analysis (similarity profile, SIMPROF test) revealed significant differences in special distribution patterns of various interstitial taxa.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Priyalakshmi Geetha is thankful to the authorities of Bharata Mata College, Thrikkakara, Kochi, for the assistance and help and N. R. Menon the EU Incolab Grant no. 295092 for the work. They are also thankful for the constructive comments and suggestions of two anonymous referees on an earlier draft of the paper.