Forensic studies of phytoplankton ecology of two water bodies of Kurukshetra area of Haryana, State in India

In forensic investigations of death cases in which the body has been recovered from a water body, there is an urgent need to prove whether the death is due to drowning or not. Moreover, another important issue arises, whether a person was alive at the time of drowning (anti-mortem drowning) or he was thrown dead into the water body (post-mortem drowning) To answer these questions, phytoplankton species of forensic importance from two water bodies (Braham Sarovar and Prithus Pools) were analyzed to develop the Phytoplankton maps (P-maps). The P-maps are basically a detailed list of categories of species which are common, site-specific, rare and seasonally occurring phytoplanktons along with their photomicrographs. In the present study, water samples have been collected from Braham Sarovar and Prithus Pools situated at Kurukshetra area of Haryana state, India. The phytoplankton species present in both ponds were extracted by mild acid digestion (5% HCl). The physico-chemical properties of these water bodies such as pH, water temperature, total dissolved solids, electrical conductivity, hardness, alkalinity, calcium, sodium, potassium, and chloride ions have been recorded to establish a correlation with the distribution of phytoplankton using principal component analysis. Total 138 species of 59 genera of phytoplanktons belonging to Cyanophyceae, Chlorophyceae, Bacillariophyceae, Euglenophyceae and Dinophyceae have been observed. In Braham Sarovar, Cylindrospermum sp., Epithemia turgida, Eunotia rhomboidea, and Westella botryoides were the site-specific phytoplanktons whereas in Prithus pool, the site specific phytoplanktons were Leptolynbya granulifera, Arthospira jenerii, Cymbella cymbiformis, Frustulia vulgaris. It is hoped that results/data obtained from the present study (P-maps) may go a long way in solving the mysteries related to the deaths due to drowning and can be usefully used as a reference data for determining the exact site of drowning.


Background
The recovery of a dead body from any water body does not always mean that the death is due to drowning. There can be other reasons also. In forensic cases, it is urgently required to investigate the real cause of death, whether it is due to drowning or not. In addition, few important questions arising in relation to it are: Whether the drowning is anti-mortem or post-mortem?
The place where actual drowning took place? To answer these questions, a diatom test has been considered to be the best marker of drowning but it confronted with many difficulties and hardships in terms of its reliability. It was shown by some of the workers that diatoms were also present in non-drowned victims due to postmortem invasion, autopsy contamination (Lunetta et al. 2013), inhalation of diatoms from the air (Otto 1961;Spitz et al. 1965), water and foodstuff (Yen and Jayaprakash 2007). However, the reliability of diatom test has been established by the majority of researchers by considering qualitative and quantitative criteria (Auer and Möttönen 1988;Ludes et al. 1999;Pollanen 1997;Singh et al. 2005Singh et al. , 2006aSingh et al. , 2006bThakar et al. 2011;Thakar and Singh 2009;Thakar and Singh 2010). In 1988, Auer and Mottonen suggested that minimum 20 number of diatom per microscope slide from lung samples is required to define the death due to drowning. Similarly Ludes et al. (1999) also suggested the limit of 20 diatoms per slide per 100 μL pellet for lungs and 5 diatoms per slide per 100 μL of a pellet sediment extracted from 2 g of tissue samples such as brain, kidney, liver and bone marrow. But under certain environmental conditions, diatoms may be scanty or absent in a particular water body, leading to a false negative diagnosis of drowning. Many studies on qualitative and quantitative assessment of diatoms during different seasons and at different time periods have reported lesser number or absent of diatom species in one season and their bloom in the other (Canini et al. 2013;Cumming and Smol 1993;Singh et al. 2010;Sushanth and Rajashekhar 2012). Pollanen (1997) also reported highly positive diatom test outcome in the body of drowned victims in the months in which diatom population in the water body was at a peak. This showed variations in the number of diatoms over a period of time in water bodies.
The forensic usefulness of phytoplankton came to light when Chardez and Lambert (1985) (cited in Lunette and Modell 2005) detected the presence of phytoplankton other than diatoms in the drowned bodies. They also found aquatic organisms belonging to division Chlorophyta, Dinoflagellates, invertebrates, protozoan ciliates in the blood of drowned victims. Later, Yoshimura et al. (1995) developed a method to recover phytoplankton from the tissue of drowned victims by using soluene-350. Various species of green algae (Staurastrum and members of Zygnemataceae) as well as diatoms (Navicula, Cymbella and Melosira) were detected in both putative water samples and tissue samples of drowned victims. In an another case a serious attempt has been made to link the suspects and victims to the crime scene by analyzing the sediment encrusted sneakers from both assailants and victims along with numerous species of diatoms and scaled Chrysophytes (planktonic algae) from the sneakers as well as reference samples of pond sediment (Siver et al. 1994).
Further in Díaz-Palma et al. 2009 suggested to extend the diatom test to detect drowning by taking into consideration other groups of phytoplanktons. They have standardized the methods for recovery of microalgae from drowned victims by using different extraction procedures and concluded that dinoflagellates and some chlorophytes have cell walls and other resistant structures similar to diatoms. As a result, these organisms can be found in tissues of drowned victims.
Phytoplanktons are the photoautotrophic organisms which prevail in all types of water bodies such as streams, lakes, rivers, and ponds and are an important component of the food chain in the aquatic environment. Phytoplankton groups may be distinguished from each other based on a combination of characteristics, including photosynthetic pigments, shape, size and cell covering. They are microscopic in nature and "invisible" to perpetrators of crimes, thus it helps in linking a suspect to the crime scene (Siver et al. 1994).
Keeping in mind the needs of environmental and forensic experts, the present study has been undertaken to study the phytoplankton diversity along with seasonal variations in two forensically important ponds (Braham Sarovar and Prithus Pools) of Kurukshetra region of Haryana, India. These two ponds are forensically important because of two reasons: Firstly, sufficiently large amount of water remains in these unguarded ponds throughout the year, and secondly, a large crowd gathers around them to take a dip during the year on many occasions. The data generated in the present study has been used to generate Phytoplankton-maps (P-maps). The characteristics Phytoplankton-maps of these two water bodies will hopefully help in unfolding the mysteries of deaths due to drowning and the data generated can be utilized as a reference data for determining the exact site of drowning.

