
Citation: | Hayato Iijima, Takuo Nagaike. Susceptible conditions for debarking by deer in subalpine coniferous forests in central Japan[J]. Forest Ecosystems, 2016, 2(1): 33-33. DOI: 10.1186/s40663-015-0059-y |
In recent years, the population densities of large ungulates, especially deer species, have increased worldwide (Stewart and Burrows 1989; Fuller and Gill 2001; Rooney 2001; Apollonio et al. 2010; Iijima et al. 2013). An increase of deer density has been shown to result in more prevalent debarking of trees (Akashi and Nakashizuka 1999; Nagaike and Hayashi 2003; Iijima and Nagaike 2015) and browsing of saplings and understory vegetation (Gill and Beardall 2001; Beguin et al. 2009; Suzuki et al. 2013) in forests. To deal with the overabundance of deer and its effect on forest ecosystems, control of nuisance deer and the construction of fences to protect vegetation have been conducted. However, such countermeasures are difficult to perform in subalpine forests because of difficult access (e.g., steep slopes, complex topography, great distances from roads). Deer seasonally migrate to higher elevations, where delayed bud flush caused by low temperature provides fresh leaves even in mid-summer and they can escape from predators and more intense hunting pressure in low-elevation areas (Mysterud et al. 2011). These factors also have caused the number of deer in many higher elevation areas to increase (Takatsuki 2009), and debarking in subalpine forests has increased (Yokoyama et al. 2001; Takeuchi et al. 2011) in recent years. Therefore, factors correlated with debarking by deer in subalpine forests should be clarified to help effectively conserve these forests.
Many studies have examined the factors influencing debarking and browsing in lower elevation forests. Debarking was shown to depend on the size of trees (e.g., Nagaike and Hayashi 2003; Koda and Fujita 2011; Borkowski and Ukalski 2012), their species (e.g., Kay 1993; Akashi and Nakashizuka 1999; Moore et al. 1999; Takeuchi et al. 2011), proportion of coniferous stands (Ligot et al. 2013), the distance from forest road (McLaren et al. 2000), and snow depth (Iijima and Nagaike 2015). In addition to these factors, spatial variation of deer impact was observed in lower elevation forests: the higher the deer density, the higher the proportion of debarked trees (Iijima and Nagaike 2015) and browsed saplings (Akashi et al. 2011) and the lower the sapling density (Beguin et al. 2009). The spatial variation of deer density across low-elevation areas was explained by the presence of attractive habitat (e.g., artificial grassland; Kamei et al. 2010; Iijima et al. 2013). However, the attractive habitat in subalpine forests, which are located on steeper slopes and have more complex topography than lower elevation areas, has not been well studied.
Takeuchi et al. (2011) reported that elevation was an important factor affecting expansion of deer density into subalpine zones. In addition to elevation, the amount of attractive habitat (e.g., subalpine grassland) may affect deer density. Deer generally prefer herbaceous plants to woody plants (Takatsuki 1986; Winnie 2012), and severe browsing of herbaceous plants in subalpine grasslands by deer has been reported (Schütz et al. 2003; Nagaike 2012; Nagaike et al. 2014). Thus, subalpine forests which were surrounded by large subalpine grasslands would be more susceptible for debarking by deer. Thus far, however, direct evaluation of the relationship between these factors and the intensity of debarking in subalpine forest is rare.
The objective of this study was to clarify the conditions which are correlated with debarking by sika deer in subalpine forests in order to improve the effective management and conservation of these subalpine forests. We hypothesized that the occurrence of debarking and sapling density would increase and decrease in low altitude and abundant grassland.
The present study was conducted in subalpine forests in Minami-Alps National Park (357.5 km2) in central Japan (Fig. 1), where sika deer (Cervus nippon) density has been increasing in recent years (Izumiyama and Mochizuki 2008; Izumiyama et al. 2009). The park is characterized by steep slopes and a concentration of many mountains within a small area (Fig. 1). Subalpine coniferous forests dominate from 2000 to 2500 m elevation, mainly composed of Abies mariesii, Abies veitchii, Picea jezoensis var. hondoensis, and Tsuga diversifolia (Fig. 2). Subalpine grasslands consisting of herbaceous species and dwarf tree community (Pinus pumila) are present at elevations higher than 2500 m, and mixed conifer-broadleaf forests are present below 2000 m. Sika deer are widely distributed in Japan. Their shoulder heights range from 90 to 190 cm and they have a highly varied diet (they eat graminoids in northern Japan and eat leaves and fruits in southern Japan; Ohdachi et al. 2009). They can reach bark and leaves up to ca. 2 m in height.
We established 24 rectangular plots (each of 10 m×40 m) in natural coniferous forests at elevations of 2000 m (low-elevation zone) and 2500 m (high-elevation zone) and surveyed in September 2011 (Fig. 1). For all trees with the diameter at breast height (DBH) > 3 cm in the plots, we noted the species and measured DBH and the presence of debarking. Furthermore, for saplings survey, 8 subplots (each of 1 m×2 m) were established at 10-m intervals along the two long sides of each plot (8 subplots per a plot). Height and the occurrence of debarking on saplings > 30 cm in height and < 3 cm in DBH were measured and recorded in each subplot. In this study, we could not determine the timing of the occurrence of debarking and did not distinguish new and old debarking.
The location of subalpine grassland in this region was obtained from the Natural Environment Information geographic information system (http://www.biodic.go.jp/trialSystem/top.html, accessed 20 April 2015) provided by the Biodiversity Center of Japan (Fig. 1). The minimum distance from each plot to subalpine grassland was calculated by QGIS ver 2.8.1 (Quantum GIS Development 2015).
Sapling density and the occurrence of debarking of saplings and trees were analyzed using a generalized linear mixed model (GLMM). For all models, we adopted each plot as random effect.
For analysis of sapling density, minimum distance from each plot to subalpine grassland, elevation zone, tree density in each plot (indicating understory light conditions), and species of each sapling were used as explanatory variables. We used the number of saplings as response variable. The means of the number of saplings and sapling density were equal because surveyed areas of the number of saplings were same for all plots. We also added basal area (BA) of conspecific trees (as an index of seed source) as an explanatory variable because the abundance of the seed source may affect sapling density. A Poisson distribution was adopted as the error distribution with a log-link function. Picea jezoensis was omitted from the species data because none of the surveyed individuals had been debarked, which made it impossible to estimate the coefficient.
For analysis of the occurrence of debarking of saplings, minimum distance from each plot to subalpine grassland, elevation zone, species, and sapling height were included as explanatory variables. A binomial distribution was adopted as the error distribution with a logit link function.
For analysis of the occurrence of debarking of trees, minimum distance from each plot to subalpine grassland, elevation zone, DBH, and species were used as explanatory variables. A binomial distribution was adopted as the error distribution with a logit link function.
All analyses were conducted using R ver 3.1.3 (R Core Team 2015) and the glmmML package (Brostrom and Holmberg 2011). Model selection was performed using the Akaike information criterion (AIC, Akaike 1973) in conjunction with a backward elimination procedure. We determined that explanatory variables used in conjunction with all models which had delta AIC compared to the lowest AIC < 2 are significant and explanatory variable used in conjunction with any models which had delta AIC compared to the lowest AIC < 2 are marginally significant. At the same time, we calculated Akaike weight (w) and the relative importance of variable (Burnham and Anderson 2002) for comparing the importance of each variable. Akaike weight is defined as
wi=exp(−Δi2)∑Rr=1exp(−Δi2) |
where Δi=AICi -AICmin is the difference between an AIC of each model and the minimum AIC among all candidate models (including null model and full model). This value, referred to as the Akaike weight, provides a relative weight of evidence for each model. The relative importance of predictor variable can be calculated as the sum of the Akaike weights over all of the models in which the parameter of interest appears. The relative importance range from 0 (least important) to 1 (most important).
The percentages of debarked live trees differed among plots (Fig. 3). The differences among sapling types was large, with dead saplings being debarked at very high rates (mean=81.6 %; Fig. 3).
Debarking of saplings was significantly related to a taller height of saplings and lower elevation zone, marginally related to minimum distance to grassland, and was not affected by species (Table 1). Sapling density significantly increased with increasing BA of conspecific trees and increasing minimum distance to grassland and marginally increased with increasing elevation and decreasing tree density (Table 2). Sapling density of A. mariesii was significantly higher than that of other species (Table 2). Debarking of trees was significantly related to small DBH and low elevation zone, and marginally decreased with long distance to grassland. Debarking of trees marginally differed by species (Table 3). Overall, Size of saplings and trees had strong effects on the occurrence of debarking (Table 4). Elevation zone, minimum distance to grassland, and species were marginal factors of the occurrence of debarking of saplings and trees and sapling density (Table 4).
AIC | ΔAIC | Elevation zonea 2500 m | Dist. grassland | Height of saplings | Speciesb, c | |
Abies veitchii | Tsuga diversifolia | |||||
398.028 | 0.000 | −2.634 | 0.037 | |||
399.976 | 1.948 | −2.619 | −1.269×10−4 | 0.037 | ||
aCoefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. bCoefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0. c Picea jezoensis was excluded from the analysis of saplings because no P. jezoensis saplings were debarked, so it was impossible to estimate the coefficient |
AIC | ΔAIC | Elevation zonea 2500 m | BAb of conspecific trees | Dist. grassland | Speciesc | Tree density | ||
Abies veitchii | Picea jezoensis | Tsuga diversifolia | ||||||
1382.409 | 0.000 | 0.048 | 9.310×10−4 | −0.555 | −2.675 | −1.237 | ||
1382.624 | 0.214 | 0.745 | 0.048 | 9.561×10−4 | −0.559 | −2.676 | −1.236 | −0.018 |
1383.631 | 1.272 | 0.210 | 0.048 | 9.372×10−4 | −0.557 | −2.680 | −1.239 | |
1384.286 | 1.877 | 0.048 | 9.311×10−4 | −0.554 | −2.672 | −1.239 | −0.002 | |
aCoefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. bBA, basal area. cCoefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0 |
AIC | Δ AIC | Elevation zonea 2500 m | DBHb of trees | Dist. grassland | Speciesc | ||
Abies veitchii | Picea jezoensis | Tsuga diversifolia | |||||
526.047 | 0.000 | −0.684 | −0.014 | −4.146×10−4 | |||
526.261 | 0.214 | −0.771 | −0.014 | −4.196×10−4 | 0.127 | −1.379 | −0.947 |
527.907 | 1.856 | −0.592 | −0.014 | ||||
527.974 | 1.926 | −0.666 | −0.014 | 0.076 | −1.453 | −0.959 | |
aCoefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. bDBH, diameter at breast height. cCoefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0 |
Debarking of saplings | Sapling density | Debarking of trees | |
Elevation zone | 1.000 | 0.506 | 0.785 |
BA of conspecific species | Not included | 1.000 | Not included |
DBH of trees | Not included | Not included | 1.000 |
Dist. grassland | 0.274 | 1.000 | 0.674 |
Height of saplings | 1.000 | Not included | Not included |
Species | 0.137 | 1.000 | 0.463 |
Tree density | Not included | 0.457 | Not included |
We showed the relationship between debarking of trees and saplings in subalpine coniferous forests and various factors. The percentages of debarked live saplings were higher than those of live trees, and the percentages for dead saplings were much higher than those for live saplings (Fig. 3). These results suggest that debarking by deer first occurs on saplings.
Tree and sapling sizes affected the occurrence of debarking but the effects of sizes showed opposite patterns for trees and saplings (Tables 1 and 3). Higher percentages of debarked small trees were also reported in many previous studies (Akashi and Nakashizuka 1999; Yokoyama et al. 2001; Nagaike and Hayashi 2003; Jiang et al. 2005; Kiffner et al. 2008; Takeuchi et al. 2011; Borkowski and Ukalski 2012), although the reasons were not clarified. One possibility for the selective debarking of small trees is the relative ease of debarking smaller trees (Ando et al. 2004). In contrast to trees, the relationship between sapling size and the occurrence of debarking has rarely been examined. Takatsuki and Gorai (1994) compared the size structure of trees between two distinct areas in Japan that differed in deer density. They showed that forests with high deer density lacked saplings > 90 cm in height, whereas there were abundant saplings shorter than 90 cm in forests both under high and low deer densities. In the present study, one reason for the low debarking intensity on small saplings (Table 1) is that they are more difficult to find than large saplings. Consequently, saplings > 30 cm in height and < 3 cm in DBH might be more susceptible to debarking by deer, suggesting that these might be a possible optimal size for debarking by sika deer.
P. jezoensis saplings appear to be less palatable for sika deer because no P. jezoensis saplings was debarked and we found marginal effect of species on debarking of trees (Table 3). Generally, Abies species were preferentially debarked over co-occurring Picea species (Heuze et al. 2005), and the different palatability among tree species suggests that deer debarking may alter the species composition of forests. However, P. jezoensis has been heavily debarked in some areas in Japan (Yokoyama et al. 2001; Ando et al. 2003), indicating that palatability of debarking may be affected by surrounding vegetation (Moore et al. 1999; McLaren et al. 2000). Therefore, P. jezoensis saplings at the sites examined in the present study may be debarked in the future if the deer density increases.
In addition to size and species as factors for debarking, elevation zones and grassland area also were associated with the damage to subalpine forests by deer. Low sapling density and high risk of debarking of trees and saplings were found in low elevation plots (Table 2). Other studies in Japan also found the intensities of debarking on trees and saplings in subalpine coniferous forests to be greater at lower elevations (Takeuchi et al. 2011). In our study area, sika deer migrate to high-elevation areas in summer and to low-elevation areas from late autumn to late spring (Izumiyama and Mochizuki 2008) although there are some deer in low-elevation area even in summer. Then, opportunities of debarking for saplings and trees were larger in low-elevation zone than in high-elevation because deer stay longer time in low elevation-area than in high elevation-area. In addition to elevation, we found that the large subalpine grassland caused low sapling density (Table 2) and high risk of the occurrence of debarking of saplings and trees (Tables 1 and 3). Deer would likely be attracted to subalpine grassland because it has numerous herbaceous plants that serve as rich food sources (Takatsuki 1986; Winnie 2012). Trees and saplings in subalpine forests near subalpine grassland would be eaten by deer before sprout of grasses in early spring and after defoliation of grasses in late autumn. Some sika deer migrate to high-elevation area before sprout of grasses and after defoliation to escape the competition with other deer in low-elevation areas (Shin-ichiro Hamasaki, personal communication). Minimum distance to subalpine grassland, however, had only a marginal effect on the occurrence of debarking of saplings and trees (Tables 1 and 3). As the result of intensive debarking on saplings (Fig. 3), low sapling density near subalpine grassland occurred (Table 2) and then it might be obscure the relationship between the distance to grassland and the occurrence of debarking of saplings. In the future, the intensity of debarking of tree would increase if deer density remains at this level or increases.
Our analysis revealed that large saplings within the range of 30 cm in height and 3 cm in DBH in subalpine forests were most susceptible for debarking by sika deer. In addition to the size, saplings in low elevation zone or near subalpine grasslands were more susceptible to debarking by sika deer. This information may be useful to improve monitoring deer impact, or determining where to construct deer fences in the subalpine region. At the same time, deer population control in the subalpine region should be conducted in forests of low elevation zone or near subalpine grasslands.
[1] |
Akaike H (1973) Information theory and an extension of the maximum likelihood principle, in 2nd International Symposium on Information Theory. Tsahkadsor, Armenian SSR, pp 267-281
|
[2] | |
[3] | |
[4] | |
[5] | |
[6] |
Apollonio M, Andersen R, Putman R (2010) European ungulates and their management in the 21st century. Cambridge University Press, NY.
|
[7] |
Beguin J, Pothier D, Prévost M (2009) Can the impact of deer browsing on tree regeneration be mitigated by shelterwood cutting and strip clearcutting? For Ecol Manage 257: 38-45
|
[8] | |
[9] | |
[10] |
Burnham KP, Anderson DR (2002) Model selection and multimodel inference: a practical information-theoretic approach, 2nd edn. Springer, New York
|
[11] | |
[12] | |
[13] | |
[14] | |
[15] |
Iijima H, Nagaike T (2015) Appropriate vegetation indices for measuring the impacts of deer on forest ecosystems. Ecol Ind 48:458-463
|
[16] |
Izumiyama S, Mochizuki T (2008) Seasonal range use of sika deer, which inhabits the subalpine zone in the southern Japan Alps. Bull Shinshu Univ Alps Field Center 6:25-32 (in Japanese with English summary)
|
[17] |
Izumiyama S, Mochizuki T, Takii A (2009) GPS tracking of Sika deer which inhabits the sub-alpine zone in the Southern Japan Alps. Bull Shinshu Univ Alps Field Center 7:63-71 (in Japanese with English summary)
|
[18] | |
[19] | |
[20] | |
[21] |
Kiffner C, Rößiger E, Trisl O, Schulz R, Rühe F (2008) Probability of recent bark stripping damage by red deer (Cervus elaphus) on Norway spruce (Picea abies) in a low mountain range in Germany: a preliminary analysis. Silva Fennica 42:125-134
|
[22] | |
[23] | |
[24] | |
[25] | |
[26] |
Mysterud A, Loe LE, Zimmermann B, Bischof R, Veiberg V, Meisingset E (2011) Partial migration in expanding red deer populations at northern latitudes: A role for density dependence? Oikos 120: 1817-1825
|
[27] | |
[28] | |
[29] |
Nagaike T, Ohkubo E, Hirose K (2014) Vegetation recovery in response to the exclusion of grazing by sika deer (Cervus nippon) in seminatural grassland on Mt. Kushigata, Japan. ISRN Biodiversity: Article ID 493495
|
[30] |
Ohdachi SD, Ishibashi Y, Iwasa MA, Saitoh T (2009) The Wild Mammals of Japan. Shoukadoh Book Sellers, Kyoto
|
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] |
Takatsuki S (1986) Food habits of sika deer on Mt. Goyo, northern Honshu. Ecol Res 1:119-128
|
[38] | |
[39] | |
[40] | |
[41] | |
[42] |
1. | Yuichiro Hiraoka, Naoyuki Nishimura, Yasuhiro Koyama, et al. Effects of Sika Deer on Sub-alpine Coniferous Forest Dynamics in Northern Yatsugatake Based on a 20-Year Investigation. Journal of the Japanese Forest Society, 2023, 105(6): 216. DOI:10.4005/jjfs.105.216 |
2. | Akiko Takii, Masaaki Ozeki, Chihiro Takahata, et al. Habitat selection of large herbivores evidenced as threats to alpine ecosystem. Acta Oecologica, 2022, 114: 103812. DOI:10.1016/j.actao.2022.103812 |
3. | Hayato Iijima, Tetsuya Maruyama, Hiroyuki Sakaniwa, et al. The Similarity and Difference of Factors that Determine the Intensity of Debarking by Sika Deer:. Journal of the Japanese Forest Society, 2021, 103(5): 344. DOI:10.4005/jjfs.103.344 |
4. | Takuo Nagaike. Effects of heavy, repeated bark stripping by Cervus nippon on survival of Abies veitchii in a subalpine coniferous forest in central Japan. Journal of Forestry Research, 2020, 31(4): 1139. DOI:10.1007/s11676-019-00940-x |
5. | Kanoko Takarabe, Hayato Iijima. Abundant artificial grasslands around forests increase the deer impact on forest vegetation. European Journal of Forest Research, 2020, 139(3): 473. DOI:10.1007/s10342-020-01262-y |
6. | Judith S. Nantongo, Brad M. Potts, Hugh Fitzgerald, et al. Quantitative Genetic Variation in Bark Stripping of Pinus radiata. Forests, 2020, 11(12): 1356. DOI:10.3390/f11121356 |
7. | Takuo Nagaike. Bark Stripping by Deer Was More Intensive on New Recruits than on Advanced Regenerants in a Subalpine Forest. Forests, 2020, 11(5): 490. DOI:10.3390/f11050490 |
8. | Kanoko Takarabe, Hayato Iijima. Contrasting Effect of Artificial Grasslands on the Intensity of Deer Browsing and Debarking in Forests. Mammal Study, 2019, 44(3): 173. DOI:10.3106/ms2018-0082 |
9. | Hiroyuki Uno, Yoshihiro Inatomi, Mayumi Ueno, et al. Effects of sika deer (Cervus nippon) and dwarf bamboo (Sasa senanensis) on tree seedlings in a cool-temperate mixed forest on Hokkaido Island, Japan. European Journal of Forest Research, 2019, 138(6): 929. DOI:10.1007/s10342-019-01214-1 |
10. | Hayato Iijima. Seasonal Change of Deer Occurrence and Damage of Plant Biomass in the Mosaic Landscape of Artificial Grasslands and Forests. Mammal Study, 2018, 43(2): 91. DOI:10.3106/ms2017-0087 |
11. | Hayato Iijima. The Effects of Landscape Components, Wildlife Behavior and Hunting Methods on Hunter Effort and Hunting Efficiency of Sika Deer. Wildlife Biology, 2017, 2017(1): 1. DOI:10.2981/wlb.00329 |
12. | Hayato Iijima, Takuo Nagaike. The factors that determine the intensities of deer browsing and debarking on broadleaf tree around artificial grasslands. Journal of Forest Research, 2017, 22(3): 199. DOI:10.1080/13416979.2017.1305262 |
13. | Hayato Iijima, Mayumi Ueno. Spatial heterogeneity in the carrying capacity of sika deer in Japan. Journal of Mammalogy, 2016, 97(3): 734. DOI:10.1093/jmammal/gyw001 |
14. | Haruka Ohashi. Sika Deer: Life History Plasticity and Management. Ecological Research Monographs, DOI:10.1007/978-981-16-9554-4_3 |
15. | Akiko Takii, Shigeyuki Izumiyama. Sika Deer: Life History Plasticity and Management. Ecological Research Monographs, DOI:10.1007/978-981-16-9554-4_10 |
AIC | ΔAIC | Elevation zonea 2500 m | Dist. grassland | Height of saplings | Speciesb, c | |
Abies veitchii | Tsuga diversifolia | |||||
398.028 | 0.000 | −2.634 | 0.037 | |||
399.976 | 1.948 | −2.619 | −1.269×10−4 | 0.037 | ||
aCoefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. bCoefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0. c Picea jezoensis was excluded from the analysis of saplings because no P. jezoensis saplings were debarked, so it was impossible to estimate the coefficient |
AIC | ΔAIC | Elevation zonea 2500 m | BAb of conspecific trees | Dist. grassland | Speciesc | Tree density | ||
Abies veitchii | Picea jezoensis | Tsuga diversifolia | ||||||
1382.409 | 0.000 | 0.048 | 9.310×10−4 | −0.555 | −2.675 | −1.237 | ||
1382.624 | 0.214 | 0.745 | 0.048 | 9.561×10−4 | −0.559 | −2.676 | −1.236 | −0.018 |
1383.631 | 1.272 | 0.210 | 0.048 | 9.372×10−4 | −0.557 | −2.680 | −1.239 | |
1384.286 | 1.877 | 0.048 | 9.311×10−4 | −0.554 | −2.672 | −1.239 | −0.002 | |
aCoefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. bBA, basal area. cCoefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0 |
AIC | Δ AIC | Elevation zonea 2500 m | DBHb of trees | Dist. grassland | Speciesc | ||
Abies veitchii | Picea jezoensis | Tsuga diversifolia | |||||
526.047 | 0.000 | −0.684 | −0.014 | −4.146×10−4 | |||
526.261 | 0.214 | −0.771 | −0.014 | −4.196×10−4 | 0.127 | −1.379 | −0.947 |
527.907 | 1.856 | −0.592 | −0.014 | ||||
527.974 | 1.926 | −0.666 | −0.014 | 0.076 | −1.453 | −0.959 | |
aCoefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. bDBH, diameter at breast height. cCoefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0 |
Debarking of saplings | Sapling density | Debarking of trees | |
Elevation zone | 1.000 | 0.506 | 0.785 |
BA of conspecific species | Not included | 1.000 | Not included |
DBH of trees | Not included | Not included | 1.000 |
Dist. grassland | 0.274 | 1.000 | 0.674 |
Height of saplings | 1.000 | Not included | Not included |
Species | 0.137 | 1.000 | 0.463 |
Tree density | Not included | 0.457 | Not included |