Susceptible conditions for debarking by deer in subalpine coniferous forests in central Japan

Recently, deer have expanded their distribution to higher altitude ranges including subalpine forests. However, culling deer and construction of deer fence in subalpine forests are difficult because of steep slopes and complex topography. Thus it is necessary to clarify the factors which are associated with debarking by deer for the effective protection of subalpine forests. In this study, we examined which factors are associated with debarking by sika deer (Cervus nippon) in subalpine coniferous forests. We conducted our survey in Minami-Alps National Park, central Japan. We established 24 10 m × 40 m plots and surveyed the occurrence of debarking on saplings >30 cm in height and <3 cm in diameter at breast height (DBH) and on trees >3 cm in DBH, as well as sapling density within each plot. Minimum distances to nearest grassland of plots were calculated (tentatively assuming grassland would attract deer and would cause high debarking pressure in the surrounding subalpine forests). The mean percentage of debarked live saplings was higher than that of live trees. The mean percentage of debarked saplings which had already died was 81.6 %. Debarking of saplings increased with lower elevation, taller sapling size, and marginally increased near grassland. Sapling density was lower in plots with low basal area of conspecific trees near grassland and differed among species. Sapling density marginally decreased with decreasing elevation and increasing stand tree density. Debarking of trees was positively related to small DBH and low elevation, and marginally increased near grassland and differed among species. Our results suggest that tall saplings in subalpine forests of low elevation or near subalpine grassland were susceptible to debarking by deer and monitoring of these areas may permit the early detection of the impacts of deer in subalpine coniferous forests.


Background
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.

Study area
The present study was conducted in subalpine forests in Minami-Alps National Park (357.5 km 2 ) 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.

Field surveys
We established 24 rectangular plots (each of 10 m × 40 m) in natural coniferous forests at elevations of 2000 m (lowelevation 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 minimum distance from each plot to subalpine grassland 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).

Statistical analysis
Sapling density and the occurrence of debarking of saplings and trees were analyzed using a generalized linear . 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

Results
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).

Discussion
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. Coefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. b Coefficient 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 Coefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. b BA, basal area. c Coefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0 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 highelevation 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.

Conclusions
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.

Competing interests
This study was funded by the Comprehensive Research Organization for Science and Technology of Yamanashi Prefectural Government and Mitsui & CO., LTD. The authors declare that they have no competing interests about non-financial aspects.
Authors' contributions HI and TN substantially contributed about setting the experimental design and collecting data. HI and TN deeply discussed and approved the content of manuscript. Both authors read and approved the final manuscript. Coefficient of elevation zone was calculated as a relative value when the coefficient of the low-elevation zone (2000 m) was set to 0. b DBH, diameter at breast height. c Coefficient of species was calculated as a relative value when the coefficient of Abies mariesii was set to 0