Comparison between phytoplankton bio-diversity and various indices for winter monsoon and inter monsoon periods in north-eastern Arabian Sea

Similar documents
Supplementary Figure 1 Surface distribution and concentration of the dinoflagellate N. scintillans

Seasonal cycle of phytoplankton community composition off Newport, Oregon, in 2009

Fluorometry Project Chlorophyll Temperature Time Series

Comparison of chlorophyll concentration in the Bay of Bengal and the Arabian Sea using IRS-P4 OCM and MODIS Aqua

2001 State of the Ocean: Chemical and Biological Oceanographic Conditions in the Newfoundland Region

Diversity of Zooplankton in some Reserviours in and around Karwar- Uttara Kannada District Karnataka

Ocean Boundary Currents Guiding Question: How do western boundary currents influence climate and ocean productivity?

FLOW: Amigos de Bolsa Chica Citizen Science Program. Plankton Collection and Identification Report

Ireland: Predictions. NMP Current closures ASP AZP DSP PSP

VIKRANT A. KULKARNI. CRZ and Marine Environment Consultant. emissions.

Ireland: Current Conditions

Ireland: Current Conditions

Turbulence and the Spring Phytoplankton Bloom

TEACHER VERSION: Suggested student responses are included. Seasonal Cycles: the North Atlantic Phytoplankton Bloom

manuscript associated with each of your comments are described below. Your comments are in bold

Multivariate time-series forecasting of the NE Arabian Sea Oil Sardine fishery using satellite covariates

Species composition and long-term dynamics of potentially toxic dinoflagellate species in benthic assemblages of Peter the Great Bay, Sea of Japan

2. What is a phytoplankton bloom and when does it generally occur in the North Atlantic?

INTRODUCTION MATERIALS AND METHODS

Recent Harmful Algal blooms (HABs) Events in Indonesia

Homework 5: Background Ocean Water Properties & Stratification

OUR COUNTRY INDIA LOCATIONAL SETTING

Management of recreational waters in relationship with harmful microalgae blooms (HAB) in the Mediterranean Sea

CHAPTER VI GENERAL CONCLUSION

Government of India India Meteorological Department (Ministry of Earth Sciences) Subject: Current weather conditions and forecast for next two weeks

SPATIAL CHARACTERISTICS OF THE SURFACE CIRCULATION AND WAVE CLIMATE USING HIGH-FREQUENCY RADAR

The bloom of the dinoflagellate (Noctiluca miliaris) in the North Eastern Arabian Sea: Ship and Satellite study

UNIT 5: ECOLOGY Chapter 15: The Biosphere

Announcements. Third problem site due November 30 Review calculations

Phytoplankton. Zooplankton. Nutrients

Ireland: Predictions. NMP Current closures ASP AZP DSP PSP

Biological survey of species diversity of sea grass beds in selected sites of southern Sri-Lanka.

Satellite-derived environmental drivers for top predator hotspots

Long Range Forecast Update for 2014 Southwest Monsoon Rainfall

Vancouver Lake Biotic Assessment

53 contributors for 35 individual reports in 2009 show 5% of figures today

Tree distribution and Habitat preference with respect to the Elevation, Rain and Soil type of Western-ghats region of Karnataka.

Honors Biology Unit 5 Chapter 34 THE BIOSPHERE: AN INTRODUCTION TO EARTH S DIVERSE ENVIRONMENTS

Remote sensing of evolution and coupling of green Noctiluca and diatom blooms in the northern Arabian Sea using value-added timeseries

The Microbial World. Microorganisms of the Sea

Climate. Annual Temperature (Last 30 Years) January Temperature. July Temperature. Average Precipitation (Last 30 Years)

Serial No. N4470 NAFO SCR Doc. 01/83 SCIENTIFIC COUNCIL MEETING SEPTEMBER 2001

Biogeochemistry of trace elements and isotopes in the Indian Ocean

Marine biodiversity - PBBT102

Impact of Alaskan Stream eddies on chlorophyll distribution in the central subarctic North Pacific* Hiromichi Ueno 1,

Forecasting inshore red tide blooms using recent past offshore conditions on the West Florida Shelf

Comparison of chlorophyll distributions in the northeastern Arabian Sea and southern Bay of Bengal using IRS-P4 Ocean Color Monitor data

Inter comparison of wave height observations from buoy and altimeter with numerical prediction

