Difference between revisions of "Measurements of biodiversity"

From Coastal Wiki
Jump to: navigation, search
Line 1: Line 1:
 
{{Review
 
|name=Job Dronkers|AuthorID=120|
 
}}
 
 
  
 
A variety of objective measures have been created in order to empirically measure biodiversity. The basic idea of a diversity index is to obtain a quantitative estimate of biological variability that can be used to compare biological entities, composed of diverse components, in space or in time. It is important to distinguish ‘richness’ from ‘diversity’. Diversity usually implies a measure of both species number and ‘equitability’ (or ‘[[evenness]]’). Three  types of indices can be distinguished:
 
A variety of objective measures have been created in order to empirically measure biodiversity. The basic idea of a diversity index is to obtain a quantitative estimate of biological variability that can be used to compare biological entities, composed of diverse components, in space or in time. It is important to distinguish ‘richness’ from ‘diversity’. Diversity usually implies a measure of both species number and ‘equitability’ (or ‘[[evenness]]’). Three  types of indices can be distinguished:
Line 43: Line 38:
 
   
 
   
 
====Species richness indices====
 
====Species richness indices====
Species richness <math>R</math> is the simplest measure of biodiversity and is simply a count of the number of different species in a given area. This measure is strongly dependent on sampling size and effort. Two species richness indices try to account for this problem:
+
Species richness <math>S</math> is the simplest measure of biodiversity and is simply a count of the number of different species in a given area. This measure is strongly dependent on sampling size and effort. Two species richness indices try to account for this problem:
  
 
*'''Margalef’s diversity index''':<ref> Clifford H.T. and Stephenson W. 1975. An introduction to numerical classification. London: Academic Express.</ref>   
 
*'''Margalef’s diversity index''':<ref> Clifford H.T. and Stephenson W. 1975. An introduction to numerical classification. London: Academic Express.</ref>   
  
<math>D_{Mg} = \Large\frac{R-1}{\ln n}</math>,
+
<math>D_{Mg} = \Large\frac{S-1}{\ln N}</math>,
  
where <math>n</math> = the total number of individuals in the sample and <math>R</math> = the number of species recorded.
+
where <math>N</math> = the total number of individuals in the sample and <math>S</math> = the number of species recorded.
  
 
*'''Menhinick’s diversity index''':<ref>Whittaker R.H. 1977. Evolution of species diversity in land communities. Evolutionary Biol. 10: 1-67 </ref>   
 
*'''Menhinick’s diversity index''':<ref>Whittaker R.H. 1977. Evolution of species diversity in land communities. Evolutionary Biol. 10: 1-67 </ref>   
  
<math>D_{Mn} = \Large\frac{R}{\sqrt{n}}</math>.
+
<math>D_{Mn} = \Large\frac{S}{\sqrt{N}}</math>.
  
 
    
 
    
Line 68: Line 63:
 
<math>\alpha\,x \, , \quad \Large\frac{\alpha\, x^2}{2} \, , \frac{\alpha\, x^3}{3} , \, ...\, , \frac{\alpha\, x^n}{n}</math> ,
 
<math>\alpha\,x \, , \quad \Large\frac{\alpha\, x^2}{2} \, , \frac{\alpha\, x^3}{3} , \, ...\, , \frac{\alpha\, x^n}{n}</math> ,
  
where <math>\alpha\,x</math> is the number of species to have one individual, <math>\Large\frac{\alpha\, x^2}{2}</math> those with two individuals, and so on. Since <math>0  < \,x < 1</math> and  <math>\alpha\,</math> and <math>\,x</math> are presumed to be constant, the expected number of species will be the highest in the first abundance class. The total number of species <math>R</math> for a large sample with many species  can therefore be approximated by 
+
where <math>\alpha\,x</math> is the number of species to have one individual, <math>\Large\frac{\alpha\, x^2}{2}</math> those with two individuals, and so on. Since <math>0  < \,x < 1</math> and  <math>\alpha\,</math> and <math>\,x</math> are presumed to be constant, the expected number of species will be the highest in the first abundance class. The value of <math>\,x</math> is calculated iteratively from:
  
<math>\approx  \alpha \sum_{k=1}^{\infty} x^k / k = - \alpha \, \ln(1-x)</math>.
+
<math>\Large\frac{S}{N} = \frac{1-x}{x}.\ln\frac{1}{1-x}</math>,
  
The total number of individual can be approximated by
+
and <math>\alpha\,</math> can be calculated from the equation:
 
 
<math>n \approx \alpha \sum_{k=1}^{\infty} x^k = \alpha \, x / (1-x)</math>.
 
 
 
The values of <math>\,x</math> and <math>\alpha</math> can be estimated from these two equations by counting the numbers <math>n</math> and <math>R</math> of the sample.
 
  
 +
<math>\alpha\, = \Large\frac{N(1-x)}{x}</math>.
  
 
====Non-parametric indices====
 
====Non-parametric indices====
The first two indices are based on information theory. These indices are based on the rationale that the diversity in a natural system can be measured in a similar way to the information contained in a code or message (see Appendix A1).
+
The first two indices are based on information theory. These indices are based on the rationale that the diversity in a natural system can be measured in a similar way to the information contained in a code or message.
 
 
  
 
* The most widely used diversity index in the ecological literature is the '''Shannon-Wiener diversity index'''.<ref>Shannon C. E. and Weaver W. 1949. The mathematical theory of communication. Urbana, IL: University of Illinois Press.</ref><ref name="Pielou">Pielou E.C. 1975. Ecological diversity. New York: Wiley Interscience.</ref>.  
 
