T H E S T R U C T U R E AND REGENERATION P A T T E R N OF OLDG R O W T H S T A N D S IN T H E M O I S T C O L D E N G E L M A N N S P R U C E - S U B A L P I N E FIR S U B Z O N E O F C E N T R A L B R I T I S H COLUMBIA by PAL V A R G A B.S.F. University of Forestry, Sopron, Hungary 1991 A THESIS S U B M I T T E D I N PARTIAL F U L F I L L M E N T OF T H E R E Q U I R E M E N T S FOR T H E D E G R E E OF M A S T E R OF S C I E N C E in F A C U L T Y OF G R A D U A T E STUDIES Department of Forest Sciences We accept this thesis as conforming to the required standard T H E U N I V E R S I T Y OF BRITISH C O L U M B I A January 1997 © Pal Varga, 1997 In presenting this degree at the thesis in partial fulfilment of University of British Columbia, I agree freely available for reference copying of department this or publication of and study. this his or her representatives. P0vL^9»T The University of British Columbia Vancouver, Canada Date 'ZM f g^C-^Ae-Y that the may be It thesis for financial gain shall not permission. Department of requirements I further agree thesis for scholarly purposes by the Sc\CErc, that advanced Library shall make it by the understood be an permission for extensive granted is for allowed head that without of my copying or my written 11 ABSTRACT The need for sound management practices of high-elevation forests is becoming more and more important in British Columbia; however, there is a lack of information on the dynamics and regeneration pattern of these forests. This pilot study investigates stand dynamics by describing the structure of oldgrowth stands, and correlating the occurrence of regeneration with site factors. Three study stands, which originated after fire and appeared to be in advanced stages of succession were selected in the Moist Cold Engelmann Spruce Subalpine Fir subzone, near Smithers, B.C. All three stands showed the same general stand structure, although the sites were slightly different. The number of seedlings, saplings, alive and standing dead trees of Engelmann spruce (Picea engelmannii Parry ex. Engelm.) was consistently very low, and, hence, excluded from the analysis. The diameter frequency distributions of subalpine fir (Abies lasiocarpa (Hook.) Nutt.) approximated a balanced distribution, but did not fit with the negative exponential function. The height frequency distributions of fir indicated the presence of at least two canopy strata. The age of the trees at breast height (1.3 m) ranged from 100 to 400 years. The age distributions did not conform to distributions expected in the old-growth stage. The age distribution of saplings indicated that in the past 100 years the number of established saplings has declined. The spatial pattern of firs changed from highly aggregated, to random, to regular, with increasing size (diameter or height). The spatial pattern of the iii oldest trees was aggregated possibly indicating remnants of the colonizing cohort. Since the diameter distributions and the spatial pattern of trees with different sizes are similar to that of a steady state they are not expected to change in the future. Canopy stratification was most pronounced i n the oldest study stand, therefore, it is possible the height distribution will never show the inverse-J shape typically expected in old-growth forests. Due to weak size-age relationships the future change i n age distribution is difficult to predict. The fir regeneration showed strong correlations with the presence of decaying wood substrate. There were more than twice as many seedlings on decaying wood than expected from the total cover of the substrate. However, survivorship from seedling to sapling stage was better on the forest floor than on decaying wood. Vegetation cover, light, and humus form showed no correlations with the number and vigour of seedlings. The occurrence of seedlings was also independent of canopy openings. This result implies that gap dynamics was not the driving force for the establishment of fir regeneration. One possible explanation is that snow melts later i n canopy openings. The number of seedlings was higher than expected in the intermediate snow melt time zones which roughly corresponds to canopy edge positions. The future of the spruce component in the study stand was difficult to predict because of the very low number of individuals. The dominant position of spruce trees and presence of recently established seedlings throughout the stands would imply that there was an adequate seed source. However, without knowledge about survival of seedlings until they reach the canopy it is impossible to say whether the amount of regeneration is adequate for maintaining the spruce population. V T A B L E OF CONTENTS ABSTRACT TABLE OF CONTENTS LIST OF TABLES ii v vii LIST OF FIGURES ACKNOWLEDGMENTS 1. INTRODUCTION x xiii 1 1.1. Overview .1 1.2. Stand Structure 6 1.3. Regeneration Pattern 9 2. METHODS 13 2.1 Study Area 13 2.2. Stand Selection and Plot Layout 14 2.3. Stand Structure 15 2.3.1. Sampling and Data Collection 15 2.3.2. Data Analysis 18 2.3.2.1. Diameter Structure 24 2.3.2.2. Height Structure 26 2.3.2.3. Age Structure 27 2.3.2.4. Spatial Pattern 28 2.4. Regeneration Pattern 30 2.4.1. Sampling and Data Collection 30 2.4.2. Data Analysis 32 3. RESULTS 3.1. Stand Structure 37 37 3.1.1. Diameter Structure 37 3.1.2. Height Structure 52 3.1.3. Age Structure 61 3.1.4. Spatial Pattern 66 3.2. Regeneration Pattern 67 VI 4. DISCUSSION 82 4.1. Stand Structure 82 4.2. Regeneration Pattern 85 5. CONCLUSIONS LITERATURE CITED APPENDTX 89 91 100 Vll LIST OF TABLES Table 1. Location and characteristics of the study stands. SNR = soil nutrient regime, SMR = actual soil moisture regime 14 Table 2. Summary of comparisons of the mean dbh of fir. The plots and mean dbh values in cm are shown. Values within a column are not significantly different, since the difference between means is less than 1 cm. The boxes indicate the groups within which all means are not significantly different at cx=0.05 (based on logarithmically transformed data) 39 Table 3. Results from the contingency table analysis comparing the dbh distributions of plots. NS = non-significant, S = significant at the overall error rate of a = 0.05 45 Table 4. Results from the Kolmogorov-Smirnov test comparing the dbh distributions of plots. NS = non-significant, S = significant at the overall error rate of a = 0.05. The shading indicates results that differ from Table 3 45 Table 5. Summary of regression analysis of the negative exponential function on the dbh distributions. R^ is coefficient of determination and S E E is the standard error of estimate (in stems/ha unit) 50 Table 6. Summary of regression coefficients for the analyses shown i n Table 5. For the meaning of symbols refer to Equations 1 and 3; CC' = -PY; andq = e" 50 Table 7. Summary of the lack of fit test results of the negative exponential function on dbh distributions showing degrees of freedom (df), sum of squares (SS), F values and the probability of the F value (Prob. F) with the null hypothesis of no lack of fit 51 Table 8. Summary of regression statistics for dbh - height relationships. R is the coefficient of determination and S E E is standard error of estimate 56 Table 9. Parameter values for dbh - height regression models 56 2b 2 Table 10. Summary of spatial pattern analysis. A denotes aggregated and R regular spatial patterns which were significant (a = 0.05). The numbers indicate the side length of the quadrats in metres. In all other cases the spatial pattern is random 67 Table 11. Number and percent (calculated column-wise) of fir seedlings under the canopy and in canopy openings in the crown transects. z values were used to test the null hypothesis of no difference between percentage of seedlings under canopy and canopy cover. Values printed in bold indicate a significant difference at (X = 0.05 69 Table 12. Number and percent of seedlings in different snow melt time zones in the crown transect of plot 3A. The z values are the statistics to test whether the ratios of number of seedlings in different snow melt time zones is different from the percent cover of each snow melt time zone. The critical values for a - 0.05 are ±1.96 and n = 534. Significant values are marked in bold 70 Table 13. Summary of average percent substrate covers and number of spruce and fir seedlings encountered in the quadrats within plots. The percent covers do not add up to 100, because of the presence of trees in the quadrats. (BI = subalpine fir, Se= Engelmann spruce.) 77 Table 14. Pearson correlation coefficients between number (noBl) and mean vigour (mVig) of fir seedlings and site variables, which were significant at a = 0.05. (See methods section for variable abbreviations.) Number of observations = 169 77 Table 15. Summary statistics and the final set of variables (all of them significant at a = 0.05) of the multiple regressions. R = coefficient of determination, SEE = standard error of the estimate and n = 169 for all plots 78 Table 16. Correlation of the variables with canonical variables. REGl and REG2 are the regeneration canonical variables, SITE1 and SITE2 are the canonical variables of site factors 78 2 IX Table 17. Coefficient of determination between regeneration variables (noBl = number of fir seedlings; mVig = mean vigour of fir seedlings) and the canonical variables (SITE1 and SITE2) of site factors 79 Table 18. Mean value and 95 % confidence intervals for regeneration efficiency coefficient of seedlings calculated for each quadrat. reFF = regeneration efficiency on FF, reDW = regeneration efficiency on DW 80 Table 19. Number and regeneration efficiency of saplings in plots 2A and 3A that were less than 100 years old 81 X LIST O F F I G U R E S Figure 1. Schematic diagram illustrating the location and size of a) plots; b) subplots; and c) quadrats and the method of mapping locations within subplots 16 Figure 2. Species composition for live seedlings, saplings and trees. (BI = subalpine fir, Se= Engelmann spruce). The number of stems was calculated for seedlings by combining the data from three crown transects; for saplings and trees by combining data from all nine plots 22 Figure 3. Box plot of dbh of fir for all study plots. The extent of the box indicates the 2 5 and 75 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles denote 5 and 95 percentiles 38 Figure 4. Dbh frequency distributions of fir for all plots using class width of 1 cm 40 Figure 5. Ratio (q) of number of individuals i n neighbouring dbh classes as a function of class midpoint. The ratios were calculated using class width of 5 cm. The two graphs present the plots with lowest (3B) and highest (1A) variation in q values. The dashed line indicates the average q value 46 Figure 6. Dbh frequency distributions on a logarithmic scale. The class width was 5 cm. The two graphs present the plots with the best (2B) and worst (IB) linear trend. The dashed line indicates the linear trend expected for a negative exponential function 47 Figure 7. Plots of the best fit of the cumulative dbh frequency distribution with the negative exponential function, (a) Cumulative frequency distribution for plot 3C (thick line) and fitted negative exponential curve (thin line), (b) Residual plot 48 th th th th XI Figure 8. Plots of the worst fit of the cumulative dbh frequency distribution with the negative exponential function, (a) Cumulative frequency distribution for plot I B (thick line) and fitted negative exponential curve (thin line), (b) Residual plot Figure 9. Box plot of tree heights. The extent of the box indicates the 2 5 and 7 5 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles note 5 and 9 5 percentiles 49 th th th th 52 Figure 10. Height frequency distributions of saplings and trees using class width of 1 m. The vertical lines indicate the canopy boundaries 53 Figure 11. Density trace diagram of heights. The vertical lines indicate the canopy boundaries 54 Figure 12. Relationship between canopy boundaries and crown heights. The hatched bars on the left side of 0 are the frequency distribution of numbers of canopy trees and the grey bars are the crown height frequency of the same trees. The hatched bars on the right side of 0 represent the frequency distribution of subcanopy trees 57 Figure 13. Relationship between the canopy boundary values and the dbh distribution. The thin continuous horizontal and vertical lines mark the canopy boundary values and the corresponding dbh values; the dashed lines mark dbh = 10 and the corresponding height values. The thick continuous line in the dbh - height graph is the fitted equation and the grey dots are the data points 58 Figure 14. Box plot of sapling ages determined at the base. The extent of the box indicates the 25 and 75 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles note 5 and 95 percentiles 62 Figure 15. Box plot of tree ages determined at 1.3 m. The extent of the box indicates the 2 5 and 75 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles note 5 and 9 5 percentiles 62 th th th th th th th th Xll Figure 16. Age (determined at base) frequency distributions of saplings using class width of 10 years 63 Figure 17. Age (determined at 1.3 m) frequency distributions of trees using class width of 10 years. The open bars represent minimum ages. The dashed lines indicate the maximum sapling age at breast height (1.3 m) 64 Figure 18. Box plot of time required to reach breast height. The extent of the box indicates the 25 and 75 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles note 5 and 95 percentiles 65 Figure 19. Crown transect maps. The lines indicate crown boundaries, the thick lines stems, the grey dots seedlings on DW and the black dots seedlings on FF. The measurements on the axes are in metres 71 Figure 20. Average snow depth, snow water equivalent and snow density measured i n canopy openings (open circle), at canopy edge (grey circle) and under canopy (black circle) 74 Figure 21. Frequency of substrate covers in quadrats for the three A plots. F F = forest floor, DW = decaying wood 76 th th th th Xlll ACKNOWLEDGMENTS I thank my supervisor, Dr. Karel Klinka for introducing me to the various forest types of B.C. as well as various aspects of research; for his great patience and for his support and care during my time at U B C , without which I would not have been able to pursue my studies. Special thanks to my supervisory committee for valuable advice and comments: to Dr. Antal Kozak, for his great help inside and outside U B C and for not using my thesis as a bad example ("a thesis in ecology ..."); to Dr. Leslie Lavkulich for changing my views on soils and helping to realize their importance; and to Dr. Gary Bradfield for providing guidance i n the "not-that-simple" world of multivariate statistics. I thank the G.F. Dodd Operational Research Chair, B.C. Ministry of Forests, Prince Rupert Region for providing funding for the research and also to Dave Coates for his hospitality and help during my stay at Smithers. For the help during the field work, which was probably exciting only for me, I thank David New, Gordon Kayahara, Bernhard Splechtna for their assistance and Sybille Haeussler (Skeena Forestry Consultants) for her help collecting the snow data. For help with ring counting I thank Marcus Rudyk, Nicolas Picard, Graham Bustard, Larissa Puis, Christine Chourmouzis, Luqiong Ling; and for the advice on core preparation special thanks goes to Dr. Jaroslav Dobry. Special thanks to all fellow "Klinkoids" for your friendship and help: to Gordon Kayahara for your help to adjust to U B C , for the use of your private library and for advice in almost every area of life; to Lori Daniels for endlessly correcting my "Hunglish" and for discussions on forest dynamics; to Bob Brett for stimulating debates and for being my English translator; to Han Chen for valuable advice in the art of paper writing and for excellent suggestions on this manuscript; to Christine Chourmouzis for patiently correcting several versions of this manuscript; to Dr. Jaume Fons for helpful discussions on computers; and to all other Klinkoids, Dr. Hong Qian, Dr. Gaofeng Wang, Dr. Qingli Wang for your help i n various things. I thank all "Ponderosans" for your friendship and especially those "night owls" who kept me company. Although my parents and family were far away I still felt their support during my studies for which I am very grateful. Finally, but most importantly I thank B.V. Mary for your encouragement, kindness without which I would not have been able to complete my studies and to you this study is dedicated. 1 1. INTRODUCTION 1.1. Overview High-elevation forests in British Columbia (B.C.) are becoming increasingly important not only for timber supply but also as potential areas for conservation. In the interior of the province, logging of high-elevation forests started i n the beginning of the century, but only became extensive in the late 1970-s (Coates et al. 1994). Most of the management activities take place i n naturally developed old-growth stands, yet little research has been done on the structure and dynamics of these stands (Jull 1990, Coates et al. 1994, Farnden 1994). Also, there is a need to improve current regeneration methods, because of many cases of regeneration failure (Caza 1991, Coates et al. 1994, Farnden 1994). High-elevation ecosystems of interior B.C. are represented by three biogeoclimatic zones - Spruce-Willow-Birch, Engelmann Spruce-Subalpine Fir, and Montane Spruce, all being influenced by a continental subalpine boreal climate. The Engelmann Spruce-Subalpine Fir zone (ESSF) (Krajina 1960, 1965; Meidinger and Pojar 1991) has the largest area and the widest distribution compared not only to the interior subalpine zones but to all the other 11 biogeoclimatic zones of B.C. The E S S F zone extends from 900 to 1700 m elevation i n northern B.C. (57° N) and from 1500 to 2300 m i n southern B.C., and occupies the area between the Coastal Mountains and Rocky Mountains. Similar ecosystems dominated by Engelmann spruce and subalpine fir extend 2 southward to Utah and Colorado where they are identified as the spruce - fir forest type (Alexander 1980, 1986). The main tree species in the ESSF zone are: Engelmann spruce (Picea engelmannii Parry ex. Engelm.), subalpine fir (Abies lasiocarpa (Hook.) Nutt.), and lodgepole pine (Pinus contorta Doug, ex Loud.). Engelmann spruce is considered to be moderately shade-tolerant (Chen 1996), has a high longevity (Alexander and Shepperd 1990) and is usually a dominant canopy species (Krajina 1969, Alexander 1990, Coates et al. 1994). Subalpine fir is considered to be very shade-tolerant (Klinka et al. 1992), has a moderately high longevity (Alexander et al. 1990), and when growing together with spruce is usually a subcanopy species (Krajina 1969, Shea 1985, Veblen 1986, Aplet et al. 1988). In the upper limits of the forested subzone and in the parkland subzone of the E S S F zone, subalpine fir is a dominant canopy species, as the occurrence of spruce decreases with increasing elevation (Meidinger and Pojar 1991). Lodgepole pine is a shade-intolerant (Chen 1996) and short-lived species which dominates forest canopies in mid-seral stages, especially i n the driest subalpine climates (Lotan and Critchfield 1990, Meidinger and Pojar 1990). Oliver (1981) developed a classification that categorizes stand development after major disturbance i n temperate forests of North America into several stages. He avoided applying the climax theory of Clements (1916) which has been disputed for some time (Connell and Slatyer 1977, Cook 1996). The classification has gained a wide acceptance and has been used i n many studies, 3 including investigations of the stand dynamics in high-elevation spruce - fir forests (e.g., Aplet et al. 1988, Jull 1990). Oliver (1981) distinguished four stand developmental stages: (1) stand initiation, (2) stem exclusion, (3) understory reinitiation, and (4) old-growth (Oliver 1981, Oliver and Larson 1990). At the stand initiation stage any regeneration that survived disturbance is released and newly established seedlings utilize the available growing space. In the stem exclusion stage the established trees fill the available growing space by crown expansion thus preventing the establishment of regeneration. With the breaking up of the canopy in time, growing space becomes available for the establishment of seedlings and understory vegetation in the forest understory which characterizes the understory reinitiation stage when trees regenerate and grow without the influence of external disturbances. The final stage of stand development is the old-growth stage "when the trees regenerate and grow without the influence of external disturbances" (Oliver and Larson 1990 page 153). Old-growth forests are usually characterized by the balanced or inverse-J shaped distribution in size (diameter and height) and age and the presence of advanced regeneration. Being under the influence of a continental climate, the most prevalent disturbance in the ESSF zone is fire, which creates a landscape mosaic consisting of stands which vary in their tree species composition and serai stages. Fire usually occurs as a stand-destroying crown event (Romme and Knight 1981, Johnson and Fryer 1989, Jull 1990, Coates et al. 1994). The mean 4 fire return interval varies with climate. In wetter subzones the mean fire return interval is about 350 years, and in the drier subzones about 150 years (Coates et al. 1994). After fire disturbance, lodgepole pine, subalpine fir and/or Engelmann spruce i n various combinations colonize the site i n the stand initiation stage (Whipple and Dix 1979, Aplet et al. 1988, Jull 1990). During secondary succession lodgepole pine is gradually replaced i n the canopy (except i n the driest subzones where it may be a persistent serai species) by shade-tolerant Engelmann spruce and subalpine fir which form the old-growth stands (Meidinger and Pojar 1990, Veblen 1986, Aplet et al. 1988, Coates et al. 1994), although the tree species composition of the final successional stage is not yet known. Based on the high shade-tolerance of subalpine fir and its abundant regeneration in the understory of the old-growth stands in the E S S F zone, it could be expected that, in the absence of disturbance, fir may eventually outcompete spruce (Veblen 1986, Aplet et al. 1988). There are other theories that hypothesize coexistence of the two species. The coexistence is suggested to be possible through (1) regeneration niche differentiation, which means that the two species require different conditions for establishment; (2) different life histories, which means the greater abundance of fir i n the understory is compensated by a higher longevity of spruce (Shea 1985); (3) non-equilibrium coexistence, which means that the disturbances are so frequent that they prevent the competitive exclusion of spruce (Peet 1981). 5 Most studies investigating dynamics of the high-elevation spruce - fir forest were carried out i n the Rocky Mountains of Colorado (Whipple and Dix 1979, Shea 1985, Veblen 1986, Aplet et al. 1988), while only a few investigations were conducted i n the E S S F zone of B.C. (Day 1972, Jull 1990). Considering the large area and environmental variability within both the E S S F zone of B.C. and the spruce - fir forest of the U S Rocky Mountains, it is conceivable that different studies may produce different results with a restricted portability (Coates et al. 1994). Therefore, there is a need to investigate the old-growth E S S F forests for a better understanding of their composition and structure i n relation to stand dynamics. These kind of investigations will help silvicultural management (such as selecting and applying appropriate silvicultural systems or managing for diversity) and conservation (such as selecting and maintaining reserves and parks) (Day 1972, Weetman et al. 1990). The present pilot study aims to improve our understanding of the dynamics of E S S F forests by describing (1) structure of old-growth forests and (2) the occurrence of regeneration based on sampling data. The study questions are: (1) What are the compositional and structural characteristics of the stands considered to be representative of the most advanced stage in stand development i n the study area? and (2) What site factors (canopy cover, vegetation, substrate, humus form, light, snow) does the occurrence of regeneration correlate with? 6 The study subzone, Moist Cold (ESSFmc), represents an intermediate average fire return interval compared to the wet and dry subzones (Coates et al. 1994) . Stand structure of only one developmental stage was sampled because establishing a reliable chronosequence would require a large sample size; and the old-growth stage was selected because of its importance i n the utilization and conservation of forests (Weetman et al. 1990, Forest Practices Code of B.C. 1995) . The study stands, which were selected using structural, field identifiable old-growth criteria, were considered to represent advanced stages of secondary succession after fire disturbance. For the regeneration pattern part, most of the factors believed to influence establishment of regeneration, with the exception of wildlife, were included. In the case of soil, only the humus layer was examined because the soils are cold (Farnden 1994) and most roots, even the larger ones do not deeply penetrate the mineral soil. 1.2. Stand Structure Structure is usually defined as size-age distribution of individuals (MuellerDombois and Ellenberg 1974, Daniel et al. 1979, Kimmins 1987, Thomas et al. 1993). In this study, the definition by Oliver and Larson (1990, page 1) will be used: "Stand structure is the physical and temporal distribution of trees i n a stand. The distribution can be described by species; by vertical or horizontal spatial patterns; by size of tree parts, including the crown volume, leaf area, stem, stem cross section, and others; by tree ages; or by combinations of the above." The structure of the study stands was described by species composition, 7 diameter at breast height (dbh), height, and the age distributions and spatial pattern of saplings and trees. The structure of high-elevation spruce-fir old-growth forest shows the following characteristics. The studies in Colorado found that generally 20-30% of the stems are Engelmann spruce and 80-70% are subalpine fir. However, i n terms of basal area (BA) the ratio is reversed: spruce generally represents about 70% of the total B A and fir 30%, indicating that the large number of firs is due to abundant regeneration (Alexander 1974, Alexander 1985, Roovers and Rebertus 1993). The dbh and height distribution of fir showed the inverse-J shape (Alexander 1985, Roovers and Rebertus 1993). Spruce had a similar dbh distribution in the study by Alexander (1985), but Roovers and Rebertus (1993) found bimodal distributions, with two canopy layer boundaries at 13 and 23 m. These boundaries are indicated by local minimums in the height frequency distribution (Roovers and Rebertus 1993). The age distribution of old-growth forests is expected to show an inverse-J shape like dbh distribution. In the case of fir, this was observed by Roovers and Rebertus (1993) and Aplet et al. (1988). Spruce on the other hand showed a bimodal distribution (Roovers and Rebertus, 1993) even at the stand age of 575 years (Aplet et al. 1988). The spatial pattern of trees appears to change with size. Leemans (1991) found the spatial pattern in old-growth Norway spruce (Picea abies) forests to be random for seedlings, aggregated up to distances of 10 m for saplings, and random again for canopy trees. 8 The objectives of the structural analsysis were to (1) compare dbh and age distributions with steady state functions, (2) evaluate the relationships between dbh, height and age, (3) describe the spatial pattern, and (4) compare the structural attributes between and within stands, when possible. The definition of steady state is that in the absence of major external disturbance the stand structure remains unchanged. The comparison of plots within and between stands can provide a good indication of the variability of structural attributes. Since all of the study stands are considered old-growth based on field identifiable criteria (see Chapter 2), between-stands comparison of structures would indicate the range or variability in structural attributes of old-growth stands. Reconnaissance indicated that the site factors are relatively homogeneous within the study stands, although there are some differences in sites between stands. Therefore, if site is a major determining factor in stand development, plots within a stand are expected to be more similar than plots in different stands. On the other hand, if there is just as much within-stand as betweenstand variation, then the previous disturbance was not homogeneous within the stand, or the initial floristic composition differed, or some other stochastic processes played an important role. It is possible that stands with very different floristic composition and structure could occur on similar sites (Rowe 1961, Lyon and Stickney 1976, Oliver 1981). As it is not known how well the steady state functions should fit the distributions, several methods were used. Since regression analysis provides the most reliable results, regression analysis was accepted as final. There is also no 9 standard method for comparing distributions. Again, several methods were used and the results were compared. 1.3. Regeneration Pattern The structural variability of the tree layer, in conjunction with a variability in surface substrates, should provide ample opportunities for the establishment of seedlings of shade-tolerant species, such as subalpine fir. In fact, the presence of regeneration in the forest understory is one of the characteristics of the oldgrowth stage. It is important to study the pattern of regeneration, particularly in relation to species and microsites (Hofgaard 1993, Frolich and Quednau 1995), to be able to predict future stand structure. Canopy openings, or gaps, were shown in many studies to be related to the presence of advanced regeneration in forest understories. For example, gaps were one of the most important factors in predicting the presence of advanced regeneration of balsam fir (Abies balsamea) and white spruce (Picea glauca) in the boreal forests of Quebec (Kneeshaw and Bergeron 1996). In the subalpine forest of Japan, Yamamoto (1993) found that fir (Abies mariesii), which was usually present before gap formation, readily released, while spruce (Picea jezoensis var. bondoensis) usually established after the gap formation. The density of Norway spruce (Picea abies) regeneration in forests of east-central Sweden increased with increasing gap size (Qinghong and Hytteborn 1991). In the Engelmann spruce-subalpine fir forest of Colorado, Shea (1985) found that 10 on wet sites 63% of spruce and 45% of fir seedlings occurred i n gaps, while on dry sites the seedling occurrence decreased to 34% for spruce and 1% for fir. Understory vegetation may create favourable conditions for germination and seedling growth by providing protection from desiccation or frost damage. In other situations, it may create unfavourable conditions by hindering regeneration and decreasing seedling survival (Zasada and Gregory 1969, Noble and Ronco 1978). The cover of woody vegetation overtopping seedlings was identified as the best growth predictor for Douglas-fir (Pseudotsuga menziesii) seedlings (Wagner and Radosevich 1991). For Engelmann spruce, a competition index, which is based on the cover and height of competing vegetation, was found to be negatively correlated with the seedling growth (Comeau et al. 1993, Kneeshaw and Bergeron 1996). Comeau et al. (1993) also demonstrated that the competition index or growth of Engelmann spruce seedlings can be predicted from available light measured as "the fraction of full sunlight penetrating the vegetation canopy" above the seedlings. However, on high-elevation sites i n the Rocky Mountains of Colorado where sun angles are low, available light is not a good growth predictor for Engelmann spruce and subalpine fir seedlings (Alexander et al. 1990, Alexander and Shepperd 1990). Snow is probably also an important factor influencing the establishment of regeneration, however little is known about its influence on regeneration. Brooke (1970) has shown that snow can cause mechanical damage by creeping on steep slopes. In the study stands, where slopes were gentle, snow is more 11 likely to influence the occurrence of regeneration by the melt pattern i n early summer. It is known that more snow accumulates i n canopy openings than under the canopy (Harestad and Bunnell 1981, Hanley and Rose 1987, Askin and Dragunas 1995) and since the growing season is short (about three months), a difference of only a few weeks i n the disappearance of snow might be significant for regeneration (Brooke 1970). The quality of forest floor or humus form, especially the thickness of litter, is an important edaphic factor determining the regeneration capacity of various microsites. Litter depth greater than about 2.5 cm was found to effectively prevent the establishment of Engelmann spruce seedlings, while subalpine fir seedlings are easily established even on excessively thick forest floors (Knapp and Smith 1982). Downed coarse woody debris (or "nurse logs") and exposed mineral soil are the substrates that support abundant regeneration of Engelmann spruce and subalpine fir (Knapp and Smith 1981, Harmon et al. 1986, Geier-Hayes 1987). In humid climates, the survival of Sitka spruce (Picea sitchensis) and western hemlock (Tsuga heterophylla) seedlings increased with increasing thickness of litter even on relatively undecayed logs (Harmon 1987). However, the long-term survival and growth performance of the seedlings established on decaying wood may be poor because of the low nutrient content of wood (Harmon 1987). The structural analysis would be incomplete and the question about the future stand development could not be fully answered without examining the 12 pattern of regeneration i n the study stands. Therefore, the objectives of this study were: (1) to examine the relationship between the location of seedlings and canopy cover and time of snow melt; (2) to investigate w h i c h site factors (vegetation, light, humus form, or substrate) correlate best w i t h the number and vigour of seedlings; and (3) to make inferences about factors affecting the s u r v i v a l of seedlings to the sapling stage, particularly i n relation to substrate materials. Since this study is based on sampling, it can only provide descriptions and correlations and cannot answer the question of why the regeneration occurs where it does. 13 2. METHODS 2.1 Study Area The study area was located in the Moist Cold Engelmann Spruce Subalpine Fir subzone (ESSFmc) (Krajina 1960, 1965; Meidinger and Pojar 1991) near Smithers, B.C. (54° 47 N 127° 10' W). The E S S F zone is under the influence of the continental subalpine boreal climate. There are no long term climatic data available for the study area, but the ESSFmc subzone's climate can be characterized by cold winters, cool, dry summers, relatively low annual precipitation (450-650 mm), light snowpack (70-150 cm) and a mean annual temperature of-0.7 °C (Banner et al. 1993, Farnden 1994). Mean monthly temperatures are below 0 °C for 6 months and above 10 °C for 1 month (Farnden 1994). The zonal soils are Humo-Ferric Podzols with Hemimor humus forms (Agriculture Canada Expert Committee on Soil Survey 1987; Green et al. 1993, Banner et al. 1993). Subalpine fir and Engelmann spruce are the major tree species i n oldgrowth stands on zonal sites i n the ESSFmc subzone. Lodgepole pine is common in early-seral stands. The most common understory species include: Vaccinium membranaceum, Vaccinium ovalifolium, Menziesia ferruginea, Cornus canadensis, Rubus pedatus, Clintonia uniflora, Pleurozium schreberi, and Barbilophozia floerki. 14 2.2. Stand Selection and Plot Layout Three unmanaged, old-growth stands were selected for the study near Smithers (Table 1). Since the stands had to be identified before sampling, selection was based on the structural definition of old-growth (Hamilton and Pojar 1991, Franklin et al. 1981) instead of criteria based on population processes (Hayward 1991). In this context old-growth is defined as those forest stands that have no pioneer species (e.g. lodgepole pine), have advanced regeneration, abundant snags, and coarse woody debris (Franklin et al. 1981, Thomas et al. 1988, B C M O F 1989). A l l stands selected (Table 1) had (1) originated after fire, which is confirmed by the presence of charcoal in the forest floor (Lorimer 1985, J u l l 1990) and (2) had an intermediate position along the soil moisture - nutrient gradient (Klinka et al. 1989). Table 1. Location and characteristics of the study stands. SNR = soil nutrient regime, SMR = actual soil moisture regime. Stand 1 Stand 2 Stand 3 Location McKendrick Pass McKendrick Pass Hudson Bay Mountain Elevation (m) 1165 1040 1157 Aspect NE NW SE Slope (%) 10 5 10 SNR Medium Poor Rich SMR Fresh Slightly dry Moist 15 The field work was conducted during the summer of 1992. In each stand (1, 2 and 3) three 52 m by 52 m plots (A, B, and C) were established using a rightangle prism and measuring tape (Figure la). The accuracy was checked by laying out three sides of the plot and measuring the fourth which should be 52 m. The corners were re-established until the measurement was within 5% (2.6 m) of the expected value. On all plots a 4 m by 4 m grid was established thus providing 169 subplots (Figure lb). The four meter grid size was selected as a compromise between having a fine resolution for measuring different attributes at the center of the subplots and being able to measure light i n all quadrats i n a short time. In the centre of the 4 m by 4 m subplots, 1 m by 1 m quadrats were established to evaluate regeneration and site factors. Also, i n each A plot a 12 m wide transect was selected so that it would represent the canopy conditions within the plot (Figure lc). 2.3. Stand Structure 2.3.1. Sampling and Data Collection Individuals of the tree species were classified into three categories based on dbh and height: trees (> 10 cm dbh), saplings (< 10 cm dbh and > 1.3 m height) and seedlings (< 1.3 m height). Dbh was selected as a distinguishing characteristic, because it is the easiest to measure. Measurements on trees and saplings are described below, measurements on seedlings and site factors are presented i n section 2.4.1. - Regeneration Pattern Sampling and Data Collection. 16 Figure 1. Schematic diagram illustrating the location and size of a) plots; b) subplots; and c) quadrats and the method of mapping locations within subplots. For all live trees species, dbh, height, crown height, crown diameter, vigour, age and location were recorded on the A plots; species, dbh and height were recorded on the B plots; and only species and dbh on the C plots. In the case of snags (standing dead trees), location was noted on the A plots; species and dbh in all plots. Dbh was measured outside the bark to the nearest 0.1 cm using a dbh tape. Height and crown height (height of the lowest live branch) was measured to the 17 nearest 0.1 m using a clinometer; and the measurements were taken from at least tree height distance from the trees. Crown diameter was measured i n two perpendicular directions with a measuring tape to the nearest 10 cm and then the two measurements were averaged. Vigour was subjectively estimated based on the following criteria: 1 - dead individuals. 2 - broken or damaged top, individuals thought to be severely damaged, and probably will not survive for the next growing season. 3 - moderately damaged individuals that are likely to survive several growing seasons, but growth is retarded. 4 - relatively healthy individuals with some growth. 5 - very healthy individuals with excellent growth. To get an estimate of age, trees were cored at breast height using methods described by Jozsa (1988). The cores were mounted on wooden holders and then sanded. The trees were cored at breast height because (1) during timber cruising the age is measured at breast height, (2) heart rot is less extensive at breast height than closer to the ground. The locations of trees within the subplots were recorded by measuring X and Y coordinates from the 4 m by 4 m gridlines to the side of trees (Figure lc). The measurement error is estimated to be about 10 cm. The coordinates were adjusted by half of dbh so the coordinates would represent the center of the trees. If any measurement was greater than 4 m, then all other measurements within the subplot in the same direction were adjusted. 18 In the case of saplings, species, dbh, height, crown diameter, vigour, location, base age and, for some saplings breast height (bh) age and rooting substrate were recorded on the A plots; species and dbh on the B and C plots. Dbh, crown diameter, vigour and location were measured i n the same way as for trees. Height was measured to the nearest 0.1 m using a measuring rod. To determine the base age of saplings, disks at the base were cut as close to the germination point as was possible. Also, in a systematically selected subsample, disks at breast height were cut as well. Breast height disks were obtained from saplings so that time to reach breast height could be estimated. The disks were sanded and rings were counted using a dissecting microscope with a maximum of 40 x magnification. The rooting substrate of saplings was determined in the richest and poorest A plots (Table 1) (plot 2A and half of plot 3A because of time constraints). The substrate was categorized as: forest floor (FF), decaying wood (DW), and mineral soil (MS). 2.3.2. Data Analysis The inverse-J distribution plays an important role i n describing balanced diameter as well as steady state age distributions; therefore, its derivation will be shown i n some detail. DeLiocourt (in Meyer 1952) observed that the ratio (q) of number of trees i n neighbouring dbh classes is constant i n balanced dbh distributions (balanced meaning that the distribution stays the same over time in the absence of major disturbance) or dbh distributions covering large areas. 19 The constancy of ratios implies that the dbh frequency distribution follows a negative exponential function (Meyer 1952, Leak 1965). This result can be derived i n the following way. If the number of trees i n a dbh class (noted by y) is proportional to the number of trees in a neighbouring class, then the difference between the neighbouring classes (noted by Ay) is also proportional to the number of trees i n either dbh class. The proportion also depends on the dbh class width (noted by Ax). Therefore, the difference in number of trees between two neighbouring dbh class is: Ay = kyAx, Ay — = ky, Ax where k is a constant. if Ax->0, then y'=ky This simple differential equation can be solved by separation. The solution is: y = Ce kx Equation 1 where C and k are constants, x is dbh and e = 2.71828... If k < 0 (and C > 0) then the curve is inverse-J shaped and is called a negative exponential function (Meyer 1952, Leak 1965). Equation 1 can be rewritten i n a more practical form: n = ae ' hD i where Equation 2 is the number of trees in the ith dbh class, Dj is the midpoint of the ith dbh class, and a and b are parameters. Although forest stands in steady state have been described with constant quotient value (Meyer 1952), it is uncertain how close the dbh frequency 20 distribution should match the negative exponential function to be able to say that the dbh structure will stay the same over time. In the case of age distribution, C is the number of individuals entering the population at age zero and | k | is the mortality rate i n Equation 1 (Hett and Loucks 1976). This means that if the mortality rate is the same for all age classes, then the age distribution will also have an inverse-J shape. For the description of structure a summary table of some simple statistics (mean, standard deviation, extreme values) was compiled for dbh, height, crown diameter, age, basal area and volume of trees and saplings. (Because of low values, basal area and volume were not calculated for saplings). Volume was calculated using the B.C. Forest Service (1976) equations for Forest Inventory Zone J , which is applicable for the study area. Before proceeding with the comparison of plots it is important to consider the species composition. The number of spruces was consistently very low i n all size categories and i n all plots (Figure 2, Appendix). Because of the low numbers, spruce was excluded from further analysis. For each structural attribute, the description and comparison of plots/stands involved the following steps. Box plots were used to visually compare plots, and histograms were used to evaluate the shape of the frequency distribution. If it seemed useful, density trace diagrams were used to further investigate the shape of the frequency distributions or to identify minimum or maximum points. The means were then compared by one-way analysis of 21 variance. If the distributions appeared to be similar in shape, then contingency tables and the Kolmogorov-Smirnov test were used for within- and betweenstand comparisons. Because of the central significance of the inverse-J distribution, graphical and regression analysis was used to evaluate the fit of the structural attribute that appeared to follow the inverse-J shape. Finally, the relationships between the various structural attributes were also evaluated using regression analysis. On some occasions, the selection of the next analysis depended on the results, so it was necessary to refer to some of the results i n the methods section. Although histograms are an excellent way of inspecting data distributions, there is a problems with them: the selection of class width is usually arbitrary, and changes i n class width can dramatically change the shape of the histogram (Chambers et al. 1983). One way of determining the optimal class width without changing the shape of the frequency distribution is by using cumulative frequencies and regression analysis (see page 25). Another way to make the class width selection less arbitrary as well as "smoothing" the graph, is to use density trace diagrams (Chambers et al. 1983). The density trace diagram is based on calculating the local density over some interval width h according to the formula: . . . number of observations in[x - h I2, x + h I 2] local density at x = — h x total number of observations 7 22 20000 18000 16000 14000 12000 H co j= 10000 -j 0 w 8000 H 6000 4000 2000 0 Seedlings I 1 BI Saplings Trees Se 2. Species composition for live seedlings, saplings and trees. (BI = subalpine fir, Se= Engelmann spruce). The number of stems was calculated for seedlings by combining the data from three crown transects; for saplings and trees by combining data from all nine plots. Figure Usually the local density is calculated at s, 2s, 3s, etc. values (where s < h/2) so two parameters, s and h would define the density trace diagram for a given set of data. Generally the higher the s and h values are, the smoother the graph becomes. The s and h values were determined by using several 23 combinations and selecting the one which provides relatively smooth curves without too much loss of detail. To further smooth the graph, a cosine weight function was used to calculate the density trace (Chambers et al. 1983). The mean values of the structural attributes were compared using one-way analysis of variance and Tukey's multiple range test (Neter et al. 1990). Only those means were compared which showed differences greater than the minimum significant difference. The minimum significant difference was set based on expected measurement errors. Homoscedasticity was checked using Bartlett's test with a = 0.01 (Walpole 1982). To statistically compare distributions, contingency table analysis (Walpole 1982) and the KolmogorovSmirnov test (Neave and Worthington 1988, Shea 1985, Veblen 1986) were used. Since the Kolmogorov-Smirnov test is sensitive to any kind of differences between two distributions (Neave and Worthington 1988), not just differences in shape, standardized distributions were compared to eliminate possible differences i n mean or variance. Both the contingency table analysis and the Kolmogorov-Smirnov test require several comparisons to hold at the same time; therefore, the a value was adjusted using Bonferoni's adjustment to ensure an a = 0.05 overall error rate (Neter et al. 1990). For both comparisons, plots within stands were compared to see whether there were any significant withinstand differences. Then, using the results from the graphical comparison, plots which appeared to be similar but were from different stands were also compared. The SYSTAT package (Wilkinson 1990) was used for the Kolmogorov-Smirnov 24 test, all other statistical analyses were done using SAS release 6.08 (SAS Institute Inc. 1989). 2.3.2.1. Diameter Structure Since dbh was intensively sampled, it was possible to compare distributions within- and between-stands. In the box plot graph, the plots were ordered according to visual similarity. A class width of 1 cm was used in the histograms, as the minimum significant difference was set to 1 cm. The dbh distributions were non-normal and the Bartlett's test indicated heterocedasticity; therefore, logarithmically transformed means were compared. Since the dbh distributions appeared to be quite similar, they were further compared using contingency tables and Kolmogorov-Smirnov tests. In all plots, the dbh distributions conformed closely to an inverse-J shape. The fit of the negative exponential function on the dbh frequency distributions was evaluated graphically. Plots were created using the ratios (q) of neighbouring dbh classes and logarithmically transformed frequency values, and the dbh class midpoints. The ^-values were calculated using a dbh class width of 5 cm — a smaller width would result in many 0 and non-defined values and a larger class width i n too few points. If the distribution is close to the negative exponential, the q-values should be the same, and frequency values should show a linear trend on logarithmic scale. The regression analysis followed the methods described i n Varga and Klinka (1996). Instead of the frequency distribution, the cumulative frequency 25 distribution was used to avoid problems with selecting class width and classifying dbh values equal to class boundaries. The following equation was fit (which is the indefinite integral of Equation 1): cf =a e^ hh +/ Equation 3 where c/is cumulative frequency, a, ft and fare parameters. To evaluate the fit, residuals were plotted and an approximate lack of fit test was used (Varga and Klinka 1996). As a comparison, the regression analysis and the lack of fit test were carried out on frequency distributions with a class width of 1 cm using the same procedure as i n the case of cumulative frequencies. The g-values were also calculated using the b values from the fit of Equation 2 on frequency distributions: q = e' (Meyer 1952). To check whether the results depend on the 2b sample size, the analysis was also carried out on plots combined within stands and then on all plots combined. The class width was selected so that it was the closest to the reconstructed negative exponential function, which best fitted the cumulative frequency distribution. The cumulative frequency distribution of plot 3C had the best fit with the negative exponential function, therefore, this plot was used to select the dbh class width. First, Equation 1 was reconstructed from the regression coefficients a and f3 resulting in a = 89.56 and b = -0.05110. Then, the class width was selected which provided the minimal mean squared error (MSE), 26 ^(sasee- ™--^) 0 2 MSE = ^ k-2 where n{is the number of trees in the ith dbh class, Z ) ; is the midpoint of the ith dbh class and k is the number of classes. 2.3.2.2. Height Structure The analysis was carried out for combined sapling and tree heights on the A plots. The class width of 1 m, representing the minimum significant difference, was chosen for histograms. The height distributions were non-normal and the Bartlett's test indicated significant heterocedasticity, therefore, logarithmically transformed means were compared. Because of the low number of distributions (three) and the apparent differences in the shape of distributions, no statistical comparison was applied. Also, regression analysis was not used because the shape of height distributions did not approximate the inverse-J distribution. However, the occurrence of bi- and trimodal height distributions, indicating a stratified canopy, warranted further investigation. The canopy strata were compared with crown heights and with the dbh distributions using the dbh height regression results. The selection criteria for the dbh - height regression model were (1) good fit with the data; (2) as few variables as possible; and (3) as simple form as possible. Since a regression model will be used within the data range, the curve behaviour for dbh values >80 cm was not a consideration. Using the 27 Tablecurve® software package by Jandel Scientific the following model was selected: height = b +b dbh + b dbh 2 0 1 25 2 + b dbh 3 3 The fit of the regression model was evaluated with an approximate lack of fit test. Since the intercept should be at 1.3 m (breast height), the intercepts were tested to see whether they were significantly different from 1.3. Using the developed regression models, the canopy boundary values were projected to the dbh distributions to compare the location of minimum/maximum values on dbh and height frequency distributions. (Since the developed dhb - height regressions had a good fit these regressions were used to predict dbh values from height, although the equations may not be the least square solution for height dbh relationship.) 2.3.2.3. Age Structure Separate analyses were used for saplings and trees, since saplings were aged at base and trees at breast height. As trees were distinguished from saplings by dbh, the minimum age for trees varied from plot to plot. A relatively high number of cores (24% of all cores for plot 1A, 16% for plot 2A and 23% for plot 3A) could only provide minimum age (number of recognizable rings) due to heart rot. No age correction was applied considering the very weak relationships between age and dbh or height. The minimum significant difference was set to 10 years resulting i n a 10 year class width i n the histograms. In the case of trees with heart rot the minimum age was plotted. The Bartlett's test indicated 28 significant heterocedasticity, therefore, logarithmically transformed means were compared. No statistical comparison was applied because of the low number of distributions and apparent differences in the shape of distributions. In addition, regression analysis was not used, as the shape of the age distributions did not approximate the inverse-J distribution. Since the trees were cored at breast height (bh), the time required to reach bh was also calculated using the saplings which were aged at both base and bh. The correlation between time to reach bh and base age, base diameter and dbh was calculated to determine whether the time to reach bh is independent from these variables. For the dbh - age and height - age regressions, simple linear models were used, as transformations or more complicated models only marginally improved the fit and made interpretations much more difficult. 2.3.2.4. Spatial Pattern To evaluate the spatial pattern, Morisita's Is index was used (Morisita 1959). This index is not sensitive to density (Morisita 1962) but does depend on quadrat size. The calculation of Is: I =a— * y N(N-l) where q = number of quadrats, m - number of trees i n the i quadrat and th q N =^?i ; ox N= total number of trees in all quadrats. 29 If the spatial pattern is random, (i.e. the quadrat counts follow the Poisson distribution (Thompson 1956, Pielou 1977)) the expected value of I§, Efo) = 1. For aggregated patterns, the expected value of I^is between 1 and q, and for regular patterns, less than 1 (but greater than 0). The Rvalue can be tested with a x test, since Is (N-l)-t-q-N has a x distribution with q-1 degrees of 2 2 freedom. Contiguous quadrat sampling was used (Greig-Smith 1983) because the Morisita's index provides more reliable information with systematic sampling than random sampling (Hairston et al. 1971). Since the trees and saplings were mapped, it was possible to re-sample them using any kind of sampling technique. To aid the calculations, a spatial analysis program (SPA) was developed using Borland C++ language. The program makes it possible to restrict the calculations to certain subsets of trees (e.g. trees with dbh 10-20 cm) or to certain parts of the plot. The randomness of sapling and tree spatial patterns was tested for each 10 cm dbh class, upper and lower canopy height classes and 50-year age classes at a = 0.05 level. For each size/age class, quadrats with side length of 1 to 26 m at 1 m increments were used. The starting point for the contiguous quadrat was set so the analysis was carried out at the central portion of the plot. 30 2.4. Regeneration Pattern 2.4.1. Sampling and Data Collection Seedlings were measured along the crown transects and in the quadrats. In the case of crown transects, the species, location, and rooting substrate of seedlings were recorded. Also the crown projection of trees and saplings were mapped and on plot 3A snow was measured within the transect. Location and rooting substrate were recorded as for the saplings. The crowns were mapped by measuring crown radius in four directions parallel with the plot boundaries. Since the crowns were often distorted, this provided a more accurate representation of crowns, than the average crown diameter. Snow depth, water equivalent and melting pattern were determined from December 23, 1992 to May 20, 1993. Snow depth and the water equivalent (determined by weighting) were measured using a snow tube once a month around 28 poles. (The poles were arranged so that 13 were in the middle of canopy openings, 7 were under the canopy edge, and 8 were under the canopy.) Water equivalent, which is the vertical depth of water that the melting of snowcover would produce, was measured by weighing the snow tube; and snow density was calculated by dividing the water equivalent by the snow depth (Gray and Male 1981). In shallow snow, (less than 50 cm deep) snow depth was measured with a ruler because the snow tube becomes relatively inaccurate (Gray and Male 1981). When the snow started to melt (April 15, 1993), areas 31 without snow were mapped once a week until the snow disappeared (May 20, 1993). In the quadrats, the following attributes were recorded or measured: species, vigour, and rooting substrate of seedlings; and understory vegetation cover, substrate, light level, and humus form. Vigour and rooting substrate were recorded as for the saplings. The percent vegetation cover was estimated according to life form (shrubs, herbs, ferns, mosses and lichens) to speed up sampling. Substrate was estimated as the percent cover of forest floor (FF), decaying wood (DW), and mineral soil (MS) within the quadrat. To characterize available light percent above canopy light (PACL) was used. Above the quadrats photosynthetically active radiation (PAR) was measured i n the morning and afternoon under clear sky conditions using a 'Sunfleck Ceptometer' (Model SF-80; Decagon Devices Inc. 1987). The measurements were taken in two perpendicular directions and averaged (Klinka and Carter 1992). At the same time open sky P A R was continuously measured using Li-Cor LI-190SA quantum sensor and LI-1000 datalogger (Li-Cor Inc. 1988). P A C L was calculated as the ratio of P A R measured above the quadrat and the open sky PAR at the corresponding time and expressed as a percent. For the humus layer, the thickness of the L , F, and H layers were recorded (Green et al. 1993). 32 2.4.2. Data Analysis Considering the very low number of Engelmann spruce seedlings observed, this species was excluded from the analysis (Table 13); thus, only the data for subalpine fir were analyzed. Using the crown transect data, each subalpine fir seedling was evaluated as to whether it was located beneath forest canopy cover or a canopy opening. The canopy of each individual tree was represented by a quarter-ellipse (Koop 1989). If there were only two diameter measures for the crown, then the crown was represented by a circle with diameter of the average crown diameter. This situation happened for all saplings i n plot 2A and for some saplings and trees i n plots 1A and 2A. The total canopy cover within the crown transects was estimated by generating a 5 cm by 5 cm grid i n the computer. Using the coordinates of trees and saplings and the four crown radii values, the centre of each grid cell was marked as being either beneath the canopy (a canopy cell) or not (a canopy opening cell). The ratio of canopy cells to the total number of cells was used to determine the total canopy cover. Although the small cell size increased the computing time, it did not contribute to the measurement and representational errors (crown shapes were approximated by quater ellipses) already present. The percent of seedlings beneath the canopy was compared to the canopy cover by a z - test (Walpole 1982). It was assumed that the cover of DW and F F are independent of canopy cover. Snow depth, density, and water equivalent values were plotted against time to check whether the general trend of snow accumulation i n openings applies to 33 the study stands. To assign a snow melt time to the location of each seedling within the canopy transect in plot 3A, the snow melt map was manually digitized by overlaying a 0.25 m by 0.25 m grid on the map and assigning a snow melt time for each cell and eventually to each seedling. The percent of seedlings within snow melt time zones were then compared to the total area of the snow melt time zone using a z - test (Walpole 1982). The analysis, examining relationships between site factors and the pattern of advance regeneration of subalpine fir, used regeneration and environmental variables. The number of seedlings per quadrat (noBl) and the mean of vigour values of seedlings within the quadrat (mVig) were used as regeneration variables. The mean vigour was set to zero if the quadrat did not contain any seedlings. The following environmental variables were used to describe surface substrates: the cover (%) of forest floor (FF), decaying wood (DW), and mineral soil (MS); thickness of the L, F, and H horizons i n the humus form (Green et al. 1993); the cover (%) of shrubs, herbs, ferns, mosses, and lichens, and the mean percent of above-canopy light (PACL) (Carter and Klinka 1992). The percent cover of M S was omitted for plot 1A because it was always zero. As an exploratory step, Pearson correlation coefficients were calculated between all variables to detect the presence of strong correlations among pairs of variables which might cause multicollinearity i n further analysis. A single high correlation between the regeneration and environmental variables could inflate the canonical correlation values (Dillon and Goldstein 1984). Multiple 34 correlation analysis with backward elimination was used to detect correlation between individual regeneration and environmental variables (Neter et al. 1990). Residual plots were inspected for trends which would indicate the inadequacy of the model. The advantage of multiple correlation is that it is relatively simple and the results are easy to interpret. The disadvantage of this method is that relationships between regeneration variables are not taken into account i n the analysis (Manly 1986). To overcome this limitation, canonical correlation analysis was carried out using the correlation matrix of the variables (Grittins 1985). Multiple correlation analysis selects a linear combination of independent variables which correlates best with a dependent variable, while canonical correlation analysis provides a linear combination of one set of variables that correlates best with a linear combination of another set of variables (Dillon and Goldstein 1984). As canonical correlation analysis was used as an exploratory method, no attempt was made to conform the variables to multivariate normality and homoscedasticity. Since the results indicated that the surface substrate is more important than the cover of understory vegetation and P A C L , the success of seedling establishment on DW versus F F was investigated further. The total cover of M S was very low on plots 2A and 3A so it was not included i n the analysis. To quantify the success of seedling establishment, a regeneration efficiency coefficient (Geier-Hayes 1987) was used: 35 ,.. . ... . % of all seedlings growing on a given substrate regeneration efficiency coefficient = — :— % cover of the substrate in the area When the given substrate is not present in an area, the regeneration efficiency coefficient is not defined. If the coefficient is greater than one then it indicates that the given substrate provides more favourable conditions for germination and early establishment than other kinds of substrates. A value of less than one indicates less favourable regeneration conditions. Prior to calculating the regeneration efficiency coefficient, the % cover of F F and DW i n different plots was compared. First, the frequency of F F and DW substrate cover within quadrats was plotted. Second, the mean cover was calculated for each plot by averaging the cover values of the individual quadrats. The minimum significant difference between the means was set to 10%. As the difference between the means was greater than 10% only for plots 1A and 3A, these two plots were compared with a one-way analysis of variance at a=0.05. Both variables had homogeneous variances at the «=0.05 level as indicated by a Bartlett's test. Third, the regeneration efficiency coefficient was calculated for each quadrat which (1) contained seedlings, and (2) had both F F and DW because the regeneration efficiency is based on a between-substrates comparison. The mean regeneration coefficient values were calculated for both F F and DW for each plot and were compared with a one-way analysis of variance. The preferred way to determine seedling survival would be to monitor seedlings for a certain growth period. However, such a monitoring could not be 36 done i n this study because it might take decades for a suppressed seedling to reach breast height. One way to resolve this problem is to assume that the saplings present i n the study stands represent the seedlings which have survived. Assumptions involved in this approach include: (1) constancy of seed source, (2) constancy of appropriate seedbed (in this case the constancy of the F F and DW cover), and (3) constancy in rates of survivorship. The two substrates types (FF and DW) were used to evaluate seedling survival because they were highly correlated with regeneration variables. Also, the between-plot comparison provided some evidence for constancy of the F F and DW cover. The cover of F F and DW was very similar i n plots 2A and 3A, but the difference i n the maximum ages of trees was at least 100 years. Therefore, only those saplings which were <100 year-old at the base were used in the survival analysis. For those saplings the regeneration efficiency coefficient was calculated i n the same way as shown earlier and compared to the results for seedlings. 37 3. RESULTS 3.1. Stand Structure Most structural attributes showed considerable variation (Appendix). The maximum dbh (84 cm), height (49 m) and volume (794 m /ha) values for fir were 3 high, considering the growth-constraining climatic conditions of the E S S F zone. Stand 3, which was the most productive stand, had the largest maximum sizes (dbh, height). However, the oldest fir trees (about 400 years old) were i n plot 2A. The age of these trees approaches the expected life span of fir (Oosting and Reed 1952). 3.1.1. Diameter Structure The ordering of plots using the box plots indicates that most plots are more similar to the plots from other stands than to the plots from the same stand (Figure 3). Stand 1 showed the most variation and stand 2 the least. Plot I B represents one of the extremes because of the very high number of saplings (Appendix). The box plots showed that the dbh distributions are skewed for all plots: the median is consistently lower than the mean and the 5 and 2 5 th th percentiles are much closer than the 7 5 and 9 5 percentiles. The histograms th th confirmed this observation, all plots were positively skewed (Figure 4). The general trend of the frequency distribution closely approximated an inverse-J shape, although the individual class frequencies varied considerably around the general trend. 38 The ranking of the mean dbh values for the plots was similar to the order the box plot (Table 2). The mean values for plots IC, 2C and 3B were similar and the means in plots 3A and I B were the extreme values. Comparing the means for plots within stands i n the case of stand 2 and 3 there were no significant differences. Also means for plots i n stand 1 and 2 were not significantly different from each other with the exception of plot I B . 60 50 40 I 30 a 20 10 3A 1A Figure 3. Box plot of dbh of fir for all study plots. The extent of the box indicates the 2 5 and 7 5 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles denote 5 and 9 5 percentiles. th th th th 39 Table 2. Summary of comparisons of the mean dbh of fir. The plots and mean dbh values i n cm are shown. Values within a column are not significantly different, since the difference between means is less than 1 cm. The boxes indicate the groups within which all means are not significantly different at «=0.05 (based on logarithmically transformed data). 3A 19:3 3C 17.2 1A 16.6 3B 15.3 2C 15.1 1C 14.9 2A 13.6 2B 12.4 1B 10.0 The contingency table analysis showed that there were significant differences within stands except in the case of stand 2 (Table 3). The results were consistent with the previous results for stand 2 (plots within stand were similar) and for plot I B (different from plots 1A and IC). Plots IC, 2C and 3B appeared to be similar based on box plots and means, therefore, these three plots were also compared using the contingency tables and Kolmogorov-Smirnov test i n addition to within-stand comparisons. According to the contingency table analysis the three selected plots (IC, 2C, and 3B) were not significantly different indicating that there can be more within-stand than between stand variation. 01? zv 43 The Kolmogorov-Smirnov test gave similar results to the contingency table analysis but i n some cases it appeared to be more sensitive (Table 4). The Kolmogorov-Smirnov test indicated significant differences for plot 2B and plots 2A and 2C; and plot 3A became significantly different from plot 3B. From the three similar plots identified before (IC, 2C, and 3B), the 1C-3B comparison was significantly different and the 2C-3B comparison was close to significant (P < 0.041). The graphical evaluation of the fit of the negative exponential function on the dbh distributions provided contrasting results. The ratio of neighbouring dbh classes (q) showed large variation around the mean value, whereas for the negative exponential function the q value is constant (Figure 5). On the other hand, the dbh frequency values generally showed a strong linear trend on a logarithmic scale, and in some cases, the fit appeared to be close to that expected for the negative exponential function (Figure 6). The regression analysis showed that the negative exponential function did not provide an adequate fit for the cumulative frequency and most of the frequency distributions. The coefficient of determination was high and the standard error of estimate was low (Table 5). However, i n the case of the cumulative frequency distributions, all plots showed lack of fit, as indicated by the residual plots (Figures 7 and 8), and the lack of fit test (Table 7). These results may seem contradictory, but, because of the high total variation i n cumulative frequencies, the regression could explain a much higher percent of the total variation than the 44 mean. Cumulative frequencies were consistently overestimated for medium dbh values and underestimated for large dbh values with the exception of plot 3C. Plot 3C showed the best fit, therefore this plot was used to identify the optimal class width which was 1 cm. The regression coefficients for the frequency distributions were quite close to the reconstructed values (Table 6). In the case of frequency distributions, only plots 3A and 3B did not show significant lack of fit (Table 7). The dbh distribution of plot 3C, which had the best fit using cumulative frequencies, showed significant lack of fit for the frequency distribution. The combined plots also showed significant lack of fit which indicated that the results do not change with increasing sample size. The q - plots showed more consistent results with the regression analyses (Figure 5) than the other plot types. Plot 3C, did not show closer fit on the q plots than the other plots in the graphical analysis even though it had the closest fit using regression on cumulative frequencies. Both plot 3A and 3B had good fit using regression on frequencies, but only plot 3B had a good fit based on q values. 45 Table 3. Results from the contingency table analysis comparing the dbh distributions of plots. NS = non-significant, S = significant at the overall error rate of a = 0.05. Stand 1 1A IB IC Stand 2 2A 2B 2C Stand 3 3A 3B 3C 1A IB IC - S NS S - 2A 2B 2C - - NS NS NS - 3A 3B 3C - - NS S - NS - Table 4. Results from the Kolmogorov-Smirnov test comparing the dbh distributions of plots. NS = non-significant, S = significant at the overall erroi rate of a = 0.05. The shading indicates results that differ from Table 3. 1A Stand 1 1A IB IC Stand 2 2A 2B 2C Stand 3 3A 3B 3C IB IC - S NS S - 2A 2B 2C - - NS S 3A 3B 33111 3C - s NS - 46 Plot 3B in in re o _c cr in U) _re u c & Plot 1A in in re u c cr in in re o c cr 0) 30 40 50 60 70 80 90 Dbh (cm) Figure 5. Ratio (q) of number of individuals in neighbouring dbh classes as a function of class midpoint. The ratios were calculated using class width of 5 cm. The two graphs present the plots with lowest (3B) and highest (1A) variation in values. The dashed line indicates the average q value. 47 Plot 2B Plot-IB 0 10 20 30 40 50 60 70 80 90 Dbh (cm) Figure 6. Dbh frequency distributions on a logarithmic scale. The class width was 5 cm. The two graphs present the plots with the best (2B) and worst (IB) linear trend. The dashed line indicates the linear trend expected for a negative exponential function. 48 Dbh (cm) 200 100 CO •g <n a> DC -100 -200 500 1000 Predicted Values Figure 7. Plots of the best fit of the cumulative dbh frequency distribution w i t h the negative exponential function, (a) Cumulative frequency distribution for plot 3 C (thick line) and fitted negative exponential curve (thin line), (b) R e s i d u a l plot. 49 Figure 8. Plots of the worst fit of the cumulative dbh frequency distribution with the negative exponential function, (a) Cumulative frequency distribution for plot IB (thick line) and fitted negative exponential curve (thin line), (b) Residual plot. 50 rH ^H C P rH CO l-H O rO a CO fl >> O • rH -r^ • C M d d cd CO • rH rd . i—1 >H CP rd rS ^ I S fl B o S tn O ft C M rH rH r^ d ^ 2 rd d rH CO CO * 00 C M o o d d C M rH C D m C M i-H rH l-H O m c ~ C D rH cn C O o o d d t> l-H l-H i-H o d o C M C D C D i-i d d d rH cn i-H cn l-H C M i-H C O CO CP o in co d • cn C O i > cn C D m C D rH o i-H C D C D C M O cn o o d cn o cn C M o o od d i—1 rd i—H -rH W 2° « CQ 3 3 -O err a) U rj-J a d Cn 2 °° 2 Cn C D ^ cn co 00 C O o d o d C O C M C O cn t> O l-H i-! C O m CM O lO H O H m m c- ^t< m <N T f l C~- 0 0 (M 0 0 O J M CO C l O ) C - CO 0 0 PS o d d o d d o d d m M> o "43 d > >rHa ^ cl) a "S a ^ ° d CP •rH d) o i—i {rt « CP C O 00 co oo co O C O cn d rH S © II d cn C M i-H rH C M o C D C O rH l> 00 1—1 S -+3 CQ a, 1 S X CM ^3 « C - •rH CO rH CO C O cn C M rf C M c o C n M co H Tj< C M r-i i n cn oo co ai Tf t> CO cn cn cn C l cn cn Cn Ci cn cn cn cn o d d cn cn cn o d d o d d <IfflO C M C M C M < PQ O CO CO CO d • d 1 00 C O t - o d i O l rH C M t> 1> C M C M 00 C D O cn t - Tt o o O d d d• d • i t C M C O 00 00 rH 1 d !>• o C O i-H C D 00 o ro be co. co l-H o rH ^* 1 in 00 d O l o 1—1 rH m o d m m C D m i > m d o d oo i—i rH cn C O C M i-H C O I—I i-H cn r H m i n cn C M in in C M cn m i n r-H d C O C O i > i n rH C M C O m 0 0 cn o C M i-H C O rH C O rH C M C M ^* C M rH C O 00 C M C M C M 00 l-H i—1 i-H t H d H CD o C M 00 00 C O d 00 C O 8 o< co CO C H >> o H C O o X! O o CO cd d c o C O -ri d •rH TJ< d OJ O m C O i—J i n d 00 CO >i CO rfi rH i—i m C M * I© -a co co co co d 0 0 cn r-H cd d cd . , CO CO CP CO >> Q) X o rH o ra C O C M C M m rd cp Cd C M C M m rH d CP C O C M E-i cd 2 • H CP r—l rQ cd a. a cd 10 id rS3 W •2 3 rQ • H H H H > rH O fa a <HH O E O rH rH C P %X > t> t> cm I—t rH -4-» « ° m CO <HH O cp £ CO rO ' 11 o >> * . a o cn 0 0 r H 0 0 0 0 I—I C O in C O 00 t > c- co cn t > o i—i o o o o o d d d di d d d i i i i i i C M 8 a «> a a cl • d tO o • rH i—i C M rO cd d H W 3 O oo ci 00 o 1—1 1 m C M C M C M C O C M m rH I rH rH I in in <! CQ o i-H i-H in C M t> C O C O rH i m cn in co o C O C D 00 t ~ d C O rH 00 rH i C M rH I C M C M 00 C M i—1 i-H C O l-H i i < PQ O < PQ C M rH co i n o o d d C M • o C O o cC M l O t> rH I O C O 51 be • i—t 2 2 rH ^ Is rH o o o ^ Tf 00 <N <N O o © Tf <N CO 00 rH © CN © SH co o - ° 13 eo ta CN $ > .2 P H o ^ cn oi CN r H t > rH r H CN (N CO m Ol 00 CN m Ol Tf TJ •2 > >> O t> in oo o oi rn in oo in CO O l IO M H CO oo o t> in oo co C N m rmH Coi N co n N m w co I N oi CN oi H a CD Tf CN IO CN CO 00 CO CO CO t> co m oi iri 2m HH r-H ! H 03 o © Tf PH •rH SH o © r0 0 nr. <D eo o o cr^ cu O CN rH Ol m 0T1f CN o rH O CD © Tf CD t C N co co in CN rH m CN 00 00 oo m t- oo co Tf CN (N rH r3 in co in co in co in co in C D in C D in co m C D in co CM m CN m CN m C N in C M m C N m C N m C N m C N m ft 5 ® 2 CD ft o3 +J CD £I rO © 13 CC * 0 CO S "5 CO rt CD rH rH +S H TJ co o o © © © Gl 05 © CD 00 01 Tf Tf as CO 01 Tf <N Ol CM CO © q co co co CD l > CO (O IO IO O l CD HH HH O cu ICS Tf CD Tf 0 01 Tf CO CO (N O H O l CD Ol 00 oo r n ci m co rH rH CN Tf CO o Tf Tf ^ H O CO O CN cq N t> H rH CN rH CO rH M CN rH Tf CN rH CN O l rH CO 00 rH t> CN rH CO rH CN CN in in rn co oo in N N co i> oi co in t- rH CN rH CN rH a B£ 2 HH 3 CO «HH O CO CD CD rH bx> CD o CO CO U U § 03 CD O r > > H^H g id rH S Tf to oo co H Tf O C N r H 00 C D oi m C N co in co o co N in co CO C O O O C O 01 co co oo m oi rH rH Ol oo C N rH m T f m co Tf rH co CN > ^ r. o CD •* Ol oi oi o Tf 3 CO © © © in t> i-j O 03 Cto idO to co d d 03 C NM toi io oO Tf CD CO rH CT Nf 0O 01 1 >> r-H CD rCH © © © " CO •H rM o © © PH 03 O O o o O > CD o o o PH «4H O CO M o ccj PH h-l PH P3 PH co p J PH n J PQ CN p PH rH co o CN o u u rH rH o rJ rH 0 +3 o> 01 p o ccj o PH CD PH < eo PQ eo O eo CN 52 3.1.2. Height Structure The mean height values (including trees and saplings) for plots 1A, 2A and 3A were 11.0, 9.5 and 13.3 m respectively. The logarithmically transformed means were not significantly different for plots 1A and 2A, but the mean for plot 3A was significantly different from both of them. The box plot of heights indicated positively skewed distributions similar to the dbh distributions (Figure 9). The histogram of height frequencies confirmed this observation (Figure 10). However, the shape of the frequency distributions did not approximate an inverse-J; rather, it showed modes. Since modal height distribution is consistent with a stratified canopy, it was assumed that the modes represent canopy strata and the minimum values represent the boundaries of the canopy strata. 40 3 5 -\ 10 H 2A Plots Figure 9. Box plot of tree heights. The extent of the box indicates the 2 5 and 75 percentiles of the data. Inside the box the solid line indicates,the median, the dashed line shows the mean. Circles note 5 and 9 5 percentiles. th th th th 89 55 Density trace diagrams with s = 2 and h = 8 were used to identify minimum values i n the frequency distribution, which were regarded as the height of canopy boundaries (Figure 11). Using different s and h values or different starting points generally does not affect the minimum values: the maximum change is 1 m. The height of canopy boundaries were also plotted on the frequency histogram (Figure 10). The lower canopy boundaries were recognizable on the frequency histogram. The second canopy boundary for plot 1A is not as clearly recognizable on the histogram as it is i n the density trace diagram. Also, there are low values around 23 m on the histogram for plot 3A. Since the density trace diagram provides a relatively objective way of assessing the shape of a distribution, the canopy boundaries identified on the density trace diagram were used in further analysis. The maximum crown height of canopy trees was close to the canopy boundary in plot 2A and to the upper canopy boundary i n plot 1A (Figure 12). For the canopy boundary i n plot 3A and the lower canopy boundary i n plot 1A, there is no apparent relationship with the crown height of canopy trees. In all plots, the most frequent height classes (1.3 to 2.3 or 3.3 m classes) were smaller than the minimum crown height of canopy trees. There was a strong relationship between dbh and height (Table 8, Figure 13). There was no significant lack of fit (details are not shown), all variables were significant, and the intercepts were not significantly different from 1.3 (breast height) (Table 9). Since the dbh - height relationship was strong, the canopy 56 boundary values were projected to the dbh frequency distribution and the 10 cm dbh value was projected to the height frequency distribution. The projection of canopy boundary values on the dbh frequency distributions generally did not coincide with the peaks or depressions of the dbh frequency distribution. However, i n the case of plot 1A the lower canopy boundary seemed to coincide to a small depression in the dbh frequency curve at dbh « 17 cm. In all three plots, the 10 cm dbh corresponded well to the point where the height distribution started to level off. This somewhat justifies the choice of 10 cm dbh as a distinguishing criterion between trees and saplings. Table 8. Summary of regression statistics for dbh - height relationships. R is the coefficient of determination and S E E is standard error of estimate. 2 Plot n R 1A 2A 3A 338 366 351 0.964 0.954 0.942 SEE 2 2.05 1.87 2.83 Table 9. Parameter values for dbh - height regression models. Plot Intercept 1A 2A 3A 1.4553 1.0190 1.3545 x 0.1044 0.1281 0.1146 2.5 -0.0212 -0.0279 -0.0234 0.001144 0.001610 0.001258 57 40 Plot 1A 30 Upper Canopy 20 '55 _"=o wer Zan o py 10 fe Subcanopy ZL 0 • I 40 20 • ' ' I • 1 20 1 ' I » • • 1 40 i 1 60 100 80 120 35 Plot 2A 30 25 20 CT> '5 Canopy S3 15 Subcanopy 10 i 5 0 1 1 60 1 I 1 40 1 1 ' I 1 20 ' 1 1 I 1 1 1 1 I 1 1 1 ' I 1 20 0 50 1 1 1 I 1 40 1 1 1 60 I 1 1 1 I 1 80 1 1 1 1 I ' 1 100 1 120 140 Plot 3A 40 30 T 20 Canopy Z3 10 /////} Subcanopy V /////'////•///////// 0 40 20 0 Frequency 20 ^ 40 60 80 (stems/plot) Figure 12. Relationship between canopy boundaries and crown heights. The hatched bars on the left side of 0 are the frequency distribution of numbers of canopy trees and the grey bars are the crown height frequency of the same trees. The hatched bars on the right side of 0 represent the frequency distribution of subcanopy trees. 58 Plot 1A Relationship between the canopy boundary values and the dbh distribution. The thin continuous horizontal and vertical lines mark the canopy boundary values and the corresponding dbh values; the dashed lines mark dbh = 10 and the corresponding height values. The thick continuous line i n the dbh - height graph is the fitted equation and the grey dots are the data points. F i g u r e 13. 59 Plot 2A Figure 13. Continued. 60 Plot 3A 0.04 H i 0 10 20 30 1 1 1 r 40 50 60 70 Dbh (cm) Figure 13. Continued. 80 61 3.1.3. Age Structure The box plots of sapling ages determined at base and tree ages determined at bh (1.3 m) indicated symmetrical, rather than skewed distributions as was the case for the size distributions (Figure 14 and 15). The frequency histograms also confirmed the symmetry (Figure 16 and 17). The age distribution of saplings was unimodal in all plots, approximating a normal distribution. Some saplings were over 200 years old and still had a dbh < 10 cm. The slow growth could be explained by extended periods of suppression. The shape of the age distribution of trees appears to be slightly bimodal with a small peak in the oldest age classes. In the case of tree ages caution should be exercised because: (1) a high number of cores could provide only minimum age; (2) the distinction between saplings and trees was based on dbh, and as a result, the frequency in the youngest age classes might be underestimated, compared to the complete census of bh ages. The time required to reach 1.3 m varied from a minimum of 11 years to a maximum of 140 years (Figure 18). The average was 50 years for plot 1A, 70 years for plot 2A, and 30 years for plot 3A. For plots 1A and 2A the time required to reach 1.3 m was significantly and positively correlated with the base age: for plot 1A the correlation coefficient was 0.64 (P< 0.01), and for plot 2A, 0.44 (P< 0.01). 62 300 200 - £ n Ol >> 150 01 < 100 - 2A Plot Figure 14. Box plot of sapling ages determined at the base. The extent of the box indicates the 25 and 75 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles note 5 and 9 5 percentiles. th th th th 350 Figure 15. Box plot of tree ages determined at 1.3 m. The extent of the box indicates the 2 5 and 75 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles note 5 and 9 5 percentiles. th th th th S9 3 ^ CD CO CD C M fD B P' ^ B CD P CD CD g 11 F° CD Frequency (stems/plot) P. p fD ft ^ P '00 CO p- ft£ p "! pj» 3 CD CD jQ C - CD P 3 B- o P p. CD co P B^ P J C CT PD O X P 3 Bo 05 > •< ro o> 3 B hg > ro CD CO P P P ». CT5 3 cn _1_ co Cfq . o O _L_ ro P <-»- CO CO pr ^ p CD HH- _ O O ^ |—1 CD V-i P CO O O CO CD O CD P cr P CO f9 O O CO > 65 140 120 H (2 CO o >- Figure 18. Box plot of time required to reach breast height. The extent of the box indicates the 2 5 and 75 percentiles of the data. Inside the box the solid line indicates the median, the dashed line shows the mean. Circles note 5 and 9 5 percentiles. th th th th The dbh - age and height - age relationships are not suitable for predictive purposes because of the high variation of the data. Although for all plots, the relationships were significant at a = 0.05 level, the coefficient of determination was around 0.2, and the standard error of estimate was large, about 50 years. This means that the 95% confidence interval for age predictions is about 200 years wide, half of the total age range. 66 3.1.4. Spatial Pattern The spatial pattern varied from highly aggregated, to random, to regular with increasing size (Table 10). The saplings represented the smallest class for dbh (0 - 10 cm) and most of the individuals in the subcanopy height class. The spatial pattern of saplings was aggregated for all plots and all quadrat sizes, which means the aggregates of saplings were also aggregated. The scale of regularity for trees in the largest dbh classes was close to the maximum crown diameter (8 m), but trees i n the upper canopy show regularity at smaller scales (2 - 4 m) as well. In relation to age, young- and medium-aged trees had aggregated and random spatial patterns respectively, as was the case for size. However, for the oldest trees, the spatial pattern was regular only i n plot 1A, while highly aggregated for plots 2A and 3A. 67 Table 10. Summary of spatial pattern analysis. A denotes aggregated and R regular spatial patterns which were significant (cx = 0.05). The numbers indicate the side length of the quadrats in metres. In all other cases the spatial pattern is random. Dbh Class 0 - 10 10 - 20 20 - 30 30 - 40 40 - 50 50 < Plot 1A A: All sizes A: 7, 10 A: 8 R: 8 Height Class Subcanopy Lower canopy Upper canopy Plot 1A A: All sizes Age Class (Trees) 0 - 100 100 - 150 150 - 200 200 < Plot 1A Plot 2A A: All sizes Plot 3A A: All sizes R: 6 R: 10, 12, 13 - Plot 2A A: All sizes Plot 3A A: All sizes - R: 2, 4, 5, 11-16 A: 7, 8, 10 R: 11 - R: 2, 4, 12, 21, 22 R: 3 Plot 2A Plot 3A - A: 3, 5-7, 9-13, 24 A: 3, 8, 13, 16 A: 5, 6, 10 A: 5, 7, 8, 10 3.2. Regeneration Pattern The canopies of the study stands were quite open, especially in plot 2A (Table 11, Figure 19). Because of the openness of the entire canopy, it was difficult to delineate individual canopy openings. The crown transect maps indicate that (1) the location of downed coarse woody debris (decaying logs) is 68 quite well delineated by the seedlings established on them, and (2) the number of subalpine fir seedlings in plot 3A is considerably less than i n the other two plots. In all plots the percent of seedlings established beneath the canopy on the forest floor was closely related to the percent of the total canopy cover even where the 2-test indicated significant differences between the observed and expected frequencies (Table 11). Compared to the percent of the total canopy cover, the percent of seedlings established beneath the canopy on decaying wood was lower i n plot 1A and 3A and the same i n plot 2A. There were significant differences between the observed and expected frequencies for these seedlings i n plots 1A and 3A. The snow accumulation pattern was similar to what was expected: the average snow water equivalent and snow depth were consistently higher i n the canopy openings than under the canopy, but there was no difference i n snow density. Generally more seedlings occurred i n the earlier snow melt time zones than was expected with the exception of 22 April, but the differences were significant only for the three snow melt time zones (Table 12). (The total number of seedlings i n Table 12 and Table 11 are different because 23 seedlings occurred on MS). 69 Table 11. Number and percent (calculated column-wise) of fir seedlings under the canopy and in canopy openings in the crown transects, z values were used to test the null hypothesis of no difference between percentage of seedlings under canopy and canopy cover. Values printed in bold indicate a significant difference at a = 0.05. Plot 1A C. opening Canopy Total z value DW FF Total 473 (50%) 474 (50%) 947 -17.73 131 (23%) 432 (77%) 563 0.95 604 (40%) 906 (60%) 1510 -13.46 DW FF Total 397 (41%) 561 (59%) 958 0.35 148 (43%) 196 (57%) 344 -0.38 545 (42%) 757 (58%) 1302 0.10 DW FF Total 115 (31%) 256 (69%) 371 -2.19 34 (24%) 106 (76%) 140 0.46 149 (29%) 362 (71%) 511 -1.63 Canopy cover 75% Plot 2A C. opening Canopy Total z value Canopy cover 58% Plot 3A C. opening Canopy Total z value Canopy cover 74% 70 Table 12. Number and percent of seedlings i n different snow melt time zones in the crown transect of plot 3A. The z values are the statistics to test whether the ratios of number of seedlings i n different snow melt time zones is different from the percent cover of each snow melt time zone. The critical values for a - 0.05 are ±1.96 and n = 534. Significant values are marked i n bold. Time of snow melt 15 Apr. 22 Apr. 29 Apr. 6 May 13 May 20 May No. of fir seedlings % cover of snow melt time zones z - value 0% 27 (5%) 8% 92 (17%) 13% 216 (40%) 36% 153 (29%) 32% 46 (9%) 11% -0.65 -2.16 2.53 2.22 -1.64 -1.77 - 71 72 74 Average Snow Water Equivalent Average Snow Density 1-Jan 1-Feb 1-Mar 1-Apr 1-May Calendar Date Figure 20. Average snow depth, snow water equivalent and snow density measured i n canopy openings (open circle), at canopy edge (grey circle) and under canopy (black circle). 75 The frequency distribution of substrate covers showed similar trends for all plots (Figure 21): DW had a low coverage (< 40%, Table 2), while F F often covered the whole quadrat. The mean F F and DW covers were significantly different between plots 1A and 3A (Table 13). The Pearson correlation between the number or mean vigour of seedlings and site factors was very low for all plots (Table 14). However, the F F , DW and moss covers were consistently correlated with the regeneration variables - F F was negatively and DW and moss cover were positively correlated. Shrub and fern cover were consistently negatively correlated, and herb and lichen cover were consistently positively correlated with the regeneration variables. The multiple correlation analysis showed weak relationships between regeneration variables and site factors (Table 15) and the residual plots also indicated poor fit (not shown). The variables selected for this analysis were those which had the highest Pearson correlations. The variables that were found significant were mainly related to vegetation layers and substrates, and in two plots to available light. As a consequence of the previous results, it was not surprising that canonical correlation analysis also indicated weak relationships between regeneration and site variables. The canonical correlation for the first canonical variable was around 0.4 for all plots. The redundancy analysis showed that about 14% of the variation in regeneration variables was explained by the first site canonical variable i n all plots. 76 140 120 - I FF V7771 DW I Plot 1A o 100 c 0) 80 3 cr 60 o 40 20 0 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% 140 120 o 100 c 0) 80 H 3 cr 60 o I Plot 2A 40 in 20 0 — 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% 140 120 o c 100 cr o 60 0) 3 Plot 3A 80 40 20 0 A - R — 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Substrate Cover Figure 21. Frequency of substrate covers in quadrats for the three A plots. F F = forest floor, DW = decaying wood. 77 Table 13. Summary of average percent substrate covers and number of spruce and fir seedlings encountered in the quadrats within plots. The percent covers do not add up to 100, because of the presence of trees i n the quadrats. (BI = subalpine fir, Se= Engelmann spruce.) FF(%) DW(%) MS(%) Total number of Se seedlings Total number of BI seedlings Plot 1A Plot 2A Plot 3A 67.8 26.3 0.0 3 58.4 37.4 0.8 5 56.4 34.2 2.0 7 369 392 169 Table 14. Pearson correlation coefficients between number (noBl) and mean vigour (mVig) of fir seedlings and site variables, which were significant at a = 0.05. (See methods section for variable abbreviations.) Number of observations = 169. Plot 1A noBl FF DW Shrub Fern Herb Moss Lichen mPACL mVig Plot 2A noBl mVig -0.19 -0.23 0.22 -0.16 -0.20 0.24 -0.16 Plot 3A noBl mVig 0.16 -0.22 0.22 -0.23 * 0.30 0.17 0.28 0.16 0.24 0.18 0.20 0.16 0.21 78 Table 15. Summary statistics and the final set of variables (all of them significant at a - 0.05) of the multiple regressions. R = coefficient of determination, S E E = standard error of the estimate and n = 169 for all plots. 2 R2 0.11 0.14 0.10 0.09 0.06 0.14 Plot 1A noBl mVig Plot 2A noBl mVig Plot 3A noBl mVig SEE 4.76 1.42 4.52 1.23 1.93 1.30 Variables FF, Moss Shrub, Herb, Moss Shrub, Moss, Lichen FF, L Moss, m P A C L DW, Fern, m P A C L Table 16. Correlation of the variables with canonical variables. REG1 and REG2 are the regeneration canonical variables, SITE1 and SITE2 are the canonical variables of site factors. noBl mVig FF DW MS L F H Shrubs Ferns Herbs Mosses Lichens mPACL Plot 1A REG1 0.835 0.844 SITE1 -0.312 0.418 REG2 -0.550 0.537 SITE 2 0.588 -0.614 - - -0.316 -0.222 -0.167 -0.266 -0.265 0.338 0.779 0.269 -0.009 -0.107 -0.026 -0.113 -0.429 -0.243 0.298 -0.085 0.336 0.170 Plot 2A REG1 -0.719 -0.923 SITE1 0.486 -0.500 -0.116 0.413 0.294 0.228 0.460 0.292 -0.163 -0.313 -0.351 0.354 REG2 0.695 -0.386 SITE 2 0.283 -0.159 -0.123 0.207 0.119 0.021 -0.204 0.107 -0.474 0.264 0.564 -0.103 Plot 3A REG1 0.627 1.000 SITE1 -0.472 0.480 0.327 -0.231 0.027 0.118 -0.256 -0.498 0.181 0.426 0.325 0.449 REG2 -0.779 -0.006 SITE 2 0.065 -0.152 0.281 -0.003 0.370 0.478 0.148 -0.016 0.206 -0.426 0.055 -0.272 79 Table 17. Coefficient of determination between regeneration variables (noBl = number of fir seedlings; mVig = mean vigour of fir seedlings) and the canonical variables (SITE1 and SITE2) of site factors. Plot 1A SITE1 noBl mVig 0.136 0.139 SITE2 0.158 0.160 Plot 2A SITE1 0.096 0.158 SITE2 0.146 0.173 Plot 3A SITE1 SITE2 0.083 0.212 0.099 0.212 Substrate variables (FF and DW) generally had the highest correlations with site canonical variables in all plots, but always with the opposite sign (Table 16). Humus form variables (L, F, and H) played a role only in 2A (L) and 3A (H) plots. Among the vegetation variables, shrub cover was highly correlated with site canonical variables in plots 1A and 2A; herb cover only in plot 2A, moss cover i n plots 1A and 3A, and lichen cover in plot 2A. Light as a variable seemed to be important only i n plot 3A. The important variables i n canonical analysis were about the same as those which were significant in the multiple correlation analysis. However, the thickness of H horizon was not important i n multiple correlation analysis and the fern cover was not important i n canonical analysis. The site canonical variables generally correlated better with the mean vigour than the number of seedlings (Table 17). This trend is most pronounced in plot 3A. The low number of seedlings i n plot 3A might have contributed to 80 the low R i n multiple correlation analysis as well as to the low correlation with 2 site canonical variables. The regeneration efficiency coefficient showed values of 2 to 4 for DW for all plots (Table 18). This means that at least twice as many seedlings occurred on DW than was the percent cover of DW. In other words, downed coarse woody debris was utilized by regeneration of subalpine fir to a much greater extent than the forest floor. Table 18. Mean value and 95 % confidence intervals for regeneration efficiency coefficient of seedlings calculated for each quadrat. reFF = regeneration efficiency on FF, reDW = regeneration efficiency on DW. Plot 1A n Mean St. error Upper limit Lower limit reFF 68 0.40 0.074 0.25 0.54 Plot 2A reDW 68 3.56 0.618 2.33 4.80 reFF 97 0.27 0.064 0.14 0.40 Plot 3A reDW 97 2.64 0.263 2.12 3.16 reFF 61 0.22 0.097 0.03 0.42 reDW 61 3.09 0.404 2.28 3.89 The regeneration efficiency of saplings less than 100 years old was different from that determined for the seedlings (Table 19). The regeneration efficiency of saplings on DW was lower than the regeneration efficiency of seedlings on DW and on F F it was higher; but saplings were still utilizing DW better than F F . As the variance of regeneration efficiency coefficients for saplings can not be calculated, the regeneration efficiency values between seedlings and samplings 81 can not be compared. However, the regeneration efficiency coefficients for saplings are outside of the 95% confidence interval of the regeneration efficiency coefficients for seedlings. This would mean that the seedling and sapling regeneration efficiency coefficients are significantly different if the variance of the sapling efficiency coefficient were known. Table 19. Number and regeneration efficiency of saplings in plots 2A and 3A that were less than 100 years old. Number of saplings on DW on FF Regeneration efficiency on DW on FF Plot 2A 17 10 1.54 0.57 Plot 3A 34 18 1.89 0.61 82 4. DISCUSSION 4.1. Stand Structure The structure of all study stands could be described as follows: marginally mixed species, at least two-storied, and uneven-aged. The low number of Engelmann spruce trees did not permit a detailed analysis, but all plots had both spruce saplings and trees (Appendix). In all A plots the spruce trees had dominant positions (belonged to the canopy), and their age ranged from 180 to 360 years, therefore, they had the potential to produce viable seed (Alexander and Shepperd 1990). The spruce sapling ages at the base ranged from 31 to 104 years i n the A plots (Appendix), indicating recent establishment of regeneration. However, the present study can not evaluate whether the established spruce regeneration could survive to reach the canopy. There was no evidence of recent spruce mortality, since there were only three dead spruce trees and no dead spruce saplings i n all of the study plots. The snag-bole decay rate (Harmon et al. 1986) for Engelmann spruce is 0.012 and the half-time of decay is 67 years; for snag-bark the decay rate is the same and the half-time is 57 years (Mielke 1950). Based on the decay rates it takes roughly 70 years for spruce snags to decompose so that the species would not be easily identifiable any more, which means there was very little mortality of spruce trees i n the past 70 years. Spruce had a much lower mortality rate compared to fir, since fir snags were abundant and their decay rate is higher (decay constant is 0.317 and half-time is 4 years) (Lyon 1977). Higher fir mortality was also observed i n other studies (Oosting and 83 Reed 1952, Day 1972, Veblen 1986) and it represents one of the arguments for coexistence of spruce and fir (Oosting and Reed 1952, Miller 1970, Whipple and Dix 1979). The structure of the study stands was generally very similar despite some differences i n site quality and, possibly in stand history. This would imply that the findings of this pilot study will likely be valid for other stands on similar sites i n the ESSFmc subzone. On the other hand, there was just as much within-stand as between-stand variation in the dbh distribution. This means that within-stand processes are just as important i n influencing dbh distributions as are site differences. The dbh distribution and spatial pattern of trees and saplings showed characteristics of the old-growth stage, but the height and age distributions did not. The dbh distributions become inverse-J shaped very early after disturbance, therefore, they are not reliable indicators of the old-growth stage (Harper 1977, Aplet et al. 1989). In this study plot 3C had the closest fit with the negative exponential function while stand 3 had the youngest saplings and trees (assuming that ages in plot 3A are representative for the whole stand). One possible explanation is that on the more productive site (stand 3) competition is more intensive resulting in greater differentiation (Piatt et al. 1988, Knox et al. 1989, Begon et al. 1990). The spatial pattern of trees with different sizes might be subject to differentiation due to competition as well; 84 therefore, spatial pattern may not be a reliable indicator of old-growth processes (Williamson 1975, Nakashizuka and Numata 1982, Cale et al. 1989). The height distributions were not consistent with the old-growth stage, because they indicated canopy strata (Roovers and Rebertus 1993). Similar height structure i n old-growth stage was observed only for the less shadetolerant Engelmann spruce (Roovers and Rebertus 1993). It was assumed that trees either make it to the upper canopy and continue growing, or remain suppressed i n the lower canopy without any further growth. This i n turn would cause a depression i n the height distribution where the canopy boundary occurs. Since the stands are quite open, the canopy strata can not be attributed to canopy closure (Oliver and Larson 1990). The modal height distribution was most prominent i n the plot (2A) with the oldest trees, which may indicate that it is a steady state distribution for the study stands. The age distributions in the old-growth stage were different from the age distributions of the study stands in (1) the unimodal distribution of sapling ages and (2) the bimodal distribution of tree ages. For subalpine fir, the expected age distribution is inverse-J in the later stages of stand development (Whipple and Dix 1979, Aplet et al. 1988, Roovers and Rebertus 1993). The sapling age distribution indicated a peak in regeneration establishment about 100 - 150 years ago and a decline since then. The lack of established regeneration might be due to lack of ingress or high mortality (Johnson et al. 1994). However, there were no apparent reasons why the ingress should decrease, since with increasing 85 age sites suitable for regeneration increase (gaps, decaying wood, exposed mineral soil from blowdown) (Knapp and Smith 1982). Also, the seedling mortality probably decreased with increasing stand age, since the time to reach bh (1.3 m) was lower for younger seedlings, which may indicate better growing conditions. The shape of the tree age distributions should be interpreted with caution as the number of incomplete cores is high. The age distributions did not show inverse-J shape which means they were not i n a steady state assuming that the mortality rate is independent of age (Hett and Loucks 1976). The future stand structure, especially in terms of size and spatial pattern, probably will only conform closer to what is expected in a steady state (Meyer 1952, Aplet 1988, Oliver and Larson 1990). The future height distribution can not be predicted, because the reasons for canopy stratification are not known, so its change with time can not be evaluated. The future age distribution is also difficult to predict, because the relationships between size and age are very weak, and growth and mortality are more related to size than to age (Harper 1977, Peet 1981, Parker 1988). 4.2. Regeneration Pattern For Engelmann spruce seedlings, a similar pattern was found as for spruce saplings, and trees: the number of spruce seedlings was very low but they were present i n all A plots. Almost all (90%) spruce saplings and seedlings encountered occurred on DW which is consistent with the findings of Knapp and Smith (1981), Harmon et al. (1986), and Geier-Hayes (1987). The seedlings were 86 spread out along the whole length of the crown transect i n all A plots, which indicates that the spruce trees in or around the plot are sufficient i n seeding the plots. In clearcuts, about 50 to 70% of spruce seeds fall within 30 m of the stand edge and about 95% within 100 m (Noble and Ronco 1978, Alexander and Edminster 1983), which compares well with the plot size of 52 m. The ages of spruce seedlings which were encountered in the quadrats varied from 3 years to 154 years, but most were younger than 40 years. This indicates that spruce seedlings have established in the recent past. However, the future of spruce regeneration can not be predicted because the mortality rate of seedlings i n the study stands is unknown. The distribution of subalpine fir seedlings was independent of canopy opening with the exception of plot 1A. This result contrasts with findings of Shea (1985), Qinghong and Hytteborn (1991), Kneeshaw (1992), Yamamoto (1993), and Kneeshaw and Bergeron (1996) who found that gaps were important in determining the location of seedlings. Although gaps are usually defined as canopy openings caused by tree mortality, i n this study no distinction was made between canopy openings caused by mortality or other factors. However, the study plots are relatively homogeneous and there are no apparent edaphic reasons for canopy openings. The percent of seedlings under canopy was only significantly different from the canopy cover for seedlings rooting in DW in plots 1A and 3A where there were more seedlings on DW in canopy openings than expected. This result is consistent with the results of Harmon and Franklin (1989) who suggested that 87 competition could be the reason for high percentage of seedlings occurring on DW. Since plots 1A and 3A are more productive than 2A, the competition is more intense especially i n canopy openings. The DW in canopy openings can provide a relatively competition free environment (from the root competition of both trees and understory vegetation) which would explain the higher percent of seedlings. But on FF, the distribution of seedlings is independent of the canopy cover. This result may question the applicability of gap dynamics (Runkle 1981) for the study stands. The relative number of seedlings i n canopy openings might have been relatively low because of the late snowmelt. The intermediate snow melt time zones had relatively the highest number of seedlings which roughly corresponds to the snow accumulation trend of canopy edges. However, the results should be interpreted with caution, since the measurements were made only i n one plot and only through one snow melt period. The correlation between site variables and the abundance of seedlings was low. Bradfield and Scagel (1984) found similarly low correlation between vegetation strata (including tree species) and environmental variables using canonical correlation analysis. Possible explanations for the weak relationship are: (1) a high variation i n site variables, or (2) non-linear relationships between the variables which the analysis could not detect (Bradfield and Scagel 1984), or (3) some other factor is determining the occurrence of regeneration, like soil temperature or seed predation (Farnden 1994). The relative importance of DW compared to F F i n seedling establishment was the only statistically 88 supported relationship. The observed mean regeneration efficiency coefficient values (2.6-3.5) for decaying wood corresponded well to the value 3.0, determined by Knapp and Smith (1982). However, Kneeshaw (1992) found that substrate was not a determining factor for the regeneration of subalpine fir i n the SubBoreal Spruce zone (the zone below ESSF) - in the same area as this research was conducted i n - although he did not provide data on the relative importance of substrates. If competition is the main reason for the high regeneration efficiency on DW, as was suggested by Harmon and Franklin (1989), then one would expect a strong negative correlation with competing vegetation, especially with shrubs i n the most productive 3A plot. In fact, the regeneration efficiency on DW was higher for 3A plot than for the least productive 2A plot. However, the vegetation cover was a weak predictor of both seedling abundance and vigour. The seedling survival, judged by the presence of saplings, indicated a higher survival on the F F than DW substrate. This means that DW might be more favourable for seedling establishment than F F , but F F may be more favourable for long-term seedling survival than DW. 89 5. CONCLUSIONS The study stands showed similar general structure despite some differences in site and possibly i n stand history. The observed structure was similar to oldgrowth structure of spruce - fir forests i n Colorado in terms of diameter distribution and spatial pattern but different i n terms of height and age distributions. The diameter distribution of trees and saplings was close to a steady state distribution as indicated by the constancy of neighbouring dbh class ratios. The spatial pattern of trees and saplings changed from aggregated to regular with increasing size which was observed i n other old-growth stands. On the other hand, the height distributions consistently indicated modal distribution which was taken as evidence for canopy layers. Since the canopy stratification is most pronounced in the oldest stand, the modal height distribution may represent the steady state condition in these forests. The most important factor determining fir seedling establishment was the presence of DW: the percent of seedlings and saplings on DW were consistently higher than what is proportional to the relative cover of DW. This can not be explained by competition because vegetation correlated poorly with either number or vigour of seedlings. Other site factors, like humus form and light also showed poor correlation with regeneration. The percent of seedlings i n the intermediate snow melt time zone was higher than the relative cover of this snow melt time zone. However, the spatial distribution of seedlings was 90 independent of canopy openings. Therefore, gap dynamics was not an important factor i n seedling establishment in the study stands, but it may play a role i n the survival of seedlings. The future dbh structure and spatial pattern are expected to remain the same since they are close to a steady state. The modal character of the height distribution is also expected to remain unchanged considering that the oldest stand exhibited it most clearly. The future age distribution is difficult to predict because (1) the age distribution is not i n a steady state; and (2) the amount of established regeneration has declined i n the last 100 years. The number of Engelmann spruce seedlings, saplings and trees was consistently low, but they were present in all plots and they had much lower mortality rates than fir. Since spruce seedlings are establishing i n the present conditions, the spruce trees i n the upper canopy are likely sufficient for seeding. However, it is unlikely that the spruce component would increase i n the future. 91 LITERATURE CITED Agriculture Canada Expert Committee on Soils. 1987. The Canadian Soil Classification. 2 ed. Agric. Can. Publ. 1646. 164 pp. nd Alexander, R.R. 1974. Silviculture of subalpine forests in the central and southern Rocky Mountains: The status of our knowledge. U S D A For. Ser. Rocky Mtn. For. Range Exp. Sta., Res. Paper RM-121. Alexander, R.R. 1980. Engelmann spruce-subalpine fir 206. p. 86-87. In F.A. Eyre (ed.). Forest cover types of the United States and Canada. Society of American Foresters, Washington, D.C. Alexander, R.R. 1985. Diameter and basal area distributions in old-growth spruce-fir stands in Colorado. U S D A For. Ser. Rocky Mtn. For. Range Exp. Sta., Res. Note RM-1. Alexander, R.R. 1986. Silvicultural systems and cutting methods for old-growth spruce-fir forests in the central and southern Rocky Mountains. U S D A For. Ser. Rocky Mtn. For. Range Exp. Sta., Gen. Tech. Rep. RM-126. Alexander, R.R., R.C. Shearer, and W.D. Shepperd. 1990. Abies lasiocarpa. p. 6070. In Burns, R.M. and B . H . Honkala (Tech. Coords.). Silvics of North America. Vol. 1, Conifers. U S D A For. Ser. Agric. Handbook 654. Alexander, R.R. and W.D. Shepperd. 1990. Picea engelmannii. p. 187-203. In Burns, R . M . and B . H . Honkala (Tech. Coords.). Silvics of North America. Vol. 1, Conifers. U S D A For. Serv. Agric. Handbook 654. Alexander, R.R., and C B . Edminster. 1977. Uneven-aged management of old growth spruce-fir forests: cutting methods and stand structure goals for the initial entry. U S D A For. Serv. Rocky Mtn. Range Exp. Sta., Res. Pap. RM-186. Askin, R.W., and V.P. Dragunas. 1995. MASS: Snow hydrology pilot study. In Arnott, J.T., W.J. Beese, A . K . Mitchell and J . Peterson (eds.) Montane Alternative Silvicultural Systems (MASS). F R D A Rep. No. 238. Aplet, G.H., R.D. Laven and F.W. Smith. 1988. Patterns of community dynamics in Colorado Engelmann spruce - subalpine fir forests. Ecology 69: 312-319. 92 Aplet, G.H., R.D. Laven and F.W. Smith. 1989. Stemwood biomass and production during spruce-fir development. J. Ecol. 77: 70-77. Banner, A., W. MacKenzie, S. Haeussler, S. Thomson, J. Pojar, and R. Trowbridge. 1993. A Field Guide to Site Indentification and Interpretation for the Prince Rupert Forest Region. B.C. Ministry of Forests, Land Management Handbook No. 26. B.C. Forest Service. 1976. Whole stem cubic metre volume equations and tables. Forest Inventory Division, B.C. Min. of For., Victoria, B.C. BCMOF (British Columbia Ministry of Forests). 1989. Towards an old growth strategy. Summary of the old growth workshop Nov. 3-5, 1989. Integ. Resour. Branch, Min. of For., B.C. 39 pp. Begon, M., J.L. Harper, and CR. Townsend. 1990. Ecology. 2 ed. Blackwell Scientific Publications, London. 945 pp. nd Bradfield, G.E., and A. Scagel. 1984. Correlations among vegetation strata and environmental variables in subalpine spruce-fir forests, southeastern British Columbia. Vegetatio 55: 105-114. Brooke, R.C., E.B. Peterson, and V.J. Krajina. 1970. The subalpine mountain hemlock zone. Ecol. Western N. Am. 2: 153-349. Cale, W.G., G.M. Henebry, and J.A. Yeakley. 1989. Inferring process from pattern: Natural communities. Bioscience 3: 600-605. Caza, CL. 1991. The ecology of planted Engelmann spruce (Picea engelmannii Parry) seedlings on subalpine forest cutovers. Ph.D. Thesis, Fac. of Forestry, Univ. of British Columbia, Vancouver, B.C. 308 pp. Chambers, J.M., W.S. Cleveland, B. Kleiner, and P.A. Tukey. 1983. Graphical Methods for Data Analysis. Wadsworth and Brooks/Cole Publishing Company. Pacific Grove, CA. 395 pp. Chen, Y.H.Y. 1996. Responses of some conifers to light availability in central British Columbia. Ph.D. Thesis Manuscript, Fac. of Forestry, Univ. of British Columbia, Vancouver, B.C. 118 pp. Clements, R.F. 1916. Plant succession: an analysis of the development of vegetation. Carnegie Inst. Pub. 242, Washington, D.C Coates, K.D., S. Haeussler, S. Lindeburgh, R. Pojar, and A.J. Stock. 1994. Ecology and silviculture of interior spruce in British Columbia. FRDAII. B.C. Min. of For. Victoria, B.C. 182 pp. 93 Comeau, P.G., T.F. Braumandl, C.-Y. Xie. 1993. Effects of overtopping vegetation in light availability and growth of Engelmann spruce (Picea engelmannii) seedlings. Can. J . For. Res. 23: 2044-2048. Connell, J . H . and R.O. Slatyer. 1977. Mechanisms of succession in natural communities and their role i n community stability and organization. Am. Nat. I l l : 1119-1144. Cook, J . E . 1996. Implications of modern successional theory for habitat typing: A review. For. Sci. 42: 67-75. Daniel, T.W., J.A. Helms, and F.S. Baker. 1979. Principles of Silviculture. 2 McGraw-Hill Series i n Forest Resources. New York. 500 pp. n d ed. Day, R. J . 1972. Stand structure, succession, and use of southern Alberta's Rocky Mountain forest. Ecol. 53: 472-478. Decagon Devices Inc. 1987. Sunfleck Ceptometer Users Manual. Decagon Devices Inc., Pullman, Wash. Dillon, W.R., and M . Goldstein. 1984. Multivariate Analysis. John Wiley and Sons, Inc. New York. 587 pp. Farnden, C. 1994. Forest regeneration in the E S S F zone of north-central British Columbia. Can. For. Serv. Pacific and Yukon Region, Information Rep. BC-X-351. 31 pp. Forest Practices Code of B.C. 1995. Biodiversity Guidebook. 99 pp. Franklin J.F., K . Cromack, Jr., W. Denison, A. McKee, C. Maser, J . Sedell, F. Swanson, and G. Juday. 1981. Ecological characteristics of oldgrowth Douglas-fir forests. U S D A For. Ser. Pac. N W Res. Sta. Gen. Tech. Rep. PNW-118. Frohlich, M . , and H.D. Quednau. 1995. Statistical analysis of the distribution pattern of natural regeneration in forests. For. Ecol. Mgt. 73: 45-57. Geier-Hayes, K . 1987. Occurrence of conifer seedlings and their microenvironments on disturbed sites i n central Idaho. U S D A For. Ser. Intermtn. Res. Sta. Res. Pap. INT-383. Good, B . J . and S.A. Whipple. 1982. Tree spatial patterns: South Carolina bottomland and swamp forests. Bull. Torrey Bot. Club 109: 529-536. Gray, D . M . and D.H. Male (eds). 1981. Handbook of Snow. Priciples, Processes, Management and Use. Pergamon Press, Toronto. 776 pp. 94 Green, R.N., R.L. Trowbridge, and K. Klinka. 1993. Towards a taxonimic classification of humus forms. For. Sci. Monogr. 29: 1-48. Greig-Smith, P. 1983. Quantitative Plant Ecology. 3 ed. Univ. of California Press, Berkeley, CA. 359 pp. rd Gittins, R. 1985. Canonical Analysis. Springer-Verlag, Berlin. 351 pp. Hairston, N.G., R.W. Hill, and U. Ritte. 1971. The interpretation of aggregation patterns, p. 337-356 In Patil, G.P., E.C. Pielou, and W.E. Walters (eds.). Statistical Distributions. 582 pp. Eclology Volume 1. Spatial Patterns and Statistical The Pennsylvania State Univ. Press. University Park. Hamilton, E. and J. Pojar. 1991. Old growth forests in British Columbia, p. 449457. In XIX IUFRO World Congress Proceedings. Montreal, Quebec. 5-11 Aug. 1990. Hanley, T.A. and CL. Rose. 1987. Influence of overstory on snow depth and density in hemlock-spruce stands: Implications for management of deer habitat in southeastern Alaska. USDA For. Ser. Pac. N.W. Res. Sta. Res. Note PNW-459 11 pp. Harestad, A.S. and F.L. Bunnell. 1981. Prediction of snow-water equivalents in coniferous forests. Can. J. For. Res. 11: 854-857. Harmon, M.E., J.F. Franklin, F.J. Swanson, P. Sollins, S.V. Gregory, J.D. Lattin, N.H. Anderson, S.P. Cline, N.G. Aumen, J.R. Sedell, G.W. Lienkaemper, K. Cromack, Jr., and K.W. Cummins. 1986. Ecology of coarse woody debris in temperate ecosystems. Advances in Ecological Research 15: 133-302. Harmon, M.E. 1987. The influence of litter and humus accumulations and canopy openness on Picea sitchensis (Bong.) Carr. and Tsuga heterophylla (Raf.) Sarg. seedlings growing on logs. Can. J. For. Res. 17: 1475-1479. Harmon, M.E. and J.F. Franklin. 1989. Tree seedlings on logs in Picea-Tsuga forests of Oregon and Washington. Ecol. 70: 48-59. Harper, J.L. 1977. Population Biology of Plants. Academic Press, London. 892 pp. Hayward, G. D. 1991. Using population biology to define old-growth forests. Wildlife Soc. Bull. 19: 111-116. 95 Hett, J.M. and O.L. Loucks 1976. Age structure models of balsam fir and Eastern hemlock. J. Ecol. 64: 1029-1044. Hofgaard, A. 1993. Structure and regeneration patterns in a virgin Picea abies forest in northern Sweden. J. Veg. Sci. 4: 601-608. Hozumi, K., K. Shinozaki, and Y. Tadaki. 1968. Studies on the frequency distribution of the weight of individual trees in a forest stand I. A new approach toward the analysis of the distribution function and the -3/2th power distribution. Jap. J. Ecol. 18: 10-20. Hutchings, M.J. and C.S.J. Budd. 1981. Plant competition and its course through time. Bioscience 31: 640-645. Johnson, E.A., and G.I. Fryer. 1989. Population dynamics in lodgepole pine Engelmann spruce forests. Ecology 70: 1335-1345. Johnson, E.A., K. Miyanishi, and H. Kleb. 1994. The hazards of interpretation of static age structures as shown by stand reconstructions in a Pinus contorta - Picea engelmannii forest. J. Ecol. 82: 923-931 Jozsa, L. 1988. Increment core sampling techniques for high quality cores. Special Pub. No. SP-30. Forintek Canada Corp., Vancouver, B.C. 26 pp. Jull, M.J. 1990. Long-term stand dynamics in high-elevation Engelmann spruce - subalpine fir forests. M.Sc. Thesis, Fac. of Forestry, Univ. of British Columbia, Vancouver, B.C. 145 pp. Kimmins, J.P. 1987. Forest Ecology. Macmillan Publishing Co., New York. 531 pp. Klinka, K., V.J. Krajina, A. Ceska, and A.M. Scagel. 1989. Indicator Plants of Coastal British Columbia. Univ. of British Columbia Press. Vancouver, B.C. 288 pp. Klinka, K., Q. Wang, R. E. Carter, H. Y. H. Chen. 1996. Height growth-elevation relationships in subalpine forests of interior British Columbia. For. Chron. 72: 193-1998. Klinka, K., Q. Wang, G.J. Kayahara, R.E. Carter, and B.A. Blackwell. 1992. Light-growth response relationships in pacific silver fir (Abies amabilis) and subalpine fir (Abies lasiocarpa). Can. J. Bot. 70: 1919- 1930. 96 Knapp, A.K. and W.K. Smith 1982. Factors influencing understory seeedling establishment of Engelmann spruce (Picea engelmannii) and subalpine fir (Abies lasiocarpa) in southeast Wyoming. Can. J. Bot. 60: 2753-2761. Kneeshaw, D.D. 1992. Tree population dynamics of some old sub-boreal spruce stands. M.Sc. Thesis. Fac. of Forestry, Univ. of British Columbia, Vancouver, B.C. 152 pp. Kneeshaw, D.D. and Y. Bergeron. 1996. Ecological factors affecting the abundance of advance regeneration in Quebec's southwestern boreal forest. Can. J. For. Res. 26: 888-898. Knox, R.G., R.K. Peet, and N.L. Christensen. 1989. Population dynamics in loblolly pine stands: changes in skewness and size inequality. Ecology 70: 1153-1166. Koop, H. 1989. Forest Dynamics. SILVI-STAR: a comprehensive monitoring system. Springer-Verlag, Berlin. 229 pp. Krajina, V.J. 1965. Biogeoclimatic zones and biogeocoenoses of British Columbia. Ecol. West. N. Am. 1: 1-17. Krajina, V.J. 1969. Ecology of forest trees in British Columbia. Ecol. West. N. Am. 2: 1-152. Leak, W.B. 1965. The J-shaped probability distribution. For. Sci. 11: 405-409. Lertzman, K.P. 1992. Patterns of gap-phase replacement in a subalpine, oldgrowth forest. Ecology 73: 657-669. Li-Cor Inc. 1988. LI-100 Datalogger Instruction Manual. Li-Cor Inc., Lincoln, Neb. Lorimer, C.G. 1985. Methodological considerations in the analysis of forest disturbance history. Can. J. For. Res. 15: 200-213. Lotan, J.E., and W.B. Critchfield. 1990. Pinus contorta. p. 302-315. In Burns, R.M. and B.H. Honkala (Tech. Coords.). 1990. Silvics of North America. Vol. 1, Conifers. USDA For. Serv. Agric. Handbook 654. Lyon, L.J. 1977. Attrition of lodgepole pine snags on the Sleeping Child Burn, Montana. USDA For. Serv. Intermt. For. and Range Exp. Stn. Res. Note INT-219. 97 Lyon, L.J., and P.F. Stickney. 1976. Early vegetational succession following large Rocky Mountain wildfires. Proc. Tall Timbers Ecol. Conf. 14. pp. 355-373. Manly, B.F.J. 1986. Multivariate Statistical Methods. A Primer. Chapman and Hall, London. 159 pp. Meidinger, D. and J. Pojar (eds). 1991. Ecosystems of British Columbia. B.C. Min. of For. Special Rep. Ser. 6. 330 pp. Meyer, A.H. 1952. Structure, growth, and drain in balanced uneven-aged forests. J. For. 50: 85-92. Mielke, J.L. 1950. Rate of deterioration of beetle-killed Engelmann spruce. J. For. 48: 882-888. Miller, P.C. 1970. Age distributions of spruce and fir in beetle-killed forests on the White River Plateau, Colorado. Am. Midi. Nat. 83: 206-212. Morisita, M. 1959. Measuring of interspecific association and similarity between communities. Mem. Fac. Sci. Kyushu. Univ. Ser. E. (Biol.) 3: 65-80. Morisita, M. 1962. Ig-index, a measure of dispersion of individuals. Res. Popul. Ecol. 4: 1-7. Nakashizuka, T. and M. Numata. 1982. Regeneration process of climax beech forests: I. Structure of a beech forest with the undergrowth of Sasa. Jpn. J. Ecol. 32: 57-67. Neave, H.R. and P.L. Worthington. 1988. Distribution-free tests. Unwin Hyman, London, 430 pp. Neter, J., W. Wasserman, and M.H. Kutner. 1990. Applied Linear Statistical Models. R.D. 3 ed. Irwin Inc., Homewood, Ilinois. 1181 pp. rd Noble, D.L., and F. Ronco, Jr. 1978. Seedfall and establishment of Engelmann srpuce and subalpine fir in clearcut openings in Colorado. USDA For. Serv., Rocky Mtn. For. Range Exp. Sta., Res. Pap. RM-200. Oliver, CD. 1981. Forest development in North America following major disturbances. For. Ecol. Manage. 3: 153-168. Oliver, CD. and B.C. Larson. 1990. Forest Stand Dynamics. McGraw-Hill Inc., New York. 467 pp. Oosting, H.J., and J.F. Reed. 1952. Virgin spruce-fir forests of the Medicine Bow Mountains, Wyoming. Ecol. Monogr. 22: 69-91. 98 Parker, A.J. 1988. Stand structure in subalpine forests of Yosemite National Park, California. For. Sci. 34: 1047-1058. Peet, R.K. 1981. Forest dynamics of the Colorado Front Range: Composition and dynamics. Vegetatio 45: 3-75. Peet, R.K., and N.L. Christensen. 1987. Competition and tree death. BioScience 37: 586-595. Pielou, E.C. 1977. Mathematical Ecology. Wiley and Sons Inc., New York. 385 pp. Qinghong, L. and H. Hytteborn. 1991. Gap structure, disturbance and regeneration in a primeval Picea abies forest. J. Veg. Sci. 2: 391-402. Romme, W.H. and D.H. Knight. 1981. Fire frequency and subalpine forest succession along a topographic gradient in Wyoming. Ecology 62: 319-326. Roovers, L.M. and A. J. Rebertus. 1993. Stand dynamics and conservation of an old-growth engelmann spruce-subalpine fir forest in Colorado. Nat. Areas J. 13: 256-267. Rowe, J.S. 1961. Critique of some vegetational concepts as applied to forests of northwestern Alberta. Can. J. For. Res. 39: 1007-1017. Runkle, J.R. 1981. Gap regeneration in some old-growth forests of the eastern United States. Ecology 62: 1041-1051. SAS Institute Inc. 1989. SAS/STAT User's Guide, Version 6. 4 ed. Volume 1-2, th SAS Institute Inc., Cary, North Carolina. 1686 pp. Shea, K.L. 1985. Demographic aspects of coexistence in Engelmann spruce and subalpine fir. Am. J. Bot. 72: 1823-1833. Sokal, R.R. and F.J. Rohlf. 1981. Biometry. 2 ed. W.H. Freeman and Company, New York, 859 pp. nd Thomas, J.W., L.F. Ruggliero, R.W. Mannan, J.W. Schoen, and R.A. Lancia. 1988. Management and conservation of old-growth in the US. Wildlife Soc. Bull. 16: 252-261. Thomas, J.W., M.G. Raphael, R.G. Anthony, E.D. Forsman, A.G. Gunderson, R.S. Holthausen, B.G. Marcot, G.H. Reeves, J.R. Sedell, and D.M. Solis. 1993. Viability assessments and management considerations for species associated with late-successional and old-growth forests of the Pacific Northwest. USDA Nat. For. System, For. Ser. Res. 530 pp. 99 Thompson, H.R. 1956. Distribution of distance to nth neighbour ina population of randomly distributed individuals. Ecology 37: 391-394. Varga, P., and K. Klinka. 1996. Diameter distribution of some subalpine fir stands in central British Columbia. Env. Mon. Assessm. 39: 601-610. Veblen, T.T. 1986a. Age and size structure of subalpine forests in the Colorado Front Range. Bull. Torrey Bot. Club 113: 225-240. Veblen, T.T. 1986b. Treefalls and the coexistence of conifers in subalpine forests of the central Rockies. Ecology 67: 644-649. Wagner, R.T. and S.R. Radosevich. 1991. Neighbourhood predictors of interspecific competition in Douglas-fir plantations. Can. J. For. Res. 21: 821-828. Walpole, R.E. 1982. Introduction to statistics. 3 ed. Macmillan Publishing Co., Inc., New York, 521 pp. rd Weetman, G.F., E. Panozzo, M. Jull, and K. Marek. 1990. An assessment of opportunities for alternative silvicultural systems in the SBS, ICH and ESSF biogeoclimatic zones of the Prince Rupert Forest Region. B.C. Min. For., Prince Rupert For. Reg., Smithers, B.C. Contract Rep. 154 pp. Whipple, S.A. and R.L. Dix. 1979. Age structure and successional dynamics of a Colorado subalpine forest. Am. Midi. Nat. 101: 142-158. Wilkinson, L. 1990. SYSTAT: The system for statistics. SYSTAT Inc., Evanston, Illinois. 676 pp. Williamson, G.B. 1975. Pattern and serai composition in an old-growth beechmaple forest. Ecology 56: 727-731. Yamamoto, S.-L 1993. Gap characteristics and gap regeneration in a subalpine coniferous forest on Mt Ontake, central Honshu, Japan. Ecol. Res. 8: 277-285. Zasada, J.C., and R.A. Gregory. 1969. Regeneration of white spruce with reference to interior Alaska: a literature review. USDA For. Ser. Pac. NW For. Exp. Sta., Res. Pap. PNW-79. 100 APPENDIX Summary of structural attributes of live trees and saplings. (BI = subalpine fir, Se= Engelmann spruce, Ht = height, Cr. diam. = average crown diameter, BA = basal area, No. of obs. = number of observations, B. diam. = diameter at base, BH Age = age at breast height) TREESPlot BI 1A No. o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha IB IC Standi N o . o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha N o . o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha N o . o f obs. Mean St. dev. Min. Max. Sum/plots Total/ha Stems Dbh (cm) Ht (m) Age (yrs) Cr. d i a m (cm) 155 31.4 15.3 10.0 73.9 155 20.7 8.8 4.9 39.5 118 155 367 86 155 595 149 26.3 13.4 10.2 60.2 149 16.7 8.1 4.7 32.7 155 573 149 551 171 26.3 12.7 10.0 61.5 475 586 Volume (m ) 3 2 155 0.096 0.090 0.008 0.429 14.855 54.94 155 1.05 1.22 0.02 6.13 163.04 603 149 0.068 0.066 0.008 0.285 10.194 37.70 149 0.66 0.81 0.02 3.50 97.88 362 171 0.067 0.061 0.008 0.297 11.469 42.41 171 632 475 28.0 14.0 10.0 73.9 157 62 62 295 BA (m ) 475 0.077 0.074 0.008 0.429 304 18.7 8.7 4.7 39.5 3 6 4 5 5 1 7 0 2 3 0 4 0 86 106 0 02 6 13 2 6 0 9 4 g 2 2 101 Appendix Continued. TREES Plot BI 2A N o . o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha 2B 2C N o . o f obs. Mean St. d e v . Min. Max. Sum/plot Total/ha N o . o f obs. Mean St. d e v . Min. Max. Sum/plot TotaUha S t a n d 2 N o . o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha Stems Dbh (cm) Ht (m) Age (yrs) 174 23.9 10.2 10.2 57.0 174 16.9 7.2 3.9 32.0 145 188 55 88 403 214 23.9 12.1 10.0 64.5 212 17.2 7.1 4.1 33.8 174 643 214 791 174 25.9 12.8 10.1 66.9 562 693 174 331 79 160 560 BA (m ) 2 Volume (m3) 174 0.053 0.045 0.008 0.255 9.195 34.01 174 0.48 0.53 0.02 3.16 83.58 309 214 0.056 0.059 0.008 0.327 12.004 44.39 212 0.54 0.70 0.02 4.20 113.73 421 174 0.066 0.063 0.008 0.352 11.436 42.29 174 643 562 24.5 11.8 10.0 66.9 Cr. d i a m (cm) 386 17.0 7.1 3.9 33.8 562 0.058 0.057 0.008 0.352 32.634 40.23 386 0.51 0.63 0.02 4.20 197.31 365 102 Appendix Continued. TREES Plot BI 3A No. o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha 3B Stems No. o f obs. Mean Age Cr. d i a m BA (m) (yrs) (cm) (m ) 195 195 150 195 195 Volume 2 195 31.4 15.7 21.5 9.5 147 328 0.097 57 104 0.093 1.30 10.0 3.7 45 114 0.008 0.02 83.5 49.3 291 725 0.548 18.840 69.67 8.57 214.74 204 204 205 204 1.10 794 26.6 19.9 0.068 0.73 12.8 8.2 0.065 0.84 0.02 Min. 10.0 5.7 0.008 Max. 61.8 37.4 0.300 3.78 13.971 147.82 Sum/plot 205 Total/ha 758 No. o f obs. Mean St. dev. 51.67 266 0.064 13.0 0.071 10.0 0.008 Max. 83.9 0.553 266 984 S t a n d 3 No. o f obs. Mean 547 266 25.3 Min. Sum/plot Total/ha 16.906 62.52 666 399 666 399 27.5 20.7 0.075 St. dev. 14.0 8.9 0.077 1.10 Min. 10.0 3.7 0.008 0.02 83.9 49.3 Max. All Stands Ht (cm) 195 721 St. dev. 3C Dbh Sum/plot 666 Total/ha 821 No. o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha 0.553 8.57 49.716 362.56 61.29 1703 1089 413 524 0.91 1703 670 1089 26.6 18.8 164 341 0.070 0.75 13.4 8.4 60 92 0.071 0.96 10.0 3.7 45 114 0.008 0.02 83.9 49.3 403 725 0.553 8.57 1703 118.868 820.79 700 48.84 506 103 Appendix Continued. TREES Se Stems Plot 1A IB IC No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plots Total/ha Dbh Ht Age Cr. d i a m BA (cm) (m) (yrs) (cm) (m ) 2 2 2 2 32.8 10.2 25.6 40.0 22.2 11.0 14.4 29.9 2 7 6 6 23.7 24.0 12.5 72.5 14.2 11.1 8.8 36.8 6 22 6 6 22 14 17 268 67 250 345 2 0.089 0.052 0.051 0.126 0.177 0.65 2 6 0.081 0.162 0.012 0.413 0.489 1.81 6 0.90 0.83 0.31 1.48 1.79 7 0.97 2.22 0.05 5.51 5.85 22 6 0.080 0.070 0.008 0.192 0.480 1.78 28.7 15.4 103.0 49.5 14 27.1 18.2 10.3 72.5 114 40 86 142 Volume 2 8 16.2 10.9 8.8 36.8 14 0.082 0.111 0.008 0.413 1.146 1.41 8 0.96 1.91 0.05 5.51 7.64 14 104 Appendix Continued. TREES Plot Se 2A No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha 2C Stems Dbh Ht BA Volume (m) Age (yrs) Cr. d i a m (cm) (cm) (m ) (mS) 2 2 2 2 2 2 31.1 13.7 22.1 233 5.7 21.4 40.8 0.083 0.81 180 303 32 0.067 0.75 18.1 105 280 0.036 0.28 26.1 360 325 0.131 0.167 1.34 1.62 0.62 6 2 7 13 2 13 21.5 14.4 13 0.041 8.0 5.7 0.032 0.21 10.7 4.4 21.4 0.009 0.02 0.115 0.84 13 0.530 3.23 48 1.96 12 38.2 N o . of obs. Mean 4 St. dev. 0.014 3.8 0.007 Min. 10.0 0.008 Max. 17.6 0.024 Sum/plot 4 0.056 Total/ha 15 0.21 19 St. dev. 0.25 4 13.0 S t a n d 2 N o . o f obs. Mean 13 15 29 1 5 20.7 15.4 0.040 0.32 9.0 6.1 0.036 0 34 Min. 10.0 4.4 0.008 0 02 Max. 40.8 26.1 0.131 1 34 Sum/plot 19 Total/ha 23 0 7 5 3 0.93 4 8 5 9 105 Appendix Continued. Se TREES Stems Dbh Plot 3A 3B No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha a 3C x - N o . of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha All Stands (cm) 3 3 3 3 27.9 5.4 21.8 32.1 193 8 186 202 225 45 180 270 3 11 1 1 67.9 - 33.5 Volume (m3) 3 0.114 0.052 0.055 0.155 0.341 1.26 3 1.31 0.71 0.50 1.79 3.94 15 1 0.362 4.44 0 V362 4' 44 - 1 4 4 15 14 27.1 18.2 10.3 72.5 4 41 27 29.3 5.2 21.8 33.5 14 17 27.4 16.8 10.0 76.5 3 4 1 6 4 0.177 0.191 0.051 0.460 0.710 2.63 43.2 23.0 25.6 76.5 41 17 BA 1 1" 4 No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha Cr. diam (yrs) 37.2 9.5 26.5 44.4 N o . of obs. Mean St. dev. Min. ^Sum/plot TotaUha Ht (m) Age (cm) 17.7 8.9 4.4 36.8 7 7 182 90 86 360 268 57 180 345 14 0.082 0.111 0.008 0.413 1.146 1.41 4 41 27 0.081 0.106 0.008 0.460 3.311 1.36 2.10 1.67 0.50 4.44 8.38 16 0.77 1.33 0.02 5.51 20.87 13 106 Appendix Continued. SAPLINGS Plot BI 1A No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha IB IC No. of obs. Mean St. dev. Min. Max. Sum/plot TotaUha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha S t a n d 1 No. of obs. Mean St. dev. Min. Max. Sum/plots Total/ha Stems Dbh (cm) B. Diam. (cm) Ht (m) 183 148 183 62 147 183 4.1 6.2 2.8 68 112 150 2.7 0.2 9.9 2.8 1.6 13.0 1.3 1.3 7.6 20 31 109 34 45 290 63 25 330 183 677 348 3.0 2.5 0.2 9.9 348 1287 174 3.7 2.7 0.2 9.9 174 643 705 3.4 2.7 0.2 9.9 705 869 B H Age Base A g e C r . d i a m (yrs) (yrs) (cm) 107 Appendix Continued. SAPLINGS BI Stems Plot 2A 2B 2C N o . o f obs. Mean St. dev. Min. Max. Sum/plot TotaUha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha N o . of obs. Mean St. dev. Min. Max. Sum/plot Total/ha S t a n d 2 N o . of obs. Mean St. dev. Min. Max. Sum/plot Total/ha Dbh (cm) B . Diam. (cm) Ht (m) 192 172 192 4.4 2.6 0.0 9.9 192 710 301 4.2 2.7 0.2 9.8 301 1113 178 4.5 2.7 0.2 9.9 178 658 705 3.4 2.7 0.2 9.9 705 869 5.7 2.3 2.1 11.7 B H Age Base Age C r . d i a m (yrs) (yrs) (cm) 63 2.8 1.4 1.3 8.3 163 92 36 15 189 192 148 45 58 266 151 55 50 290 108 Appendix Continued. SAPLINGS BI Stems Plot 3A 3B 3C No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha S t a n d 3 N o . of obs. Mean St. dev. Min. Max. Sum/plot Total/ha All Stands B. Diam. (cm) Ht (m) 156 128 156 B H Age Base Age C r . d i a m (yrs) (yrs) (cm) 44 124 156 4.1 5.4 3.0 61 75 3.0 2.8 1.6 24 28 172 62 0.2 1.2 1.3 23 31 65 9.9 12.3 9.9 108 137 335 156 577 212 4.4 2.7 0.2 9.9 212 784 173 4.8 2.5 0.2 9.7 173 640 541 4.4 2.7 0.2 9.9 541 667 No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha Dbh (cm) 1917 1917 788 448 531 169 461 531 4.0 5.8 2.9 75 114 157 2.7 2.6 31 46 60 0.0 1.2 1.4 1.3 15 31 25 9.9 13.0 9.9 189 290 335 109 Appendix Continued. SAPLINGS Plot Se 1A N o . o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha IB IC Standi N o . o f obs. Mean St. dev. Min. Max. Sum/plot TotaUha N o . o f obs. Mean St. dev. Min. Max. Sum/plot Total/ha N o . o f obs. Mean St. dev. Min. Max. Sum/plots Total/ha Stems Dbh (cm) B. Diam. (cm) 7 B H Age Base Age Cr. d i a m (yrs) (yrs) (cm) 7 2.5 2.5 0.4 6.7 7 26 24 3.5 2.9 0.2 8.5 24 89 6 0.6 0.5 0.3 1.5 6 22 37 2.8 2.8 0.2 8.5 37 46 Ht (m) 1 2.3 1.2 1.3 4.3 5 44 - 7 74 30 31 104 104 51 40 190 110 Appendix Continued. SAPLINGS Plot Se 2A No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha 2B 2C No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha S t a n d 2 No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha Stems Dbh (cm) B. Diam. (cm) 8 B H Age Base Age Cr. d i a m (yrs) (yrs) (cm) 8 2.4 1.5 1.0 5.4 8 30 9 3.9 3.4 0.2 9.0 9 33 6 1.8 1.4 0.2 3.6 6 22 23 2.8 2.5 0.2 9.0 23 28 Ht (m) 3 2.0 0.4 1.6 2.8 6 30 19 8 44 8 58 20 36 83 124 53 75 235 Ill Appendix Continued. SAPLINGS Plot Se 3A No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha 3B 3C No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha Stand 3 No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha All Stands Stems Dbh (cm) 14 Ht (m) B H Age Base Age C r . d i a m (yrs) (yrs) (cm) 14 13 2.9 2.0 0.9 7.2 2.6 1.1 1.4 5.2 14 47 12 30 69 138 48 70 255 14 52 2.0 1.2 0.4 3.3 4 15 1 3.3 3.3 3.3 1 4 19 2.8 1.8 0.4 7.2 19 23 No. of obs. Mean St. dev. Min. Max. Sum/plot Total/ha B. Diam. (cm) 29 24 2.7 1.9 0.4 7.2 29 12 29 4.8 3.1 1.8 15.6 7 2.4 1.0 1.3 5.2 24 32 13 8 44 29 56 20 30 104 126 50 40 255
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- The structure and regeneration pattern of old-growth...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
The structure and regeneration pattern of old-growth stands in the Moist Cold Engelmann Spruce - Subalpine.. Varga, Pál 1997-12-31
pdf
Page Metadata
Item Metadata
Title | The structure and regeneration pattern of old-growth stands in the Moist Cold Engelmann Spruce - Subalpine Fir subzone of Central British Columbia |
Creator |
Varga, Pál |
Date | 1997 |
Date Issued | 2009-03-09T21:19:22Z |
Description | The need for sound management practices of high-elevation forests is becoming more and more important in British Columbia; however, there is a lack of information on the dynamics and regeneration pattern of these forests. This pilot study investigates stand dynamics by describing the structure of oldgrowth stands, and correlating the occurrence of regeneration with site factors. Three study stands, which originated after fire and appeared to be in advanced stages of succession were selected in the Moist Cold Engelmann Spruce - Subalpine Fir subzone, near Smithers, B.C. All three stands showed the same general stand structure, although the sites were slightly different. The number of seedlings, saplings, alive and standing dead trees of Engelmann spruce (Picea engelmannii Parry ex. Engelm.) was consistently very low, and, hence, excluded from the analysis. The diameter frequency distributions of subalpine fir (Abies lasiocarpa (Hook.) Nutt.) approximated a balanced distribution, but did not fit with the negative exponential function. The height frequency distributions of fir indicated the presence of at least two canopy strata. The age of the trees at breast height (1.3 m) ranged from 100 to 400 years. The age distributions did not conform to distributions expected in the old-growth stage. The age distribution of saplings indicated that in the past 100 years the number of established saplings has declined. The spatial pattern of firs changed from highly aggregated, to random, to regular, with increasing size (diameter or height). The spatial pattern of the oldest trees was aggregated possibly indicating remnants of the colonizing cohort. Since the diameter distributions and the spatial pattern of trees with different sizes are similar to that of a steady state they are not expected to change in the future. Canopy stratification was most pronounced in the oldest study stand, therefore, it is possible the height distribution will never show the inverse-J shape typically expected in old-growth forests. Due to weak size-age relationships the future change in age distribution is difficult to predict. The fir regeneration showed strong correlations with the presence of decaying wood substrate. There were more than twice as many seedlings on decaying wood than expected from the total cover of the substrate. However, survivorship from seedling to sapling stage was better on the forest floor than on decaying wood. Vegetation cover, light, and humus form showed no correlations with the number and vigour of seedlings. The occurrence of seedlings was also independent of canopy openings. This result implies that gap dynamics was not the driving force for the establishment of fir regeneration. One possible explanation is that snow melts later in canopy openings. The number of seedlings was higher than expected in the intermediate snow melt time zones which roughly corresponds to canopy edge positions. The future of the spruce component in the study stand was difficult to predict because of the very low number of individuals. The dominant position of spruce trees and presence of recently established seedlings throughout the stands would imply that there was an adequate seed source. However, without knowledge about survival of seedlings until they reach the canopy it is impossible to say whether the amount of regeneration is adequate for maintaining the spruce population. |
Extent | 4790495 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2009-03-09 |
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.0075198 |
URI | http://hdl.handle.net/2429/5796 |
Degree |
Master of Science - MSc |
Program |
Forest Sciences |
Affiliation |
Forestry, Faculty of Forest and Conservation Sciences, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 1997-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- ubc_1997-0113.pdf [ 4.57MB ]
- Metadata
- JSON: 1.0075198.json
- JSON-LD: 1.0075198+ld.json
- RDF/XML (Pretty): 1.0075198.xml
- RDF/JSON: 1.0075198+rdf.json
- Turtle: 1.0075198+rdf-turtle.txt
- N-Triples: 1.0075198+rdf-ntriples.txt
- Original Record: 1.0075198 +original-record.json
- Full Text
- 1.0075198.txt
- Citation
- 1.0075198.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 12 | 0 |
China | 6 | 19 |
Russia | 4 | 0 |
France | 1 | 0 |
Canada | 1 | 0 |
City | Views | Downloads |
---|---|---|
Ashburn | 10 | 0 |
Saint Petersburg | 4 | 0 |
Shenzhen | 4 | 19 |
Unknown | 2 | 46 |
Absecon | 2 | 0 |
Beijing | 2 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
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-0075198/manifest