Upper Ocean Circulation

Name Period Part I: INVESTIGATING OCEAN CURRENTS: PLOTTING BUOY DATA

Parting the Red Seas: The Optics of Red Tides

Modeling Indian Ocean Biogeochemistry Iron Limitation and Dipole-Zonal Mode Impacts

Bugs in JRA-55 snow depth analysis

The Earth s Ecosystems

Ireland: Predictions. NMP Current closures ASP AZP DSP PSP

BIOLOGICAL OCEANOGRAPHY

Carbon Exchanges between the Continental Margins and the Open Ocean

C M E M S O c e a n C o l o u r S a t e l l i t e P r o d u c t s

DENSITY AND BIODIVERSITY OF BGA IN RICE FIELDS OF GOA

Ireland: Current Conditions

Biogeographic Approach to Coastal Assessments & Spatial Planning

Effects of environmental variables on midsummer dinoflagellate community in the Neva Estuary. Mikhail Golubkov, Vera Nikulina, Sergey Golubkov

Understanding Oceans, Gulfs & Tides

Life on Earth

Summary and Conclusions

Figure 14 p.385 5/11/2016. Plankton Production and Food Webs. Plankton Production and Food Webs

Primary Productivity (Phytoplankton) Lab

Fig Available seismic reflection, refraction, and magnetic profiles from 107 the Offshore Indus Basin close to the representative profile GCDH,

The seasonal and interannual variability of circulation in the eastern and western Okhotsk Sea and its impact on plankton biomass

Time-series observations in the Northern Indian Ocean V.V.S.S. Sarma National Institute of Oceanography Visakhapatnam, India

APPENDIX B PHYSICAL BASELINE STUDY: NORTHEAST BAFFIN BAY 1

Wave Climate Variations in Indonesia Based on ERA-Interim Reanalysis Data from 1980 to 2014

Lesson: Primary Production

SEAWIFS VALIDATION AT THE CARIBBEAN TIME SERIES STATION (CATS)

Advanced Hydrology. (Web course)

Factors impacting the formation & modification of sinking oil snow : Processes and Pathways

Lesson 2. Antarctic Oceanography: Component I - Ice/Glaciers Component II - Marine Snow

Marine Ecology I: Phytoplankton and Primary production

H U Solanlki, R M Dwivedi & S R Nayak

Project 1.3.1: Improved knowledge of biota, habitats and risks. Project Leader: Dr Mick Haywood, CSIRO Division of Marine and Atmospheric Research

Ireland: Current Conditions

THE ARABIAN SEA: A NATURAL EXPERIMENT IN PHYTOPLANKTON BIOGEOGRAPHY

STUDY AREA AND METHODOLOGY

Size scaling deviation in phytoplankton photosynthesis and the energy flow through a

First bloom record of toxic dinoflagellate Prorocentrum lima and climate change interactions in the Dardanelles (Turkish Straits Sistem)

Chiang Rai Province CC Threat overview AAS1109 Mekong ARCC

EPA Region 3 Mid-Atlantic State s Algae Identification Workshop

3. HYDROMETEROLOGY. 3.1 Introduction. 3.2 Hydro-meteorological Aspect. 3.3 Rain Gauge Stations

UPDATE OF REGIONAL WEATHER AND SMOKE HAZE (February 2018)

Marine Life. and Ecology. 2. From phytoplanktons to invertebates

Government of Sultanate of Oman Public Authority of Civil Aviation Directorate General of Meteorology. National Report To

Key Finding: Long Term Trend During 2014: Rain in Indian Tradition Measuring Rain

GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L14610, doi: /2007gl029633, 2007

Seasonal Climate Outlook for South Asia (June to September) Issued in May 2014

Abiotic Dominant Dominant Factors Plants Animals

Ocean Circulation. In partnership with Dr. Zafer Top

Trends and Variability of Climatic Parameters in Vadodara District

Seasonal and mesoscale variability of phytoplankton in the Arabian Sea. from satellite observations to models

Comparison of drop size distribution between stations in Eastern and Western coasts of India

Parallel Measurements of Light Scattering and Characterization of Marine Particles in Water: An Evaluation of Methodology

Transcription:

