QUANTIFYING FOREST STAND DIVERSITY USING STAND STRUCTURE by Christina Lynn Staudhammer B.Sc, University of California (Davis), 1990 A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE in THE FACULTY OF GRADUATE STUDIES DEPARTMENT OF FORESTRY We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA December 1999 © Christina Lynn Staudhammer, 1999 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department The University of British Columbia Vancouver, Canada Date DE-6 (2/88) Abstract Stand structure is an important component of the overall description of biological diversity in a stand. The diversity of tree sizes within a stand affects tree growth and yield, and is highly correlated with the biodiversity of a stand. Four measures of stand diversity based on the stand distribution of basal area by height, diameter, and species were proposed, assuming a baseline of maximum diversity corresponding to a uniform distribution. First, the Extended Shannon Index (ESI), a measure based on Shannon's Index, was derived. Second, a measure based on fitting the univariate and bivariate distributions of diameter and height was investigated. The third measure, STRI, was derived as a modified R-squared, based on the fit of the distribution to that of a uniform distribution. The fourth measure, STVI, was derived based on comparing the variance of a distribution to that of a uniform distribution. The four measures were evaluated with simulated data and inventory data from the Malcolm Knapp Research Forest (MKRF). The ESI and STVI resulted in logical orderings of both the simulated data and the MKRF data. The advantages of the ESI include its known sampling distribution and accepted use in forestry; however, the ESI depends on allocating data to arbitrary classes. The STVI does not require classifying data; however, its sampling distribution is unknown. ii Since the measures derived require only normally collected inventory plot data, they are relatively inexpensive to use. However, a thorough investigation of the properties of each should be undertaken. Since evaluating structural diversity is becoming an increasingly important part of forest assessment, a quantitative measure is needed to measure diversity. These measures provide a starting point toward finding an inexpensive, practical measure of structural diversity which gives reasonable results. iii TABLE OF CONTENTS ABSTRACT ii TABLE OF CONTENTS iv LIST OF TABLES vi LIST OF FIGURES vii ACKNOWLEDGEMENT viii CHAPTER 1: INTRODUCTION 1 CHAPTER 2: COMMON DIVERSITY MEASURES 5 2.1 SHANNON'S INDEX 6 2.2 SIMPSON'S INDEX 7 2.3 T H E Q STATISTIC 8 2.4 RELATIONSHIPS A M O N G INDICES 9 2.5 USING INDICES TO MAKE INFERENCES ABOUT POPULATIONS 11 2.6 DISCUSSION OF COMMON DIVERSITY MEASURES IN U S E 12 CHAPTER 3: STRUCTURAL DIVERSITY MEASURES 16 3.1 MEASURES OF VERTICAL DIVERSITY 18 3.2 MEASURES OF HORIZONTAL DIVERSITY 20 3.3 OTHER STRUCTURAL DIVERSITY MEASUREMENTS 22 3.4 DISCUSSION OF MEASURES 24 CHAPTER 4: POSSIBLE MEASURES OF STAND DIVERSITY USING STAND STRUCTURE 26 4.1 DESIRABLE CHARACTERISTICS OF A MEASURES OF STAND DIVERSITY USING STAND STRUCTURE 27 4.2 POSSIBLE ALTERNATIVE MEASURES OF STAND DIVERSITY USING STAND STRUCTURE 29 4.2.1 Shannon's Index Extended to Diameter, Height, and Species 29 4.2.1.1 Post-hoc Method 30 4.2.1.2 Combination Method 31 4.2.1.3 Properties of the Extended Shannon's Index 32 4.2.2 Structural Index based on Curve Fitting 33 4.2.2.1 Univariate Curve Fitting 35 4.2.2.2 Bivariate Curve Fitting 36 4.2.2.3 Properties of a Structural Index Based on Curve Fitting 37 4.2.3 Structural Index Based on R2 (STRI) 38 4.2.3.1 STRI Univariate Method 39 4.2.3.2 STRI Bivariate Method 41 iv 4.2.3.3 Properties of STRI 42 4.2.4 Structural Index based on Variance Comparisons (STVI) 43 4.2.4.1 STVI Univariate Method 44 4.2.4.2 STVI Bivariate Method '. 49 4.2.4.3 Properties of STVI 52 CHAPTER 5: EVALUATION OF ALTERNATIVE MEASURES OF STRUCTURAL DIVERSITY 54 5.1 SIMULATED TREE DATA 54 5.1.1 Methods '. 54 5.1.1.1 Baseline Maximal Diversity 54 5.1.1.2 Simulated Stands of Varying Diversity 55 5.1.2 Results and Discussion 62 5.1.2.1 Shannon's Index Extended to Diameter, Height, and Species 62 5.1.2.2 Structural Index based on Curve Fitting 65 5.1.2.3 Structural Index based on R 2 (STRI) 71 5.1.2.4 Structural Index based on Variance Comparisons (STVI) 72 5.1.3 Overall Discussion of Proposed Indices Computed with Simulated Data 76 5.2 MALCOLM KNAPP RESEARCH FOREST DATA 76 5.2.1 Methods 77 5.2.2 Results and Discussion 80 5.2.2.1 Shannon's Index Extended to Diameter, Height, and Species 80 5.2.2.2 Structural Index based on Variance Comparisons (STVI) 84 5.2.3 Discussion of Extended Shannon Index and Univariate STVI Computed with MKRF Data 88 CHAPTER 6: CONCLUSION : 89 LITERATURE CITED 95 APPENDIX I 101 v List of Tables Table 2.1. Summary of performance and characteristics of selected diversity measures. 14 Table 5.1. Plot type abbreviations for simulated cases 60 Table 5.2. Extended Shannon indices computed using simulated stand data 64 Table 5.3. Resulting coefficients and RMSE for the univariate Weibull fitted to dbh using simulated stands 66 Table 5.4. Resulting coefficients and RMSE for the univariate Weibull fitted to height using simulated stands 67 Table 5.5. RMSE and resulting values of coefficient as for the bivariate Weibull fitted to simulated stands 69 Table 5.6. Resulting coefficients bi and bj for the bivariate Weibull fitted to simulated stands 69 Table 5.7. Resulting coefficients c/ and C2 for the bivariate Weibull fitted to simulated stands 70 Table 5.8. Univariate (STRId+h) and bivariate (STRIdh) STRI using simulated stands. 71 Table 5.9. Univariate (STVIa+h) and bivariate (STVIdh) STVI using simulated stands. 73 Table 5.10. Average sample variances for simulated cases compared to the variance of a theoretical uniform distribution 75 Table 5.11. Data for height-dbh equation development by species 78 Table 5.12. Descriptive statistics of selected MKRF plots 80 Table 5.13. Extended Shannon indices computed with MKRF data 82 Table 5.14. Univariate (STVI d b h, STVIh eight, and STVI d + h) and bivariate (STVId h) STVI using MKRF data 85 Table 5.15. Sample variances and average sample variances by species for chosen MKRF stands compared to the variance of a theoretical uniform distribution... 87 v i List of Figures Figure 4.1. STRI and three equal dbh ranges 44 Figure 4.2. Comparison of STVI values with m=l and varying values of pi and p2.47 Figure 4.3. Comparison of STVI values with pt=p2=2 and varying values of m 48 Figure 5.1. Possible combinations of dbh versus height 57 Figure 5.2. Graphical description of simulated dbh and height values 59 vii Acknowledgement I would like to extend thanks to the many people that helped me to successfully complete this thesis, especially my supervisor Dr. Valerie LeMay, and my thesis committee, Drs. Tony Kozak, Guillaume Therien, Rick Routledge (Department of Statistics, Simon Fraser University), and Peter Marshall. Each and everyone one of them contributed significantly to my education and my experience at U.B.C., and I will not forget their generosity. I would also like to thank Dr. LeMay for helping me to find funding for this project and securing the data necessary to complete this project. Last, I would like to thank my family and friends for their amazing support during this time in my life, and extend an especially warm thanks to William Kovacs, whose love of the forest was my first inspiration. Thanks for letting me grow. viii Chapter 1 Introduction Over the past two decades, the focus of the sustainable forestry debate in North America has been increasingly drawn toward the concept of diversity. The maintenance of forest diversity has been linked to robust bird and insect populations, and is seen as an indicator of the wellbeing of a biological system (Magurran 1988). Interest in forest diversity has prompted both the Canadian Council of Forest Ministers (Canadian Standards Association 1996) and the United States federal government (U.S. federal law 36 Code of Federal Regulations Sec. 217.27) to require forest managers to conserve biological diversity at natural levels. The term 'biodiversity' has gained wide acceptance in describing biological diversity, encompassing the total variety of a biological system, including species diversity of plants, animals, fungi and bacteria, as well as the genetic and individual variation that exists within those species (Dobson 1996). Biodiversity may be evaluated at many levels: over an entire ecosystem, in a forest stand, or within a species at the genetic level. At the stand level, biodiversity has been defined to include not only species diversity, but also the variation of size, age, and genetic composition of the stand (Lahde et al. 1999). Thus, the definition of stand diversity may include not only the composition of species, but also the size distribution of individuals of those species. 1 Along with governmental and non-govemmental directives to conserve biological diversity comes the necessity to measure biological diversity. For example, U.S. federal law (36 C.F.R. Sec. 219.26) requires that inventories include "quantitative data making possible the evaluation of diversity". The evaluation of diversity, however, has long been a subject of heated debate. Many measures of diversity, in the form of indices, are in use in forestry and ecology; however, none are without critics (Magurran 1988). It has even been asserted that since there is no precise accepted definition of biological diversity, there can be no effective quantitative measure (Hurlbert 1971, Silbaugh and Betters 1995). Several traditional biodiversity indices, such as Simpson's index (Simpson 1949) and Shannon's index (Shannon and Weaver 1949) have gained wide acceptance in forestry as quantitative measures of species diversity (Swindel et al. 1984, Magurran 1988, Lewis et al. 1988, McMinn 1992, Silbaugh and Betters 1995). These indices rely on the distribution of individuals by species to arrive at a quantitative value of biological diversity. However, because of the variation that exists in tree size, traditional diversity indices are not entirely suited to the measurement of forest stand diversity (Lahde et al. 1999). Furthermore, constraints of time and money limit the forest manager to measuring only a subset of species represented in a single sample plot of forest. Surrogate measures for ecological composition are unavoidable, as the number of existing species in forested ecosystems is vast, and there are inherent difficulties in measuring them. To help remedy this, the use of indicator species has been proposed 2 (Noss 1990, Schneider 1997). Desirable criteria for species used as biodiversity indicators include (Noss 1990): • Wide distribution across regions • Potential to act as surrogate for other species » Inexpensive and easy to monitor • Pre-existing data available • Biology of species reasonably known • Sensitive to forest management disturbances • Sufficient sample size attainable • Population dynamics stable Since trees meet these criteria, a useful measure might be derived using data from forest inventories or other regular forest surveys (LeMay et al. 1997, Lahde et al. 1999). Stand structure has been defined in terms of the mixture of tree sizes and species, without explicit recognition of other flora and fauna (Buongiomo 1995). The diversity of stand structure indicated by the mixture of tree size and species has been shown to be an important element of total biological diversity (MacArthur and MacArthur 1961, Willson 1974, Franzreb 1978, Temple et al. 1979, Aber 1979, Ambuel and Temple 1983, Freemark 1986), with the highest diversity found in stands where there are multiple species and multiple tree sizes (Buongiomo et al. 1994). It has been argued that managing forests for biodiversity can be accomplished by managing for structural diversity (Onal 1997). Thus, a useful measure of the total diversity of a forest stand may be found by measuring the structural diversity of a stand. 3 Stand structure can be measured based on the distribution of tree ages and sizes (Lahde et al. 1999). Although age is difficult to measure quickly and cheaply, tree size may be readily measured by diameter at breast height (dbh) and height. An index of stand structure as a function of species composition, dbh, and height should be correlated with overall forest diversity, and would be a function of tree variables found in forest inventory samples. The objective of this study is to present possible measures of stand diversity that quantify the diversity in a stand based on both species and size diversity, using commonly measured tree attributes. Several measures are suggested and discussed. The measures were then evaluated using simulated stand data, and applied to the University of British Columbia's Malcolm Knapp Research Forest (MKRF), located in coastal British Columbia. The measures were evaluated using objective criteria and conclusions were drawn as to the performance of the suggested measures. 4 Chapter 2 Common Diversity Measures The most common diversity measures in use are species diversity measures. Species diversity measures can be applied at various spatial and temporal scales, incorporating landscape patterns, as well as stand level variation (LeMay et al. 1997). Biodiversity is often broken up into three spatial components. The a-component of the diversity of organisms inhabiting a variety of habitats is loosely defined as the average diversity within a single habitat (Routledge 1980). The (^ -diversity can be defined as the amount of turnover between habitat types (Routledge 1979). T-diversity describes the variation in a and P -diversity as one moves across regional landscapes and along major environmental gradients (Kimrnins 1997). This paper will focus on a-diversity, that is, diversity at the stand level. Species diversity can be broken into two elements. Species richness refers to the number of species. Species evenness, or species equitability, refers to how equally abundant the species are. From these definitions, species diversity can be classified into three categories (Magurran 1988): 1. Species richness indices are a measure of the number of species in a defined sampling unit. 2. Species abundance models describe the distribution of species as a function, or parameter of an associated function. 3. "Diversity indices" or "heterogeneity indices" are indices based on proportional abundances of species, combining richness and evenness into a single index. There are many measures of species diversity (see Magurran 1988). The diversity measures chosen for discussion in this thesis are Simpson's Index of Concentration, Shannon's Information Index (which is often incorrectly referred to as the Shannon-Weaver or Shannon-Wiener index (Magurran 1988)), and the Q statistic. Based on information theory, Simpson's Index and Shannon's Index are perhaps the most popular diversity measures. Both are nonparametric indices; that is, neither relies on any assumptions about the underlying species abundance distribution. The Q statistic is a diversity measure derived from an abundance model. 2.1 Shannon's Index Shannon's index is a measure of how uncertain we are that an individual picked at random from an infinitely large community will be of a certain species. The more uncertainty one has about the identity of an individual, the higher the diversity of the community. As well as assuming that the population is 'indefinitely large' and that individuals are selected at random, the calculation of this index requires that all species be represented in the sample. The failure to represent all species can be a source of considerable error (Magurran 1988). Shannon's index, H', is defined as follows (Shannon and Weaver 1949): H' = tpihPi [2-1] i=i where: p, = proportion of individuals in the t'th species; and S is the number of species in the population. 6 When the population proportions are unknown, Shannon's Index is calculated with the maximum likelihood estimator (MLE) of Pi = ntIN [2.2] where: «, = the number of individuals sampled of species /, and N - total number of individuals sampled. In practice, the terms pu nu and N are not necessarily computed with the number of individuals, as individual trees, plants, or animals can vary greatly in size. The proportion of a species has been measured in a variety of ways, including: number of individuals (Franzreb 1978, Swindel et al. 1991, Niese and Strong 1992, Condit et al. 1996), basal area (McMinn 1992, Harrington and Edwards 1995), stems per ha (McMinn 1992, Harrington and Edwards 1995), foliar cover (Swindel et al. 1984, Lewis et al. 1988, Qinghong 1994, Corona and Pignatti 1996), crown cover (Corona and Pignatti 1996), and biomass (Swindel et al. 1984, Swindel etal. 1991). Shannon's Index is often applied in the form of a relative index. This relative index is the ratio of absolute diversity to the maximum diversity possible. This index can be taken as a measure of evenness (Pielou 1977): E = H'IH'mx=H'l\xiS [2.3] 2.2 Simpson's Index Simpson's index of concentration, A, is the probability that two individuals selected at random from an infinite population are of the same species (Simpson 1949). Simpson's 7 index differs from Shannon's in that Simpson's index decreases as diversity increases. Thus, Simpson's Index is often expressed as l/A or exp(A). Simpson's index is defined as follows (Simpson 1949): where: /?, = proportion of individuals in the ith species. When a sample is chosen at random from a population, Simpson's index can be computed with the following unbiased estimator (Simpson 1949): where: n, is the number of individuals of the ith species; and N = total number of individuals sampled. As with Shannon's Index, the terms p„ n„ and N are not necessarily computed with the number of individuals. The proportion of a species has been measured in a variety of ways, including: number of individuals (Swindel et al. 1984, Condit et al. 1996), basal area (Harrington and Edwards 1995), stems per ha (Harrington and Edwards 1995), foliar cover (Swindel et al. 1984, Lewis et al. 1988), crown cover (Swindel et al. 1991), and biomass (Swindel et al. 1984, Swindel etal. 1991). 2.3 The Q Statistic The Q statistic can be considered to be a species richness index, but it is a statistic derived from the species abundance distribution of a particular population. The Q statistic [2.4] 8 measures the inter-quartile slope of the cumulative species abundance curve, providing an indication of the diversity of the community without weighting toward the very abundant or very rare species (Kempton and Wedderburn 1978): 0.5n„, +yR2 '« +0.5nB1 Q = — [2.6] log(/?2//fl) where: nr is the total number of species with abundance r; Rl and R2 are the 25% and 75% quartiles, respectively; and nm and n«2 are the number of individuals in the classes where Rl and R2 fall, respectively. This approach takes into account the distribution of species abundances, but does not actually entail fitting a model. 2.4 Relationships Among Indices Hill's (1973) indices provide an elegant relationship between S, A, and H'. Hill's indices were derived from the observation that these indices can be arranged by their propensity to weight toward rare species or common species (i.e., toward species richness or evenness). Hill (1973) defined a diversity index as "reciprocal mean proportional abundance", with the ath order of diversity given by: A l = (5>r I"""" . for s o m e « > 0 [2-7] where: pi = proportion of individuals in the /th species. It has been proven that if an index of diversity, N, satisfies the following properties, it must be one of Hill's indices (Routledge 1979): 9 HI: N(l/S, 1/S,. . . , 1/5) = S . For a population with a given number of species, the measure of diversity will be at a maximum when all the species are present in equal proportions (or with maximum evenness). U2:N(pi, p2,...,ps)<S H3: is a continuous function in each pt, for pi in [0, 1]. H4: Species diversity = genus diversity times average species diversity within genera, where the average should satisfy: 1. Consistency: If all the genera are equally diverse, the average is N. 2. Strict monotonicity: If the diversity of any genus is increased, the average increases. 3. Quasi-linearity: Averages can be performed in stages. 4. Homogeneity: If the diversity of each genus was increased by a factor c, then the overall diversity should be increased by c. 5. Symmetry: The average will be independent of the ordering of the genera. S, Simpson's Index, and Shannon's Index satisfy the above properties, as: N 0 = S N , =exp(#') Several authors (Routledge 1980, Kobayashi 1981) advocate the use of Shannon's and Simpson's indices in Hill's form, as they are more stable and less sensitive to sample size. 10 2.5 Using Indices to Make Inferences about Populations Since diversity indices attempt to compress a multidimensional description into a single number, they can be ambiguous when making comparisons between populations (Hurlbert 1971). However, the effects of such ambiguity can be greatly reduced by restricting the use of indices to comparisons of similar collections of species (Zahl 1977). The comparability of two communities will depend on: 1) taxonomic definition, 2) temporal scale, and 3) the region or space that the community occupies (Pielou 1975). Differences in estimated diversity from samples of two different populations can result from three effects (Smith and Grassle 1977): 1. true differences in population diversity; 2. sample size dependent bias of the diversity estimator; or 3. random sampling error. It may be difficult to separate these effects, particularly when the sample size is small and the bias is large (Smith and Grassle 1977). The traditional method for testing for significant differences in diversity indices is a t-test, here shown with Shannon's Index (Hutcheson 1970): t, H' - H'2 (d.f.) (VaxH' + VarH'2) [2.9] d.f.= (Var/^+Var/ /^) 2 [2.10] (Var//', ) 2/N, +(Var H\ )2/N2 11 This test, although conservative, may be used only if the populations are normally (or nearly normally) distributed (Hutcheson 1970), the population is randomly sampled, and the total number of species is known (Heltshe and Bitz 1979, Routledge 1980, Pielou 1966). 2.6 Discussion of Common Diversity Measures in Use Diversity indices attempt to describe both the richness and evenness of a community. A major criticism of diversity indices is that they confound a number of variables that characterize community structure in combining richness and evenness (Ludwig and Reynolds 1988): 1. the number of species; 2. the relative species abundances; and 3. the homogeneity and size of the area sampled. Despite these criticisms, they are widely used in forestry (e.g., Swindel et al. 1984, Lewis et al. 1988, McMinn 1992, Qinghong 1994, Corona and Pignatti 1996). Pielou (1975) gives two reasons to be cautious in the use of Shannon's index. First, it has been shown that Shannon's Index, as it is commonly computed with equation [2.1], is biased. Although this bias can be removed, it requires that S*, the number of species in the population be known; S calculated from a sample is merely a lower bound for S*. Secondly, using the MLE of pt makes an implicit assumption that the sample units are a random sample from a community (Pielou 1975). Although random locations of plots can be assured, the patterns of the organisms' locations may not be random. In natural communities, organisms are more likely to be located in clumps or patches, and thus, this assumption will not be met. 12 One may adjust for these shortcomings by using a re-sampling technique, such as Pielou's pooled quadrat method (Pielou 1966) or a jackknife estimator (Routledge 1980). However, there is a potential for severe bias using the pooled quadrat method (Routledge 1980). When a small number of quadrats is sampled, the jackknife has been found to give better estimates of H' than traditional computational methods or bootstrapping methods (Smith 1984). The relativized Shannon index, E, is often used as a measure of evenness in a community. However, E contains the term 5, which often is not known. Estimates of S from a sample can vary greatly as a result of sample size or stochastic variation, greatly affecting the computation of E. Evenness will be overestimated when S is underestimated and thus, this relative index can be misleading (Peet 1974). Species evenness should only be compared between communities when collection sizes are equal (Kobayashi 1981). The major criticism levelled at Simpson's Index is its insensitivity to rare species (Hurlbert 1971, Pielou 1975, Magurran 1988). However, this criticism does not recognize that diversity is an aspect of community structure, and that structurally, rare species are minor components of their community (Hurlbert 1971). If this viewpoint is taken, Simpson's index may be a better representation of community composition than other diversity indices. The choice of a diversity index may depend on the characteristic of a community that one wishes to measure and how proportions of each species (pi) are computed. Simpson's index 13 is most affected by changes in the abundances of the commonest species (dominance), whereas Shannon's index is more affected by rare species (richness). However, this may be compensated for at least in part by choosing basal area or stems per ha to compute p,. In forest stands, the use of basal area has been recommended, as it gives a better measure of the occupation of a stand (McMinn 1992). Furthermore, one may prefer to use an index with a biological interpretation. Although Simpson's index can be described as the probability of inter-specific encounter, Shannon's index lacks such an interpretability. Kempton and Wedderburn (1978) recommend the Q statistic in preference to Simpson's or Shannon's index for two reasons. First, X and to a lesser extent H', are strongly influenced by the abundance of the commonest species, while Q depends on a wide range of species with medium abundance. Second, Q shows a smaller variability than X or H' when testing replicate samples. There is much debate on the discriminant abilities of diversity indices. Magurran (1988) offers a summary of the performance and characteristics of some selected diversity statistics. Table 2.1 contains a partial list. Table 2.1. Summary of performance and characteristics of selected diversity measures (Magurran 1988). Index Type of index Discriminant Ability Sensitivity to Sample Size Q statistic Richness Good Low S (species count) Richness Good High Shannon Richness Moderate Moderate Simpson Dominance Moderate Low Shannon Evenness Evenness Poor Moderate 14 Authors vary on their preference of indices. Simpson's is preferred over Shannon's by Peet (1974) and Routledge (1979), while Kempton and Wedderburn (1978) favour the Q statistic. Magnussen and Boyle (1995) indicate that Shannon's index is more efficient than the Simpson Index in terms of number of sample plots required; however, this is contrary to results found by Routledge (1979). Species diversity measures can be applied at various spatial and temporal scales; however, this review focused on diversity at the stand level. Stand level species diversity can measure species richness, species evenness, or both. Diversity indices may be more practical, as they combine richness and evenness into a single index, and there is no need to fit a model to abundance data. Regardless of the diversity index chosen, a consistent approach should be used in order to facilitate comparisons and provide meaningful diversity measures. 15 Chapter 3 Structural Diversity Measures Structural diversity can be as important to overall species diversity as is the diversity of plant species in the forest plant community (Kimmins 1997). This viewpoint is supported by research into avian and insect diversity. Stand structure can also have a profound impact on stand growth efficiency, and measurements of structure can play important roles in predicting the growth of trees. Birds have long been observed to frequent more structurally diverse habitats, depending on certain structures or serai stages for "niches" or "ways of life" (Whittaker 1972). Spatial complexity was cited as a key element of the maintenance of bird species diversity in northeastern and north central mixed hardwood-conifer forests (Temple et al. 1979), while the vegetation profile was found to be more important to bird species diversity than the number of plant species present in the southwestern United States (Franzreb 1978). In northern California, Moen and Gutierrez (1997) reported that spotted owls were found to use areas of high structural diversity as roosts and nest sites. Shrubs, snags, and downed wood, which occur in a wide range of forest ages, were found to be an important part of western North American bird habitat, and bird species abundance was found to be highly related to the surrounding habitat pattern (Vernier et al. 1997). In Australian Eucalypt forests, more species of arthropods and a greater abundance of individuals were found on the most productive plants and where the physical and structural 16 complexity was high (Recher et al. 1996). Furthermore, data also showed that bird species diversity was correlated with the structural attributes of the vegetation and the number of plant genera in the understory. These data led Recher et al. to conclude that the high structural and floristic complexity of the forests which are rich in this regard were correlated with abundant and diverse arthropod fauna. Bird diversity, as well as arthropod diversity, has been shown to be significantly correlated with mature (and more structurally diverse) forests and site productivity (Recher et al. 1996). Since both structural and floristic diversity determine the energy rates, nutrient flux, water, and hence, available resources in the forest (Aber 1979), species' abilities to coexist will be affected by structural and/or floristic diversity and will determine species richness and diversity (Recher et al. 1996). Stand structure can have an impact on stand growth. Gaps in the forest cover can allow for rapid regeneration, while large leaf areas and tree crowns can intercept sunlight and precipitation (Spies 1998). There are numerous studies that show a variety of patterns of growth relating to structural complexity, depending on the species involved (see Oliver and Larson 1996). In general, mixed species stands are considered to have a positive effect on growth, as they tend to slow or prevent the spread of pathogens (Hunter 1990). Multi-cohort stands can have higher growth rates than that of single-cohort stands in situations where lower-strata trees do not die readily (Oliver and Larson 1996). On the other hand, some studies show that yield can be less in these types of stands, but tree size and quality can be higher (Oliver and Larson 1996). 17 Many stand structure diversity measures have been developed. Most measurements of stand structure are 'ad hoc' methods; that is, they are derived from modifying an existing diversity measurement. However, structural indices can be functions of many different stand variables. Since most use either height or diameter size densities, the indices have been classed as vertical, horizontal, or "other". 3.1 Measures of Vertical Diversity Pioneering research on the relationship of structural diversity to songbird populations led to the first attempts at a structural diversity index. MacArthur and MacArthur (1961) constructed foliage height profiles by measuring the amount of vegetation at sixteen heights, from six inches to 60 feet above ground. The foliage height diversity index, FHD, has the same form as Shannon's Index and is calculated as follows: FHD = -JJiPi\npi [3.11] • i where: pt = the proportion of total foliage in the ith layer. MacArthur and MacArthur (1961) used this index in various analyses and found that knowledge of plant species diversity did not provide additional knowledge of bird species diversity. Linear regression of bird species diversity (calculated with Shannon's Index) versus FHD yielded a significant relationship. Moreover, the variability in bird species diversity not accounted for by the variation in foliage height diversity was not accounted for by variations in plant species diversity or latitude (MacArthur and MacArthur 1961). 18 Subsequent researchers into avian community structure in the United States, Germany, and Britain have continued to use MacArthur and MacArthur's (1961) approach, with slight modifications to the calculation of pt. In a variety of midwestern habitats, ranging from grasslands to forest, Willson (1974) calculated pt using presence/absence inthree defined vertical layers. In an analysis of bird communities in Germany, Erdelen (1984) calculated P i using foliage area in twelve defined vertical layers. The FHD of British forests were evaluated by Ferris-Kaan (1998) with P i calculated as the percent foliar cover in four defined vertical height bands. In other research, MacArthur and MacArthur's (1961) original approach has been modified. In a study of northern hardwood forests, Aber (1979) used that original approach, as well made a second parallel analysis, dividing the forest canopy into four strata and calculating pt as the relative distribution of foliage by strata. The diversity of vegetation in southern Wisconsin forests was evaluated by Ambuel and Temple (1983) in two ways. First, FHD1 was evaluated similarly to MacArthur and MacArthur's (1961) original approach, but using Simpson's Index as a base: FHDAT=-Zl/P: [3.12] where: /?, = the proportion of total foliage in vertical layer i. Second, a measure of forest profile diversity, FPD (Ambuel and Temple 1983) was also introduced: FPD = -JJiP\nPi [3.13] 1 To avoid confusion between the definition of FHD used in MacArthur and MacArthur (1961) and that of Ambuel and Temple (1983), the latter will be referred to as FHDAT. 19 where: />, = the proportion of profiles of type /. Ambuel and Temple (1983) designated sixteen profile types, using the presence/absence of vegetation in each possible combination of four height zones. Willson (1974) and Erdelen (1984) confirmed MacArthur and MacArthur's (1961) results, finding that FHD was highly correlated to bird species diversity. However, these relationships did not hold when only forested plots were included in the analyses. When Aber (1979) used MacArthur and MacArthur's (1961) layer definitions, FHD correlated strongly with stand age. On the other hand, there was not a significant relationship with stand age when FHD was calculated with four strata. Results of Ambuel and Temple's (1983) study revealed a strong relationship between FHDAT and the presence of long-distance migrants. f 3.2 Measures of Hor izonta l Diversity The first, measures of horizontal diversity emerged from forest planning models. Distributions of diameters are used as inputs to some natural resources mathematical programming models, and so can be used to juxtapose the evaluation of diversity with economic returns in a forest. Gove et al. (1991) measured diameter diversity using the dichotomous rarity measures presented by Patil and Taillie (1979, 1982), with nine two-inch diameter classes in place of species. In this approach, rarity is measured as: 1 —Vs 18+1 A,(g)= y , P > - \ [3.14] where: K.t = the proportion of individuals in diameter class i; and 20 s = the number of diameter classes. These rarity indices are synonymous with Hill's indices: when /?= -1, Afi (n) is the species count; when /?= 0, A^(^;)is Shannon's index; when /?= 1, A^(^)is Simpson's index (Gove etal. 1991). Later researchers used Shannon's Index to predict the effects of some management regimes on the diversity of tree sizes. Buongiomo et al. (1994) used the diameter distribution as a measurement of stand structure for several reasons. First, because it plays a key role in the structure of forest stands, it is a key element in the biological diversity of stands. Second, since diameter - height relationships are well defined and tree height is positively correlated with FHD, the diameter distribution is a good proxy for FHD. Third, a measure using diameters is more economically feasible than FHD. Buongiomo et al. (1994) measured diameter diversity in northern hardwoods in Wisconsin using Shannon's Index with p, = the proportion of trees in the ith two-inch diameter classes. In the French Jura Mountains, Buongiomo et al. (1995) used a similar technique with five-cm diameter classes. A "normalized absolute deviation" (NAD) index was used by Onal (1997) to measure divergence of a diameter distribution from an arbitrarily specified target distribution: [3.15] n 21 where: y, is the number of individuals in species i; Pi = yi I Zi yi; {Pi }is a specified distribution with /?, >0 and S, /?,; = 1 ; n is the number of classes; and £ = m i n ( # ) . This measure exhibits properties similar to that of Shannon's index when the target distribution is uniform (i.e., /?,= l/h). Like Shannon's index, the maximum value of NAD is reached when the proportion in each class is equal, NAD is scale neutral, and the value of NAD is independent of class ordering2 (Onal 1997). Onal (1997) used simulations to show that there is a strong significant quadratic relationship between NAD and Shannon's Index. 3.3 Other Structural Diversity Measurements Various other measures of structural diversity have been developed in the literature. However, none has gained widespread acceptance. The following methods are only a sampling of the various methods developed. Freemark and Merriam (1986) introduced an index of spatial variability. Their habitat heterogeneity index (HH) is a measure of plant species and forest structure diversity derived from Shannon's index (after Orloci 1970)3: HH=-±±XAn(Xij/Xi) • [3.16] /=i y=i 2 See section 2.4: HI, point 4 in H4, and point 5 in H4, respectively. 22 where: c = number of plots; r = number of classes; and Xij is! the proportion of individuals in the /th class of they'th plot. Equation [3.16] differs from Shannon's index in that the denominator contains the average value of the class, not the total value of the class. HH was computed separately for eight components: tree density, tree dbh class, canopy closure, foliar cover in five defined vertical bands, average canopy height, herb height, percent litter, and percent bare ground. Since components varied in units, magnitude, and number of classes, HH was scaled to facilitate comparability and combinability. To arrive at a total value of forest HH, the components were summed across forest type. Using this index, the authors concluded that the more hetrogeneous habitats in the agricultural landscape near Ottawa had more bird species (Freemark and Merriam 1986). Some researchers have included a measure of structure in Shannon's and Simpson's indices by modifying the way /?, is calculated. Cain and Shelton (1995) computed the p,'s with importance values: IV = relative density + relative height + relative frequency where: "relative" refers to each species' contribution to the total for all species. Then, /?, = proportion of TV in species /. The resulting values for both Shannon's and Simpson's indices were then used to compare the diversity of woody understory growth after disturbance in the southeastern United States. 3 Although this measure of diversity is presented as a forest-level measure (whereas the other measures presented are stand level measures), the measure could easily be modified to a stand level one by removing the second summation sign. 23 A variety of tree and stand variables have been used to arrive at an index of structural diversity. One such recent example (Lahde et al. 1999) uses seven indicator variables by species: stem distribution by diameter size, basal area of growing stock, volume of standing dead trees, volume of fallen dead trees, undergrowth density, occurrence of "special trees", and volume of charred wood. All variables were delineated into classes, e.g., trees were classed into three diameters groups and two basal area groups. Based on the values,of the indicator variables, stands were given diversity scores by species, which were then combined into an index for the entire stand. The index was found to differentiate between forest development types and site types fairly well (Lahde et al. 1999). 3.4 Discussion of Measures Reviews of bird communities and vegetation stand structure subsequent to MacArthur and MacArthur (1961) emphasized that the location of plots greatly affects the discriminant abilities of the foliage height diversity index. However, all but one of the studies cited was not directly comparable to MacArthur and MacArthur's original work. Willson (1974) and Erdelen (1984) assert that the mere existence of a tree layer in the canopy determines bird diversity, regardless of the foliage volume or distribution among layers; however, their results are not truly comparable as classes were not defined in exactly the same manner. Aber (1979) used four strata in the calculation of FHD, finding it uncorrelated with stand age. Again, the calculation of this FHD was quite different from MacArthur and MacArthur's since Aber's strata were of different sizes and at different heights, depending on 24 the sample site. This result is only comparable to MacArthur and MacArthur's results if stand age is seen as a proxy for bird diversity. A variety of techniques were used to assess vegetation cover and FHD, and results depended, at least in part, on the selection of class boundaries. Height and diameter are continuous variables, and as such, the division of data into classes will necessarily depend on the distribution of the data. Since many studies are incomparable, caution should be used in drawing conclusions (Erdelen 1984). Further, no study compared the use of vertical and horizontal stand structure indices, and thus, no comparison may be drawn. Many stand structure measures have been developed to characterize both vertical and horizontal stand structure. However, most measurements of stand structure are derived by modifying existing diversity measurements. Unfortunately, no one method has been used consistently, which makes comparisons and evaluations difficult. Furthermore, results depend on splitting continuous variables into classes. This class division could have a marked effect on results. Some researchers have also used variables such as volume of coarse woody debris or crown ratio, which are not always measured in common forest surveys, and thus, these measures would be difficult to use from an economic standpoint. A more practical measure of stand structure will be one that uses readily available data, does not depend on subjective division of continuous variables into classes, and directly correlates to overall stand biological diversity. 25 Chapter 4 Possible Measures of Stand Diversity Using Stand Structure An index of diversity is needed which will characterize the diversity of stand structure. Stand structure diversity has been defined to include the richness of tree species and sizes (Lande et al. 1999). Further, structural diversity increases as the number of species and number of sizes increases (Buongiomo et al. 1994). This situation may be described using a probability distribution where every value over a specified range is represented equally. That is, every possible value of dbh and height may be found in a stand with equal probability. The probability distribution of maximal diversity under these circumstances is described mathematically by the uniform distribution. Low diversity stands would be defined as those with few species and few tree sizes, and minimum diversity would occur where all trees in a stand are of single tree species and single size. This distribution would be a spike at a single point. In an effort to incorporate stand structure into a diversity measure, several possible indices were derived. For all of these indices, the standard for maximum diversity was a uniform distribution. However, this assumption may be challenged since there is no quantitative definition of high structural diversity. Although a highly diverse stand is well characterized with this distribution, proposed measures could be modified with relative ease for use with another distribution. 26 There are several characteristics which a stand diversity index should possess. The alternative measures will be evaluated with a goal of meeting all or most of these criteria. 4.1 Desirable Characteristics of Measures of Stand Diversity Using Stand Structure An index of forest stand structure should have the following characteristics: 1) The index should measure both richness and evenness of structural and species diversity. A diversity index of this type would take into account the number of individuals of different species and size, as well as the distribution of individuals over those species and sizes. 2) The index should use commonly measured tree attributes only. This would enable forest managers to use the index with inventory data or other common forest sample data without additional expense. 3) The index should not rely on combining data into diameter or height classes. Many of the indices presented in Chapter 3 require that data be segregated into arbitrary classes. This practice is undesirable, as a change in class definitions could cause a change in index. 4) The index should not depend on sample size or plot size. Sample plots from forest inventories may consist of one plot or a cluster of several plots. In order to make comparisons between populations, the index should be insensitive to cluster size and plot size. 5) Horizontal and vertical diversity should contribute equally to the index. It has been shown that both tree diameter and height contribute to diversity by creating a variety of habitats. Thus, given two populations with equal horizontal diversity, the population with greater vertical diversity would have a higher index. 27 6) Basal area should be used to measure proportions of trees in species and size distributions. Although some authors (for example, Lewis et al. 1988, McMinn 1992) have used stems per ha, LeMay et al. (1997) recommended basal area, as it better represented the proportion of resource use. 7) The index should take on values between zero and S. An index value of S should be associated with stands of high diversity, and an index value of zero should be associated with stands of low diversity. 8) The index should measure the deviation between a particular community and that of a maximally diverse community. a) Given a set range, the index should be tree size insensitive, or location invariant. For example, a population with small trees evenly distributed over a ten cm diameter range would have the same index value as a sample with large trees evenly distributed over a ten cm diameter range (assuming all other things being equal). This would ensure that larger individuals would not drive the index. b) The index should be related to the range of values in diameter or height for continuously occupied ranges. In other words, a population with a larger diameter and/or height range that occurs evenly over that range would yield a larger diversity index than a population with values occurring over a smaller range, or over only parts of the range (as in a bimodal distribution). This would ensure that stands that exhibit greater size diversity would have higher index values. c) The index associated with a bimodally distributed population with two separate and distinct canopy layers with ranges of Ranget and Range2, should be at least as great as that of a unimodal distribution with a continuous range = Range/ + Range2. For 28 instance, an old growth stand with regenerating trees in the understory would be considered at least as diverse as a stand with a wide single-canopy layer. If these characteristics are met, the index should be easy to use and interpret, and the maximum value of the index will result for uniform stands with the greatest range. Some of the criteria may not be met, and evaluations with data analysis should reveal whether these criteria are truly necessary. 4.2 Possible Alternative Measures of Stand Diversity Using Stand Structure Several alternative measures were derived based on existing measures, such as Shannon's Index, or based on techniques familiar to foresters and biologists, such as curve fitting and variance measurements. 4.2.1 Shannon's Index Extended to Diameter, Height, and Species Shannon's index (Shannon and Weaver 1949) has been computed in the forestry and ecology literature with a species distribution, a diameter distribution, or a height distribution. However, Shannon's index computed with diameter, height, and species was not found in the literature. The index could be used separately for each of these variables, producing H'cl, H'h, and H's, respectively, and then these indices could be combined to produce a stand index, H'd+h+s. Alternatively, data could be classed simultaneously into diameter, height and species classes, producing one index, H'ilhs. The former method will be referred to as the 'Post-hoc Method', and the latter will be referred to as the 'Combination Method'. 29 4.2.1.1 Post-hoc Method The procedure for computing the post-hoc Shannon's index for a sample plot is: 1) Tree data are combined into diameter classes using dbh. The proportion of trees in the ith diameter class, /?,, is the percentage of basal area in that diameter class. Then, Shannon's Index is computed with P l and equation [2.1] to yield H'(l. 2) Tree data are combined into height classes. The proportion of trees in the jth height class, Pj, is calculated as the percentage basal area in that height class. Then, Shannon's Index is computed with pj and equation [2.1] to yield H'h . 3) Tree data are combined by species. The proportion of trees in the kxh species, pk, is calculated as the percentage basal area in that species. Then, Shannon's Index is computed withpk and equation [2.1] to yield H's. 4) The final index is computed as the average of the diameter, height, and species indices, as it is a convenient way to keep the index on a similar scale of values as the original Shannon's Index: H'd+h¥, = average (H'd,H'h,H's)= H ' " + H > + H * [4.17] The minimum value of the Extended Shannon Index under the post-hoc method is zero. This minimum value results when only one species, one diameter class, and one height class are present in a sample. The maximum value of the Extended Shannon Index under the post-hoc method is reached when H'd, H'h , and H's are at their maximum values. The maximum value of the H's is ln(5) (Magurran 1988). This maximum is achieved when each species is represented equally within a plot (uniformly distributed). Similarly, if all diameter and 30 height classes are represented equally within a plot, the maximum values of H'd and H'h are \r\{ni) and ln(n2), respectively, where ni is the number of possible dbh classes and n2 is the number of possible height classes. The overall maximum of the Extended Shannon Index under these circumstances is the average of these maximum values: m a x ( „ - ^ ) = M ^ % l ± i ! ^ [ 4 . 1 8 ] 4.2.1.2 Combination Method The procedure for computing the combination index for a sample plot representing one forest stand is to first combine the data simultaneously by species, diameter class, and height class. The proportion of basal area in the ith diameter, 7th height, and klh species class, />//*, is then calculated. Then, Shannon's Index is computed with and summed over all classes to yield H'dhs . As with the post-hoc method, the minimum value of the Extended Shannon Index under the combination method is zero. This value results when there is only one species in a plot, and this species occurs in only one diameter-height class combination. The maximum value of the Extended Shannon Index under the combination method will be reached when all diameter and height classes, as well as species are represented equally within a plot (uniformly distributed). This maximum value is: max(if ^ ) = ln(5 x n, x n2) = 3 x max(H'd+h+s) [4.19] 31 4.2.1.3 Properties of the Extended Shannon's Index The Extended Shannon Index, in both the post-hoc and combination versions, measures richness and evenness of stand structure, giving equal weight to horizontal and vertical diversity. The indices use basal area and require only commonly measured tree attributes in their computation. Shannon's Index also has been shown to be relatively insensitive to sample size (Magurran 1988), and attains its highest values when distributions are perfectly even (i.e., uniform). The Extended Shannon indices should be insensitive to tree size, as indices are weighted only on the proportion of basal area occurring in a particular class. Both forms of the Extended Shannon Index possess the desirable characteristics listed as items 1, 2, 4, 5 and 6 (Section 4.1). Item 8 cannot be considered without further using data. Item 7 is not satisfied because the maximum values of these indices are functions of the number of possible diameter and height classes. This dependence results in indices that have a lower bound of zero, but are not bounded by S. Item 3 is not met because the indices rely on pre-defined classes. This is undesirable, as this creates issues such as where to put class boundaries. A change in class boundaries or an increase in the number of classes would invariably change the value of the indices. Even if set classes were used, there would be pronounced differences simply due to geography. For instance, in Coastal British Columbia, trees are larger on average and have a wider range of sizes than in the interior of the province. More classes would be represented in Coastal areas, resulting in larger indices. On the other hand, more species may be represented in sample plots from Interior areas of the province, which would also result in larger index values. For this reason, care should be taken in comparing diversity indices (of any type) from different geographic areas. 32 To alleviate problems associated with classes, an integral form of the Extended Shannon Index could be presented. However, this option would rely on finding functional forms for species, diameter, and height distributions. While a functional form may be found, it may not have a closed form, or may not be flexible enough to fit the variety of distributions that may be found. This issue will be addressed in the next section. 4.2.2 Structural Index based on Curve Fitting When a continuous variable, such as height or dbh, is translated into a class variable, information about the distribution of that variable is invariably lost. In order to avoid using diameter and height classes, which may affect the index, an index based on a continuous distribution is desirable. Foresters and biologists have used continuous distributions to describe the frequency distributions of diameter and height both univariately and bivariately. These distributions have included the Beta, Johnson's SB, and Weibull, with varying degrees of success (Hafley and Schreuder 1977). In fitting such a distribution, the resulting functional form may be used in an integral form of Shannon's Index, or the Q statistic could be used to measure diversity based on the derived curve. However, if data are sparse, a derived curve may cover a range much greater than that of the data. The use of the Extended Shannon Index or the Q statistic in this type of situation would involve making an inference about the population using the fit of a distribution with a restricted range. This may be inappropriate, and could produce misleading results. Hence, this option was not further investigated. 33 Alternatively, in fitting a regression curve, the regression coefficients may indicate the centre, shape, or scale of a distribution. For instance, the two coefficients in the function describing a univariate normal distribution indicate the mean and variance of the distribution. A known distribution could be fit to the data and its regression coefficients could be analyzed as indicators of diversity. Of the candidate distributions, the Beta distribution offers the most flexibility in terms of skewness, kurtosis and shape of the curve (Hafley and Schreuder 1977). However, the Beta is a complex function, which involves the computation of Gamma function values and may be extremely difficult to fit. Johnson's SB and its bivariate form, Johnson's S B B have been shown to perform consistently better than the Weibull or Beta distributions in fitting diameter and height data (Schreuder and Hafley 1977). However, this probability density function does not have a closed-form cumulative density function (Schreuder and Hafley 1977), and thus, cannot be used without allocating data to classes or making an approximation with numerical integration. The Weibull distribution has been used extensively in forestry applications (Clutter et al. 1983). It has advantages in that it is fairly flexible in shape and much simpler than some other functions. Although the univariate and bivariate Weibull distributions can closely approximate a normal curve under certain conditions (Bailey and Dell 1973), approximations 34 to univariate and bivariate uniform distributions are only adequate over a portion of the range of the data. A desirable distribution for investigation would be one that is flexible in shape, but relatively easy to fit. Considering the advantages and disadvantages, the Weibull distribution (in both univariate and bivariate forms) was selected as the base curve for fitting. 4.2.2.1 Univariate Curve Fitting The cumulative univariate Weibull probability functions for diameter and height based on basal area are as follows: where: b/, b2, c/, and c2 are coefficients to be estimated; and dbh and height are the independent regression variables. Equations [4.20] and [4.21] would be fit for dbh and height separately by species using non-linear fitting methods. A good fit would be indicated by a small root mean square error (RMSE), computed as: V n — m where: y, is the dependent variable, the cumulative distribution of the sample (Fi(dbh), or F2( height)); [4.20] [4.21] [4.22] 35 y. = is the predicted value of the cumulative distribution under the assumed theoretical distribution; n = the number of observations; and m = the number of coefficients in the regression model. The estimated coefficients could be analyzed for patterns relating to diversity and a coefficient chosen as a representative diversity index. By species, the diversity index would be the average of the dbh and height diversity indices. The total diversity index for a plot would be computed as the sum of all species diversity indices. Thus, if two plots have the same structure, the plot with more species would have a higher diversity index. 4.2.2.2 Bivariate Curve Fitting The bivariate Weibull has been derived in many forms (see Lu and Bhattacharyya 1990); however, based on pilot investigations, the best performing form was given by Hafley and Schreuder (1976): G(dbh, height) = F, (dbh) x F 2 (height) x [l + a 3 x (l - F, (dbh))x (l - F 2 (height))] [4.23] where: Fi(dbh) and F2(height) are defined as above; and a3, along with bi, b2, Cj, and c2 from [4.20] and [4.21], are coefficients to be estimated. F(dbh, height) is computed as the cumulative basal area at each dbh - height combination represented in the plot. Equation [4.23] would be fit for dbh and height simultaneously, by 36 species, using non-linear fitting methods. Starting values for bj, b2, cj, and c2 could be obtained from the univariate regressions. The resulting estimated coefficients would then be analyzed for patterns relating to diversity. If possible, a coefficient would be chosen as a representative diversity index. As in the univariate case, the total diversity index for a plot would then be computed as the sum of all species diversities. Thus, all things being equal, if two plots have the same structure, one with more species would have a higher diversity index. 4.2.2.3 Properties of a Structural Index Based on Curve Fitting Both the univariate and bivariate methods possess the desirable characteristics listed as items 1, 2, 3, 5, and 6 (Section 4.1). An index derived from the Weibull distribution would measure richness and evenness using basal area and commonly measured tree attributes. The index would not rely on combining data into classes, and would account for vertical and horizontal diversity equally. However, items 4, 7, and 8 cannot be considered without further evaluation using data. Both the univariate and bivariate methods presented above rely on the fit of the data to a known distribution. To assure that no inappropriate inferences are made, this method should be used only if samples yield similar fit statistics, such as small residual error and bias. A more significant problem may arise in that no coefficient may be adequate as a measure of diversity. 37 4.2.3 Structural Index Based on R 2 (STRI) The coefficient of determination, or R , is the measure of how much of the variation in a regression line (or curve) is explained by the independent variable(s). Thus, a function representative of maximum diversity, such as the uniform, could be put into the form of a cumulative distribution function (cdf), and compared to a cumulative empirical sample distribution. If the cumulative empirical sample is regressed on the function of the theoretical cdf, the resulting R 2 will give a measure of the deviation of the sample from maximum diversity. A high R 2 (close to one) would indicate that the empirical sample came from a distribution that was similar to that of a high diversity population; a low R 2 (close to zero) would indicate that the sample came from a population with a distribution that deviated substantially from that of a high diversity population. R is traditionally computed as follows: R 2 _ 1 residual sum of squares _ 1 2X1 (y i~& ) 2 [424] total corrected sum of squares (y -y)2 where: y,- is the dependent variable, the cumulative distribution of the sample (F(dbh), F(height), or F(dbh, height)); y,. = the predicted value of the cumulative basal area distribution under the assumed theoretical distribution; y = the mean value of the cumulative basal area distribution of the sample; and n - the number of observations. If the theoretical distribution is taken to be a uniform distribution, then the regression form would necessarily go through the origin. This changes the formula for R to: 38 total uncorrected sum of squares residual sum of squares = 1- SL,(yj-y,-)2 i ; = , ( ^ , ) 2 [4.25] Although, [4.25] is algebraically correct, it no longer has the same interpretation as [4.24] (Kozak and Kozak 1995). However, this R 2 could be used as a measure of diversity, as it describes the amount of dispersion between some theoretical distribution of maximal though it may fall below zero in cases where the shape of the empirical distribution is much steeper than that of the thoretical. Since this is a variation on the R value, this measure will be known as STRI. 4.2.3.1 STRI Univariate Method For the univariate case, dbh and height are regressed separately by species. The cumulative empirical basal area distributions of dbh and height for species k are labelled as: F, (dbh)), and F2 (height)k. If it is assumed that a uniform distribution represents the highest possible diversity in a community, then the theoretical predicted values of the cumulative distribution functions of dbh and height for any species are represented by: diversity and a sample. R 2 gives a numerical value between zero and one in most cases, a4 < dbh < bA [4.26] F2 (height) = height — a , a5 < height < b, [4.27] where: a4 and as are the minimum values of dbh and height, respectively; and b4 and bs are the maximum values of dbh and height, respectively. Then, S T R L j b h and STRIheight for species k are: 39 S T R I - i a f a ^ ^ - ^ ^ J R 4 ? 8 1 S T R I d b h t - l [4.28] Z , - = i ( * i ( ^ , - ) i t ) i ( F 2 ( ^ g f a j ) t ~ ^2(heightj)f S T R I H E I G H T T -1 v"2 /j? ,u • , V [ 4 - 2 9 ] L^Wiiheight^) where: dbhj corresponds to the ith value of dbh; heightj corresponds to they'th value of height; ni is the number of observations with measured dbh; and n2 is the number of observations with measured height. To arrive at a measure of structural diversity for species k, STRLjbh^ and STRIheight^ are averaged to produce STRI ( d + h ) t : S T R I j k h "t" STRI h- i ( > h t STRI ( d + h )^ = A ^ - 2 ^ [4.30] An overall measure of diversity for a sample plot will beSTRI d + h , the sum of the values of all STRI ( d + h ) j over a plot: S T R I d + h = l t i S T R I ( d + h ) t [4-31] The maximum value of the index is S, the number of species. This maximum would occur if both of the cumulative distributions for dbh and height match that of the uniform for each species. The minimum value of this index is negative, and would occur if dbh or height were a single value at either the minimum or maximum value of the possible range. 40 4.2.3.2 STRI Bivariate Method For the bivariate case, dbh and height are considered simultaneously by species. The cumulative empirical sample distribution of dbh and height for species k is labelled as: G(dbh, height)k. If the assumption is made that a bivariate uniform distribution represents the highest possible diversity in a community, then the theoretical predicted cumulative basal area distribution function of dbh and height for species k is: height - a7 G(dbh, height) k dbh - a6 b 6 ~ a 6 x a6 < dbh < b6 [4.32] a7 < height < b7 b7 -a7 where: and 0,7 are the minimum possible values of dbh and height, respectively; and bf, and by are the maximum possible values of dbh and height, respectively. Then, STRIdh for species k is: Q T D T 1 (G^dbhJ • j )k - G(dbht, height j )f O I K I J U — 1 7 r - 4.33 llUllMdbh^height;),)2 1 where: (dbhi, height/) corresponds to the ith value of dbh andy'th value of height. An overall measure of diversity for a sample plot is STRIdh, the sum of the values of all STRLjh^ over a plot: S T R I d h = X L S T R I d h t [4.34] The maximum value of STRI^ is S. The maximum would occur if the distribution of the sample was a bivariate uniform for each species in the sample. The minimum value of this index is negative, and would occur in cases where the surface of the empirical joint distribution is much steeper than that of the bivariate uniform distribution. 41 4.2.3.3 Properties of STRI Both the univariate and bivariate STRI possess the characteristics 1, 2, 3, 5 and 6 listed in Section 4.1. An index derived from the R 2 would measure richness and evenness using basal area and commonly measured tree attributes. The index would not rely on combining data into classes, and would account for vertical and horizontal diversity equally. Item 4 will not be met, as R is dependent on sample size. If n is equal to the number of coefficients in the regression equation, then R 2=l; as n gets large (greater than 30), R 2 decreases, regardless of the fit of the data (Neter et al. 1996). Items 7 and 8 will only partially be met. While the index would be correlated to the range of values in dbh and/or height and would adequately describe bimodal distributions, the index will be sensitive to tree size. For example, consider a forest type where the maximum range of dbh values is 0 to 100 cm, and where there are three stands which have three univariate uniform dbh distributions, covering three equal portions of this 100 cm dbh range: 0.0 to 33.3 cm, 33.3 to 66.7 cm, and 66.7 to 100 cm. If the three distributions have the same proportions of basal area equally spaced over their respective ranges, the denominators of each of the three STRI's derived from the distributions will be equal. The numerator of the univariate STRI is dependent on the area between the empirical and theoretical cdf's; it is the squared distance between the two curves. Graphically, this is represented by the shaded areas in Figure 4.1. 42 Figure 4.1. STRI and three equal dbh ranges. 0 10 2 0 3 0 4 0 5 0 6 0 7 0 80 90 100 EUi(on) The grey region and the blue region (corresponding to the small stand and the large stand) have equal areas. However, it can be proven mathematically that the yellow region (corresponding to the mid-range stand) has a smaller area than that of the large and small stands. Thus, the STRI associated with a uniform distribution centred in the middle of the theoretical range will be higher than that of uniform distributions centred either at the low end or high end of the range. 4.2.4 Structural Index based on Variance Comparisons (STVI) The distribution of diameters and/or heights in a structurally diverse forest may be described in terms of variances. The uniform distribution yields a very high variance, though not the highest possible variance for any distribution over a given range. The diversity of a stand may be evaluated on the measured stand variance relative to that of a uniform distribution. 43 4.2.4.1 STVI Univariate Method The variance of an empirical dbh or height distribution can be given by: S2 = ^ ' - l L ' — [4.35] S,-=iw« where: xt = dbhi or heightf, x is the mean of dbh or height; Wj is the weight of the ith cell, defined as proportion of basal area of the plot represented by the ith dbh or height point; and n is the number of dbh or height points. This estimate of the variance differs from that of the usual estimator used in forestry, which is weighted by trees per hectare. This weighting is consistent with the other proposed measures. The variances of univariate uniform distributions of dbh and height weighted by basal area are given by: ^ = ^ f ^ [4.36] where: a8 and ag are the minimum values of dbh and height, respectively; and bs and b9 are the maximum values of dbh and height, respectively. The maximum possible variance of a distribution occurs when the distribution is maximally bimodal. For the univariate case with dbh alone, this occurs when half of the basal area is at as and half the basal area is at bs. In this situation, the variance estimator weighted by basal area is given by: 44 'dbhm — X 2 +—x 2 [4.38] Similar results can be obtained for height by substituting coefficients a9 and b9 for as and respectively. Assuming that a uniform distribution represents a maximally diverse population, the difference between the variance of the dbh distribution from an extremely diverse population (Sdbhaiv) a n d t n a t ° f t n e uniform distribution (Sdhhu ) will be close to zero. However, the difference between the variance of a low diversity dbh distribution (Sdhh] ) and that of a uniform is not as predictable. If the distribution consists of trees that are all of the same dbh sdbhkm w i l 1 b e c lose to zero and Sdhh]m -Sdbhu will be large and negative (at most, -Sdhhv ). If the distribution is maximally bimodal, then Sdhh^ will be very large and Sdhhi -Sdhhv will be large and positive (at most, Sdhhm^ - Sdhhu ). In order to use these results to yield an index that measures the difference between the variance of a sample plot and that of the maximally diverse case, both instances must be considered. The following formula can be used to define a diversity index, STVLjbh for a species k: STVI d b h t = < 1-f ^2 _ n Jdbh,, °dbk 2 y 'k V 'dbh. ,when5 2 <S2 1 -U J 2 ( s2 - s Jdbhk ^dbhu m x S d b h m a ~Sdbhu \A>2 [4.39] where: Sdhh - is the variance of dbh for species k; 45 pi and p2 are constants > 0; and m is a constant > 1.0 . The constants pi and p2, define the shape of the curve relating the value of the index to the sample variance: when pt (or p2)< 1, the curve is concave upward; when pt (or pi) = 1, the curve is segmented linear; when pt (or p2) > 1, the curve is concave downward (Figure 4.2). If pi =p2 > 1, then a smooth, continuous function results. Figure 4.2. Comparison of STVI values with m= \ and varying values of pi and p2. 1.000 "j 7 1 - ^ 1 1 The coefficient m controls the value of the index when the distribution is maximally bimodal (i.e., equation [4.38], for dbh). If m= I, then the index will be zero for a maximally bimodal distribution; as m gets larger, the index value increases for the maximally bimodal case (Figure 4.3). 46 Figure 4.3. Comparison of STVI values with p1-p2=2 and varying values of ra. 1.000 0.000 -m=l •m=l . l •m=1.2 • uniform • maximum Variance The values for pi, p2, and ra may be chosen by placing three constraints on the index to yield certain index values under defined conditions. First, when a stand is uniform over half the maximum possible range, the value of the index was constrained to equal 0.5. Second, when a stand is distributed bimodally, such that half of its values are uniformly distributed over the lower quartile of the maximum possible range and half of its values are uniformly distributed over the upper quartile of the maximum possible range, the index was constrained to equal 0.5. Third, when a stand is distributed such that its variance is at the maximum value, the index was constrained to equal 0.1. Algebraically, these constraints are: 0.5 = STVI d b h i =1-( K1 - ?2 °dbhu °dbhil}u 'dbh,. [4.40] 0.5 = STVI d b h t = 1-f s2 -s2 ^dbtl^ ^dbhB m X S d b h m M ~Sdbhu \P2 [4.41] 47 where: Sdhhn is the variance associated with the half-range uniform distribution described in the first constraint; and Sdhh is the variance associated with the bimodal distribution described in the second constraint. After algebraic manipulations, [4.40], [4.41], and [4.42] simplify to: ln(l/2) „ A n n A " - l ^ 3 2 - 4 0 9 4 [ 4 ; 4 3 ] „2= 1 ^ = 0.5993 [4.44] ln(3/8) and m — — 3 1 + e x p f ln(0.9)xln(0.75)-(ln(2))2 ^ V ln(9/5) 1.1281 [4.45] With these choices of/?/, p2, and m, plots with variances close to that of the uniform would have index values close to one. Plots that are of low diversity and unimodal (e.g., the variance is close to zero) would have index values close to zero. Plots that are of low diversity and bimodal (e.g., variance is close to the maximally bimodal case) would have index values close to 0.1. 48 Similarly, a diversity index, STVIheight, based on height could be developed. To arrive at a measure of structural diversity for species k, S T V L j b h ^ and S T V I h e i g h ^ are then averaged to produce S T V I ( d + h ) t : , STVI d b h +STVI h e i e h t S T V I ( d + h ) t = [4.46] An overall measure of diversity for a sample plot is labelled as S T V I d + h , the sum of the values of all S T V I ( d + h ) t over a plot: S T V I d + h = 2:L S T V I ( d + h ) , [4-47] The maximum value of STVI^+h is S. The maximum occurs when the univariate distributions of dbh and height are uniform for each species in the plot. The minimum value is zero, and results in cases where each of dbh and height only exist in a single value for each species represented in the plot. 4.2.4.2 STVI Bivariate Method The variance of an empirical bivariate distribution is usually described with the generalized variance (Johnson and Wichern 1998), which is equal to the determinant of the variance-covariance matrix: ' dbh, height = det 'dbh cov(dbh, height) cov(dbh, height) ' height [4.48] where: Sdhh and Sleight, are given by equation [4.36]; and cov(dbh, height) = • S"=i lwi x (dbfy - dbh) x (height, - height) Ti=\wi [4.49] 49 where: wt is the weight of the ith cell, defined as proportion of basal area of the plot represented by the ith (dbh, height) pair; and n is the number of (dbh, height) pairs. The variance of the bivariate uniform distribution of dbh and height is given by the determinant of its variance-covariance matrix: det 'dbh,. cov(dbh, height) v cov(dbh, height) v heightu [4.50] where: Sdhh{j and S%eightu are given in [4.36] and [4.37]; and cov(dbh, height) v - 0. Thus, [4.50] simplifies to: (b8-a^)2 det 12 0 0 (b9 -a9)2 12 (h - Q 8 > 2 > < ( f e 9 ~ a 9 ) 2 144 [4.51] As in the univariate case, the maximum possible variance occurs when the population is maximally bimodal. This occurs when one quarter of the basal area is, at (as, b$), one quarter is at (as, bg), one quarter is at (a% b8), and one quarter is at (ag, b9). In this situation, the variance of dbh is as in [4.37] and the variance of height can be computed similarly. The covariance of dbh and height is identical to that of the uniform: cov(dbh, height)^ = 0 [4.52] The maximum variance of a joint distribution of dbh and height is given by the determinant of its variance-covariance matrix: 50 SL, = det Q>%~a*Y 4 0 0 (b9-a9)7 _ (fr 8 -a 8 ) 2 x(fr 9 -a 9 Y 16 [4.53] The formula given as [4.39] for STVIdbh can now be modified to use the generalized variance to arrive at an index of diversity for a species k using a bivariate distribution: STVI d , = \ 1-f ° U ^k s2 - s2 , when S; < 5 0 \P2 [4.54] , when S2k > S2, where: S% is the generalized variance for species k; pi and p2 are constants > 0; and m is a constant > 1.0 . As in the univariate case, the powers, pi and p2, define the shape of the curve, and the coefficient m controls the value of the index when the distribution has the maximum variance. Imposing the same constraints on the index as in the univariate case, the same values of pi, p2, and m would be used in the bivariate case. Using these values of pi, p2, and m, sample plots with variances close to that of the uniform would have index values close to one. Sample plots that are of low diversity would have index values close to zero if they were unimodal or 0.1 if they are bimodal. An overall measure of diversity for a sample plot is labelled as STVLjh, the sum of the values of all STVI d h ; over a plot: S T V I d h = l f = 1 S T V I d h A [4.55] 51 The maximum value of STVLjt, is S. The maximum occurs when the joint distribution of dbh and height is uniform for each species in the sample. The minimum value is zero, and results when there are data at only one dbh - height combination for each species represented in the sample. 4.2.4.3 Properties of STVI The STVI in both the univariate and bivariate forms possesses the characteristics presented as items 1, 2, 3, 5, 6, and 7 (Section 4.1). Both forms of the STVI would measure richness and evenness using basal area and commonly measured tree attributes. The indices would not rely on combining data into classes, and would account for vertical and horizontal diversity equally. In order to satisfy criteria 4, the indices should not depend on sample size; however, this criteria would need to be evaluated using data. Since the indices are based on the variances of the dbh and height distributions, weighted by basal area, they should be insensitive to nominal tree size with respect to the distribution of basal area . However, since basal area is used as a weight, larger trees are given more weight. The indices should also be j correlated to range, as the variance of distributions with wider ranges tends to be larger. Furthermore, providing that pi,p2, and m are well chosen, the index value for a bimodal distribution with two distinct ranges should be the same as that of a unimodal distribution with a continuous range equal to the sum of the bimodal ranges. Thus, the desirable characteristics listed as items 8a, 8b, and 8c (Section 4.1) should be met; however, further evaluation is necessary. 52 The STVI may be the most flexible of the proposed indices presented, as the shape of the curve and its response to bimodal distributions may be modified to fit the users' expectations. Furthermore, this index does not require fitting a distribution to the data, making its application more simple. However, a disadvantage of these indices are their lack of theoretical background. Existing indices, such as Shannon's Index, have been well studied and have well-defined sampling distributions. The sampling distributions of these proposed indices will not be known unless further research is conducted. 53 Chapter 5 Evaluation of Alternative Measures of Structural Diversity Evaluations of the alternative measures with respect to item 8 (Section 4.1) were performed using simulated stands and using data collected during the 1995 Malcolm Knapp Research Forest (MKRF) inventory. Methods for generating the simulated data and collecting the MKRF data are presented, followed by results in using the proposed indices. 5.1 Simulated Tree Data Proposed measures were evaluated using simulated stands of varying structural diversity. Simulated stands are an invaluable tool in testing indices in that the shape of the distribution of dbh and height can be controlled, allowing for the evaluation of well-defined cases, which may be difficult to find in natural data sets. 5.1.1 Methods A baseline of maximally diverse structure was first chosen, then simulated stands were compared to this baseline. 5.1.1.1 Baseline Maximal Diversity High diversity stands were defined as those stands that have a large number of tree sizes and species present, and maximum diversity would occur when the maximum number of tree sizes and species are present. Mathematically, this is described with a uniform distribution: fW = , * • ,min(X)<X<max(X) [5.56] max(X)-min(X) 54 where: X represents the variable of interest, which may be dbh or height; min(X) is the minimum value of X; and max(X) is the maximum value of X. The uniform distribution may also be extended to the bivariate case, for simultaneous consideration of dbh and height: where: X and Y represent the variables of interest, which may be dbh and height; min(X) and min(F) are the minimum values of X and Y, respectively; and max(X) and max(T) are the maximum values of X and Y, respectively. By convention, if X is distributed as a univariate uniform variate over a range r, to r2, it is denoted as X ~ U[o, r2]. If X and Y are distributed jointly as bivariate uniform variates with ranges for each of X and Y of ru to rt2 and r2] to r22, respectively, their joint distribution is denoted as (X, Y) ~ \J[(rn, r,2), (r21, r22)]. Using a bivariate uniform distribution over the possible range of diameter and height as the standard for maximum diversity, (dbh, height) ~ U[(min(dbh), max(dbh)), (mm(height), max(height))]. 5.1.1.2 Simulated Stands of Varying Diversity Simulated stands were used to represent the most diverse case (the baseline), as well as other cases representing various degrees of diversity. Trees of various sizes can be generated 8(X,Y) = max(X)-min(A r) max(F) - min(F) 1 1 min(X)<X <max(Z) min(F) < Y < max(F) [5.57] 55 assuming a variety of distributions. There is evidence to suggest that the distribution of sizes of real tree data may approximate known distributions. For instance, the distribution of stems per hectare of single-cohort, single-species stands have been characterized as following the normal distribution, whereas well-differentiated stands can be bimodal or skewed (Oliver and Larson 1996). A skewed distribution may result when intermediate or suppressed trees die readily, or in very old or very young stands (Oliver and Larson 1996). In mixed species stands where tolerant suppressed trees form lower strata, a bimodal distribution may result (Oliver and Larson 1996). In order to evaluate the possible structural diversity measures, stands were simulated that represent different degrees of diversity in terms of the spread and shape of the distribution of tree sizes by basal area. To represent these distributions, notation similar to that introduced for the uniform will be used. If X is distributed as a univariate normal variate with mean ju and standard deviation cr, it is denoted as X ~ N[//, d\. The bivariate normal notation for X ~ N[jU/, 0 7 ] and Y ~ N[//2, cr?] is (X, Y) ~ N[(///, 0 7 ) , {fi2, oi)]. Given set ranges of dbh and height, Figure 5.1 shows the possible combinations of values of dbh versus height that may exist in a simulated sample plot. 56 Figure 5.1. Possible combinations of dbh versus height. r 22 r2\ + 2r 22 2r 2, + r22 I } rn | | | r\\ 2r I, + r l 2 rn + 2r l 2 r, 2 3 3 Dbh Simulated plots were generated under ten descriptions: 1. Very diverse: (dbh, height) ~ U[(r//, r!2), (r2i, r22)] 2. Moderately diverse and nominally large, with a range two-thirds that of case one and centred at the upper end: (dbh, height) ~ \J[((2ru+r12)/3, r]2), ((2r2!+r22)/3, r22)] 3. Moderately diverse and nominally small, with a range two-thirds that of case one and centred at the lower end: (dbh, height) ~ U[(r//, (2r/2+r//)/3), (r2i, (2r22+r2i)/3)] 4. Not very diverse and nominally large, with a range one-third that of case one and centred at the upper end: (dbh, height) ~ \J[((2r12+ru)/3, r / 2), ((2r22+r2i)/3, r22)] 5. Not very diverse and nominally mid-range, with a range one-third that of case one and centred at the mean of case 1: (dbh, height) ~ U[((2r//+A-/2)/3, (2r/2+r//)/3), ((2r2;+r22)/3, (2r22+r2/)/3)] 6. Not very diverse and nominally small, with a range one-third that of case one and centred at the lower end: (dbh, height) ~ U[(o/, (2rn+rI2)/3), (r2i, (2r2/+r22)/3)] 57 7. Moderately diverse and normally distributed, with a mean equal to that of case one and a standard deviation of 20% of the mean: (dbh, height) ~ N[((rn+r12)/2, (r„+r / 2)/10), ((r2/+r22)/2, (r2/+r22)/10)] 8. Not very diverse and normally distributed, with a mean equal to that of case one and a standard deviation of 10% of the mean: (dbh, height) ~ N[((rn+r/2)/2, (r//+r72)/20), ((r2/+r22)/2, (r2/+r22)/20)] 9. Moderately diverse and bimodally distributed, with half the basal area distributed as in case four and half the points distributed as in case six 10. Not very diverse and maximally bimodal, with half the basal area at (rn, r2i), and half the basal area at (r!2, r22). Theoretically, all of the above simulated plots are possible; however, there are many values that are biologically impossible. To avoid unreasonable dbh-height combinations, and to ensure that the marginal distributions of dbh and height were uniform (where applicable), restrictions were applied. For instance, in simulated case one, univariate uniform distributions of dbh and height were desired. In order to accomplish this without producing unreasonable (dbh, height) pairs, the uniform random variables for dbh and height were generated in twelve very small ranges. This resulted in dbh ~ U[ry/, rl2] and height ~ U[r 2/, r 2 2]. Although the joint distribution of (dbh, height) is not a bivariate uniform distribution, the resulting simulated plots were more biologically reasonable. Also, for the normal cases (seven and eight), randomly generated dbh/height ratios might be outside of a reasonable range. Thus, dbh and height values were generated with restrictions on the dbh/height ratio. 58 Figure 5.2 shows the possible values of dbh and height for each of the numbered simulated plots in graphical form. Figure 5.2. Graphical description of simulated dbh and height values. Simulated case 1 Dbh Simulated case 2 Dhh Simulated case 3 Dbh Simulated case 4 D b h Simulated case 7 Dbh Simulated case 10 Dbh Simulated case 5 Simulated case 8 Dbh Simulated case 6 Dbh Simulated case 9 Dbh The dbh and height ranges for the simulations were selected to roughly correspond with the ranges found in the MKRF (presented in Section 5.2). Trees grown in this region can be extremely large, reaching heights over 60 m and dbh over 120 cm in some species. Thus, the 59 simulated ranges were wide, with rn = 0.1 cm, rl2 = 120 cm, r2\ = 0.1 m, and r 2 2 = 60.0 m. Species distributions were not simulated in these plots; however, species distributions were evaluated for the MKRF inventory presented later (Section 5.2). For ease of presentation, the simulated cases will be referred to with plot type abbreviations (Table 5.1): Table 5.1. Plot type abbreviations for simulated cases. Simulated Plot Case No. Type Plot Description 1 High High diversity, with basal area uniformly distributed 2 Med-Lrg Medium diversity, with basal area skewed to large trees 3 Med-Sml Medium diversity, with basal area skewed to small trees 4 Low-Lrg Low diversity, with basal area skewed to large trees 5 Low-Mid Low diversity, with basal area skewed to mid-sized trees 6 Low-Sml Low diversity, with basal area skewed to small trees 7 Med-Norm Medium diversity, with basal area normally distributed 8 Low-Norm Low diversity, with basal area normally distributed 9 Med-Bim Medium diversity, with basal area bimodally distributed 10 Low-Bim Low diversity, with basal area bimodally distributed The "High diversity" case represents the baseline of maximum diversity as closely as possible, where the entire range of both dbh and height is represented equally in terms of basal area. In terms of stems per hectare, this would be similar to a reverse-J distribution. Because smaller trees represent larger numbers of trees per hectare than larger trees, they are inflated more than large trees when converting from basal area to stems per hectare, resulting in a distribution that is reverse-J shaped. Cases that are of "Medium diversity" have approximately two-thirds the range of data as the maximum diversity case. The "Low diversity" cases have very small ranges, where basal area represents one-third or less of the possible dbh and height ranges. 60 Simulated dbh and height values were generated using random number generators in SAS (SAS Institute Inc. 1988). For the uniform distributions (simulated cases one through six and nine), the SAS function RANUNI was used to generate random numbers. First, a random number was generated and scaled to conform to the dbh range for the plot type. Then, a second random number was generated and scaled to conform to the corresponding height range. This process was repeated over all dbh-height ranges within a particular plot. For example, there were twelve dbh-height ranges for simulated case 1, and four dbh-height ranges for simulated case four. Equal amounts of basal area were generated separately over all dbh-height ranges for a particular plot. The SAS function RANNOR was used to generate normally distributed numbers. For the two normally distributed cases (seven and eight), first, the dbh was generated with an underlying mean and variance. Then, the height was generated with an underlying mean and variance. The (dbh, height) pair was then tested for reasonableness by computing the dbh/height ratio. If the ratio did not fall in the range (0.8, 2.4), additional heights were generated until the ratio did fall in the range. This range was again based on the MKRF data; approximately 90% of dbh/height ratios fell in the range from 0.8 to 2.4. To test the normality of the simulated distributions, the SAS procedure PROC UNIVARIATE (SAS Institute Inc 1989) was used. Dbh and height distributions were not significantly different from the normal at a = 0.05. Because of the restrictions in generating the data, the joint distribution of dbh and height likely only approximated a bivariate normal distribution. 61 Each simulated plot consisted of 60 cells (dbh-height combinations), with each cell contributing equally to the total basal area of the plot. This ensured a uniform distribution of basal area over the plot. Each simulated plot type was generated 100 times. Since there are ten plot types, there are 60,000 simulated cells in all. This method of simulation produced tree size distributions which were univariately uniform by basal area (or normal); however, because of the restrictions made to ensure biologically reasonable tree sizes, the stands are not jointly uniform. The method of simulating normal stands was slightly less restrictive and produced stands that may be approximately jointly normal. 5.1.2 Results and Discussion The indices presented in Section 4.2 were computed with the simulated data. The results were examined and discussed using criteria 8 listed in Section 4.1. 5.1.2.1 Shannon's Index Extended to Diameter, Height, and Species The procedure for computing the Extended Shannon Index, using the post4ioc method was as follows: 1. Tree data were classed in ten-centimetre diameter classes, using the diameter outside bark at breast height (1.3 m above ground; dbh): 0 to 10 cm through 110 to 120 cm. Then, Shannon's Index was computed with /?,-, the proportion of basal area in the ith diameter class, and equation [2.1] to yield H'd. 2. Tree data were classed into five-metre height classes: 0 to 5 m through 55 to 60 m. Then, Shannon's Index was computed with pj, the proportion of basal area in the jth height class, and equation [2.1] to yield H'h. 62 3. The final index, H'd+h, was computed as the average of the diameter and height indices. The procedure for computing the Extended Shannon Index under the combination method was slightly different. The data were classed simultaneously by diameter and height, using the same classes as with the post-hoc method. Then, Shannon's Index was computed with pij, the proportion of basal area in the ith diameter, y'th height class, and summed over all classes to yield H'dh. The results of calculating the two types of Shannon's indices with the simulated data are shown in Table 5.2. The average, minimum, maximum, and standard deviation show the spread of the index over the sample plots of the same type. The range of values represents the variability in the generation of the data. The simulated stands were ranked by average index values, with the highest index receiving a rank of one and the lowest receiving a rank of ten. Index values within 0.01 received the same rank. 63 Table 5.2. Extended Shannon indices computed using simulated stand data. Standard Plot Average Minimum Maximum Deviation Rank Type Hd+h H'dh Hd+h Hdh H'd+h H'dh Hd+h Hdh Hd+h Hdh High 2.485 2.522 2.476 2.485 2.506 2.652 0.006 0.038 1 2 Med-Lrg 2.079 2.115 2.073 2.077 2.127. 2.273 0.008 0.039 2 4 Med-Sml 2.079 2.119 2.068 2.077 2.124 2.271 0.009 0.046 2 4 Low-Lrg 1.391 1.436 1.381 1.386 1.466 1.606 0.014 0.051 6 7 Low-Mid 1.392 1.438 1.385 1.386 1.448 1.570 0.014 0.051 6 7 Low-Sml 1.392 1.441 1.382 1.386 1.447 1.629 0.014 0.060 6 7 Med-Norm 1.494 2.676 1.339 2.418 1.643 2.895 0.069 0.113 5 1 Low-Norm 0.950 1.812 0.796 1.542 1.079 2.067 0.068 0.121 9 6 Med-Bim 2.082 2.121 2.073 2.077 2.112 2.226 0.010 0.045 2 3 Low-Bim 0.693 0.693 0.693 0.693 0.693 0.693 0.000 0.000 10 10 The indices associated with the combination index tended to be larger than that of the post-hoc index, as did the standard deviations. This was not unexpected, as the simultaneous relegation of data to dbh and height classes necessarily produced more dbh-height combinations, and hence produced a larger index. Both types of indices were ranked in the same manner, with two exceptions. Whereas the post-hoc index (H'd + h) ranked the medium diversity normal stand (Med-Norm) with the three low-diversity stands (Low-Lrg, Low-Med, Low-Sml), the combination index (H'dh) ranked Med-Norm higher than that of the high diversity stand (High). The low-diversity normal stand (Low-Norm) also ranked higher when the combination method was used. The normal stands may have ranked higher under the combination method because they were generated in such a way as to be more diverse in a bivariate sense. That is, under the normal simulation method, stands had larger bivariate ranges than those that were generated as uniform distributions, and this could have resulted in illogical orderings. 64 With respect to the criteria listed in Section 4.1, both methods of computing the Extended Shannon Index produced indices that possess the characteristics listed as 8a, 8b, and 8c. The indices are tree size insensitive, are correlated with the range of diameter and height values, and rank bimodally distributed stands in a logical way. The post-hoc method may have performed better than the combination method in that it produced a more logical ordering of normally distributed stands; however, this may have been a consequence of the method of random number generation. 5.1.2.2 Structural Index based on Curve Fitting The Weibull function was first fit univariately and separately to the cumulative distribution functions of dbh and height by plot. The regressions were performed using the SAS procedure PROC NLIN and the Marquardt method (SAS Institute Inc. 1989). The fit statistics of the equations were first examined to ensure that the Weibull provided a good relationship. As shown in Table 5.3, the RMSE was low (below 0.05) for every case except the medium-diversity bimodal (Med-Bim), where it was nearly twice as high. The poor fit of the Weibull function to a bimodal stand was expected, since it has been shown that the Weibull function has only one maximum (Bailey and Dell 1973). The perfect fit (RMSE equal to zero) of the extremely bimodal case was also expected, since the two-parameter regression was fit with only two unique points, which reduced the regression to the solution of two equation in two unknowns. 65 Table 5.3. Resulting coefficients and RMSE for the univariate Weibull fitted to dbh using simulated stands. Standard Plot Average Average Minimum Maximum Deviation Rank Type RMSE b, Cl b, Cl bi Cl bi Cl bi Cl High 0.041 70.43 1.576 68.86 1.515 71.26 1.663 0.475 0.030 3 1 Med-Lrg 0.030 89.15 3.338 87.64 3.189 90.46 3.488 0.451 0.064 2 5 Med-Sml 0.043 46.49 1.565 45.32 1.430 47.67 1.671 0.487 0.050 8 8 Low-Lrg 0.034 104.67 8.589 103.25 7.698 105.42 9.512 0.412 0.348 1 2 Low-Mid 0.032 64.52 5.051 63.35 4.574 65.48 5.604 0.450 0.188 5 4 Low-Sml 0.045 23.52 1.581 22.30 1.342 24.48 1.808 0.448 0.094 9 6 Med-Norm 0.027 64.05 5.803 60.53 4.116 67.76 7.508 1.608 0.696 6 3 Low-Norm 0.030 62.08 11.603 59.85 9.255 63.87 15.136 0.841 1.324 7 1 Med-Bim 0.079 67.76 0.977 66.13 0.935 69.09 1.027 0.610 0.020 4 9 Low-Bim 0.000 0.23 0.436 0.23 0.436 0.23 0.436 0.000 0.000 10 10 The resulting average coefficients bj and cj were then examined (Table 5.3). The value of the coefficient bi was roughly correlated to the centre of the distribution. The low-diversity large tree (Low-Lrg), medium-diversity large tree (Med-Lrg), and very diverse (High) cases had the largest values of bj. However, the lowest value of b} Occurred for the low-diversity bimodal distribution (Low-Bim), whereas the distribution centred at the smallest dbh value was the low-diversity small tree case (Low-Sml). There was some overlap in the range of values of bj, as the most diverse case (High) and the medium-diversity bimodal case (Med-Bim) produced similar values of b/, as did the low-diversity mid-sized tree stand (Low-Mid), and the two normally distributed stands (Med-Norm and Low-Norm). The value of the coefficient cj was between 1.58 and 3.4 for the more diverse stands (High, Med-Lrg, and Med-Sml). For the most part, a much larger ci resulted for the 66 nominally larger and less diverse stands (Low-Lrg and Low-Norm), and a smaller c} resulted for the nominally smaller and less diverse stand (Low-Sml). However, the coefficient c} was very small for the two bimodal stands (Med-Bim and Low-Bim). There was considerable overlap in the values of c}, with three groups (High, Med-Low, and Low-Sml; Low-Mid and Med-Norm; and Low-Lrg and Low-Norm) that produced similar values. The resulting average coefficients, b2 and c2 for height are listed in Table 5.4, along with the average RMSE. The RMSE and resulting coefficient b2 tracked with the results for dbh, as did the cases with overlapping coefficients. The Low-Lrg, Med-Lrg, and High cases had the largest values of b2\ the lowest values of b2 occurred in the Low-Bim case! The value of the coefficient c2 was between 1.56 and 3.34 for the more diverse stands (High, Med-Lrg, and Med-Sml); however, the value of c2 also fell in this range for the Low-Sml case. As in the univariate dbh case, for the most part, a larger c2 resulted for the nominally larger and less diverse stands (Low-Lrg and Low-Norm), and the coefficient c2 was very small for the bimodal cases. Table 5.4. Resulting coefficients and RMSE for the univariate Weibull fitted to height using simulated stands. Standard Plot Average Average Minimum Maximum Deviation Rank Type RMSE b2 bi Cl b2 Cl bi Cl b2 Cl High 0.041 35.20 1.564 34.61 1.481 35.72 1.648 0.229 0.031 3 8 Med-Lrg 0.030 44.56 3.339 43.98 3.068 45.10 3.489 0.218 0.071 2 5 Med-Sml 0.042 23.17 1.571 22.71 1.437 23.62 1.733 0.213 0.048 8 7 Low-Lrg 0.035 52.34 8.494 51.70 7.700 52.91 9.234 0.238 0.319 1 2 Low-Mid 0.033 32.26 5.087 31.65 4.585 32.82 5.681 0.221 0.208 6 4 Low-Sml 0.044 11.77 1.596 11.27 1.391 12.41 1.813 0.234 0.072 9 6 Med-Norm 0.029 34.29 7.425 32.67 5.612 35.79 9.681 0.713 0.860 4 3 Low-Norm 0.032 31.51 12.970 30.74 10.19 32.46 19.616 0.342 1.655 7 1 Med-Bim 0.079 33.83 0.982 32.96 0.934 34.89 1.040 0.428 0.024 5 9 Low-Bim 0.000 0.21 0.490 0.21 0.490 0.21 0.490 0.000 0.000 10 10 67 The investigation of the univariate Weibull distributions did not lead to a useful index of diversity. No coefficient was found to possess the desirable characteristics presented as item 8 (Section 4.1), in that the coefficients produced were tree size sensitive and did not differentiate between stands of low, medium, and high diversity. The results of bivariate Weibull curve fitting follow in Tables 5.5 to 5.7. The bivariate Weibull models were fit to the joint cumulative distribution function of dbh and height. As in the univariate case, the regressions were performed using the SAS procedure PROC NLFN and the Marquardt method (SAS Institute Inc. 1989). Starting values for bj, b2, c/, and c2 were obtained from the univariate regressions. The RMSE was first analyzed to ensure that the Weibull provided a good fit to the data. As shown in Table 5.5, RMSE values were low (less than 0.05) for all but the bimodal stands. The upper limiting value for the coefficient aj (1.00) resulted for all stands except those based on the normal distribution (Table 5.5). Convergence to this maximum value may have been a result of lack of fit. Because data were generated under univariate uniform conditions, but not uniform bivariate conditions, the bivariate Weibull function may not have been appropriate for the uniformly generated stands. 68 Table 5.5. RMSE and resulting values of coefficient a3 for the bivariate Weibull fitted to simulated stands. Plot Standard Type RMSE Average Minimum Maximum Deviation High 0.045 1.000 1.000 1.000 0.000 Med-Lrg 0.047 1.000 1.000 1.000 0.000 Med-Sml 0.046 1.000 1.000 1.000 0.000 Low-Lrg 0.053 1.000 1.000 1.000 0.000 Low-Mid 0.049 1.000 1.000 1.000 0.000 Low-Sml 0.043 1.000 1.000 1.000 0.000 Med-Norm 0.029 0.923 -0.704 1.000 0.217 Low-Norm 0.029 0.645 -0.237 1.000 0.375 Med-Bim 0.059 1.000 1.000 1.000 0.000 Low-Bim 0.076 1.000 1.000 1.000 0.000 As in the univariate case, values of bi and b2 were roughly correlated to size (Table 5.6); the simulated stands with the largest average tree size (Med-Lrg and Low-Lrg) had the largest values for bj and b2. As in the univariate results, the coefficient ranges for the High, Med-Norm, Low-Norm, Low-Mid, and Med-Bim cases overlapped. Table 5.6. Resulting coefficients b\ and b2 for the bivariate Weibull fitted to simulated stands. Standard Plot Average Minimum Maximum Deviation Rank Type b, b2 b, b2 h b2 bi bi b. b2 High 61.52 30.78 58.08 29.71 63.51 32.08 0.97 0.49 6 6 Med-Lrg 84.10 41.90 81.39 40.46 86.87 43.41 1.37 0.67 2 2 Med-Sml 40.87 20.16 38.86 18.99 43.01 21.10 0.92 0.42 8 8 Low-Lrg 102.46 51.16 99.27 49.57 104.94 52.48 1.51 0.72 1 1 Low-Mid 62.20 31.08 58.69 29.51 64.97 32.49 1.38 0.69 5 5 Low-Sml 20.85 10.39 17.91 9.40 23.42 11.61 1.16 0.52 9 9 Med-Norm 64.00 33.99 60.21 31.80 67.86 35.81 1.78 0.89 3 . 3 Low-Norm 62.25 31.52 59.56 30.73 64.26 32.51 0.96 0.40 4 4 Med-Bim 56.16 27.92 46.29 21.65 68.64 33.71 5.02 2.53 7 7. Low-Bim 0.15 0.15 0.15 0.15 0.15 0.15 0.00 0.00 10 10 69 Values of the coefficients c\ and c2 are shown in Table 5.7. The more diverse cases had average values of c ; and c2 in the range 1.3 to 3.2. However, the Low-Sml case also had values of c/ and c2 in this range. As in the univariate case, values of ct and c2 were affected by the range and centre of the distribution, with higher values of c/ and c2 associated with larger ranges and mean values. The bivariate values of a and c2 also were prone to large ranges, and many of the values of these coefficients overlapped from case to case. Table 5.7. Resulting coefficients cy and c2 for the bivariate Weibull fitted to simulated stands. Standard Plot Average Minimum Maximum Deviation Rank Type Cl c2 Cl Cl Cl Cl Cl Cl Cl Cl High 1.371 1.356 1.308 1.254 1.483 1.452 0.038 0.041 6 8 Med-Lrg 3.181 3.214 2.923 2.880 3.541 3.585 0.139 0.154 5 5 Med-Sml 1.341 1.366 1.207 1.206 1.494 1.512 0.059 0.061 8 7 Low-Lrg 8.462 8.472 7.016 6.711 10.142 10.621 0.796 0.888 2 2 Low-Mid 4.900 4.989 4.054 3.879 5.948 6.350 0.446 0.506 4 4 Low-Sml 1.350 1.378 1.055 1.127 1.639 1.767 0.115 0.129 7 6 Med-Norm 5.999 7.285 4.036 5.298 8.994 10.855 0.795 1.079 3 3 Low-Norm 12.235 12.614 9.126 9.755 19.431 19.850 1.912 1.782 1 1 Med-Bim 0.905 0.912 0.791 0.790 1.073 1.071 0.060 0.065 9 9 Low-Bim 0.436 0.490 0.436 0.490 0.436 0.490 0.000 0.000 10 10 The investigation of the bivariate Weibull distribution did not lead to a useful index of diversity. No coefficient was found to possess the desirable characteristics presented as item 8 (Section 4.1), in that the coefficients produced were tree size sensitive and did not differentiate between stands of low, medium, and high diversity. 70 5.1.2.3 Structural Index based on R 2 (STRI) Both the univariate and bivariate forms of the STRI were computed with the maximum and minimum dbh and height values used to generate the simulated stands. Therefore, the theoretical distributions used in calculating the univariate STRI were: F (dbh) = d b k ° A , 0.1 < dbh < 120.0 120.0-0.1 [5.58] F2 (height) = h e i g k t 0 , 1 , 0.1 < height < 60.0 2 60.0-0.1 5 The theoretical distribution used in calculating the bivariate STRI was: G(dbh,height) — dbh-OA 120.0-0.1 x height -0.1 60.0-0.1 0.1 <dbh< 120,0 0.1 < height < 60.0 [5.59] [5.60] Results of the STRI methods follow in Table 5.8. Index values within 0.01 received the same rank. Table 5.8. Univariate (STRId+h) and bivariate (STRIdh) STRI using simulated stands. Standard Plot Average Minimum Maximum Deviation Rank Type STRI d + h STRIdh STRId+h STRIdh STRI d + h STRIdh STRI d + h STRIdh STRI d + h STRId h High 0.999 0.927 0.999 0.918 0.999 0.935 0.000 0.003 1 1 Med-Lrg 0.895 0.882 0.881 0.866 0.903 0.896 0.004 0.006 3 3 Med-Sml 0.877 0.562 0.871 0.547 0.883 0.579 0.003 0.006 5 8 Low-Lrg 0.582 0.121 0.529 -0.018 0.640 0.259 0.022 0.054 9 10 Low-Mid 0.892 0.778 0.881 0.761 0.905 0.794 0.005 0.005 3 5 Low-Sml 0.548 0.168 0.539 0.159 0.561 0.179 0.005 0.004 10 9 Med-Norm 0.877 0.794 0.855 0.726 0.899 0.830 0.009 0.017 5 4 Low-Norm 0.821 0.685 0.792 0.638 0.846 0.721 0.009 0.020 7 6 Med-Bim 0.972 0.894 0.967 0.874 0.976 0.908 0.002 0.005 2 2 Low-Bim 0.800 0.571 0.800 0.571 0.800 0.571 0.000 0.000 8 7 71 The very diverse and medium-diverse bimodal stands (High and M e d - B i m , respectively) ranked the highest under both the univariate and bivariate versions of this index. The lowest rank under both types of STRI occurred in the low diversity large and small tree stands (Low-Lrg and Low-Sml , respectively). A s expected, higher values of STRI appeared to be correlated with distributions centred at the middle of the range. There was much overlap in STRLj+n values for most of the middle-centred stands (Med-Sml, Med-Lrg , L o w - M i d , and Low-Norm); however, the STRLj h index values overlapped less. The Med-Lrg and Med-Sml cases produced very different STRLjh index values, whereas the STRLj+h values for these two stands were very close. Due to the nature of the R 2 calculation, the middle-centred distributions, such as L o w -M i d , Med-Norm, and Low-Norm, w i l l always receive higher STRI values than skewed or bimodal distributions. Both types of STRI indices are not desirable in that they are tree size sensitive, are not well-correlated to range, and do not rank bimodal stands logically. Thus, the STRI does not possess the desirable characteristic presented as item 8 (Section 4.1). 5.1.2.4 Structural Index based on Variance Comparisons (STVI) Both the univariate and bivariate S T V I were computed with the maximum and minimum dbh and height values used in generating the simulated data. Thus, the theoretical variances for univariate uniform distributions of dbh, height, and the joint distribution of (dbh, height) were: '2 _ dhhu ~ (120.0-0.1)2 712 = 1198 [5.61] = (60 .0 -0 .1 ) 2 /12 = 299 [5.62] 72 S 2 = (120.0 - 0 . 1 ) 2 x (60.0 - O . I ) 2 / 1 4 4 =358,203 [5.63] The maximum variances used in the STVI calculations were: S l m a x = ( 120 .0 -0 .1 ) 2 /4=3,594 [5.64] Slight = ( 6 0 . 0 - 0 . 1 ) 2 / 4 = 8 9 7 [5.65] S 2 a x = (120.0 -0 .1 ) 2 x (60 .0 -0 . 1 ) 2 /16=3,223,829 [5.66] Results of calculating indices with the STVI follow in Table 5.9. Index values within 0.01 received the same rank. The index values associated with the univariate method were much larger than that of the bivariate method, as were the standard deviations. Table 5.9. Univariate (STVI^ +h) and bivariate (STVIdh) STVI using simulated stands. Standard Plot Average Minimum Maximum Deviation Rank Type STVI d + h STVI d h STVI d + h STVIdh STVI d + h STVI d h STVI d + h STVIdh STVI d + h STVIdh High 0.946 0.032 0.894 0.018 1.000 0.044 0.024 0.005 1 2 Med-Lrg 0.772 0.014 0.747 0.010 0.798 0.020 0.011 0.002 2 4 Med-Sml 0.761 0.014 0.730 0.009 0.787 0.019 0.011 0.002 3 4 Low-Lrg 0.248 0.003 0.223 0.002 0.271 0.005 0.010 0.001 5 6 Low-Mid 0.252 0.004 0.226 0.002 0.279 0.005 0.011 0.001 5 6 Low-Sml 0.251 0.004 0.231 0.002 0.274 0.005 0.009 0.001 5 6 Med-Norm 0.230 0.020 0.170 0.009 0.325 0.042 0.031 0.006 8 3 Low-Norm 0.065 0.002 0.047 0.001 0.082 0.003 0.009 0.000 10 9 Med-Bim 0.620 0.047 0.600 0.033 0.646 0.065 0.009 0.008 4 1 Low-Bim 0.086 0.000 0.086 0.000 0.086 0.000 0.000 0.000 9 10 The highest values of the univariate STVI occurred for the very diverse case (High), followed by the medium-diverse uniform stands (Med-Lrg and Med-Sml). Under the 73 bivariate method, however, the highest rank occurred for the medium-diverse bimodal case (Med-Bim), followed by the very diverse case. As was the case with the Extended Shannon indices, the normally distributed cases ranked higher using the bivariate method. Under the univariate method, the indices for the medium-diverse normal stand (Med-Norm) were ranked far below the other medium-diverse stands. However, under the bivariate method, Med-Norm was ranked above the other medium-diverse stands. Both methods ranked the three low-diversity uniform stands (Low-Lrg, Low-Mid, and Low-Sml) in the same range, with the low-diversity normal stand (Low-Norm) slightly lower. Higher values of STVI were correlated with widely distributed stands. The overlap that occurred was mostly between stands of similar diversity (e.g., Low-Lrg, Low-Mid, and Low-Sml; and Med-Lrg, Med-Sml, and Med-Norm). However, the high diversity and medium-diverse bimodal cases also overlapped in range. The univariate STVI method produced indices that possess the desirable characteristics listed under item 8 (Section 4.1). The index is tree size insensitive, is correlated to continuously occupied ranges, and gives bimodal stands logical index values. The bivariate index did not produce indices that possessed these characteristics, in that the indices did not correlate well with continuously occupied ranges and made illogical rankings with respect to the high and medium diversity uniform simulations. However, this may have been a function of the methods used to generate the simulated stands. Simulated uniform stands were generated such that they followed a univariate uniform distribution, and were thus well described by the univariate STVI. However, simulated uniform stands were not generated to 74 follow a bivariate distribution. One consequence of this is that the sample variance of the stands never reached the value of the theoretical uniform variance (Table 5.10). Table 5.10. Average sample variances for simulated cases compared to the variance of a theoretical uniform distribution. Dbh Height Joint (Dbh, Height) Plot Type S2 ^ dbh Rank S2 ^ height Rank S2 ° dbh, height Rank High 1,219 3 306 3 4,815 2 Med-Lrg 549 4 137 4 2,141 4 Med-Sml 539 5 133 5 2,075 5 Low-Lrg 133 9 34 8 520 8 Low-Mid 137 7 34 7 523 7 Low-Sml 136 8 34 6 527 6 Med-Norm 143 6 26 9 2,959 3 Low-Norm 36 10 8 10 256 9 Med-Bim 1,768 2 441 2 7,052 1 Low-Bim 3,655 1 912 1 0 10 S2 1,198 299 358,203 When dbh and height are considered univariately, the high-diversity case has a variance near that of the uniform, and the two bimodal cases have variances greater than that of the uniform. However, when dbh and height are considered jointly (bivariately), all simulations have variances much less than that of the uniform. The medium-diverse bimodal case had the highest generalized variance of all the simulated cases, but this variance was not high enough to use the right-hand side of the STVI curve. Thus, the univariate STVI may be preferred over the bivariate STVI, in that it may better meet the criteria presented in Section 4.1. 75 5.1.3 Overall Discussion of Proposed Indices Computed with Simulated Data The proposed indices were evaluated with respect to the criteria presented in Section 4.1. Although the Extended Shannon indices rely on combining data into classes, they ranked the simulated stands in a logical manner and possess the desirable characteristics presented as item 8 in Section 4.1. Curve fitting techniques may provide a method to derive indices without relying on class data; however, further analysis is required with more flexible curves, such as the Beta distribution, before deciding whether this approach shows promise. The STRI was derived to measure the deviation between empirical data and the uniform distribution; however, due to the properties of R 2, this index did not discriminate well between simulated cases of varying diversity. All of the desirable characteristics presented in Section 4.1 were met by the univariate STVI. The STVI in its bivariate form did not perform as well; however, this may have been a consequence of the methods used to generate the simulated data. Further investigation with real data is necessary to evaluate this measure. The most desirable indices of those presented are the Extended Shannon indices and the STVI. The univariate STVI meets all of the criteria laid out in Section 4.1; however, it lacks a theoretical basis. While Shannon's Index has a solid theoretical background, it may be undesirable because it relies on combining data into classes. 5.2 Malcolm Knapp Research Forest Data Evaluation of the Extended Shannon indices and the STVI was also performed with data collected in a recent forest inventory of the University of British Columbia Malcolm Knapp Research Forest (MKRF). The approximately 5000 ha MKRF is located in the Coastal 76 Western Hemlock (CWH) biogeoclimatic (BEC) zone of southwestern British Columbia (BC), where the maximum tree size for most species is quite large, reaching over 1 m in dbh and over 50 m in height. The region encompasses low to middle elevations west of the Coastal Mountains and is the wettest BEC zone in BC, receiving 1000 to 4000 mm of precipitation annually (Meidinger and Pojar 1991). The climate is cool mesothermal, with a mean annual temperature of 8° C and mild winters. The most common species in the forest cover are Douglas-fir (Pseudotsuga menzeisii (Mirb.) Franco), western hemlock (Tsuga heterophylla (Raf.) Sarg.), western red cedar (Thuja plicata Donn), amabilis fir (Abies amabilis (Dougl.) Forbes), and Sitka spruce (Picea sitchensis (Bong.) Carr.). 5.2.1 Methods In 1995, sample plot data were collected using the Ministry of Forests, Vegetation Inventory Sampling Procedures (Resources Inventory Branch 1994) as the basis for plot establishment, plot size, and measurement procedures. Eighty-two clusters of plots were systematically located over the research forest using square spacing. At each grid intersection, a central plot was located, and then additional plot centres were located in each cardinal direction, 25 m from the main plot centre. If the total basal area of a stand was estimated to be less than 15 m /ha, fixed radius plots were used; otherwise, variable radius plots were used. Because plots occurred throughout the whole of the MKRF, the data represent a variety of microsites and growing conditions. For all trees above two cm dbh, the species, dbh (cm) and tree class (live or dead, standing or fallen) were obtained. For a subset of trees, the total height (m), crown class 77 (dominant, codominant, intermediate, or suppressed), and height to live crown (m) were also measured. The 82 clusters yielded 327 plots containing 1,832 live trees. Height was measured on only 763 trees (44%). Height prediction models were developed to estimate the remaining heights The major species in the MKRF are western red cedar, Douglas-fir, western hemlock, and alder (Alnus rubra Bong, and Alnus tenuifolia Nutt.). Height prediction equations were fit to these four species groups separately. Models were not fit separately for other species because of few observations. An equation was developed for all other broad-leaved species by grouping them with the two species of alder. Heights for the remaining coniferous species were obtained by using the prediction equations for species with similar growth patterns, as shown in Table 5.11. Table 5.11. Data for height-dbh equation development by species. Number Grouped With Species of Trees (If Applicable) True Fir (Abies amabilis (Dougl.) Forbes) 1 Douglas-fir Sitka Spruce (Picea sitchensis (Bong.) Carr.) 1 Douglas-fir Red Cedar (Thuja plicata Donn) 210 -Alder (Alnus rubra Bong, and A. tenuifolia Nutt.) 28 Douglas-fir (Pseudotsuga menzeisii (Mirb.) Franco) 230 -Other Broad-leaved Trees 22 Alder Western Hemlock (Tsuga heterophylla (Raf.) Sarg.) 261 Lodgepole Pine (Pinus contorta Dougl.) 8 W. Hemlock Yew (Taxus Brevifolia Nutt.) 3 Red cedar Twelve models were evaluated by species (see Kovacs and LeMay 1999a, and Kovacs and LeMay 1999b). Based on these evaluations, the final height prediction equations by 78, species were selected. Where trees were not measured for height, the following height prediction equations were used: Alder: height = 1.3 + (68.0529 + 0.019287 x TSTEMSHA - 0.77315 x TBAHA) x [l-exp((-0.00877 + 0.00000044 xTSTEMSHA- 0.0006 XxTBAHA) [5.67] x dbh O M 9 9 ] - ° M m m x T S T E M S H A + o m 5 4 x T B A H A yj Cedar and yew: height = 1.3 + (59.5789 + 12.5373 x REL_dbh - 36.98 x REL_dbh 2)x [1 - exp((0.00071 - 0.03818 x REL_dbh + 0.01551 x REL_dbh 2) ! [ 5 - 6 8 l x dbh 1 6 5 7 7 ~ 2 0 2 1 9 9 x R E L - d h h + l 4 1 5 9 9 xREL-dhhl)yj Douglas-fir, true fir, and Sitka spruce: height = 1.3 + (15.759 + 767.419 x REL_dbh - 699.761 x REL_dbh 2) x [1 - exp((0.00025 - 0.00332 x REL_dbh - 0.0118 x REL_dbh 2) [5-69] X dbh 2 0 7 5 4 " 3- 3 8 7 3 x S £ z - -dbh+2.1785 xREL _dbh1 Hemlock and lodgepole pine: height = 1.3 + (-1.25839 +163.146 x REL_dbh -126.716 x REL_dbh 2 )x [1 - exp((0.00021 - 0.06091 x REL_DBH + 0.05552 x REL_dbh 2) x t 5 - 7 ° ] dbh 2 3 4 1 3 5 - 4 1 ' 2 3 9 x R £ Z . _ d f o A +3.31119 xREL_dbh 2 All other broad-leaved species: height = 1.3 + (1.07897 -1.69747 x REL_dbh + 0.772644 x REL_dbh 2)x (0.868644 +0.276184 xREL_DBH -0.198905 xREL_dbh2)x [5.71] dbh 1 75194 ~1 •17265 x REL - dbh + 1 3 2 6 8 x REL - dhh 2))] Since the analysis and presentation of all 82 cluster plots would be cumbersome, nine plots were selected to match as closely as possible the dbh and height distributions in simulated cases one through nine; a plot similar to plot ten was not found in the MKRF data. Descriptive statistics for the selected plots are given in Table 5.12. 79 Table 5.12. Descriptive statistics of selected MKRF plots. Plot MKRF Species Average Minimum Maximum Std. Dev. Type Plot# Count Dbh Height Dbh Height Dbh Height Dbh Height High . 57 3 52.80 34.67 8.4 5.9 101.0 51.5 25.32 11.70 Med-Lrg 09 3 84.80 47.15 31.7 26.8 180.5 66.2 42.99 11.13 Med-Sml 46 3 34.55 25.90 6.5 4.6 61.9 42.9 14.04 8.82 Low-Lrg 58 3 58.24 32.00 26.3 21.0 99.3 46.9 17.47 . 7.33 Low-Mid 49 3 56.17 28.58 4.6 4.2 93.5 40.8 19.66 6.97 Low-Sml 52 6 12.47 10.07 2.5 3.2 24.7 15.9 6.31 3.79 Med-Norm 15 3 44.84 29.56 17.8 13.0 73.6 42.2 12.23 6.07 Low-Norm 77 3 33.74 23.46 21.3 16.5 44.9 30.2 7.27 4.47 Med-Bim 45 3 95.42 42.50 3.7 1.3 174.4 68.3 54.97 16.74 Graphs of the distribution of basal area by plot and species are found in Appendix I. Plots were chosen on the basis of dbh distribution and height distribution for all species combined, and thus, the distributions by species may differ substantially from the chosen plot type. 5.2.2 Results and Discussion The Extended Shannon Index and the STVI were computed with the MKRF data by five-point cluster plot, using all live trees. The indices were examined using criteria 8 listed in Section 4.1. 5.2.2.1 Shannon's Index Extended to Diameter, Height, and Species The steps in computing the post-hoc Extended Shannon Index with the MKRF data were modified from Section 5.1.2.1: 1. Tree data were classed in ten-centimetre diameter classes, using the diameter outside bark at breast height (1.3 m above ground; dbh): 0 to 10 cm through 140 to 150 cm, 80 and 150+ cm. Then, Shannon's Index was computed with p„ the proportion of basal area in the ith diameter class, and equation [2.1] to yield H'd. 2. Tree data were classed into five-metre height classes: 0 to 5 m through 55 to 60 m, and 60+ m. Then, Shannon's Index was computed with pj, the proportion of basal area in the jth height class, and equation [2.1] to yield H'h. 3. The data were combined by species, and Shannon's Index was computed with p^ the proportion of basal area in the kth species, and equation [2.1] to yield H's. 4. The final index was computed as the average of the diameter, height, and species indices and labelled H'd+h+s. The procedure for computing the Extended Shannon index under the combination method was also modified to include species. The data were classed simultaneously by species, diameter, and height, using the same classes as with the post-hoc method. Shannon's Index was then computed with pijk, the proportion of basal area in the ith diameter, y'th height, and kih species class, and summed over all classes to yield Hdhs. The results of computing both types of Extended Shannon index follow in Table 5.13, along with Shannon's Index computed individually for the diameter, height, and species distributions. As was the case with the simulated data, the indices associated with the combination index tended to be larger than that of the post-hoc index. This was not unexpected, as the simultaneous relegation of data to species, dbh, and height classes necessarily produced more combinations, and hence produced a larger index. 81 Table 5.13. Extended Shannon indices computed with MKRF data. Plot MKRF Rank Type Plot* K K K Hd+h+s Hdhs Hd+h+s Hdhs High 57 2.186 2.091 0.979 1.752 2.860 1 3 Med-Lrg 09 2.173 1.935 1.067 1.725 3.103 2 1 Med-Sml 46 1.748 1.915 0.997 1.553 2.781 4 6 Low-Lrg 58 1.732 1.736 0.791 1.420 2.651 5 7 Low-Mid 49 1.859 1.480 0.876 1.405 2.811 6 5 Low-Sml 52 0.958 1.117 1.659 1.245 2.399 8 8 Med-Norm 15 1.616 1.562 1.011 1.396 2.974 7 2 Low-Norm 77 0.995 1.149 0.981 1.042 2.062 9 9 Med-Bim 45 2.141 1.851 1.098 1.697 2.827 3 4 In general, both types of index ranked the plots in the same manner. Exceptions occurred for the medium-diversity small tree plot, the low-diversity large tree plot, and the medium-diversity normal plot. Although the plot chosen to represent high diversity was not the highest ranked under the combination method, the index values for this plot type were close to that of the highest ranked plot (2.86 versus 3.10). In general, the indices ranked the high and medium-diversity stands above the low-diversity stands, with one exception: under the combination method, the stand with low-diversity and mid-sized trees (Low-Mid) was ranked as more diverse than the stand with medium-diversity and small trees. The low-diversity normal (Low-Norm) and the low-diversity small tree (Low-Sml) stands always ranked last; however, this was no doubt due to their very limited ranges in comparison to the other low-diversity stands. Although the dbh and height distributions were similar for the two stands, Low-Sml ranked higher due to its higher species diversity. Al l of the selected stands except Low-Sml have only three species and H's ranged from 0.791 to 82 1.011. Low-Sml, however, has six species and H's was 1.659. Species diversity doubtlessly played a larger role in the final value of the index for this stand than in the other selected stands. The computation of the Extended Shannon indices with real data brought to light some compelling practical problems. Whereas dbh and height for the high-diversity stand were well distributed over most of the possible range, other stands (e.g., Med-Lrg and Med-Bim) had taller and larger dbh trees. Thus, there were more classes represented in some of the less equally distributed stands, and this may have had a compensating effect. Had the classes been chosen such that there was one class for dbh greater than 120 cm, indices for the plots labelled as Med-Lrg and Med-Bim would have been lower. As was the case for the simulated data, the post-hoc method produced more logical rankings of plots. For the simulated data, this may have been a function of how the random plots were generated. For the MKRF data, this may be a result of how the plots were chosen. Both the simulated plot generation and the MKRF plot selection were performed by dbh and height dimension separately, and the post-hoc method relies on the dbh and height dimensions separately. The combination method depends on the three dimensional dbh-height-species distribution, and this multi-dimensional distribution was not controlled in the simulations or selected for in the MKRF data. 83 5.2.2.2 Structural Index based on Variance Comparisons (STVI) The first step in computing both the univariate and bivariate STVI was to select the maximum and minimum dbh and height values. One option was to use the actual maxima and minima present in the MKRF. However, in the forests of Coastal British Columbia, the maximum dbh may be 2 metres or more. The practical use question is whether to use the maxima "commonly found" or to use the actual maxima found in the sample data. When computing the STRI, the actual minima and maxima must be used because failure to do so could produce unreasonable indices. For example, if a sample had dbh and height values present over a larger range than had been selected for the theoretical range, the computed STRI could be negative, even if the slope of the empirical data matched that of the theoretical uniform distribution. On the other hand, when computing the STVI, negative values could never result. However, if a sample was perfectly uniform with values outside of the range chosen for computing the STVI, the resulting index would be less than one, indicating a less diverse stand. Considering these arguments, the maxima and minima were chosen such that the trees that fell outside of the range were rare: dbh in the range 2 to 150 cm, height in the range 1.3 to 60 m. Thus, the theoretical variances for univariate uniform distributions of dbh, height, and the joint distribution of (dbh, height) used in the calculation of the STVI were: i2 _ dbhu -(150.0-2.0)2 712=1825 [5.72] S .2 heightu = (60.0 -1.3)2 712 =287 [5.73] (150.0-2.0)2 x(60.0-1.3)2 7144=524,128 [5.74] The maximum variances used in the STVI calculations were: 84 dbh, [max = (150.0 -2.0) 2 74=5,476 [5.75] height\ max = (60.0 -1.3) 2 74=861 [5.76] 5, 2 (150.0- 2.0)2 x (60.0 - 2.0)2 716 =4,717,150 [5.77] max Results of calculating indices with the STVI follow in Table 5.14. As was the case in analyzing the simulated data, the bivariate STVI values were smaller than that of the univariate. The maximum value of STVI is S. Hence, the maximum value for all the selected plots except Low-Sml is 3.00; the maximum possible value for Low-Sml is 6.00. Table 5.14. Univariate (STVIdbh, STVI h e igh,, and STVI<,+h) and bivariate (STVId h) STVI using MKRF data. Plot MKRF Rank Type' Plot# STVI d b h STVI h e i K h t STVI d + h STVI d h STVI d + h STVI d h High 57 2.061 1.126 1.594 0.085 1 3 Med-Lrg 09 1.090 1.984 1.537 0.498 3 1 Med-Sml 46 0.733 0.391 0.562 0.009 7 7 Low-Lrg 58 0.582 0.678 0.630 0.031 6 5 Low-Mid 49 0.890 0.551 0.721 0.042 5 4 Low-Sml 52 0.321 0.295 0.308 0.000 9 9 Med-Norm 15 1.028 0.697 0.863 0.013 4 6 Low-Norm 77 0.482 0.172 0.327 0.001 8 8 Med-Bim 45 1.298 1.853 1.575 0.215 2 2 The highest values of the univariate STVI occurred for the very diverse case (High), followed by the medium-diverse bimodal and medium-diverse large tree stands (Med-Bim and Med-Lrg). Under the bivariate method, however, the highest values occurred for the medium-diverse large tree stand (Med-Lrg), followed by the medium-diverse bimodal 85 stand(Med-Bim), and then the very diverse case (High). Both methods ranked the low-diversity large tree and mid-sized tree stands (Low-Lrg and Low-Mid) higher than that of the medium-diverse small tree stand (Med-Sml). Under the bivariate method, the medium-diverse normal tree stand ranked lower than that of Low-Lrg and Low-Mid. Higher values of STVI were roughly correlated with more widely distributed stands. The univariate STVI may possess the desirable characteristics listed under item 8 (Section 4.1). With the exception of the medium-diverse small tree stand, the index appeared tree size insensitive, was correlated to continuously occupied ranges, and gave bimodal stands logical index values. The bivariate index does not possess these characteristics, in that it did not correlate well with continuously occupied ranges and made illogical rankings with respect to many of the stands. However, this may be a function of the methods used to choose the representative stands. Stands were chosen by looking at the univariate distributions of dbh and height, and were thus well described by the univariate STVI. Stands were not, however, chosen by looking at the bivariate distributions of dbh and height by species. One consequence of this is that the sample variance of the stands never reached the value of the theoretical uniform variance (Table 5.15). 86 Table 5.15. Sample variances and average sample variances by species for chosen MKRF stands compared to the variance of a theoretical uniform distribution. Plot Type MKRF Plot# Species Count Dbh S2 J dbh ~ 2 *dhhk Heij ° height jht ^height, Joint (Dbh, Height) s2 si High 57 3 641 351 137 113 17,667 78 Med-Lrg 09 3 1,848 1,280 124 50 95,103 196 Med-Sml 46 3 197 104 78 32 3,047 26 Low-Lrg 58 3 305 293 54 38 10,680 47 Low-Mid 49 3 386 154 49 39 10,061 55 Low-Sml 52 6 40 47 14 8 225 4 Med-Norm • 15 3 149 197 37 48 1,700 30 Low-Norm 77 3 53 45 20 20 725 9 Med-Bim 45 3 3,022 991 280 86 271,129 126. 1,825 287 524,128 When dbh and height are considered univariately for all species combined, the high-diversity (High) and medium-diverse large (Med-Lrg) stands had variances near that of the theoretical uniform distribution, and the medium-diverse bimodal (Med-Bim) stand had a variance greater than that of the uniform. However, when dbh and height were considered univariately by species, all stands except Med-Lrg had average variances much less than that of the uniform. The influence of the distribution of dbh and height by species is apparent, in that the stand with the highest variance by species (Med-Lrg) has a high value for STVI d + h , and by far the highest value of STVIdh. When dbh and height were considered jointly (bivariately), all stands had variances much less than that of the uniform. The Med-Lrg stand had the highest average generalized variance by species of all the selected stands, but this variance was not high enough to use the right-hand side of the STVI curve. Nonetheless, the STVI for Med-Lrg was by far the 87 highest of the selected stands. The univariate STVI may be preferred over the bivariate STVI, in that it may better meet the criteria presented in Section 4.1, and provides a more logical ranking of the stands. 5.2.3 Discussion of Extended Shannon Index and Univariate STVI Computed with MKRF Data The Extended Shannon indices and the univariate STVI gave adequate results with respect to the criteria presented in Section 4.1. The indices were fairly tree size insensitive, gave logical values for bimodal stands, and were roughly correlated to the range of tree sizes. Although some illogical rankings occurred, these rankings may have resulted from the way plots were selected. Practical problems may occur with either type of index. The computation of the Extended Shannon indices requires the selection of dbh and height class boundaries. The use of classes leads to loss of information about the distribution of values within classes, and the method of class selection may have a large impact on the final values of the indices. The computation of STVI requires the selection of maximum and minimum values of dbh and height; however, this selection may be done by choosing ranges for dbh and height such that only rare values occur outside the range. While the STVI may be favoured in that there is no necessity to choose classes, an index based on Shannon's Index may be preferable because it has a solid basis in theory and a history of use in biology and forestry. 88 Chapter 6 Conclusion Tree diversity plays an important role in forested ecosystems (Lande et al. 1999). Previous research has shown that the diversity of tree species and size is correlated with the presence of avian and insect species (e.g., see Franzreb 1978, MacArthur and MacArthur 1961, and Willson 1974), which are in turn correlated with a variety of other indicators of overall ecosystem function (Recher et al. 1996). Stand structure can also impact stand growth and yield, as structure is an important factor in predicting the growth of trees (Oliver and Larson 1996). Previous attempts have been made to incorporate tree size diversity into existing diversity indices; however, no commonly used index is desirable from both a practical and biological standpoint. The objective of this thesis was to develop a quantitative measure of diversity that included both the diversity of tree species and tree size. The index should measure richness and evenness using basal area and commonly measured tree attributes. The index should not rely on combining data into classes, and should account for vertical and horizontal diversity equally. The size of samples should not affect the index, and the index should be insensitive to tree size and correlated to the range of diameters and heights. Four alternative measures were proposed. The Extended Shannon indices were derived based on Shannon's Index, measuring the diversity of tree heights, tree 89 diameters, and tree species simultaneously. Curve fitting techniques were used to explore the use of the regression coefficients from fitted diameter and height distributions based on the Weibull function as indices of diversity. The STRI was derived to measure the diversity of diameter and height distributions by measuring the fit of the distribution to that of a uniform distribution. The variance of the i distribution of diameters and heights was used to measure tree size diversity in the STVI. In order to evaluate the proposed measures, a baseline of maximum diversity was chosen. This baseline was a uniform distribution of dbh and height with respect to basal area. Many patterns have been assumed to describe the distribution of tree sizes and ages. All-aged multi-cohort stands can be described with a reverse-J shape distribution with respect to stems per hectare (Oliver and Larson 1996), which would correspond to a uniform distribution of basal area per hectare. However, several studies have shown that while the stand in total may be reverse-J, the individual species may be distributed differently, e.g., normally (Oliver and Larson 1996). Therefore, the distribution of maximum diversity, may be a much more complex mixture of several distributions. However, allowing for this complexity would require a much more complex analysis, and would result in indices with: difficult interpretations. In order to create an index which may be useful to forest managers, a simple approach seemed more appropriate. 90 All four proposed measures are based on the distribution of basal area per hectare of a stand. Although the use of stems per hectare is common in diversity indices, basal area was chosen, as it has been found to be a better representation of stand occupation, and may provide a more accurate assessment of diversity than stems (Solomon and Gove 1999). If the analysis were repeated with distributions based on stems per hectare, the indices would be weighted toward smaller trees, which represent a greater proportion of the stems per hectare in a plot. The indices were evaluated with simulated data and inventory data. Simulated data were created to represent, as closely as possible, several different levels of structural diversity. To keep this phase of the analysis simple, species distributions were not simulated. Because simulated tree sizes were limited to biologically possible values, the baseline maximum case, with dbh and height distributed as a bivariate uniform, could not be simulated. Resulting simulated distributions were controlled in a univariate sense, but not in a bivariate sense. Nine plots were selected from the 82 MKRF inventory plots to test the indices in a less controlled sense. Because of the volume of data, these plots were chosen to match the simulated cases with respect to the univariate diameter and height distributions of all species combined. Thus, the selected plots may have had radically different bivariate distributions by species. The selection criteria for the MKRF plots and the simulation definitions may have impacted the results of testing the indices. 91 The four proposed measures were initially evaluated with the simulated data. The curve fitting technique did not yield an adequate index of diversity, perhaps because of the choice of the Weibull function as the base. Further analysis with a more flexible function, such as the Beta may yield better results. Although the concept of deviation from the uniform distribution provided a good starting point for a measure of diversity, the STRI proved to be too sensitive to tree size. A more appropriate avenue for investigation may come from the area of goodness of fit testing, whereby distributions of dbh and height (univariate or bivariate) may be evaluated against that of the uniform. Again, further analysis may yield an index based on this concept. The remaining two measures gave better results, and were investigated with both the simulated data and the MKRF data. The Extended Shannon Indices were easily applied and adequately described the simulated stands and the MKRF data both in its univariate and bivariate forms. However, their computation requires that data be combined into classes. Classes may be chosen arbitrarily, and changing class boundaries or range sizes will invariably change the resulting indices. Application of the Extended Shannon Indices to data in different geographical regions may necessitate re-defining classes, making comparisons inappropriate. The idea of applying Shannon's index to a continuous distribution via an integral form may hold promise; however, one must ensure that any functional description of the distribution of dbh and height provides an adequate fit to the data. This would certainly make the application of the index more complicated, but would provide a good alternative to arbitrary classes. Because 92 Shannon's Index has been researched extensively, and the sampling properties of this index are well known, it may be worthwhile investigating this option further. The STVI was formulated to compare the variance of a sample plot to that of a distribution representing a baseline of maximum diversity. It is a very flexible index, as it may be derived assuming almost any baseline distribution. It may also be modified to fit user defined criteria. While both forms of the STVI were easily applied and performed adequately using the simulated data, only the univariate form of the STVI was adequate with the MKRF inventory data. This may have been due at least in part to the way in which plots were selected from the MKRF, as bivariate distributions were not evaluated. However, it may be inferred from the results of the MKRF, that the bivariate standard of maximal diversity may never be met with real data. No stand in the MKRF had a generalized variance close to that of the uniform distribution, even when all species are combined. Given this practical reality, the univariate STVI may be a better measure of diversity than the bivariate STVI, unless the standard for maximum diversity is modified. The STVI does not rely on combining data into arbitrary classes, making it an attractive option compared to the Extended Shannon Indices. However, applying the index to other geographical regions may involve changing the ranges needed as inputs for deriving the index. This may make comparisons between regions difficult. : Furthermore, STVI lacks theoretical background, whereas Shannon's Index benefits from a long history of use in forestry and ecology. Clearly, analysis of the sampling 93 distribution of the index and the sensitivity of the index to the input ranges is required to adequately evaluate the index. As the concept of diversity is gaining importance in forestry, a quantitative measure of diversity is needed to provide objective criteria to evaluate changes in diversity. The alternative measures presented are only a step toward deriving an index that will adequately give a quantitative measure of stand diversity incorporating structure. However, to properly evaluate the measures presented, further analysis is necessary. Before either index can be used, more information is necessary to describe the sampling properties of the indices, and describe the relationships of these indices to other stand variables. Once these properties are known, standards should be determined to make the indices adaptable to different geographic areas and standards of maximum diversity. 94 References Ci ted Aber, J. D. 1979. Foliage-height profiles and succession in northern hardwood forests. Ecology 60: 18-23. Ambuel, B., and S. A. Temple. 1983. Area-dependent changes in the bird communities and vegetation of southern Wisconsin forests. Ecology 64: 1057-1068. Bailey, R. L., and T. R. Dell. 1973. Quantifying diameter distributions with the Weibull function. For. Sci 19: 97-104. Buongiorno, J., S. Dahir, H. C. Lu, and C. R. Lin. 1994. Tree size diversity and economic returns in uneven-aged forest stands. For. Sci. 40: 83-103. Buongiorno, J., J. L. Peyron, F. Houllier, and M . Bruciamacchie. 1995. Growth and management of mixed-species, uneven-aged forests in the French Jura: implications for economic returns and tree diversity. For. Sci. 41: 397-429. Cain, M. D., and M . G. Shelton. 1995. Thirty-eight years of autogenic, woody understory dynamics in a mature, temperate pine-oak forest. Canadian Journal of Forest Resources 25: 1997-2009. Canadian Standards Association. 1996. A Sustainable Forest Management System: Specifications Document - Z809. Canadian Standards Association, Toronto, Ontario. 23 Pp. Clutter, J. L., F. C. Fortson, L. V. Peinaar, and R. L. Bailey. 1983. Timber Management -a quantitative approach. John Wiley and Sons, New York. 333 Pp. Condit, R., S. P. Hubbell, J. V. LaFrankie, R. Sukumar, N. Manokaran, R. B. Foster, and P. S. Ashton. 1996. Species-area and species-individual relationships for tropical trees: a comparison of three 50-ha plots. Journal of Ecology 84: 549-62. Corona, P., and G. Pignatti. 1996. Assessing and comparing forest plantations proximity to natural conditions. J. of Sustainable Forestry 3: 37-46. Dobson, A. 1996. Conservation and Biodiversity. Scientific American Library, New York, New York. 264 Pp. Erdelen, M . 1984. Bird communities and vegetation structure: Correlations and comparisons of simple and diversity indices. Oecologia 61: 277-284. Ferris-Kaan, R., A. J. Peace, and J. W. Humphrey. 1998. Assessing structural diversity in managed forests. Pages 331-342. Assessment of biodiversity for improved forest planning. Proceedings of the conference on assessment of biodiversity of improved forest planning ,7-11 October 1996, Monte Verita, Switzerland. Kluwer Academic 95 Publisher, Dordrecht, Netherlands. Franzreb, K. E. 1978. Tree species used by birds in logged and un-logged mixed-coniferous forests. Wilson Bull. 90: 221-238. Freemark, K. E., and H. G. Merriam. 1986. Importance of area and habitat heterogeneity to bird assemblages in temperate forest fragments. Biol. Conservation 36: 115-141. Gove, J. FL, S. E. Solomon, S. E. Fairweather, and J. P. Patil. 1991. Maximizing the diameter class diversity of uneven-aged northern hardwood stands. P. 39-45. Pages 423 in M. A. Burford, ed. Proc. 1991 Symposium on Systems Analysis in Forest Resources. USDA For. Serv. Gen. Tech. Rpt. Hafley, W. L., and H. T. Schreuder. 1976. Some non-normal bivariate distributions and their potential for forest application. Pages 104-114. International Union Forest Research Organizations. XVI World Congress Proceedings, Div. VI, Oslo, Norway. Hafley, W. L., and H. T. Schreuder. 1977. Statistical distributions for fitting diameter and height data in even-aged stands. Can. J. For. Res. 7: 481-487. Harrington, T. B., and M . B. Edwards. 1995. Structure of mixed pine and hardwood stands 12 years after various methods and intensities of site preparation in the Georgia Piedmont. Can. J. of For. Res. 26: 1490-1500. Heltshe, J. F., and D. W. Bitz. 1979. Comparing diversity measures in sampled communities. Pages 133-144 in J. F. Grassle, G. P. Patil, W. Smith, and C. Taillie, eds. Ecological diversity in theory and practice. International Co-operative Publishing House, Fairland, MD. Hill, M . O. 1973. Diversity and evenness: a unifying notation and its consequences. Ecology 54: 427-431. Hunter, M . L. 1990. Wildlife, forests, and forestry. Prentice Hall Career & Technology, Englewood Cliffs, NJ. 370 Pp. Hurlbert, S. H. 1971. The nonconcept of species diversity: a critique and alternative parameters. Ecology 52: 577-586. Hutcheson, K. 1970. A test for comparing diversities based on the Shannon formula. J. of Theo. Biology 29: 151-154. Johnson, R. A., and D. W. Wichem. 1998. Applied multivariate statistical analysis. Prentice Hall, Upper Saddle River, NJ. 816 Pp. Kempton, R. A., and R. W. M . Wedderburn. 1978. A Comparison of three measures of species diversity. Biometrics 34: 25-37. Kimmins, J. P. 1997. Biodiversity and its relationship to ecosystem health & integrity. 96 Forest Chronicle 73: 229-232. Kobayashi, S. 1981. Diversity indices: relations to sample size and spatial distribution. Japanese Journal of Ecology 31: 231-236. Kovacs, C , and V. LeMay. 1999a. Height prediction equations using diameter and stand density measures for species in coastal British Columbia. Submitted to: Forestry Chronicle . Kovacs, C , and V. LeMay. 1999b. Height-Diameter Equations for the Malcolm Knapp Research Forest. Prepared for: U.B.C. Malcolm Knapp Research Forest: 17 pp + Appendices. Lahde, E., O. Laiho, Y. Norokorpi, and T. Saksa. 1999. Stand structure as the basis of diversity index. Forest Ecology and Mgmt. 115: 213-220. LeMay, V., D. Morgan, and U. Soderberg. 1997. Biodiversity measures: What could be obtained from forest inventory data. (Unpublished). Lewis, C. E., B. F. Swindel, and G. W. Tanner. 1988. Species diversity and diversity profiles: concept, measurement, and application to timber and range management. J. of Range Mgmt. 41: 466-469. Lu, J. C , and G. K. Bhattacharyya. 1990. Some new constructions of bivariate Weibull models. Ann. Inst. Statist. Math. 42: 543-559. Ludwig, J. A., and J. F. Reynolds. 1988. Statistical ecology: a primer on methods and computing. John Wiley and Sons, New York Pp. MacArthur, R. H., and J. W. MacArthur. 1961. On bird species diversity. Ecology 42: 594-8. Magnussen, S., and T. J. Boyle. 1995. Estimating sample size for inference about the Shannon-Weaver and the Simpson indices of species diversity. Forest Ecology and Mgmt. 78:71-84. Magurran, A. E. 1988. Ecological diversity and its measurement. Princeton University Press, Princeton, NJ. 179 Pp. McMinn, J. W. 1992. Diversity of woody species 10 years after four harvesting treatments in the oak-pine type. Can. J. For. Res. 22: 1179-1183. Meidinger, D., and J. Pojar, eds. 1991. Ecosystems of British Columbia. Province of British Columbia, Victoria, BC. Moen, C. A., and R. J. Gutierrez. 1997. California spotted owl habitat selection in the central Sierra Nevada. J. of Wildlife Management 61: 1281-1287. 97 Neter, J., W. Wasseman, and M . H. Kutner. 1996. Applied Linear Statistical Models. McGraw-Hill, Boston. 1408 Pp. Niese, J. N. , and T. F. Strong. 1992. Economic and tree diversity trade-offs in managed northern hardwoods. Can. J. For. Res. 22: 1807-1813. Noss, R. F. 1990. Indicators for monitoring biodiversity: a hierarchical approach. Cons. Bio 4: 355-364. Oliver, C. D., and B. C. Larson. 1996. Forest Stand Dynamics. John Wiley and Sons, Inc., New York. 520 Pp. Onal, H. 1997. Trade-off between structural diversity and economic objectives in forest management. Amer. J. Agr. Econ. 79: 1001-1012. Orloci, L. 1970. Automatic classification of plants based on information content. Can. J. Botany 48: 793-802. Patil, G. P., and C. Taillie, eds. 1979. An overview of diversity. International Cooperative Publishing House, Fairland, MD. Patil, G. P., and C. Taillie. 1982. Diversity as a concept and its measurement. J.of the American Statistical Association 77: 548-563. Peet, R. K. 1974. The measurement of species diversity. Ann. Rev. Ecol. System 5: 285-307. Pielou, E. C. 1966. The Measurement of diversity in different types of biological collections. Journal of Theoretical Biology 13: 131-144. Pielou, E. C. 1975. Ecological Diversity. John Wiley and Sons, New York. 165 Pp. Pielou, E. C. 1977. Mathematical Ecology. John Wiley & Sons, New York. 385 Pp. Qinghong, L. 1994. A model for species diversity monitoring at the community level and its applications. Environmental Monitoring and Assessment 34: 271-287. Recher, H. F., J. D. Majer, and S. Ganesh. 1996. Eucalypts, arthropods and birds: on the relation between foliar nutrients and species richness. Forest Ecology & Management 85: 177-195. Resources Inventory Branch. 1994. Vegetation Inventory: Draft procedures for the 1994 pilot projects. Ministry of Forests, Resources Inventory Branch, Victoria, BC. 300+ pp. Plus Appendices Pp. Routledge, R. D. 1979. Diversity indices: Which ones are admissible? Journal of Theoretical Biology 76: 503-15. Routledge, R. D. 1980. Bias in estimating the diversity of large, uncensused 98 communities. Ecology 61: 276-281. SAS Institute Inc. 1988. SAS Language Guide for Personal Computers, Release 6.03 Edition. SAS Institute Inc., Cary, NC. 558 Pp SAS Institute Inc. 1989. SAS/STAT User Guide, Version 6. SAS Institute Inc., Cary, N.C. 943 Pp. Schneider, R. 1997. Ecological Diversity Monitoring Framework (draft). Biodiversity Monitoring Working Group. Unpublished paper . Schreuder, H. T., and W. L. Hafley. 1977. A useful bivariate distribution for describing stand structure of tree heights and diameters. Biometrics 33: 471-478. Shannon, C. E., and W. Weaver. 1949. The mathematical theory of communication. University of Illinois Press, Urbana, IL. 117 Pp. Silbaugh, J. M. , and D. R. Betters. 1995. Quantitative biodiversity measures applied to forest management. Environmental Reviews 3: 277-285. Simpson, E. H. 1949. The Measurement of diversity. Nature 163: 688. Smith, E. P. 1984. Non-parametric estimation of species richness. Biometrics 40: 119-129. Smith, W., and F. J. Grassle. 1977. Sampling properties of a family of diversity measures. Biometrics 33: 283-292. Solomon, D. S., and J. H. Gove. 1999. Effects of uneven-age management intensity on structural diversity in two major forest types in New England. Forest Ecology and Management 114: 265-274. Spies, T. A. 1998. Forest Structure: a key to the ecosystem. Northwest Science 72: 34-39 Swindel, B. F., J. E. Smith, and R. C. Abt. 1991. Methodology for predicting species diversity in managed forests. Forest Ecology and Management 40: 75-85. Swindel, B. F., L. F. Conde, and J. E. Smith. 1984. Species diversity: concept, measurement, and response to clearcutting and site-preparation. For. Ecol. and Mgmt. 8: 11-22. Temple, S. A., M . J. Mossman, and B. Ambuel. 1979. The ecology and management of avian communities in mixed hardwood-coniferous forests, P.132-153. Pages 268. Management of North Central and Northeastern forests for nongame birds. Workshop Proc. USDA For. Serv. Gen. Tech. Rpt. Vernier, P., A. C. Chan-McLeod, and F. L. Bunnell. 1997. Effects of neighbouring landscape structure on songbirds in the Montane forests of British Columbia. 99 Unpublished paper.. Whittaker, R. H. 1972. Evolution and measurement of species diversity. Taxon 21: 213-251. Willson, M . F. 1974. Avian community organization and habitat structure. Ecology 55: 1017-1029. Zahl, S. 1977. Jackknifing an index of diversity. Ecology 58: 907-913. 100 DBH Distribution Appendix I Plot: High Height Distribution 0.15 -0.1 -0.05 -Species: All Lb 0.35 • 1 0 , 2 5 ^ -I 0.2-0.15 • 0.1 -0.05 • CQ Species: All 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Red Cedar 0.6 -j X a. 0.5 0.4 < 0.3 • 0.2 • tS :rcen 0.1 • Q. n 0.45 X 0.4 S. 0-35 g 0.3 < 0.25 1 °-2 ™ 0.15 | 0.1 £ 0.05 0 Species: Red Cedar ----M Wi s il "" si tm 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Douglas-fir o 0 3 " X S 0.25 -c. 3 0.2 -< •3 0 1 5 -2 o.i -a 0.05 -a. 0 -Species: Douglas-fir 0.35 I 0.3 » 0.25 I 0.2 1 0.15 CQ g °'' g 0.05 a. 0 11 $5 5?* 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Western Hemlock „ 0.6 -X fe 0.5-c g 0.4 -< •3 0.3 -« 0.2 -| 0.1 -0-Species: Western Hemlock t 0.5 K. g 0.4 < •a 0.3 ™ 0.2 c 1 0-1 a. 0 t i l 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 101 DBH Distribution Plot: Med-Lrg Height Distribution Species: All 0.3 -S 0.25 -c. s 0 2 -< S 0.15 -| 0 . 1 -I 0 0 5 " 0 -1. Species: All X 0.25 a. 3 0-2 < 13 0.15 CQ - 0.1 o I 0.05 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Red Cedar 0.25 -I X I 02 ~ I 0.15 -I 0.1 -CQ g 0.05 -a 0 -D B H class Species: Red Cedar 0.7 X 0.6 £ • 0 5 < 0.4 1 0.3 CQ S 0.2 S 0.1 a. 0 tig *>: II 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 0.25 0.2 0.15 S 0-' 0.05 0 CQ Species: Douglas-fir 3$ - — D B H class Species: Douglas-fir 0.3 0.25 0.2 0.15 0.1 c I 0-05 a. 0 CQ 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Western Hemlock S °- 5 " CL g 0.4 -< 13 0.3 -2 0.2 -a o.i -0 -• T i - R A! <W VlK, Species: Western Hemlock 0.6 = 0.5 § 0.4 < a 0 3 ™ 0.2 1 o-1 CL 0 i l f l 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 102 D B H Distribution Plot: Med-Sml Height Distribution Species: All Species: All 0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0 -0.25 • 0.2 — -0.15 -0.1 -g 0.05 -0 -CQ >5£ •pi 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Red Cedar 0.45 -X 0.4 -H 0-35 -| 0.3 -< 0.25 -| 0.2 -? 0.15 -I a i 1 | 0.05 -0 -i D B H class Species: Red Cedar 0.5 £ 0.45 h 0.4 « 0.35 « 0.25 I 0.2 S 0.15 § 0.1 0.05 0 ?€' C 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Dougias-fir Z 0.6 -| 0.5-< 0.4 -1 0.3 -CQ I 0 2 ' y o.i -0 -Species: Douglas-fir 0.7 X 0.6 « ° ' 5 < 0.4 I 0.3 CQ I o.i 0 2.5 7.5 12.5 17.5 : 2.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Western Hemlock 0.35 -? 0.3 -« 0-25 -< 0.2 -1 0.15 -CQ s 01 -I 0.05 -0 -1 it Am Species: Western Hemlock 0.4 X 0.35 8. 0.3 a 0.25 •5 0.2 1 0.15 I °-' cS 0.05 0 i t - - - r " 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 103 DBH Distribution Plot: Low-Lrg Height Distribution Species: All = 0.5 &. g °-4 < •3 0.3 = 0.2 £ 0.1 Q. (I s . — ^ r - T - n D B H class Species: All 0.5 = 0.4 J 035 | 0.3 •3 0-25 | 0.2 I 0 1 5 S o.i °" 0.05 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Red Cedar 0.5 X 0.45 S 0.4 « 0.35 o °-25 | 0.2 ~ 0.15 3 0.1 £ 0.05 0 : = : h = : «v r r n — - • Species: Red Cedar 0.35 X 0.3 « 0 2 5 I 0.2 I 0.15 CQ S 0.1 S3 0.05 a. 0 '4% it -22.5 27.5 32.5 Height class JL 0.6 -0.5 • | 0.4 -< 8 0 3 " • 0.2 -a o.i -0 -Species: Douglas-fir Species: Douglas-fir = 0.3 £ 0.25 < 0.2 1 0.15 CQ S °'' S 0.05 0 13 l i 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class <* 1-2 " X 15 1 -c K 0-8 " < •3 0.6 -CQ 0.4 -g 0.2 -°" 0 -Species: Western Hemlock 1.2 1 ' 0.8 0.6 0.4 0.2 0 Species: Western Hemlock 2.5 7.5 12.5 17.5 : !.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 104 V Plot: Low-Med DBH Distribution Height Distribution 0.3 • * 0.25 f -3 °- 2 • < •3 0.15 -I 0.1-S 0.05 • Species: All o m-[,.i">.. •» 0.45 cd 0.4 1 0.35 It 0.3 < 0.25 1 0.2 aa s 0.15 | 0.1 £ 0.05 0 Species: All II 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Red Cedar S 0.25 c | 0.2 < •3 0.15 £ o.i § 0.05 0 M h Species: Red Cedar 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 0.7 -= 0.6 4 - • I 0-5 -< 0.4 -1 0.3 -co S ( U " I o.i -°~ 0 -Species: Douglas-fir Species: Douglas-fir 0.35 = 0.3 « O-25 K < 0.2 1 0.15 CQ y o.os a. 0 ts --2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Western Hemlock 0.7 = 0.6 | « < 0.4 1 0.3 CQ S 0.2 u S °-' a. 0 lip A' kj> oS> N * N * Nfe*> s * Species: Western Hemlock 0.7 K 0.6 It °'5 < 0.4 1 0.3 CO 1 ° ' 2 S3 0.1 a. 0 *>* > 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 105 Plot: Low-Sml DBH Distribution Species: All 0.6 0.5 | 0.1 0 1 D B H class Height Distribution Species: All 0.45 0 0.4 | 0.35 s °-3 < 0.25 1 0.2 CQ g 0.15 I 0.1 £ 0.05 0 _ * ft y 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5'52.5 57.5 62.5 67.5 Height class 0.6 X 0.5 CH 0.4 0.3 CQ 0.2 s 0.1 0 Species: Red Cedar p 4> DBH class t 0.5 8. S 0.4 | 0.3 2 0.2 c a. 0 Species: Red Cedar if,1 t " S( 7.5 12.5 Height class Species: Red Alder •a 0.6 4- -I 3 1 1 ? 0.4 a 0.2 • ct 0 J DBH class 0.7 * 0.6 0 0-5 < 0.4 1 0.3 CQ 1 °-2 | 0.! 0 Species: Red Alder A' S 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class a L 2 " X 6 1 " c B 0.8 -1 M " ™ 0.4 -c S 0.2 -Species: Paper Birch * 4> DBH class 8. |3 0.8 | 0.6 ™ 0.4 I 0.2 CL 0 Species: Paper Birch 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Heighi class 106 X 0.6 s. 0.5 < 0.4 1 0.3 m 0.2 (3 0.1 0 Spccics: Douglas-fir AY S T L 4» ^ DBH class Species: Douglas-fir 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Western Hemlock = 0.6-8. g 0-5 -< 0.4 • 1 0.3 • cc S 0.2 -| a i i o -DBH class Species: Western Hemlock 0.9 X O.X S. 0.7 g 0.6 < 0.5 1 0.4 2 0.3 § 0.2 £ 0.1 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class s. I 0.8 -a 0.6 -0.4 i S 0.2 -0 -CQ Species: Vine Maple D B H class 8 M ll,.J i0-6 ™ 0.4 S 0.2 a. 0 Species: Vine Maple 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 107 DBH Distribution Plot: Med-Norm Height Distribution Species: All 0.35 E 0.3 a 0.25 < 0.2 1 0.15 ca S 0.1 1 0.05 CL 0 Species: All t 0.25 4- -8. g 0.2 -< 1 0''5 " 2 o.i 4- -g 0.05 -a. 0 -2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Red Cedar 0.6 -8 05 +• -g 0.4 -< 1 "'3 " 2 0.2 -| 0.1 4- - - -* scJ> DBH class 0 Species: Red Cedar 0.4 I 0.35 | 0.3 4- -_=Ji £ 0.25 4- • 1 0.2 • I 0.15 -I 0.1 f -1 0.05 • 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Douglas-fir g 0.2 4--< •a 0.15 -3 ? 0.1 • S 0.05 -0 -4> stf DBH class Species: Douglas-fir 0.45 X 0.4 4-K. 0.35 g 0.3 < 0.25 1 0.2 4-0.15 g 0.1 £ 0.05 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Western Hemlock a 0.4 -? 0.35 4- -i l ° ' 3 " | 0.25 f -•a 0.2 -1 0.15 -1 0.1 -| 0.05 -0 -DBH class Species: Western Hemlock 0.3 * 0.25 g. S 0.2 < t 0.1 c I 0.05 a. 0 --ft 3 -M l tt e IKS i f vm 8 Wits m 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 108 DBH Distribution Plot: Low-Norm Height Distribution Species: Al l 0.45 • = 0.4 f -8. 0.35 J g 0 .3-< 0.25 -1 0.2-2 0.15 -5 0.1 -£ 0.05 -0 -4> 4> D B H class g 0.4 | 0.3 t 0.2 c 1 0.1 a. 0 Species: A l l 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Amabalis Fir 0.7 -0.6 -' 0.5 -0.4 -0.3 -0.2 -0.1 -oJ ft DBH class 0.45 £ 0.4 8. 0.35 2 0.3 < 0.25 1 0.2 2 0.15 c g 0.1 I 0.05 0 Species: Amabalis Fir • -2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Douglas-ftr Z 0.6 8 S 0 5 < 0.4 1 0.3 CQ S 0.2 I 0.1 CL 0 [5*1 ' l' Species: Douglas-ftr 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Western Hemlock 0.6 0.4 0.2 0 D B H class Species: Western Hemlock 0.8 3 0.7 8 0.6 | °'5 1 0.4 co 0.3 I 0.2 | 0.1 0 -.8 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class 109 Plot: Med-Bim DBH Distribution Height Distribution Species: A l l « 0 2 5 " x 8. 0.2-< 0.15 -1 0 .1 -| 0.05 -CL 0 -I'M Species: Al l t 0.25 8. g 0.2 < I 0.15 ? 0.1 B § 0.05 a. 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class Species: Red Cedar t 0.3 j 8. < 0.2 J a 0.15 -CO I 0.05 -0 J . DBH class Species: Red Cedar 0.45 X 0.4 8. 0.35 g 0.3 < 0.25 1 0.2 * 0.15 I 0.1 I 0.05 0 I l l p 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Height class n 0.3 -X S 0.25 -c g 0.2 4 < 0.15 -0.1 -S 0.05 -o -a Species: Douglas-fir DBH class Species: Douglas-fir t 0.5 8. g 0.4 | 0.3 " 0.2 I 0.1 a. 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Heighi class 0.6 -in-g 0.4 4- -< •3 0.3 -a m 0.2 -m Species: Western Hemlock **- V & DBH class Species: Western Hemlock 0.6 ! "-5 CL g 0.4 < | 0.3 * 0.2 C g 0.1 CU 0 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 62.5 67.5 Heighi class no
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Quantifying forest stand diversity using stand structure
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Quantifying forest stand diversity using stand structure Staudhammer, Christina Lynn 1999
pdf
Page Metadata
Item Metadata
Title | Quantifying forest stand diversity using stand structure |
Creator |
Staudhammer, Christina Lynn |
Date Issued | 1999 |
Description | Stand structure is an important component of the overall description of biological diversity in a stand. The diversity of tree sizes within a stand affects tree growth and yield, and is highly correlated with the biodiversity of a stand. Four measures of stand diversity based on the stand distribution of basal area by height, diameter, and species were proposed, assuming a baseline of maximum diversity corresponding to a uniform distribution. First, the Extended Shannon Index (ESI), a measure based on Shannon's Index, was derived. Second, a measure based on fitting the univariate and bivariate distributions of diameter and height was investigated. The third measure, STRI, was derived as a modified R-squared, based on the fit of the distribution to that of a uniform distribution. The fourth measure, STVI, was derived based on comparing the variance of a distribution to that of a uniform distribution. The four measures were evaluated with simulated data and inventory data from the Malcolm Knapp Research Forest (MKRF). The ESI and STVI resulted in logical orderings of both the simulated data and the MKRF data. The advantages of the ESI include its known sampling distribution and accepted use in forestry; however, the ESI depends on allocating data to arbitrary classes. The STVI does not require classifying data; however, its sampling distribution is unknown. Since the measures derived require only normally collected inventory plot data, they are relatively inexpensive to use. However, a thorough investigation of the properties of each should be undertaken. Since evaluating structural diversity is becoming an increasingly important part of forest assessment, a quantitative measure is needed to measure diversity. These measures provide a starting point toward finding an inexpensive, practical measure of structural diversity which gives reasonable results. |
Extent | 5030504 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
FileFormat | application/pdf |
Language | eng |
Date Available | 2009-07-07 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0075007 |
URI | http://hdl.handle.net/2429/10338 |
Degree |
Master of Science - MSc |
Program |
Forestry |
Affiliation |
Forestry, Faculty of |
Degree Grantor | University of British Columbia |
GraduationDate | 2000-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
AggregatedSourceRepository | DSpace |
Download
- Media
- 831-ubc_2000-0154.pdf [ 4.8MB ]
- Metadata
- JSON: 831-1.0075007.json
- JSON-LD: 831-1.0075007-ld.json
- RDF/XML (Pretty): 831-1.0075007-rdf.xml
- RDF/JSON: 831-1.0075007-rdf.json
- Turtle: 831-1.0075007-turtle.txt
- N-Triples: 831-1.0075007-rdf-ntriples.txt
- Original Record: 831-1.0075007-source.json
- Full Text
- 831-1.0075007-fulltext.txt
- Citation
- 831-1.0075007.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0075007/manifest