* The most widely used diversity index in the ecological literature is the '''Shannon-Wiener diversity index'''.<ref>Shannon C. E. and Weaver W. 1949. The mathematical theory of communication. Urbana, IL: University of Illinois Press.</ref><ref name="Pielou">Pielou E.C. 1975. Ecological diversity. New York: Wiley Interscience.</ref>.  
  
It assumed that individuals are randomly sampled from a very large community, and that all species are represented in the sample. The Shannon index is calculated from the equation
+
It assumed that individuals are randomly sampled from an infinitely large community, and that all species are represented in the sample. The Shannon index is calculated from the equation:
  
<math>H' = -\sum_{i=1}^R p_i \, \ln p_i</math> ,  
+
<math>H' = -\sum_{i=1}^S p_i \ln p_i</math> ,  
 
 
where <math>p_i</math> is the proportion of individuals found in the <math>i</math>-th species (<math>\sum_{i=1}^R p_i = 1</math>).
 
  
 +
where <math>p_i</math> is the proportion of individuals found in the <math>i</math>-th species (<math>\sum_{i=1}^S p_i = 1</math>).
  
 
* Where the randomness cannot be guaranteed, for example when certain species are preferentially sampled, the '''Brillouin index''' <ref> Pielou E.C. 1969. An introduction to mathematical ecology. New York: Wiley.</ref><ref name="Pielou"/>  is the appropriated form of the information index. It is calculated as follows:
 
* Where the randomness cannot be guaranteed, for example when certain species are preferentially sampled, the '''Brillouin index''' <ref> Pielou E.C. 1969. An introduction to mathematical ecology. New York: Wiley.</ref><ref name="Pielou"/>  is the appropriated form of the information index. It is calculated as follows:
  
<math>H = \Large\frac{1}{n}\normalsize [ \ln(n!) - \sum_{i=1}^R \ln(n_i !)]  , </math>
+
<math>H = \Large\frac{1}{N}\normalsize.\ln \large\frac{N!}{ \prod_{i=1}^S N_i !} </math>.
 
 
in which <math> n_i != 1 \times 2 \times 3 \times ... \times n_i</math> and  <math>n_i </math>  = the number of individuals in species <math>i</math> and <math>n</math> is the total number of individuals in the community.
 
  
 +
in which <math> N_i != 1 \times 2 \times 3 \times ... \times N_i</math> and  <math>N_i </math>  = the number of individuals in species <math>i</math> and <math>N</math> is the total number of individuals in the community.
  
 
* One of the best known and earliest evenness measures is the '''Simpson ’s index'''<ref>Simpson E.H. 1949. Measurement of diversity. Nature 163: 688.</ref> which is given by:
 
* One of the best known and earliest evenness measures is the '''Simpson ’s index'''<ref>Simpson E.H. 1949. Measurement of diversity. Nature 163: 688.</ref> which is given by:
  
<math>\gamma\, = \sum_{i=1}^R p_i^2</math>.
+
<math>\gamma\, = \sum_{i=1}^S p_i^2</math>.
 
 
This index is used for large, sampled communities.  Simpson’s index expresses the probability that any two individuals drawn at random from an infinitely large community belong to the same species. If all species are equally represented in the sample <math>p_i=1/R</math> or <math>\gamma=1/R</math>.
 
 
 
 
 
* Another evenness index was proposed by Pielou (1966<ref>Pielou, E.C. 1966. The measurement of diversity in different types of biological collection. Journal of Theoretical Biology 13: 131–144</ref>). The '''Pielou index''' <math>J'</math> is defined as
 
 
 
<math>J' \, = \, H' / \ln(R) </math>.
 
 
 
If all species are represented in equal numbers in the sample <math>J'=1</math>. If one species strongly dominates <math>J'</math> is close to zero.
 
  
 +
This index is used for large, sampled communities.  Simpson’s index expresses the probability that any two individuals drawn at random from an infinitely large community belong to the same species.
  
 
* '''The Hill numbers'''<ref> Hill, M.O. 1973. Diversity and evenness: a unifying notation and its consequences. Ecology 54: 427–473</ref> show the relation between the species-richness indices and the evenness-indices.  Hill defined a set of diversity number of different order. The diversity number of order <math>a</math> is defined as:
 
* '''The Hill numbers'''<ref> Hill, M.O. 1973. Diversity and evenness: a unifying notation and its consequences. Ecology 54: 427–473</ref> show the relation between the species-richness indices and the evenness-indices.  Hill defined a set of diversity number of different order. The diversity number of order <math>a</math> is defined as:
  
<math>H_a = (\sum_{i=1}^R p_i^a)^{\large\frac{1}{1-a}}</math> ,
+
<math>H_a = (\sum_{i=1}^S p_i^a)^{\large\frac{1}{1-a}}</math> ,
  
 
where <math>p_i</math>  = the proportional abundance of species <math>i</math> in the sample and <math>a</math> = the order in which the index is dependent of rare species.  
 
where <math>p_i</math>  = the proportional abundance of species <math>i</math> in the sample and <math>a</math> = the order in which the index is dependent of rare species.  
Line 121: Line 102:
 