Indian Journal of Geo-Marine Sciences Vol. 43 (8), August 2014, pp. 1513-1518 Comparison between phytoplankton bio-diversity and various indices for winter monsoon and inter monsoon periods in north-eastern Arabian Sea 1 Gunjan Motwani, 2 Mini Raman, 3 Prabhu Matondkar, 3 Sushma Parab, 3 Suraksha Pednekar & 1 Hitesh Solanki 1 Department of Botany, Gujarat University, Ahmedabad, India 2 Space Application Centre, ISRO, Ahmedabad, India 3 National Institute of Oceanography, Goa - 403004, India [E-mail: gunjan_motwani@yahoo.co.in] Received 12 December 2012; revised 25 February 2014 Phytoplankton samples of north-eastern Arabian Sea were collected during ocean colour satellite validation cruises from 2003-2007. Phytoplankton community organization and distribution was analyzed using various diversity indices like Shannon Diversity Index, Simpson Diversity Index, Margalef Diversity Index, McIntosh Diversity Index, Pielou Evenness Index and dominance index. Results showed that Navicula, Thalasiothrix and Rhizosolenia were most abundant among the diatoms. Trichodesmium, a cyanobacteria, dominated the shallow coastal waters. Resultant indices were correlated with phytoplankton cell counts and it was found that Shannons index better represents the diversity than other indices. Cell counts were also correlated with in situ chlorophyll-a values, which showed that this correlation does not stand good for the bloom conditions. [Keywords: Phytoplankton diversity, diversity indices] Introduction Although 70% of the Earth's surface is occupied by the oceans, our knowledge of biodiversity patterns in marine phytoplankton is very limited in comparison to that of the biodiversity of plants on the land 1. It is well established that diversity enhances productivity and stability in communities of higher organisms; however, knowledge of such relationships between unicellular organisms like phytoplankton, which contribute to about 50% to the global primary productivity, is still lacking 2. A diversity index is a measure of species diversity within a community that consists of co-occurring populations of several (two or more) different species. It includes two components: richness and evenness. Richness is the measure of the number of different species within a sample showing that more the types of species in a community, the higher is the diversity or greater is the richness. Evenness is the measure of relative abundance of the different species with in a community. Present study is an attempt to analyze phytoplankton community structure (evenness and richness) using various diversity indices like Shannon Diversity Index, Simpsons Diversity Index, Marglef Diversity Index, McIntosh Diversity Index, Pielou Evenness Index and dominance index. The use of these indices is well established in ecological literature however choice of a suitable index depends on how well the richness and evenness is represented by each index. Materials and Methods The Site selected for this study was Northern and Eastern section of Arabian Sea, occupying an area 6.225 x 10 6 km Sq. and extends from 0º to 25ºN and 50º to 80ºE 3. It is bordered by Oman at the west, Iran at the North-West and the India at the east. Periyar, Bharathapuzha and Pamba rivers from Kerala; Kali, Netravati and Sharavati rivers from Karnataka, Tiracol, Chapora, Baga, Mandovi and Zuari rivers from Goa; Shastri, Gad, Vashishti, Savitri, Patalganga, Ulhas and Vaitarna rivers from Maharashtra; Tapti, Narmada, Mahi and Sabarmati rivers from Gujarat and Indus from Pakistan bring fresh water into Arabian Sea 4. Phytoplankton samples of North Eastern Arabian Sea were collected during ship cruises organized for ocean colour satellite validation from 2003-2007 as shown in figure 1. The details of samples collected and analyzed for various cruises are summarized in table 1.

