Habitat selection and their interspecific interactions for mammal assemblage in the Greater Khingan Mountains, northeastern China

Patterns of habitat use profoundly influence interactions among wildlife species and ecological communities, the sustainability of species and the stability of populations. The Eurasian lynx Lynx lynx, sable Martes zibellina, moose Alces alces, roe deer Capreolus pygargus and mountain hare Lepus timidus are sympatric carnivore and herbivore species in Hanma Nature Reserve in the Greater Khingan Mountains, northeastern China. We conducted snow track sample plot surveys of these mammals during the winter of 2012–2013 and analyzed habitat selection characteristics and interspecific interactions. We screened the preferred habitat variables of each species and predicted their potential occurrence probability by generalized linear mixed model (GLMM) and further, selected the best model of each species that incorporated other species' potential occurrence probability in GLMM. Our results showed that when we only considered habitat factors in the model: 1) river valley was the most favored habitat factor in winter that lynx, moose and mountain hare showed strong selection on; 2) lynx preferred areas with steep slope and coniferous forest; 3) sable only avoided the forest edge; 4) roe deer avoided sparse shrub forest but preferred coniferous forest. When we incorporated other species' presence probability in GLMM: 1) the presence probability of moose and roe deer had no significant relationship with other species; 2) lynx preferred areas with higher mountain hare potential occurrence probability; 3) sable preferred areas with higher mountain hare but less lynx occurrence probability; 4) mountain hare preferred to inhabit areas with sympatric roe deer and lynx. These results will provide guidance for species specific habitat conservation and restoration, and wildlife population management based on interspecific interactions.