Selection of sampling sites
In the present study the following two sites which are situated in Kurukshetra city and Pehowa city of Haryana, India have been selected to collect samples for the analysis of phytoplanktons:

Brahma Sarovar
It is the largest man-made water body situated in historically important Kurukshetra city, Haryana state, India (29°57′37″ N 76°50′42.2″ E). The water body is rectangular in shape, measuring 3600 ft long, 1500 ft wide and 15 ft deep. It has religious significance, every year on 'Somavati Amavasya'(a special day in Hindu religion) and on the occasion of solar eclipses; millions of people from different parts of India visit this water body to take a holy bath (Fig. 1).

Prithus pool
It is located in another religiously important city, Pehowa city, Haryana, India (29°58′51.6″N 76°34′46.4″E) about 25 Kms from Braham Sarovar. It was constructed by king Prithu after the death of his father. This place is known for funeral ceremonies (pind daan in Hindu religion), therefore, thousands of pilgrims from different places visit this pool throughout the year (Fig. 2).

Collection of samples
For phytoplankton analysis, the water samples have been collected during four seasons such as summer (May-June), autumn (August-September), winter (November-December) and spring (February-March) from two selected sites from the year 2014 to 2016. The water samples were collected in a sterilized plastic bottle of 1 l capacity and preserved in 4% formaldehyde solution.
The procedure adopted to prepare slides for the analysis of phytoplankton Approximately, 50 ml of water sample was taken in a beaker and digested with 10 ml of 5% HCl for 24 h. Next day, after thorough mixing of the water sample and HCl, 10 ml of water sample was taken in a centrifuge tube and centrifuged at 3000 rpm for 10 min. The supernatant was collected out, leaving behind only a residual material at the bottom of the tube. An excess of water was again mixed with residual material and centrifuged again as above. This process was repeated twice. After discarding the supernatant, the residual material (pellet) containing phytoplankton was poured over two sets of five serially marked (I-V) microscopic slides and smear was prepared. Then, these slides of two sets were permanently mounted using a high refractive index mounting material DPX to study diatoms and Glycerin jelly for phytoplankton analysis under the low and high power of the microscope.