The most common Hill numbers are
 
The most common Hill numbers are
 
    
 
    
<math>H_0 = \, R </math>,
+
<math>H_0 = \ S </math>,
  
<math>H_1 = \exp{H'} </math> (the limit of <math>H_a</math> for <math>a \to 1</math> corresponds to the exponential of the Shannon-Wiener diversity index, see appendix A2),
+
<math>H_1 = \exp{H'} </math> (the limit of <math>H_a</math> for <math>a \to 1</math> corresponds to the exponential of the Shannon-Wiener diversity index, see appendix A),
 
   
 
   
 
<math>H_2 =\Large\frac{1}{\gamma}</math> (the reciprocal of Simpson’s <math>\gamma\, </math>) .
 
<math>H_2 =\Large\frac{1}{\gamma}</math> (the reciprocal of Simpson’s <math>\gamma\, </math>) .
 
If all species are represented in equal numbers in the sample <math>H_0=H_1=H_2=R</math>.
 
 
  
 
===Taxonomic indices===
 
===Taxonomic indices===
Line 138: Line 116:
 
• Taxonomic diversity (<math>\Delta</math>) reflects the average taxonomic distance between any two organisms, chosen at random from a sample. The distance can be seen as the length of the path connecting these two organisms through a phylogenetic tree or a Linnean classification<ref name=WC/>. This index includes aspects of taxonomic relatedness and evenness,
 
• Taxonomic diversity (<math>\Delta</math>) reflects the average taxonomic distance between any two organisms, chosen at random from a sample. The distance can be seen as the length of the path connecting these two organisms through a phylogenetic tree or a Linnean classification<ref name=WC/>. This index includes aspects of taxonomic relatedness and evenness,
  
<math>\Delta =  \Large\frac{ \sum\sum_{i<j} \, \omega_{ij} \, x_i x_j}{n(n-1)/2} </math>,
+
<math>\Delta =  \Large\frac{ \sum\sum_{i<j} \, \omega_{ij} \, x_i x_j}{N(N-1)/2} </math>,
  
where <math>x_i=n p_i</math> denotes the abundance of the <math>i</math>-th species in the sample (<math>\sum_{i=1}^R x_i =n</math> = the number of individuals in the sample) and <math>\omega_{ij}</math> is the 'distinctness weight' given to the path length (number of steps) linking species <math>i</math> and <math>j</math> in the hierarchical classification. (The path length for two individuals of the same species is zero).
+
where <math>x_i=N p_i</math> denotes the abundance of the <math>i</math>-th species in the sample (<math>\sum_{i=1}^S x_i =N</math> = the number of individuals in the sample) and <math>\omega_{ij}</math> is the 'distinctness weight' given to the path length (number of steps) linking species <math>i</math> and <math>j</math> in the hierarchical classification. (The path length for two individuals of the same species is zero).
  
 
• Taxonomic distinctness (<math>\Delta^*</math>)  is the average path length between two randomly chosen but taxonomically different organisms. This measure is measure of pure taxonomic relatedness,
 
• Taxonomic distinctness (<math>\Delta^*</math>)  is the average path length between two randomly chosen but taxonomically different organisms. This measure is measure of pure taxonomic relatedness,
Line 150: Line 128:
 
Journal of Applied Ecology 35 (4): 523-531</ref>,  
 
Journal of Applied Ecology 35 (4): 523-531</ref>,  
  
<math>\Delta^+= \Large\frac{\sum\sum_{i<j} \, \omega_{ij}}{R(R-1)/2}</math>.
+
<math>\Delta^+= \Large\frac{\sum\sum_{i<j} \, \omega_{ij}}{S(S-1)/2}</math>.
  
 
===Functional diversity===
 
===Functional diversity===
Line 159: Line 137:
 
==Species-Abundance distributions<ref>Heip, C.H.R., Herman, P.M.J. and Soetaert, K. 1998. Indices de diversité et régularité. [Indices of diversity and evenness]. Océanis (Doc. Océanogr.) 24(4): 67-87</ref>==
 
==Species-Abundance distributions<ref>Heip, C.H.R., Herman, P.M.J. and Soetaert, K. 1998. Indices de diversité et régularité. [Indices of diversity and evenness]. Océanis (Doc. Océanogr.) 24(4): 67-87</ref>==
  
Nearly all diversity and evenness indices are based on the relative abundance of species, thus on estimates of <math>p_i</math> in which  <math>p_i = n_i / n </math>,  with <math>n_i</math> the abundance of the <math>i</math>-th species in the sample and <math>n = \sum _{i = 1}^R \, n_i</math> with <math>R</math> the total number of species in the sample.
+
Nearly all diversity and evenness indices are based on the relative abundance of species, thus on estimates of <math>p_i</math> in which  <math>p_i = \large\frac{N_i}{N} </math>,  with <math>N_i</math> the abundance of the <math>i</math>-th species in the sample and <math>N = \sum _{i = 1}^S \, N_i</math> with <math>S</math> the total number of species in the sample.
 
   
 
   
 
When recording the abundance of different species in a sample, it is invariably found that some species are rare, whereas others are more abundant. This feature of ecological communities is found independently of the taxonomic group or the area investigated.  An important goal of ecology is to describe these consistent patterns in different communities, and explain them in terms of interactions with the biotic and abiotic environment.
 
When recording the abundance of different species in a sample, it is invariably found that some species are rare, whereas others are more abundant. This feature of ecological communities is found independently of the taxonomic group or the area investigated.  An important goal of ecology is to describe these consistent patterns in different communities, and explain them in terms of interactions with the biotic and abiotic environment.
Line 169: Line 147:
 
2. '''The <math>k</math>-dominance plot''' shows the cumulative percentage (the percentage of  the <math>k</math>-th most dominant plus all more dominant species) in relation to species (<math>k</math>) rank or log species (<math>k</math>) rank.  
 
2. '''The <math>k</math>-dominance plot''' shows the cumulative percentage (the percentage of  the <math>k</math>-th most dominant plus all more dominant species) in relation to species (<math>k</math>) rank or log species (<math>k</math>) rank.  
  
3. '''The Lorenzen curve''' is based on the <math>k</math>-dominance plot but the species rank <math>k</math> is transformed to <math> (k/R) \times 100</math> to facilitate comparison between communities with different numbers of species.
+
3. '''The Lorenzen curve''' is based on the <math>k</math>-dominance plot but the species rank <math>k</math> is transformed to <math> (k/S) \times 100</math> to facilitate comparison between communities with different numbers of species.
  
4. '''The collector’s curve''' addresses a different problem.  When the sampling effort is increased, and thus the number of animals <math>n</math> caught, new species will appear in the collection. A collector’s curve expresses the number of species as a function of the number of specimens caught.  As more specimens are caught, a collector’s curve can reach an asymptotic value. However, this hardly does occur in practice due to the vague boundaries of ecological communities: as the sampling effort increases, the number of different patches increases as well.
+
4. '''The collector’s curve''' addresses a different problem.  When the sampling effort is increased, and thus the number of animals <math>N</math> caught, new species will appear in the collection. A collector’s curve expresses the number of species as a function of the number of specimens caught.  As more specimens are caught, a collector’s curve can reach an asymptotic value. However, this hardly does occur in practice due to the vague boundaries of ecological communities: as the sampling effort increases, the number of different patches increases as well.
 
    
 
    
 
5. '''The species-abundance distribution''' plots the number of species that are represented by <math>1,2, … , n </math> individuals against the corresponding abundance.  This plot can only be drawn if the collection is large and contains many species.  The species are generally grouped in logarithmic density classes.
 
5. '''The species-abundance distribution''' plots the number of species that are represented by <math>1,2, … , n </math> individuals against the corresponding abundance.  This plot can only be drawn if the collection is large and contains many species.  The species are generally grouped in logarithmic density classes.
Line 192: Line 170:
 
It assumes that a species preempts a fraction <math>k</math> of a limiting resource, a second species the same fraction <math>k</math> of the remainder and so on. If the abundances are proportional to their share of the resource, the ranked abundances list is given by geometric series:
 
It assumes that a species preempts a fraction <math>k</math> of a limiting resource, a second species the same fraction <math>k</math> of the remainder and so on. If the abundances are proportional to their share of the resource, the ranked abundances list is given by geometric series:
  
<math>k, \, k(1-k), \, …, \, k(1-k)(R-2), \, k(1-k)(R-1) </math>
+
<math>k, \, k(1-k), \, …, \, k(1-k)(S-2), \, k(1-k)(S-1) </math>
  
where <math>R</math> is the number of the species in the community.  
+
where <math>S</math> is the number of the species in the community.  
  
 
The geometric model gives a straight line on a plot of log abundance against rank (species sequence).  It is not very often found in nature, only in early successional stages or in species poor environments.
 
The geometric model gives a straight line on a plot of log abundance against rank (species sequence).  It is not very often found in nature, only in early successional stages or in species poor environments.
Line 203: Line 181:
 
'''2. The broken stick model or the negative exponential distribution.'''
 
'''2. The broken stick model or the negative exponential distribution.'''
 
   
 
   
In this model a limiting resource is compared with a stick, broken in <math>R</math> parts at <math>R-1</math> randomly located points.  The length of the parts is taken as representative for the density of the <math>R</math> species subdividing the limiting resource.  If the species are ranked according to abundance, the expected abundance <math>n_i</math> of species <math>i</math> is given by:
+
In this model a limiting resource is compared with a stick, broken in <math>S</math> parts at <math>S-1</math> randomly located points.  The length of the parts is taken as representative for the density of the <math>S</math> species subdividing the limiting resource.  If the species are ranked according to abundance, the expected abundance of species <math>i</math>, <math>N_i</math> is given by:
 
   
 
   
<math> E(n_i) = \Large\frac{1}{R} \sum_{k = i}^R \frac{1}{k}  </math>.
+
<math> E(N_i) = \Large\frac{1}{S} \sum_{k = i}^S \frac{1}{k}  </math>.
  
 
The negative exponential distribution is not often found in nature.  It describes a too even distribution of individuals over species to be a good representation of natural communities.
 
The negative exponential distribution is not often found in nature.  It describes a too even distribution of individuals over species to be a good representation of natural communities.
Line 218: Line 196:
 
through the 'octaves method' (grouping of species into octaves, such that the center of each group is double that of the previous group, by analogy with musical octaves). The distribution is traditionally written in the form:
 
through the 'octaves method' (grouping of species into octaves, such that the center of each group is double that of the previous group, by analogy with musical octaves). The distribution is traditionally written in the form:
 
   
 
   
<math> R(r) = \ R_0 \exp\large(-a^2 r^2)  </math>
+
<math> S(R) = \ S_0 \exp\large(-a^2 R^2)  </math>
  
With <math>R(r)</math> = the number of the species in the <math>r</math>-the octave to the right, and to the left of the symmetric curve; <math>R_0</math> = the number of the species in the modal octave; and  
+
With <math>S(R)</math> = the number of the species in the <math>R</math>-the octave to the right, and to the left of the symmetric curve; <math>S_0</math> = the number of the species in the modal octave; and  
  
 
<math>a = 2 (\sigma^2)^{-1/2}</math> = the inverse width of the distribution.
 
<math>a = 2 (\sigma^2)^{-1/2}</math> = the inverse width of the distribution.
  
  
 +
==Appendix A==
 +
Take <math>a=1+\epsilon</math> and <math>\epsilon \to 0</math>, then <math>p_i^a \to p_i+\epsilon \large\frac{d}{d \epsilon} p_i^{1+\epsilon} \approx p_i + \epsilon p_i \ln p_i</math>.
 +
Because <math>\sum_{i=1}^S p_i = 1</math>, we have
  
==Appendix A1==
+
<math>H_1 = \lim_{\epsilon \to 0} (1 + \epsilon \sum_{i=1}^S p_i \ln p_i)^{-1/\epsilon} = \exp(-\sum_{i=1}^S p_i \ln p_i)</math>.
The Shannon-Wiener diversity index is a measure of the amount of information (inverse of entropy) contained in a sample, where information is defined as the minimum length of a string of digits necessary to describe the sample. The length of a string of digits to describe a number is proportional to the logarithm of this number. With this definition, a measure of the amount of information in a large sample is given by the logarithm of the number <math>N</math> of all different distinct ways in which the sample can be re-ordered by permutation of individuals. If all <math>n</math> individuals in a sample belong to different species the number of distinct permutations is highest and equal to <math>n!</math> (the first individual in a ordered sample can be chosen in <math>n</math> ways, the second in <math>n-1</math> ways and so on). If there are <math>n_1</math> individuals of species 1, <math>n_2</math> individuals of species 2, and so on, then the number of distinct permutations is  <math>N = n!/(n_1! \times n_2! \times . . . .  \times n_R!)</math>. Assuming that the numbers <math>n, n_1, n_2, ….</math> are very large, one can approximate <math>\ln(n!) \approx n\, \ln(n) </math>. Taking the logarithm then gives  <math>\ln(N) \approx n \, \ln(n) – n_1 \, \ln(n_1) – n_2\, \ln(n_2) - ….. n_R \, \ln(n_R) </math>, which is <math>n \times</math> the Brillouin index <math>H</math>. For very large representative samples the probability of occurrence of species <math>i</math> is given by <math>p_i=n_i/n</math>. We further have <math>\sum_{i=1}^{R} p_i=1</math>.  Substitution gives <math>\ln(N) \approx -n \sum_{i=1}^{R}  p_i \, \ln(p_i)</math>, which is <math>n \times</math> the Shannon-Wiener index <math>H'</math>. Division by <math>n</math> makes the index independent of the sample size.
 
 
 
==Appendix A2==
 
Take <math>a=1+\epsilon</math> and <math>\epsilon \to 0</math>, then <math>p_i^a \to p_i+\epsilon \large\frac{d}{d \epsilon} p_i^{1+\epsilon}\normalsize \approx p_i + \epsilon p_i \ln p_i</math>.
 
Because <math>\sum_{i=1}^R p_i = 1</math>, we have
 
 
 
<math>H_1 = \lim_{\epsilon \to 0} (1 + \epsilon \sum_{i=1}^R p_i \ln p_i)^{-1/\epsilon} = \exp(-\sum_{i=1}^R p_i \ln p_i)</math>.
 
 
 
  
  

Revision as of 22:09, 26 September 2020

A variety of objective measures have been created in order to empirically measure biodiversity. The basic idea of a diversity index is to obtain a quantitative estimate of biological variability that can be used to compare biological entities, composed of diverse components, in space or in time. It is important to distinguish ‘richness’ from ‘diversity’. Diversity usually implies a measure of both species number and ‘equitability’ (or ‘evenness’). Three types of indices can be distinguished:

1. Species richness indices: Species richness is a measure for the total number of the species in a community (examples Fig. 1). However, complete inventories of all species present at a certain location, is an almost unattainable goal in practical applications.

2. Evenness indices: Evenness expresses how evenly the individuals in a community are distributed among the different species (examples Fig.2).


Fig. 1. Schematic visualization of the species richness: with respectively 5 and 10 species.
Fig. 2. Schematic visualization of differences in evenness for 5 species.


3. Taxonomic indices: These indices take into account the taxonomic relation between different organisms in a community. Taxonomic diversity, for example, reflects the average taxonomic distance between any two organisms, chosen at random from a sample. The distance can be seen as the length of the path connecting these two organisms along the branches of a phylogenetic tree.

These three types of indices can be used on different spatial scales [1]:

  • Alpha diversity refers to diversity within a particular area, community or ecosystem, and is usually measured by counting the number of taxa within the ecosystem (usually species level);
  • Beta diversity is species diversity between ecosystems; this involves comparing the number of taxa that are unique to each of the ecosystems. For example, the diversity of mangroves versus the diversity of seagrass beds;
  • Gamma diversity is a measure of the overall diversity for different ecosystems within a region. For example, the diversity within the coastal region of Gazi Bay in Kenia.


Diversity measurement is based on three assumptions[2]

1. All species are equal: this means that richness measurement makes no distinctions amongst species and threat the species that are exceptionally abundant in the same way as those that are extremely rare species. The relative abundance of species in an assemblage is the only factor that determines its importance in a diversity measure.

2. All individuals are equal: this means that there is no distinction between the largest and the smallest individual, in practice however the smallest animals can often escape for example by sampling with nets.

Taxonomic and functional diversity measures, however, do not necessarily treat all species and individuals as equal.

3. Species abundance has been recorded in using appropriate and comparable units. It is clearly unwise to use different types of abundance measure, such as the number of individuals and the biomass, in the same investigation. Diversity estimates based on different units are not directly comparable.


Diversity measures

Species richness indices

Species richness [math]S[/math] is the simplest measure of biodiversity and is simply a count of the number of different species in a given area. This measure is strongly dependent on sampling size and effort. Two species richness indices try to account for this problem:

  • Margalef’s diversity index:[3]

[math]D_{Mg} = \Large\frac{S-1}{\ln N}[/math],

where [math]N[/math] = the total number of individuals in the sample and [math]S[/math] = the number of species recorded.

  • Menhinick’s diversity index:[4]

[math]D_{Mn} = \Large\frac{S}{\sqrt{N}}[/math].


Despite the attempt to correct for sample size, both measures remain strongly influenced by sampling effort. Nonetheless they are intuitively meaningful indices and can play a useful role in investigations of biological diversity.

Heterogeneity measures

Heterogeneity measures are those that combine the richness and the evenness component of diversity. Heterogeneity measures fall into two categories: parametric indices, which are based on a parameter of a species abundance model, and nonparametric indices, that make no assumptions about the underlying distributions of species abundances.

Parametric indices

  • The log series index [math]\alpha\,[/math] [5] (see also log-series distributions) is a parameter of the log series model. The parameter [math]\alpha\,[/math] is independent of sample size. It describes the way in which the individuals are divided among the species, which is a measure of diversity. The attractive properties of this diversity index are: it provides a good discrimination between sites, it is not very sensitive to density fluctuations and it is normally distributed such that confidence limits can be attached to [math]\alpha\,[/math].

The log series takes the form:

[math]\alpha\,x \, , \quad \Large\frac{\alpha\, x^2}{2} \, , \frac{\alpha\, x^3}{3} , \, ...\, , \frac{\alpha\, x^n}{n}[/math] ,

where [math]\alpha\,x[/math] is the number of species to have one individual, [math]\Large\frac{\alpha\, x^2}{2}[/math] those with two individuals, and so on. Since [math]0 \lt \,x \lt 1[/math] and [math]\alpha\,[/math] and [math]\,x[/math] are presumed to be constant, the expected number of species will be the highest in the first abundance class. The value of [math]\,x[/math] is calculated iteratively from:

[math]\Large\frac{S}{N} = \frac{1-x}{x}.\ln\frac{1}{1-x}[/math],

and [math]\alpha\,[/math] can be calculated from the equation:

[math]\alpha\, = \Large\frac{N(1-x)}{x}[/math].

Non-parametric indices

The first two indices are based on information theory. These indices are based on the rationale that the diversity in a natural system can be measured in a similar way to the information contained in a code or message.

  • The most widely used diversity index in the ecological literature is the Shannon-Wiener diversity index.[6][7].

It assumed that individuals are randomly sampled from an infinitely large community, and that all species are represented in the sample. The Shannon index is calculated from the equation:

[math]H' = -\sum_{i=1}^S p_i \ln p_i[/math] ,

where [math]p_i[/math] is the proportion of individuals found in the [math]i[/math]-th species ([math]\sum_{i=1}^S p_i = 1[/math]).

  • Where the randomness cannot be guaranteed, for example when certain species are preferentially sampled, the Brillouin index [8][7] is the appropriated form of the information index. It is calculated as follows:

[math]H = \Large\frac{1}{N}\normalsize.\ln \large\frac{N!}{ \prod_{i=1}^S N_i !} [/math].

in which [math] N_i != 1 \times 2 \times 3 \times ... \times N_i[/math] and [math]N_i [/math] = the number of individuals in species [math]i[/math] and [math]N[/math] is the total number of individuals in the community.

  • One of the best known and earliest evenness measures is the Simpson ’s index[9] which is given by:

[math]\gamma\, = \sum_{i=1}^S p_i^2[/math].

This index is used for large, sampled communities. Simpson’s index expresses the probability that any two individuals drawn at random from an infinitely large community belong to the same species.

  • The Hill numbers[10] show the relation between the species-richness indices and the evenness-indices. Hill defined a set of diversity number of different order. The diversity number of order [math]a[/math] is defined as:

[math]H_a = (\sum_{i=1}^S p_i^a)^{\large\frac{1}{1-a}}[/math] ,

where [math]p_i[/math] = the proportional abundance of species [math]i[/math] in the sample and [math]a[/math] = the order in which the index is dependent of rare species.

The most common Hill numbers are

[math]H_0 = \ S [/math],

[math]H_1 = \exp{H'} [/math] (the limit of [math]H_a[/math] for [math]a \to 1[/math] corresponds to the exponential of the Shannon-Wiener diversity index, see appendix A),

[math]H_2 =\Large\frac{1}{\gamma}[/math] (the reciprocal of Simpson’s [math]\gamma\, [/math]) .

Taxonomic indices

If two data-sets have identical numbers of species and equivalent patterns of species abundance, but differ in the diversity of taxa to which the species belong, it seems intuitively appropriate that the most taxonomically varied data-set is the more diverse. As long as the phylogeny of the data-set of interest is reasonably well resolved, measures of taxonomic diversity are possible.

  • Clarke and Warwick’s taxonomic distinctness index[11] which describes the average taxonomic distance – simply the “path length” between two randomly chosen organisms through the phylogeny of all the species in a data-set – has different forms: taxonomic diversity and taxonomic distinctness.

• Taxonomic diversity ([math]\Delta[/math]) reflects the average taxonomic distance between any two organisms, chosen at random from a sample. The distance can be seen as the length of the path connecting these two organisms through a phylogenetic tree or a Linnean classification[11]. This index includes aspects of taxonomic relatedness and evenness,

[math]\Delta = \Large\frac{ \sum\sum_{i\lt j} \, \omega_{ij} \, x_i x_j}{N(N-1)/2} [/math],

where [math]x_i=N p_i[/math] denotes the abundance of the [math]i[/math]-th species in the sample ([math]\sum_{i=1}^S x_i =N[/math] = the number of individuals in the sample) and [math]\omega_{ij}[/math] is the 'distinctness weight' given to the path length (number of steps) linking species [math]i[/math] and [math]j[/math] in the hierarchical classification. (The path length for two individuals of the same species is zero).

• Taxonomic distinctness ([math]\Delta^*[/math]) is the average path length between two randomly chosen but taxonomically different organisms. This measure is measure of pure taxonomic relatedness,

[math]\Delta^* = \Large\frac{ \sum\sum_{i\lt j} \, \omega_{ij} x_i x_j}{\sum\sum_{i\lt j} \, x_i x_j} [/math].

• When only presence/absence data is considered, both [math]\Delta[/math] and [math]\Delta^*[/math] converge to the same statistic [math]\Delta^+[/math], which can be seen as the average taxonomic path length between any two randomly chosen species[12],

[math]\Delta^+= \Large\frac{\sum\sum_{i\lt j} \, \omega_{ij}}{S(S-1)/2}[/math].

Functional diversity

The positive relationship between ecosystem functioning and species richness is often attributed to the greater number of functional groups found in richer assemblages. Petchey and Gaston [13] proposed a method for quantifying functional diversity. It is based on total branch length of a dendrogram, which is constructed from species trait values. One important consideration is that only those traits linked to the ecosystem process of interest are used. Thus a study focusing on bird-mediated seed dispersal would exclude traits such as plumage color that are not related to this function, but traits such as beak size and shape should be included. With standard clustering algorithms a dendrogram is then constructed. For example, a community with five species with different traits will have a higher functional diversity than a community of equal richness but where the species are functional similar.


Species-Abundance distributions[14]

Nearly all diversity and evenness indices are based on the relative abundance of species, thus on estimates of [math]p_i[/math] in which [math]p_i = \large\frac{N_i}{N} [/math], with [math]N_i[/math] the abundance of the [math]i[/math]-th species in the sample and [math]N = \sum _{i = 1}^S \, N_i[/math] with [math]S[/math] the total number of species in the sample.

When recording the abundance of different species in a sample, it is invariably found that some species are rare, whereas others are more abundant. This feature of ecological communities is found independently of the taxonomic group or the area investigated. An important goal of ecology is to describe these consistent patterns in different communities, and explain them in terms of interactions with the biotic and abiotic environment.

Different investigators have visualized the species-abundance distribution in different ways (Fig. 3).

1. The rank/abundance plot is one of the best known and most informative method. In this plot species are ranked in sequence from most to least abundant along the horizontal (or [math]x[/math]) axis. Their abundances are typically displayed in a log10 format on the [math]y[/math] axis, so that species whose abundances span several orders of magnitude can be easily accommodated on the same graph. In addition, proportional and or percentage abundances are often used.

2. The [math]k[/math]-dominance plot shows the cumulative percentage (the percentage of the [math]k[/math]-th most dominant plus all more dominant species) in relation to species ([math]k[/math]) rank or log species ([math]k[/math]) rank.

3. The Lorenzen curve is based on the [math]k[/math]-dominance plot but the species rank [math]k[/math] is transformed to [math] (k/S) \times 100[/math] to facilitate comparison between communities with different numbers of species.

4. The collector’s curve addresses a different problem. When the sampling effort is increased, and thus the number of animals [math]N[/math] caught, new species will appear in the collection. A collector’s curve expresses the number of species as a function of the number of specimens caught. As more specimens are caught, a collector’s curve can reach an asymptotic value. However, this hardly does occur in practice due to the vague boundaries of ecological communities: as the sampling effort increases, the number of different patches increases as well.

5. The species-abundance distribution plots the number of species that are represented by [math]1,2, … , n [/math] individuals against the corresponding abundance. This plot can only be drawn if the collection is large and contains many species. The species are generally grouped in logarithmic density classes.


Species-Abundance models[2]

A diverse range of models has also been developed to describe species abundance data (Fig. 4). The fitting of a model to field data is meaningful if the parameter estimates are to be used in further analysis.

1. The niche preemption model or the geometric model.

It assumes that a species preempts a fraction [math]k[/math] of a limiting resource, a second species the same fraction [math]k[/math] of the remainder and so on. If the abundances are proportional to their share of the resource, the ranked abundances list is given by geometric series:

[math]k, \, k(1-k), \, …, \, k(1-k)(S-2), \, k(1-k)(S-1) [/math]

where [math]S[/math] is the number of the species in the community.

The geometric model gives a straight line on a plot of log abundance against rank (species sequence). It is not very often found in nature, only in early successional stages or in species poor environments.


Fig. 4. The four main species- abundance models

2. The broken stick model or the negative exponential distribution.

In this model a limiting resource is compared with a stick, broken in [math]S[/math] parts at [math]S-1[/math] randomly located points. The length of the parts is taken as representative for the density of the [math]S[/math] species subdividing the limiting resource. If the species are ranked according to abundance, the expected abundance of species [math]i[/math], [math]N_i[/math] is given by:

[math] E(N_i) = \Large\frac{1}{S} \sum_{k = i}^S \frac{1}{k} [/math].

The negative exponential distribution is not often found in nature. It describes a too even distribution of individuals over species to be a good representation of natural communities.

3. The log-series distribution.

Fisher’s logarithmic series model [5] (see also the log series index [math]\alpha\,[/math]) describes the relationship between the number of species and the number of individuals in those species.

4. The log-normal distribution.

Preston [15] first suggested to use a log-normal distribution for the description of species-abundances distributions. A lognormal distribution was fitted to study the species abundance distribution through the 'octaves method' (grouping of species into octaves, such that the center of each group is double that of the previous group, by analogy with musical octaves). The distribution is traditionally written in the form:

[math] S(R) = \ S_0 \exp\large(-a^2 R^2) [/math]

With [math]S(R)[/math] = the number of the species in the [math]R[/math]-the octave to the right, and to the left of the symmetric curve; [math]S_0[/math] = the number of the species in the modal octave; and

[math]a = 2 (\sigma^2)^{-1/2}[/math] = the inverse width of the distribution.


Appendix A

Take [math]a=1+\epsilon[/math] and [math]\epsilon \to 0[/math], then [math]p_i^a \to p_i+\epsilon \large\frac{d}{d \epsilon} p_i^{1+\epsilon} \approx p_i + \epsilon p_i \ln p_i[/math]. Because [math]\sum_{i=1}^S p_i = 1[/math], we have

[math]H_1 = \lim_{\epsilon \to 0} (1 + \epsilon \sum_{i=1}^S p_i \ln p_i)^{-1/\epsilon} = \exp(-\sum_{i=1}^S p_i \ln p_i)[/math].


See also

Marine Biodiversity
Biodiversity and Ecosystem function
Wikipedia article Diversity index


Further reading

Magurran, A. E. 2004. Measuring biological diversity, Blackwell Publishing: Oxford, UK. 256 p


References

  1. [\http://en.wikipedia.org/wiki/Biodiversity#Measurement_of_biodiversity Wikipedia article Measurement of biodiversity
  2. 2.0 2.1 Magurran, A. E. 2004. Measuring biological diversity, Blackwell Publishing: Oxford, UK. 256 p
  3. Clifford H.T. and Stephenson W. 1975. An introduction to numerical classification. London: Academic Express.
  4. Whittaker R.H. 1977. Evolution of species diversity in land communities. Evolutionary Biol. 10: 1-67
  5. 5.0 5.1 Fisher, R. A., Corbet, A. S. and Williams, C. B. 1943. The relation between the number of species and the number of individuals in a random sample of an animal population. Journal of Animal Ecology 12: 42-58
  6. Shannon C. E. and Weaver W. 1949. The mathematical theory of communication. Urbana, IL: University of Illinois Press.
  7. 7.0 7.1 Pielou E.C. 1975. Ecological diversity. New York: Wiley Interscience.
  8. Pielou E.C. 1969. An introduction to mathematical ecology. New York: Wiley.
  9. Simpson E.H. 1949. Measurement of diversity. Nature 163: 688.
  10. Hill, M.O. 1973. Diversity and evenness: a unifying notation and its consequences. Ecology 54: 427–473
  11. 11.0 11.1 Warwick R.M. and Clarke K.R. 2001. Practical measures of marine biodiversity based on relateness of species. Oceanogr. Mar. Biol. Ann. Rev. 39: 207-231
  12. Clarke K.R. and Warwick R.M. 1998. A taxonomic distinctness index and its statistical properties Journal of Applied Ecology 35 (4): 523-531
  13. Petchey O.L. and Gaston K.J. 2002. Functional diversity (FD), species richness and community composition. Ecology letters Vol. 5 (3): 402-411
  14. Heip, C.H.R., Herman, P.M.J. and Soetaert, K. 1998. Indices de diversité et régularité. [Indices of diversity and evenness]. Océanis (Doc. Océanogr.) 24(4): 67-87
  15. Preston, F.W. 1948. The commonness and rarity of species. Ecology 29: 254–283
The main author of this article is Sohier, Charlotte
Please note that others may also have edited the contents of this article.

Citation: Sohier, Charlotte (2020): Measurements of biodiversity. Available from http://www.coastalwiki.org/wiki/Measurements_of_biodiversity [accessed on 29-03-2024]