Animals make habitat selection based on the tradeoff between the costs and benefits perceived by the animals Dill 1990, Lima andBednekoff 1999). Wildlife habitat selection is driven by habitat-related factors, such as forage quality and availability, shelter, predation pressure and mating system (Creel et al. 2005, Godvik et al. 2009). The most common tradeoff occur on many large sized herbivores is that closed habitats provide shelter against harsh weather and/or predators, while exposed habitats provide better forage (Godvik et al. 2009). What's more, prey density is another important factor affecting habitat selection for carnivores. Differences in the body sizes of sympatric herbivores may mean they prefer forage of different heights and quality (Kuiters et al. 2005) and for predators, body sizes were important factors deciding their prey species, and their different roles in regulating ecosystems at the top of food net (Beschta andRipple 2009, Ritchie et al. 2012).
Much work has been done on wildlife habitat selection, correlations between mammals and climate, resource availability, habitat diversity, environment factors, sympatric species, predators and human interference (Batcheler 1960, Menéndez et al. 2007, Bjørneraas et al. 2011, Nicholson et al. 2014. Researches showed that harsh winter pushed animals to choose places with low snow cover and more abundant forage (Dussault et al. 2005, Street et al. 2015, such as river basins (Nicholson et al. 1997). When wildlife use habitat, it is critical that they adjusts their behaviors to adapt to the environment to ensure viability. The coexistence mechanisms of multiple sympatric species, how animals adjust to intraspecific and interspecific competition, and how environment resources are utilized to avoid strong competition are well studied (Kuiters et al. 2005, Jiang et al. 2010), but geographic gaps remain. The Eurasian lynx Lynx lynx, sable Martes zibellina, moose Alces alces, roe deer Capreolus pygargus and mountain hare Lepus timidus are sympatric predator and herbivore species in Hanma Nature Reserve in the Greater Khingan Mountains, northeastern China. Until recently, few researches have been carried out in this area, we tried to find factors affecting habitat selection of these species and explore possible explanations for these effects.
Environmental factors can directly regulate species richness by providing caves, nests and shelter for animals (Wright 1983). Water-energy aspects of the environment can constrain species richness via a trophic cascade: water-energy availability determines plant richness, which determines herbivore richness and in turn carnivore richness (Huston 1994, Mittelbach et al. 2001. Vegetation has direct effects on herbivore richness and abundance, so herbivore species diversity is often closely correlated with plant diversity (Siemann et al. 1998), and grazing can alter vegetation configuration and patch heterogeneity (Adler et al. 2001). Exploring driving factors for habitat selection of sympatric predators and herbivores, reactions of herbivores on habitat utilization under predation pressure, and vegetation-herbivore-predator relationships can help to deepen our understanding of herbivore and predator coexistence mechanism. Hanma Nature Reserve in the Greater Khingan Mountains, northeastern China is such an optimum area to stress this problem, where lynx L. lynx, sable M. zibellina, moose A. alces, roe deer C. pygargus and mountain hare L. timidus coexist. In this research, snow track sampling method and generalized linear mixed model (GLMM) were used to examine the characteristics of habitat selection of both herbivores and carnivores and possible interspecific interactions in this area to test the hypotheses that 1) different environment factors determine habitat selection and distribution of carnivores and herbivores 2) predators prefer areas where herbivores occur more frequently.

Study area
Hanma National Nature Reserve (51°20′02″-51°49′48″N, 122°23′34″-122°52′46″E) is located in the Greater Khingan Mountains, Inner Mongolia, northeastern China ( Fig. 1). The reserve covers an area of 1073 km 2 with a subfrigid climate: a cold dry season from October to June and a warm wet season from July to September. The annual mean temperature is -5.3°C; the highest temperature in summer is 35.4°C and the lowest temperature in winter is -49.6°C. Average rainfall of 450 mm is concentrated between July and September (70% of yearly precipitation) and the snow depth is usually around 30 cm in harsh winter. At higher elevations, the dominant tree species is Siberian dwarf pine Pinus pumila, while lower place are characterized by Dahuarian larch Larix gmelinii, and deciduous trees are white birch Betula platyphylla, willow Salix spp., hazel Corylus avellane and Siberian alder Alnus hirsuta var. sibirica, while shrub berry, reed grass and moor grass dominate the forest floor of the thickets. We divided vegetation into four types: coniferous, deciduous, sparse shrub and swamp patches. Coniferous forest is the main habitat type and accounts for 82% of total vegetation cover. Hanma Nature Reserve is far away from towns and villages and landscape keeps quite natural because no logging has ever been carried out before, hence, there are no road-net and few human activities exist in this area.

Data collection
To collect species presence data, we carried out sample plot surveys (Qi et al. 2015) during the winter of 2012-2013. Vegetation type depended stratified sampling methods were used to design transect lines, and 55 transects with a total length of 275 km were set up. We used ArcGIS 9.3 to design the line transect distribution: five transect lines composed of one sample plot, and each transect was 5 km in length and the interval between two transects was 500 m. Each transect line was surveyed by two experienced local workers and they recorded animal tracks in the snow left in 24 h and corresponding location coordinates. When field work was completed, we used ArcGIS 9.3 to acquire 25 points on each transect line with an interval of 200 m, and extracted slope, aspect, elevation, distance to the edge of the nearest coniferous patch (hereafter conifer), distance to the center of the nearest swamp patch (hereafter swamp), distance to the center of the nearest sparse shrub patch (hereafter shrub), distance to the center of the nearest deciduous patch (hereafter deciduous), distance to the nearest river (hereafter river), and distance to the nearest patch edge (Table 1) as environmental indicators for species habitat selection modeling processing. In addition, we attributed the presence information of mammal species (i.e. tracks) to the nearest points.

Data analysis
We normalized all the habitat variables in order to avoid the condition that some variables were too skewed, and for the species presence data, we used binomial distribution (presence  1 and absence  0). Pearson's correlation matrix was used to identify whether problematic collinearities exist among covariates (Loyn et al. 2001), which suggested that the deciduous variable had problematic collinearities with swamp (r  0.57), so we deleted it to avoid autocorrelation problem in our analysis.
For each species, we used logistic regression to build a linear combination of terrain and vegetation variables, through stepwise regression method; we selected the most parsimonious model that best explained their occurrence and distribution. We used Akaike's information criterion difference adjusted for small sample sizes (AICc) and Akaike model weights (w i ) to get the most parsimonious model (i.e. fewest variables to explain the most variation). After we knew each species' optimal model, we used the model to predict their potential presence probability (P i ) considering only the terrain and vegetation variables.
When examining the interactions among the five species, we defined each species a second model: for carnivores, to incorporate prey species' potential presence probability (P i ) as a variable to its first habitat selection model, and for herbivores, to incorporate potential presence probability (P i ) of predator and other herbivore competitors. And once again, we used Akaike's information criterion (AICc) and Akaike model weights (w i ) to evaluate and identify models and get the most parsimonious model for each species, so that, we got the second most parsimonious model to predict their secondary presence probability (P i ).Generalized linear mixed model analysis was carried out in R ver. 3.2.2 using the 'lme4' package (< www.r-project.org >) and the 11 sample plots (with five line transects each) were included as a random factor in the analyses.

Results
We extracted 1300 points from all transect lines of 11 sample plots for analysis. For lynx, we found 35 occurrence points (one point had two occurrences); 116 points for moose (two occurrences at two sites); 86 points for roe deer (two occurrences at seven sites); 118 points for mountain hare (two occurrences at seven sites); and sable were recorded at 200 points (two occurrences at 14 sites).

Model selection
We found the five top generalized linear mixed models when only considering the effect of terrain and vegetation variables, and found three top models (without moose and roe deer) when incorporating other species' presence probability (P i ) in the second model, and calculated w i for each model. We presented the best supported generalized linear mixed models (i.e. those with ΔAIC  2) for the five species (Table 2).
On the first condition, when only environmental factors considered, we found five top generalized linear mixed models for the five species separately. For the moose, the most parsimonious resource model included two variables (river and random covariates; no. of model parameters K  2, w i  0.205, Table 2). Roe deer's most parsimonious model consisted of three variables (shrub, coniferous and random covariate; K  3, w i  0.305, Table 2). Mountain hare's most parsimonious resource model included two variables (river and random covariate; K  2, w i  0.246, Table 2). Lynx's most parsimonious model had four variables (river, coniferous, Note: Data was processed in ArcGIS 9.3 in grid format. Data extracted from the stock map was provided by Hanma Nature Reserve.

Effect of habitat factors on animal presence
Based on the top generalized linear mixed models selected for each species, we found that the most parsimonious resource selection model for moose included only one variable, river (Table 2), which meant that moose showed strong preference for areas near the rivers. More previous researches showed that moose were typical browsers and their habitat selection was affected by seasonal forage availability (Cederlund and Okarma 1988, Olsson et al. 2011, Van Beest et al. 2010, time of day, sex and reproductive status of the cow (Bjørneraas et al. 2011(Bjørneraas et al. , 2012. In harsh winter, moose strongly selected young forest and wetter areas to obtain sufficient food resources, such as willow Salix spp., hazel Corylus avellane, and Siberian alder Alnus hirsuta var. sibirica that provided shoots to moose as winter food (Jiang et al. 2009). Moose also used watersheds in the summer as escape habitat and corridors for crossing between habitat patches in winter (Nicholson et al. 2014). In our study area, river valley is the most heterogeneous area, where plants such as willow Salix spp., Asian white birch Betula platyphylla, hazel Corylus avellane and Siberian alder A. hirsuta var. sibirica as well as many other broadleaf and shrub that provide sufficient annual new shoots as winter forage for moose. Moreover, in winter, the Siberian dwarf pine Pinus pumila, which is the dominant species on higher elevation, is buried by heavy snow, making it unavailable for moose and other lager body sized animal. And according to our surveys, the snow depth of areas with Siberian dwarf pine was usually more than 60 cm, so the dense branches and deep snow cover made it almost impossible for large and medium body sized animals to walking through it, and might be an important reason for the moose to choose river basin with high quality forage and suitable habitat complexity as its major winter home range.
The best resource selection model of roe deer included coniferous and shrub, both positively significant, which meant roe deer avoided coniferous patch edge and sparse shrub habitat (Table 2). Roe deer are medium body sized and have large mass-specific energy requirements to defend themselves from the harsh winter. They usually preferred sites with more abundant annual growth of browse plant species (Hofmann and Stewart 1972), dense forests that provided protective cover (Henry 1981), mature forest habitat of high quality, and avoided edge habitats (Mysterud et al. 1999). In our research area, the coniferous forest is characterized as mature forest, good canopy and no logging for several decades; which could provide enough food resource, sufficient bedding sites and shelters of high quality. The sparse shrub patches in the nature reserve are mostly composed of bare rock and under shrub habitat, which can hardly provide enough shelters and forage and bedding sites.
The resource selection model result showed that mountain hare also preferred areas near the river (Table 2). In Europe, the mountain hare was found predominantly in mixed forest (Naumov 1947), and it reached the highest densities in transition zone with open habitat (Lindlof et al. 1974), and in moorlands (Watson and Hewson 1973). As small herbivores, two main factors, namely food and shelter against predators decide the dense shrub layer to be the best habitat slope and random covariate; K  4, w i  0.183, Table 2). The sable's most parsimonious model consisted of three variables (coniferous, shrub and random covariate; K  3, w i  0.342, Tables 2-3) and we used it as its best resource selection model.
When considering the interspecies interactions, we found three top generalized linear mixed models for mountain hare, lynx and sable. But we didn't find the presence of moose and roe deer was significantly affected by other species in logistic regression model. For the mountain hare the most parsimonious interspecies interaction model included three variables (river, P r , and random covariate; K  3, w i  0.369, Table 4). For the lynx, the most parsimonious interspecies interaction model consisted of four variables (coniferous, slope, P h and random covariate; K  4, w i  0.378, Table 4). For the sable, the most interspecies interaction parsimonious model included three variables (shrub, P h and random covariate; K  3, w i  0.321, Table 4) interspecies interaction model.  Notes: river, the distance to the nearest river; coniferous, the distance to the nearest coniferous patch edge; shrub, the distance to the nearest coniferous patch center; slope, the slope of the point; P r ,the presence probability of roe deer; P h , the presence probability of mountain hare.
moose (Swenson et al. 1999, Odden et al. 2006), but the main predation pressure came from bear, wolf and wolverine (Haglund 1974, Valdmann et al. 2005. Even though many studies revealed that roe deer were the dominating prey species for lynx in the winter (Okarma et al. 1997, Odden et al. 2006, we failed to detected such relationship, maybe because of the limited lynx occurrences (only 35 times) and different habitat requirements between lynx and roe deer in this area. The most parsimonious interspecies interactions model of mountain hare contained two variables, river and P [r] ( Table 4). Roe deer presence probability had a much more strong positive effect on mountain hare presence compared the effect of river. It was demonstrated theoretically that sympatric prey species which share predators usually had positive effects on each other's density (Abrams and Matsuda 1996). Roe deer and mountain hare are two most common distributed herbivores browsing on similar herbs in winter, but differ in body size and digestive systems, such as roe deer are ruminants but mountain hares are hindgut fermentors (Cederlund et al. 1980, Angerbjörn 1981, Pulliainen and Tunkkari 1987. As noted by Illius and Gordon (1992), smaller hindgut fermentors are more tolerant of poor-quality diets when competing directly with larger ruminants which have a more efficient digestive system, especially at times when quality and abundance of food resource reached its nadir in winter. In our study area, mountain hare and roe deer are sympatric species, and food resource utilization partitioning maybe one important reason for mitigating competing pressure in the harsh winter. Mountain hare should browse on the vegetation near the ground when compared to roe deer. At the beginning season, the forage resource is abundance and both species prefer twigs or smaller plants with smaller diameters, the easy digestible parts of plants (Palo et al. 1992). However, as the browsing season continues and vegetable resource becomes more limited at ground level, roe deer will be able to utilize the vegetation beyond the reach of the hare, and another effective way of acquiring grass plants for roe deer is from the bedding sites where they keep warm during the freezing nights (Mysterud et al. 1999) Mountain hares, on the other hand have no such refuge, will turns on the poor-quality forage, namely those twigs and plants with larger diameters (Hulbert and Andersen 2001).
The most parsimonious interspecies interactions model of lynx was best explained by coniferous, slope and P [h] ( Table 4). Lynx showed preference to places with steep slope and avoided coniferous patch edge. In addition, lynx presence showed much more significant correlation with mountain hare presence other than coniferous and slope. The dietary composition and species preyed on by lynx depends on the configuration of herbivore community in the region (Okarma et al. 1997, Pedersen et al. 1999, Weber and Weissbrodt 1999, Koubek et al. 2001, and they tend to be an opportunistic generalist predator that has specialized on small ungulates (Jedrzejewski et al. 1993). In Europe, roe deer is the most widespread ungulates and the most common prey of lynx when they occur together (Jedrzejewski et al. 1993, Nowicki 1997, Belotti et al. 2015. The survey result showed that, in our study area, mountain hare (118 occurrence points) occurred more frequently than roe deer (86 occurrence points), as they were sympatric and mountain hare selected for mountain hare (Hiltunen et al. 2004, Hiltunen andKauhala 2006). Because in winter, the snow concealed lower herbs and shrubs, forcing hares to browse twig, and deciduous bushes (mainly birches Betula spp., rowan Sorbus spp. and willows Salix spp.), which dominated the hares' diet due to their availability above the snow cover (Angerbjörn and Flux 1995). The characteristics of habitat that mountain hare required were well met by areas near the river where there were abundant deciduous shrubs that could provide enough forage and shelter.
The lynx distribution was best explained by the combination of river, coniferous and slope (Table 2). In Norway, a research found lynx would select steep resting sites for safe hiding place in daytime (Sunde et al. 1998), which could be considered as a by-product of selection for less disturbed areas (Basille et al. 2008), meanwhile, researches and field observation had found moderately rugged terrain might be helpful for lynx to get higher prey catch ability (Krofel et al. 2007). The result suggested that lynx preferred the areas near river and with higher slope, but avoided coniferous patch edges. Lynx habitat selection had been studied in many countries and lynx were known to occur in areas with forest cover of greater than 40% (Niedziałkowska et al. 2006, Schadt et al. 2002, Mikusin´ski and Angelstam 2004. Podgórski et al. (2008) studied lynx's microhabitat selection and found that during both summer and winter, lynx selected sites characterized by high complexity, but large-scale approaches couldn't identify variation between different forest types in terms of the availability of specific environmental structures suitable as stalking cover or resting sites. In our study area, there is almost no human disturbance and either coniferous forest or deciduous forest could provide enough shelter, whereas shrub and swamp patch may be too open for lynx's predatory behavior. Lynx often used frozen rivers as long distance migration corridors, and herbivores prefer river areas for foraging (Jiang et al. 2009). Therefore, we predicted that lynx selected areas adjacent to rivers to improve stalking opportunities and used the river as a crossing path between habitat patches.
The resource selection result showed that sable preferred areas near the coniferous patch edge and avoided sparse shrub (Table 2). Sable are small sized forest carnivores selecting sites with a dense tree canopy, mixed dominant tree species, and high densities and diameters of larch tree boles (Buskirk et al. 1996). Certainly, sparse shrub could not provide enough cover in winter, whereas areas near coniferous patch edge usually were characterized with more tree species and densities, which well met the requirements of the survival of sable.

Effect of interspecies interactions
The most parsimonious interspecies interactions model of moose and roe deer that added the presence probability of other animals revealed that the presence of both species had no significant correlation with other animals so that the result was the same as their resource selection models (Table 2, 4). Indeed, lynx was one of the top predators in our study area, we once found an unsuccessful predation site, where a lynx failed to prey on an adult moose, and studies reported that lynx preyed on calves and disabled we failed to find the direct correlation between lynx and roe deer, the occurrence of roe deer was positively correlated with mountain hare, which may enhance the probability of predation of roe deer by the lynx. preferred to use place where roe deer existed. Although we didn't find the direct linkage between lynx and roe deer, but the presence probability of lynx was significantly affected by the presence of mountain hare, which meant lynx preferred to use places with a higher presence probability of mountain hare, however, the presence of mountain hare was correlated with roe deer, indicating more chances for the lynx to meet with roe deer.
The most parsimonious interspecies interaction model of sable was best explained by shrub and P [h] (Table 2). Sable avoided shrub patches and selected place where mountain hare presence probability was higher. According to Miyoshi and Higashi (2005), the sable are small sized carnivores that prefer resting in dense forests to avoid predators (e.g. red fox) and strong wind and stalking in forests with abundant prey such as voles and mice. Sable mainly hunt for small animals, but it also preys on larger mammals such as chipmunks Tamias sibiricus, pikas Ochotona, squirrels Sciurus and Pteromys, muskrats Ondatra, marmots Marmota, mountain hares, and musk deer Moschus moschiferus (Moskov 1973, Khlebnikov 1977, Monakhov and Bakeyev 1981, Zirjanov et al. 2001, Chernikin 2006. Our results suggested sable preferred areas with higher mountain hare occurrence probability, and its impact was much stronger than shrub. As a small sized carnivore, sable are preyed on by many larger predator, such as brown bear Ursus arctos, wolf Canis lupus, red fox Vulpes vulpes, wolverine Gulo gulo, yellow-throated marten Martes flavigula, lynx Lynx lynx and some big birds, like eagle-owl Bubo bubo, golden eagle Aquila chrysaetus, raven Corvus corax, hawks: goshawk Accipiter gentilis and sparrow-hawk Accipiter nisus, and two owls: great gray owl Strix nebulosa and northern hawk owl Surnia ulula (Gusev 1966, Bakeyev et al. 2003. In our study area, lynx, wolverine and some other big birds might hunt sable, shelter was very important for sable and sparse shrub was too open for sable to escape from predation.

Conclusion
Our analysis result showed that food availability and shelter were most important factor deciding habitat use by the herbivores, and also the distribution of prey abundance and cover had a great impact on habitat selection of the carnivores. Different mammals chose different habitats to meet their needs for forage, shelter or escape. Areas near river was mostly favored by moose, mountain hare and lynx, while moose and mountain hare preferred river areas for browsing shoots, lynx's strong selection on river maybe the result of using river as a quick path to other habitat patches and there were more predation opportunities. Sparse shrub habitat was too open for mammals which could not provide enough shelter or cover to prevent them from predation or danger, so roe deer and sable showed strong avoidance for sparse shrub habitat. The interspecies interaction was studied in this research, and we found that the distribution of mountain hare was correlated with roe deer maybe because of the spatial separation of food resource during the harsh winter. Both lynx and sable preyed on mountain hare, and the strong relationship was detected by lynx and sable model when incorporating the occurrence probability into the models, which showed the importance of mountain hare in this area. Even though