Analysis of slides prepared
All the slides prepared for both the water bodies season wise have been examined with the help of an optical compound microscope (Nikon Eclipse microscope) fitted with light source under low (10× X 10×), high power (10× X 40×) and oil immersion (10× X 10×) magnifications. Nikon E200 series digital camera which has been attached to the microscope was used to prepare Photomicrographs for each and every phytoplankton species. These Photomicrographs were used for taxonomic identifications of phytoplankton with the help of publications given by Prescott (1962); Komárek et al. (1983); Desikachary (1959); Philipose (1967) and Wehr et al. (2015).
Data so obtained have been used to construct a phytoplankton map (p-map) for the both water bodies. Enumeration of phytoplankton was done by using Haemocytometer. The colonial and filamentous microalgae such as Microcystis sp., Merismopedia sp. and Oscillatoria sp. have been counted by using plating method (Stanier et al. 1971). Average of two years (2014-2016) data has been given in Tables 1, 2 and 3.

Physico-chemical properties
In the both selected water bodies, the parameters such as total alkalinity, chloride and calcium levels, and total hardness have been determined according to APHA (American Public Health Association) manual (APHA 2005). Sodium and potassium contents in the water samples have also been determined with the help of flame photometer. Electrical conductivity, pH, and temperature have also been determined by using the pH meter cum conductivity meter (Aquapro digital water tester model AP-2). TDS has been determined by using TDS meter (Aquapro digital water tester model AP-1). The correlation between physico-chemical parameters and phytoplankton species was carried out by Principal Component Analysis (PCA) using the Microsoft Excel Add-Ins software XLSTAT.  Diversity indices: In the present study, the diversity indices of phytoplankton species have been calculated in terms of Species diversity method given by Shannon (2001), Species richness by Menhinick (1964) and Evenness by Pielou (1975) as follows: Shannon-Wiener diversity index: Shannon diversity index is a quantitative measure of species diversity in a community and can be calculated as follows: Where, H = Shannon and Weaver diversity index. pi = ni/N, ni being the number of individual of ith species and. N = total number of individuals in the sample. Species Richness: Species richness is a count of a number of species in a water body. It is calculated using the following formula: Where D = Species Richness, n = total number of species in a sample. N = Total number of an individual organism in a sample.
Species Evenness: Species evenness is the measures of how equals the relative abundance of different species  are in a community. It is calculated using formula as follows: Where, J = Evenness H = Shannon and Weaver diversity index. S = Total number of species in a sample.

Phytoplankton-maps (P-maps)
The data generated in the present study were used to prepare the P-maps (Tables 4 and 5) of two selected water bodies (Braham Sarovar and Prithus Pools). The P-maps depicts the detailed list of phytoplankton genera and species present in the two selected water bodies. P-maps have been divided into three parts-A, B, and C as shown in Tables 4 and 5. In part -A, phytoplanktons species have been shown, which have been divided into the following four categories: ▪ Commonly found phytoplanktons ▪ Site-specific phytoplanktons, ▪ Seasonally occurring species of phytoplanktons ▪ Rarely occurring phytoplankton The detailed information about these categories has been given below and is helpful in generating P-map for a specific water body: ❖ Commonly found species: Commonly found species are those which are common in more than one water body. They have significance in detecting the death due to drowning and are of little significance in diagnosing drowning site. ❖ Site-specific species: Site-specific species are those which are restricted to a particular water body, hence they have great forensic significance in locating the exact site of drowning, particularly when a body is found away from the actual site of drowning and also they help in linking the suspect and victim to the crime scene. ❖ Seasonally occurring species: This category includes those species which are dominant in one season and may become extinct or dominated by other species during another season in a particular water body.
These season specific species of a particular water body may help in determining the season of occurrence of crime. ❖ Rarely occurring species: These are those species which are found occasionally or by chance in a particular water body. These species are not important from the forensic point of view. In part-B of P-map for both the water bodies, the photomicrographs of Site-specific species have been shown (Tables 4 and 5), whereas part -C shows Distribution of Phytoplankton species during the different season in both water bodies (Tables 4 and 5).