1514 INDIAN J MAR SCI VOL 43 NO 8, AUGUST 2014 The water samples were collected from various sampling stations, which were decided on the basis of percent light intensity/penetration with reference to the surface irradiance in the water column. Satlantic under water Hyper-spectral radiometer was used to measure the light levels at the sampling sites. For microscopic identification and cell counts 500 ml of sea water was fixed with 1% lugol s Iodine and preserved in 3% buffered formaldehyde solution and stored under dark and cool conditions till analysis. Samples were concentrated approximately to 5-10 ml by siphoning the top layer of the sample carefully with a tube, 1 ml of sample were transferred to a Sedgwick-Rafter slide and identified and counted using an Olympus Inverted Microscope (Model IX 50) at 200% magnification. Standard taxonomic keys 5 were used for identification. The cruises were temporally categorized into Fig.I Map of stations at FORV-244, FORV-253, FORV-212 and FORV-222 winter-monsoon (Dec-March) and inter-monsoon (April-May) to better understand the seasonal variation in phytoplankton type and concentration. The species diversity indices such as Shannon Diversity Index, Simpson Diversity Index, Margalef Diversity Index, McIntosh Diversity Index, Pielou Evenness Index and dominance index were computed. Surface diversity (at 0 m depth) was evaluated for FORV-212 and FORV-22 whereas surface as well as depth wise (vertical profile) evaluation of diversity was carried out for FORV-244 and FORV-253. For analyzing the community structure of phytoplankton in the North-Eastern Arabian Sea these indices were calculated as follows: Shannon Diversity Index: This index is applied to biological systems very commonly for calculating diversity. It was derived from a mathematical formula by Shannon in 1948 6. H = -Σ [(ni / N) x ln (ni / N)] H : Shannon Diversity Index ni: Number of individuals belonging to i species Simpson Diversity Index: This diversity index is derived by Simpson in 1949 7,6. Simpson index values ( ) are between 0 1. But while calculating, final result is subtracted from 1 to correct the inverse proportion. 1 - = [Σ ni (ni -1)] / N (N-1) : Simpson Diversity Index ni: Number of individuals belonging to i species Margalef Diversity Index: It has no limit value and shows a variation depending upon the number of species. Thus, it is used for comparison between various sites 8. d = (S-1) / ln N d: Margalef Diversity Index S: Total number of species Cruise Id Table I Details of the cruise, season, sampling stations and depths Period of Study/ season No. of stations analysed Sampling Depths FORV-212 27 th February to 5 th March 2003, Winter-Monsoon 5 100% light level (Surface) only. FORV-222 21 th February to 11 th March 2004, Winter-Monsoon 3 100% light level (Surface) only. FORV-244 15 th April to 28 th April 2006, Inter-Monsoon 15 100% light level (Surface), to 1 % light level. FORV-253 28 th February to 11 th March 2007, Winter-Monsoon 8 100% light level (Surface), to 1 % light level.