Results and discussion
Phytoplanktons are the diverse group of photosynthetic organisms present in all the aquatic habitats such as lakes, rivers, ponds, canals, sea etc. Moreover, they have different shapes, sizes, colors and forms, which provide them individuality. So their presence in water bodies and potential to identify the site of drowning may help in the forensic investigation. Because of their microscopic nature and invisibility phytoplankton may adhere to the suspect, victim or their clothes or weapon of offense. In addition, different species of phytoplankton present during different seasons which aids in linking the suspect to a particular crime scene during the specific season (Hall 1997). Therefore, seasonal variation of phytoplankton has been studied in two selected water bodies of Kurukshetra region, continuously for two years to generate a baseline database of phytoplankton (Phytoplankton maps) for each water body.
In the present study, remarkable changes in diversity and abundance of phytoplankton have been observed during different seasons. In Braham Sarovar and Prithus Pools, a total of 138 species belonging to 59 phytoplankton members have been found (Tables 1 and 6). Individually in Braham Sarovar alone 47 genera belonging to 5 classes i.e. Cyanophyceae (12), Chlorophyceae (18), Bacillariophyceae (9), Euglenophyceae (4) and Dinophyceae (4) representing 70 species have been found. Similarly, in Prithus pool, 36 genera belonging to 4 classes' i.e. Chlorophyceae (15), Cyanophyceae (12), Euglenophyceae (4), and Bacillariophyceae (5)     green algae, 6 species of Cyanobacteria and 6 species of other algae in three aquatic locations of Gwang-ju area (South Korea) where drowned victims were found. The detailed list of species recovered in the present study has been given in Table 1. Figures 3 and 4 show photomicrographs of different species of phytoplankton observed in Braham Sarovar and Prithus pools water bodies.

Class -Chlorophyceae
The members of Chlorophyceae formed the most dominant group among others in both water bodies which accounted for 46% of total phytoplankton diversity with 32 species in the Braham Sarovar (Table 4c) and 54.4% of total diversity was represented by 37 species of Chlorophyceae in the Prithus Pool (Table 5c). Panigrahi and Patra (2013) also reported dominance of Chlorophyceae over other phytoplankton groups in terms of numerical value and percentage composition in a water body of Cuttack city, West Bengal, India. The dominant species of Braham Sarovar were Ankistrodesmus fusiformis, Ankistrodesmus flactum, Crucigenia tetrapedia, Oocystis marsonii, Scenedesmus bijugatus, Scenedesmus abundance, Tetrastrum minimum, Coelastrum Microporum.  Similarly, in Prithus pool, Cyanophycaeae also ranked as the second most abundant group accounting for 21% of the total population with 14 species (Tables 4c and  5c). The dominating species were found to be Merismopedia tenuiasima, Spirulina sp., Oscillatoria agardhii and, Plantothrix mougeotti.

Class-Bacillariophyceae
The class Bacillariophyceae has been ranked third in Braham Sarovar and accounted for 17.1% but the same class ranked fourth in Prithus Pool and accounted for 10.2% of total phytoplankton diversity found in this water body (Tables 4c and 5c). Total 12 species of diatoms have been found from the Braham Sarovar with Aulacoseira granulata, Nitzschia acicularis, Nitzschia palecea and Navicula cryptocephala being the dominant ones, whereas seven species of diatoms reported from Prithus pool Frustulia vulgaris, Melosira varians Frustulia vulgaris and Nitzschia palea were dominant.

Class-Euglenophyceae
Euglenophyceae class constituted around 15% (10 species) of the total phytoplankton population and formed the third most dominant group in Prithus pool, whereas it ranked fifth in Braham Sarovar accounting for 6% (5 species). Among them, Phacus longicauda, Phacus tortus and Lepocinclis fusiformis were the most abundant species in the Braham Sarovar, while Lepcinclis ovum, Euglena acus and Euglena virdis were dominant in Prithus pool (Tables 4c, 5c).

Class-Dinophyceae
Members of Dinophyceae have been only observed in the Braham Sarovar water body and constituted 10% (6 species) of species diversity with Peridinium cinctum being the dominant one.

Seasonal variation in phytoplankton population density
Mean density of 1.662 × 10 3 ± 227 phytoplankton individuals mL − 1 of water was observed at Prithus pool while in Braham Sarovar mean density has been calculated to be 1.388 × 10 3 ± 248 individuals mL − 1 . In both water bodies, Chlorophyceae was found to be the most abundant group in summer which declined in autumn and winter. A similar observation has been made by Shinde et al. (2012) in Harsool-savang dam, Aurangabad; where maximum density of Chlorophyceae was observed in summer season. Nowrouzi and Valavi (2011) while studying phytoplankton of Lake Kaftar in Iran also reported highest density of Chlorophyceae in summer and spring and the lowest in winter. Cyanophyceae in Prithus pool exhibited maximum growth in summer and minimum in winter while, in Braham Sarovar, Cyanophyceae exhibited maximum growth in winter and minimum in the autumn. The maximum density of Cyanophyceae in Braham Sarovar during winter was due to the overgrowth of Microcystis aeruginosa. Khuantrairong and Traichaiyaporn (2008) also observed Cyanophyta as the most abundant group in winter and summer in Doi Tao Lake, Thailand. On the contrary, Venkateswarlu (1969) reported the dominance of Cyanophyceae in summer season in Mossi River, Hyderabad. Singh and Swarup (1980) also reported maximum growth of cyanophyceae during summer in Lake Suraha, Uttar Pradesh. The maximum abundance of Bacillariophyceae was recorded in winter and minimum in autumn. Singh et al. (2010) observed highest density of Bacillariophyceae in winter in Mawatha lake Rajasthan. A similar trend has also been observed by Baba and Pandit (2014) while studying the diversity of phytoplankton in Wular Lake Kashmir. Members of Euglenophyceae were observed both in autumn and winter in Braham Sarovar, while in Prithus pool, members of Euglenophyceae were maximum during winter and spring and minimum in summer. The members of Dinophyceae were detected only in Braham Sarovar, that too during winter and spring season. Mean density distribution of different classes of phytoplankton in Braham Sarovar and Prithus Pools is Part-C of P-map and has been shown in Tables 4 and 5.
Various environmental factors such as light, temperature and rainfall and physico-chemical factors (pH, temperature, salinity, Nitrogen, Phosphate, Calcium, Sodium, Phosphate, Chloride, etc) of a water body get changed with seasons which in turn bring about change in composition and abundance of phytoplankton (Figueredo and Giani 2009). In the present study, the maximum density of Chlorophyceae and Cynaophyceae has been recorded during the summer season and minimum during autumn and winter season. The increase in the density of phytoplankton may be attributed to the increase in temperature during summer, which not only increases the rate of photosynthesis, but also increases the concentration of available nutrients due to increase in water evaporation. Whereas a decrease in the density during autumn and winter may be due to the increase in the water level of the water body by rainwater which decreases the photosynthetic efficiency of phytoplanktons (Shinde et al. 2012). The density of Bacillariophyceae was observed maximum during winter which may be due to a higher rate of multiplication of Bacillariophyceae at low temperature (Raibole and Singh 2011). The members of Euglenophyceae were observed during autumn and winter season only, may be because of two reasons i.e. either they cannot tolerate higher temperature or increased accumulation of organic matter in the water bodies due to heavy rainfall, favors the growth of Euglenophyceae (Tiwari and Chauhan 2006).

Diversity index
High values of Shannon's index (H′) were recorded during summer and lower during winter at both water bodies. Shannon's index has been calculated to be high (3.16 ± 0.013) in Braham Sarovar as compared to Prithus pool (3.03 ± 0.26) ( Table 7). The value of Species richness of Braham Sarovar was comparatively high during summer (4.55) and minimum during winters (3.87) and autumn (3.96) while in Prithus pool, species richness was maximum during spring (3.86) and minimum during autumn (3.42). Species evenness in Prithus pool showed little variation during different seasons, indicating that all the species were equally distributed throughout the year. In Braham Sarovar, the highest value of species evenness was 0.93 during summer and lower during winter (0.86).