MOTWANI et al.: COMPARISON BETWEEN PHYTOPLANKTON BIO-DIVERSITY AND VARIOUS INDICES 1515 McIntosh Diversity Index: It was suggested by McIntosh in 1967. The value ranges between 0 1. The closer is the value to 1; the more homogeneous is the distribution of phytoplankton 9. Mc = [N - (Σni 2 )] / [N - N] Mc: McIntosh Diversity Index ni: Number of individuals belonging to i species Pielou Evenness Index: It was derived from Shannon index by Pielou in 1966. The ratio of Shannon index to the maximum value gives the Pielou Evenness Index result. The values range between 0 1. The closer is the value to 1; the more even is the distribution of phytoplankton 10. J = H / H max J : Pielou evenness index H : The observed value of Shannon index H max: lns S: Total number of species Dominance Index: Dominance index ranges from 0 to 1. 0 represents that all taxa are equally distributed and 1 means one taxon dominates the community completely 7. D= Σ ((n i /n) 2 ) n i : Number of individuals of taxon i. Results and Discussions Phytoplankton Richness The phytoplankton community structure of the Arabian Sea was highly diverse with 274 species identified. Diatoms (Bacillariophyceae) exhibited the greatest diversity with 142 species followed by dinoflagellates (Dinophyceae) 129 species; other Algae (Cyanophyceae) with 3 species. Diatoms and dinoflagellates were the most diverse groups. Out of 142 species of diatoms 28% was contributed by three genera: Chaetoceros (18 species), Navicula (12 species), Rhizosolenia (11 species). Of the 129 species of dinoflagellates, 30% were represented by two genera: Ceratium (21 species) and Protoperidinium (17 species). As a whole, a pronounced prevalence of diatoms was typical for the phytoplankton community in the Arabian Sea during the period of analysis. On an average, diatoms contributed 52% to the total species diversity as shown in figure 2. Their prevalence was at a maximum (70.2%) during the winter-monsoon period, in February and March, and reduced to 65% during the Inter-monsoon period (April to May). Dinoflagellates contributed only 47% to the total species diversity, with 25.4% during the winter-monsoon period, and reduced to 19.24% during the Inter-monsoon period. The remaining 1% was contributed by other algae. 2 Species of Trichodesmium represented this category. Temporal variation in phytoplankton concentration at surface Total phytoplankton cells observed at the surface (0 m depth) of the Arabian Sea ranged from (11 cells/lit to 10440 cells/lit) during 2003 to 2007. This shows that the waters of Arabian Sea are highly productive. Their concentrations in winter monsoon and inter monsoon periods is summarized in table 2 and figures 3 and 4 illustrate their group wise concentration in the two periods. More than half (60%) of the total diatoms was contributed by Rhizosolenia alata (10779 cells/lit), Rhizosolenia shrubsolei (6552 cells/lit), Navicula sp. (4500 Fig II Richness of various Phytoplankton groupss Table II Phytoplankton cell concentration and percent contribution, in winter monsoon and inter monsoon periods and their group wise concentration and percent contribution in the two periods Total cell Concentration (116525 cells/lit) Winter monsoon period (60867 cells/lit) 52% Inter monsoon period (55658 cells/lit) 48% Diatoms (42568cells/lit) 70% Dinoflagellates (15407cells/lit) 25% Other Algae (2892cells/lit) 5% Diatoms (36170cells/lit) 65% Dinoflagellates (10714cells/lit) 19% Other Algae (8774cells/lit) 16%

1516 INDIAN J MAR SCI VOL 43 NO 8, AUGUST 2014 cells/lit) and Rhizosolenia hebatata (3710 cells/lit). Noctiluca scintillans (10440 cells/lit) alone contributed to 68% of the total dinoflagellates in the winter monsoon period. Trichodesmium erythraeum (6316 cells/lit) was the greatest contributor (72%) among other algae (Cyanophyceae), in the inter monsoon period. Among dinoflagellates observed in the inter monsoon period, 21% was formed by Scripsiella trachoidea (2220 cells/lit) and 14% was formed by Prorocentrum minimus (1462 cells/lit). Among diatoms that occurred in the inter monsoon period, Navicula sp. (5782 cells/lit); Thalassiothrix frauenfeldii (2504 cells/lit) and Rhizosolenia fragilissima (2220 cells/lit) were the major ones. Spatial variation in phytoplankton concentration at surface To study the spatial distribution of phytoplankton Arabian Sea was categorized as Coastal (< 50 m depth), shelf (50-200 m depth), slope (200-500 m depth) and Open Ocean (>500 m depth). Noctiluca scintillans formed massive blooms in the open ocean of northern Arabian Sea covering a large area from 17º19.40 N and 70º11.95 E to 20º28.72 N and 67º30.51 E during winter monsoon period as shown in figure 5. Whereas Trichodesmium erythraeum formed bloom in the coastal waters at 20º31.87 N and 70º34.77 E during inter monsoon period as shown in figure 6. Phytoplankton Concentration and its relation with light and chlorophyll Phytoplankton concentration in cell/lit was found to be directly proportional/linearly related to the chlorophyll concentration. Vertical distribution of phytoplankton cells was related to chlorophyll concentrations at surface, Deep Chlorophyll maxima (DCM) and 1% light level depth (euphotic depth). This correlation stands well in the inter monsoon period with r 2 = 0.968 at surface, r 2 = 0.895 at DCM and r 2 = 0.807 at 1% light level as show in figure 7. Phytoplankton concentration does not show a very good correlation with chlorophyll, in the winter monsoon period. It has r 2 = 0.195 at surface, r 2 = 0.546 at DCM and r 2 = 0.063 at 1% light level as shown in figure 8. This is because Noctiluca scintillans formed bloom at the surface in the winter monsoon period. Though the number of cells/lit was found to be high in this period but chlorophyll concentration remained comparatively low. Trichodesmium erythraeum also formed bloom in the inter monsoon period, but this bloom was restricted to the coastal waters only and coastal waters have a well mixed vertical profile with respect to phytoplankton cells, chlorophyll and light. So the bloom condition did not influence the correlation between cells and chlorophyll very significantly. Range of diversity indices The results of diversity calculation show that the coastal waters have more diversity than the open ocean waters. Also diversity of the open ocean is greatly affected by the Noctiluca scintillans bloom in Fig III and IV: Percent contribution of diatoms, dinoflagellates and other algae in winter and inter monsoon periods Fig V and VI: Spatial distribution of phytoplankton cells over the Arabian Sea during winter and inter monsoon periods. The colour bar shows cell concentration (cells/lit); green to red colour in the map shows the region covered by the bloom. Fig VII and VIII: Correlation between phytoplankton cell counts and Chlorophyll-a at surface, DCM and1% during winter and inter monsoon periods. The dotted lines are the linear fits applied to the data points.

MOTWANI et al.: COMPARISON BETWEEN PHYTOPLANKTON BIO-DIVERSITY AND VARIOUS INDICES 1517 Diversity index H d Mc J D Table III Range of diversity indices in all the cruises Minimum value 0.098 (FORV-222, station 10) 1.45 (FORV-253, station 4) 0.44 (FORV-222, station 10) 0.17 (FORV-253, station 4) 0.08 (FORV-253, station 4) 0.01 (FORV-222, station 7) Maximum value 3.28 (FORV-244, station 15) 20.91 (FORV-244 station 15) 4.64 (FORV-244 station 15) 0.78 (FORV-244 station 15) 0.43 (FORV-244, station 10) 0.69 (FORV-253 station 4) Fig 9 Correlation between phytoplankton cell counts and various Diversity indices Fig 10 Correlation between Dominance Index and Pielou Evenness Index the winter monsoon period. The minimum and maximum values of various diversity indices are tabulated in table 3. Comparison between various diversity indices The diversity indices were calculated using the cell counts obtained from the samples collected during the cruise. All the diversity indices represent the phytoplankton diversity over the Arabian Sea. To determine which index best represents the diversity, the diversity indices (Shannon index, Simpson s index, Margalef index and McIntosh index) were correlated with the phytoplankton cell concentration. The results of this correlation, as shown in figure 9, clearly indicate that the Shannon index has better correlation with the cell counts. Thus it best represents the phytoplankton diversity as compared to other indices. If one or two species dominate the phytoplankton community, the phytoplankton distribution will be uneven. Such an inverse relation between dominance and evenness was computed and their negative correlation has r 2 =0.73 and RMSE=0.003 as shown in figure 10. Analyses of these six indices show that none of these could formulate the richness component of diversity. In particular, the most preferred index among the other diversity indices, the Shannon index, cannot express the evenness component. Even other scaling diversity indices failed to formulate both the richness aspect and the evenness aspect of diversity. Acknowledgements This work is supported under the Meteorology and Oceanography (MOP) II program of Indian Space Research Organization (ISRO). In this valuable advice and encouragement of Dr. Ajai, Group Director, MPSG and Dr. Prakash Chauhan, Head of the Division, MPD of Space Applications Center (ISRO), Ahmedabad is greatfully acknowledged. References 1 Irigoien, X., Huisman, J. and Harris, R. P., Global biodiversity patterns of marine phytoplankton and zooplankton, Nature 429 (2004) 863-867. 2 Ptacnik, R., Solimini, A. G., Andersen, T., Tamminen, T., Brettum, P., Lepisto, L., Wille n E. and Rekolainen S. Diversity predicts stability and resource use efficiency in natural phytoplankton communities. Proceedings of the national academy of Sciences 105 (13) (2008) 5134-5138. 3 Qasim, S. Z., Biological productivity of the Indian Ocean. Indian Journal of Marine Sciences 6 (1977) 122-137. 4 http://www.museumstuff.com/learn/topics/list_of_rivers_of _India::sub::Rivers_Flowing_Into_Arabian_Sea. Accessed on 01/11/2010, 8:15 am.

1518 INDIAN J MAR SCI VOL 43 NO 8, AUGUST 2014 5 Tomas, C. R., Identified Marine phytoplankton. (Academic Press, New York), 1997, pp. 858. 6 Mandaville S M, Benthic Macro-invertebrates in Freshwater Taxa Tolerance Values, Metrics and Protocols, Project H - 1. (Nova Scotia: Soil & Water Conservation Society of Metro Halifax) 2002. 7 Simpson, E. H., Measurement of Diversity. Nature 163 (1949) 688. 8 Turkmen G & Kazanci N. Applications of various diversity indices to benthic macro-invertebrate assemblages in streams of a natural park in Turkey. Fourth International Scientific Conference: BALWOIS 2010 - Conference on Water Observation and Information System for decision Support organized in Ohrid, Republic of Macedonia 2010. 9 McIntosh, R. P. An index of diversity and the relation of certain concepts of diversity. Ecology 48 (3) (1967) 392 404. 10 Pielou, E.C., The Measurement of Diversity in Different Types of Biological Collections. Journal sof Theoretical Biology 13 (1966) 131-144.