Physico-chemical properties
Physico-chemical parameters of Braham Sarovar and Prithus pools have been assessed to determine the tropical status of the water body and their co-relation with the phytoplankton species diversity. The physico-chemical parameters of Prithus pool were numerically high as compared to Braham Sarovar. The high value of physico-chemical parameters was probably due to the anthropogenic activity at this pond. The physic-chemical parameters of both aquatic bodies during the different seasons are shown in Tables 8 and 9.
Correlation between physico-chemical parameters and phytoplankton community Principal component analysis (PCA) Principal Component Analysis (PCA) is a multivariate data analysis technique which reduces the multi-dimensional data into two principal components. In present study Principal Component Analysis has been used to establish the degree of correlation between physico-chemical parameters (as mentioned above) and the phytoplankton species prevailing in the water body. It helps to understand whether phytoplankton has a positive or negative relationship with different physico-chemical parameters and what are the key factors responsible for their growth and inhibition?
As in forensic investigations, it is common to find a dead body away from the actual site of drowning due to water current, tides, etc. Thus, determination of the actual site of drowning becomes important for the medico-legal expert. The abundance and diversity of phytoplankton vary between different water bodies due to the difference in various parameters such as pH, temperature, salinity, electrical conductivity, hardness, alkalinity etc. of that water body (Nowrouzi and Valavi 2011). Some of the physico-chemical parameters are favorable for the growth of one species while these may not support growth of others. The different species of phytoplankton get established in a particular water body based on their preferences for particular physico-chemical components of water. Therefore, a positive test for drowning should not only be based upon the presence or absence of phytoplankton in the body of a drowned victim, rather it should also exhibit the same type of phytoplankton in similar concentrations of the drowning medium (Sasidharan and Resmi 2014). A detailed study of phytoplanktons and their relationship with characteristics of the water body can play a significant role in connecting a drowned dead body to a specific water body source particularly when the body is found outside or lying adjacent to the water body.
In the present the study, PCA method has been employed to 10 physico-chemical parameters for 20 dominant species of phytoplankton for Prithus pool. The PCA axis 1 and 2 explained 42.81 and 36.70% variation in phytoplankton and physico-chemical parameter. Crucigeniella quadrata, Crucigenia tetrapedia, Scenedesmus dimorphus, Scenedesmus quadricauda, Planktothrix mougeotti, and Melosira varians are positively co-related to pH, hardness, alkalinity, chloride, electrical conductivity and temperature of water sample while Crucigenia crucifera is negatively co-related with all these factors. Dictyochloropsis reticulata, Tetrastrum punctatum, and Frustulia vulgaris are positively co-related to TDS (total dissolved solids) and sodium while Nitzschia palea is negatively co-related to TDS. Lepocinclis ovum and Nitzschia palea are positively correlated to potassium while Coelastrum microporum and Coelastrum astroideum are positively correlated to calcium (Fig. 5).
PCA for Braham Sarovar has been drawn by taking in to account 10 physico-chemical parameters and 14 dominant species of phytoplankton. The PCA axis 1 explained 50.14% and axis 2 explained 29.03% variations in phytoplankton and physico-chemical parameters. Aulacoseira granulata, Tetraedron minimum, Gloeocapsa minimum, Nitzschia palecea, and Coelastrum microporum, is positively correlated to temperature whereas Microcystis aeruginosa and Synechococcus sp. are negatively correlated. Chroococcus dispersus and Aphanothceae sp. have been positively correlated with calcium and Ankistrodesmus fusiformis is positively correlated with chloride. Phacus longicauda, Phacus tortus, Lepocinclis fusiformis, and Synechococcus sp. are positively co-related with pH, potassium, hardness and electrical conductivity whereas they all are negatively correlated with TDS and calcium. Similarly, Synechococcus sp. and Microcystis aeruginosa are positively co-related with pH and alkalinity (Fig. 6).

Phytoplankton maps (P-maps) of Braham Sarovar and Prithus pools
Among the total 138 species, 13 species were common in both the water bodies, therefore, these are categorized as common species and are significant in diagnosing the drowning death, but not significant in determining the exact site of drowning. Similarly, rarely occurring species in Braham Sarovar (25 species) and Prithus pools (20 species) are not forensically significant. The site-specific species in Braham Sarovar (4 species) and Prithus pools (4 species) are forensically significant in determining the exact/putative site of drowning. The P-maps with a list of common and site-specific phytoplanktons have been shown in Tables 4 and 5. The P-map generated of phytoplankton species of two water bodies can help in determining not only cause of death, but also helps in linking the suspects to a crime scene. Sarovar and Prithus pools of the historically important cities (Kurukshetra and Pehowa) of Haryana State, for two years continuously. The data analysis differentiated phytoplankton into common, rare, and site-specific species in order to construct phytoplankton-map. The P-maps generated for the two water bodies of Haryana state, India, are expected to be useful to forensic experts and forensic medicine experts in solving drowning cases and unfolding mysteries related to the putative site of drowning, especially when diatoms are present in very less number during a particular season. In addition, the data has been successfully used to find the degree of correlation between different environmental parameters and phytoplanktons prevailing in a water body with the help of Principal Component Analysis method.