Sika deer Cervus nippon in Hokkaido, Japan, have recovered from a population bottleneck about 120 years ago and their distribution has expanded rapidly in the last three decades. We tracked 53 radio-collared female sika deer, and obtained 4,430 locations during the 25-month study period from April 1997 to April 1999. We examined the seasonal distribution of female sika deer in relation to spatial landscape features (snow depth, vegetation, bamboo grass and roads) with a logistic regression model using a geographic information system database. We presented a population-landscape scale evaluation of sika deer habitat for summer and winter within the telemetry study area (TSA) using resource selection functions. We then extrapolated the model to the rest of eastern Hokkaido to discuss the seasonal migration for an expanding population. Most radio-collared sika deer (71%) moved between high-elevation summer and low-elevation winter ranges, whereas some (29%) moved between low-elevation summer and similar or high-elevation winter ranges. During winter, sika deer selected middle elevation habitats (200–400 m a.s.l.) with both a relatively low snow depth and the presence of coniferous and mixed forests. On the other hand, sika deer were widely distributed regardless of elevation during summer, although they were further from roads and less often in agricultural lands. Within the TSA, the suitable habitat was very limited during winter compared with during summer. Although migration from summer to winter ranges may depend on the abundance and distribution of suitable winter habitat at a landscape scale, migration from winter to summer ranges could not be explained from this study. Our approach is useful for understanding the relationships among seasonal habitat selection, seasonal migration and the expansion of the population.
Effective management of wildlife populations largely depends upon understanding and predicting their habitat needs (Clark, Dunn & Smith 1993). Conservation biologists and resource managers face complex challenges in managing large animals, especially migratory mammals, because they often move or migrate long distances crossing both good and poor habitat. Management across disjunctive areas of suitable habitat is necessary (Mlandenoff, Sickley & Wydeven 1999). Therefore, management can benefit from large-scale evaluation and mapping of seasonally favourable habitat distribution of these species on a landscape scale (Mlandenoff et al. 1999). From a viewpoint of wildlife management or conservation, some studies have presented the prediction of potential or favourable habitat of a target species (Mlandenoff, Sickley, Haight & Wydeven 1995, Kaji, Miyaki, Saitoh, Ono & Kaneko 2000), or the relationship between the distribution of a species and human activities (Pereira & Itami 1991, Mace & Waller 1996, Mace, Waller, Manley, Lyon & Zuuring 1996, Mace, Waller, Manley, Ake & Wittinger 1999).
Many populations of large herbivores such as wildebeest Connochaetes taurinus, zebra Equus burchelli and Thompson' s gazelle Gazella thompsoni in Africa (Maddock 1979), and mule deer Odocoileus hemionus hemionus (Garrott, White, Bartmann, Carpenter & Alldredge 1987), moose Alces alces (Edwards & Ritchey 1956), black-tailed deer O. hemionus (Loft, Menke & Burton 1984), and white-tailed deer O. virginianus (Larson, Rongstad & Terbilcox 1978) in North America migrate between seasonal ranges. Migration of these animals has probably evolved to take advantage of spatial and temporal variations in the environment (French, Reed, Calambokidis & Cubbage 1989). Fryxell & Sinclair (1988) demonstrated that seasonal migration has important implications for the structure and dynamics of large herbivore communities.
Northern ungulates migrate between their summer and winter ranges, and a great deal of data on seasonal movement and factors affecting seasonal habitat selection of ungulate species in temperate ecosystems have been reported (Loft et al. 1984, Schoen & Kirchhoff 1985, Kufeld, Bowden & Schrupp 1989, Van Deelen, Champa, Hamady & Haufler 1998). Snow depth (Klein & Olson 1960, Garrott et al. 1987), vegetation cover (Bloom 1978, Wallmo & Schoen 1980, Kufeld, Bowden & Schrupp 1988), forage (Hobbs & Spowart 1984, Garrott et al. 1987, Kufeld et al. 1989), and road systems (Rost & Bailey 1979) have been demonstrated as factors affecting winter habitat selection. On the other hand, habitat quality (Loft et al. 1984, Bowyer 1986, Kufeld et al. 1989, Nicholson, Bowyer & Kie 1997), forage (Garrott et al. 1987), vegetation cover (Loft et al. 1984), human disturbance (Nicholson et al. 1997) and predation (Riley & Dood 1984, Festa-Bianchet 1988) have been implicated as factors affecting summer habitat selection. Previous studies suggested that snow depth and bamboo grass Sasa spp. variety (Kaji et al. 2000) and vegetation cover (Hokkaido Institute of Environmental Sciences 1994) were important variables limiting sika deer Cervus nippon distribution in Hokkaido, Japan. The importance of bamboo grass as winter forage has been recognised for sika deer in Japan (Takatsuki 1986, 1993, Yokoyama, Kaji & Suzuki 2000). Most of these studies do not consider the population status or historical background of the population dynamics (but see Mlandenoff et al. 1999, Kaji et al. 2000). Therefore, studies of the relationship between seasonal habitat selection and seasonal migration in relation to expansion of population size and distribution are needed.
Sika deer historically occurred throughout Hokkaido, the northernmost island in Japan, and they migrated seasonally between the west, where snows are heavy, and east, where wintering grounds are more open (Inukai 1952). Sika deer on Hokkaido have a reduced level of genetic variation because of a population bottleneck (Nagata, Masuda, Kaji, Kaneko & Yoshida 1998) about 120 years ago, caused by overharvest and heavy snowfall (Inukai 1952). Through the genetic differences among the populations, Nagata et al. (1998) described three subpopulations (Daisetsu, Akan and Hidaka) consisting of individuals from the main three coniferous forests in Hokkaido (Fig. 1).
Under government protection, the number of sika deer gradually recovered during the last three decades (Kaji et al. 2000). Elimination of wolves Canis lupus, replacement of native mixed hardwood forests with conifer plantations, and increased pasture acreage likely contributed to the expansion in deer distribution (Kaji et al. 2000). By the mid-1970s, sika deer occupied most available habitat in the eastern half of the island (Kaji et al. 2000). Damage to agriculture and forestry by sika deer remained at low levels from the mid-1950s through the mid-1970s, but increased dramatically to nearly USD 15 million in 1990 and to over USD 30 million in 1996 (Kaji et al. 2000). As the range of sika deer expanded, agricultural and forest damage increased, resulting in relaxed hunting regulations. Total harvest fluctuated between 2,000 and 3,000 animals during 1970–1989 and increased to 16,134 animals in 1990 and 46,634 animals in 1996 as more areas were opened for hunting and antlerless harvest was allowed (Kaji et al. 2000). The population of sika deer reached a peak about 1996 and tended to decrease gradually due to heavier hunting (Hokkaido Government 1998).
Figure 1.
Location of the telemetry study area (TSA) in eastern Hokkaido, Japan, where seasonal resource selection function models were developed. All radio-collared deer were captured on the Shiranuka Hills (). On the inserted map of Hokkaido three genetically differential populations of sika deer are indicated (Nagata et al. 1998).

The deer population recovery and expansion in eastern Hokkaido may have originated from the Akan subpopulation (Nagata et al. 1998). Eastern Hokkaido is a largely continuous suitable habitat for sika deer (Hokkaido Institute of Environmental Sciences 1994, Kaneko, Kaji & Ono 1998, Kaji et al. 2000). Based on the background of the sika deer population in eastern Hokkaido, we predict that seasonal habitat selection in the expansion of the population is closely related to present patterns of seasonal migration in eastern Hokkaido. However, the relationship between seasonal habitat selection and seasonal migration remains poorly understood.
Bowyer (1986), Scarbrough & Krausman (1988) and Bowyer, Kie & Ballenberghe (1996) have reported sexual differences in the use of forages, habitats and space by deer. As female deer may be able to increase their reproductive success by producing more and healthier young, habitat selection will probably be more important for females than for males. Therefore, a possible relationship between migration and habitat selection may be clearer for females than for males, which is why we focused our study on females. In this paper, we examine the seasonal distribution of radio-collared female sika deer in relation to spatial landscape features using a logistic regression model based on a geographic information system (GIS) database, and examine the factors affecting seasonal habitat selection of an expanding population. We then present a population-landscape scale evaluation of sika deer habitat for summer and winter using resource selection functions (RSF; Manly, McDonald & Thomasl993), and demonstrate how seasonal migration is related to distribution of suitable habitat on a landscape scale. Further, we extrapolate the RSF model to the remainder of eastern Hokkaido to discuss the seasonal migration in relation to the expansion of the population.
Material and methods
We conducted the telemetry study in 7,466 km2 of eastern Hokkaido (see Fig. 1). The Shiranuka Hills is one of the largest and most important wintering areas for sika deer, which have been expanding their distributional range in eastern Hokkaido during the last three decades (Kaneko et al. 1998). Within the Shiranuka Hills wintering population, both resident and migratory individuals occur (H. Igota, unpubl. data). Resident deer remain in the Shiranuka Hills throughout the year while migratory deer migrate seasonally and have their summer ranges to the north and east of the Shiranuka Hills. The telemetry study area (TSA) boundary was defined by the minimum convex polygon encompassing all seasonal movements of the radio-collared deer. The Shoro River and several roads ran through the Shiranuka Hills from north to south. Many hunters have access to this area using the roads to hunt deer during winter (Hokkaido Institute of Environmental Sciences 1995). Plantations of Sakhalin fir Abies sachalinensis and Sakhalin spruce Picea glehnii occur in parts of the Shiranuka Hills. The climate of the study area is cool and humid with an average temperature of-10.4°C in February and 17.4°C in August, the mean annual temperature is 3.6°C and annual precipitation is 1,200 mm at Akan (Igarashi 1986).
The high elevation regions (300–800 m a.s.l.) of the TSA are primarily composed of coniferous and mixed forests. Conifers, Yezo-spruce P. jezoensis, Sakhalin fir and Sakhalin spruce, are dominant and mixed with Japanese linden Tilia japonica and Japanese oak Quercus crispula in mountainous areas over 400 m. Alpine species such as Erman' s birch Betula ermanii and Japanese stone pine Pinus pumila dominate above 800 m. Most natural forests are classified as ‘pan mixed forest’ with conifers and broad-leaved trees typical of the transition between temperate and subarctic zones (Tatewaki 1958). The low elevation regions (below 300 m) are dominated by urban development and agriculture, but deciduous forests and plantations of Japanese larch Larix leptolepis persist in this area. Deciduous broad-leaved forest is comprised of Japanese oak, Japanese white birch B. platyphylla, painted maple Acer mono, Japanese linden, Japanese elm Ulmus davidiana, Japanese ash Fraxinus mandshurica, Amur cork tree Phellodendron amurense, caster aralia Kalopanax pictus, and willow Salix spp. (Igarashi 1986). Dominant undergrowth is bamboo grass Sasa spp., which covers 90% of the Hokkaido forests (Toyooka, Sato & Isizuka 1983). There are four species of bamboo grasses. The largest, S. kurilensis (Sk), is found primarily in the west where snowfall is heavy. S. senanensis (Ss) is a medium-sized bamboo grass which is found in broad-leaved and mixed forests with intermediate snow accumulation. S. nipponica (Sn), the smallest bamboo grass, is found in the southeastern region with little snowfall. S. borealis (Sb) has a very limited distribution within the Sn range (Matsui 1963). Sn and Ss dominate over a large part of the TSA.
Table 1.
Climatic data for the northern and southern portions of the study area in eastern Hokkaido. Japan (Japan Weather Association Hokkaido Regional Head Office 1991).

The extrapolated study area (ESA), the extent to which the Akan subpopulation potentially could expand its range (Nagata et al. 1998), encompassed about 15,847 km2 of sika deer habitat in eastern Hokkaido, bounded to the west by the boundary of the Shiranuka Hills and the Daisetsu mountains (see Fig. 1). The vegetation structure of the ESA is similar to that of the TSA. High elevation regions are covered with forests, low elevation regions are dominated by urban development and agriculture, and the undergrowth is composed primarily of Sn and Ss. The Shiretoko Hills and Akan mountain ranges are continuous (see Fig. 1), and divide the study area into northern and southern portions based on climatic conditions (Table 1). The Sea of Okhotsk affects the climatic characters of the northern, and the Pacific Ocean the southern portions of the mountains, respectively (Japan Weather Association Hokkaido Regional Head Office 1991).
A total of 53 female sika deer were captured within the Shiranuka Hills with anesthetic guns in April 1997 (19 individuals), March 1998 (31 individuals), and March 1999 (three individuals). Captured deer (37 adults, six yearlings, 10 fawns) were fitted with VHF radio collars (144 MHz, ATS, Isanti, Minesota, USA or Lotek Engineering, Inc., Newmarket, Ontario, Canada) with an expected battery life of three years, and ear tags for visual identification. From April 1997 to April 1999, deer were located by triangulation procedures (White & Garrott 1990) during daylight hours from the ground with radio-equipped trucks and on foot every two or three days. Triangulation bearings were taken with a 107-cm, 3-element Yagi antennae and receivers (FT-290mk II, YAESU MUSEN, Tokyo, Japan). Locations were recorded on 1:25,000 scale topographic maps. We also recorded the location of each marked individual that was visually observed.
To account for inherent errors in triangulated deer locations, we estimated the circular error associated with telemetry by calculating the 99% confidence limits (Murakami & Mano 1998) for the average distance of our estimates from actual locations of radio-collars placed at 60 locations. The average error in locations was 145 m ± 149 (SD) and the upper 99% confidence limits was 196 m. We considered all habitats that occurred within a 196-m circle as an estimate of the habitat at that location. Moreover, radio-telemetry location from the ground has limitations regarding the area that we could search by roads. All telemetry locations (N = 4,837) from 1997 to 1999 were distributed within 5 km from road networks (99.7% of TSA), and 80% of those were within 1 km of road networks (68.9% of TSA).
Deer locations were grouped into seasons based on the movement behaviour of sika deer. We defined the seasonal ranges of migratory deer during the period when movements stabilised over time (i.e. localised movements within areas over extended periods that were separated from another area of localised movement in time). Winter was defined as November to April for resident deer, and the period when migratory deer stabilised in their winter ranges. Summer was defined as May to October for resident deer, and the period when migratory deer stabilised in their summer ranges. Spring and fall were defined as the periods when locations fell sequentially between the summer and winter ranges for migratory deer, and data for these periods were not needed in determining habitats for summer and winter ranges in this study.
To define habitat characteristics, we developed a 1-km block database map containing habitat variables with a GIS. We used the GIS software Arc View 3.1 (Environmental Systems Research Institute, Inc., Redlands, California, USA). The TSA consisted of 7,466 blocks, and the ESA consisted of 16,582 blocks excluding 120 blocks that fell within lakes. We used the snow depth and elevation data of national land numerical information (National Land Agency 1992). Although the snow depth was the mean during winter for the 30-year period (1955–1984) and was not real time data during our study period, it represented the geographical distribution of relative snow depth throughout Hokkaido. We used the digital road map (Pasco digital map, Tokyo, Japan), and distance to road of a block was defined as the distance from the center of the block to the nearest road. We generated a 30-m grid vegetation database from Natural Environment Information GIS (Environment Agency 1997). Vegetation types were classified into six types (coniferous forest, mixed forest, deciduous forest, agricultural land, Japanese larch forest and alpine vegetation) and were defined as the area percentage of each vegetation type within the block. We used a 1:200,000 distribution map of bamboo grass species in Hokkaido (Toyooka, Sato & Isizuka 1984). Bamboo grass forage value was ranked in the order Sn > Ss > Sb > Sk by Kaji et al. (2000) based on availability (Hokkaido Government 1987, Takatsuki 1992a), grazing tolerance (Matsui 1963, Agata, Kubota & Kamata 1979) and nutritional value (Ohara 1948, Kawai, Juni, Yasue, Ogawa, Hata, Kondo, Okubo & Asahida 1995, Yamane & Furubayashi 1997). We assigned the blocks the relative forage value of each species according to Kaji et al. (2000). Bamboo grass forage value (scores in parentheses) were none (0), Sk (1), Sk+Ss (2), Sb or Sn + Sb or Ss + Sb (3), Ss (4), Ss + Sn (5), Sn (6) (see Kaji et al. 2000).
Table 2.
Categorisation of habitat variables for logistic regression analysis.

Sika deer telemetry coordinates from 1997 to 1999 were pooled to assess habitat selection at a population-landscape scale (Johnson 1980). Telemetry sample sizes were 2,384 locations for summer and 2,046 locations for winter. The blocks containing the 196-m circle of the telemetry locations were defined as used blocks, and we compared habitat variables of used blocks between summer and winter using a Mann-Whitney' s U test to examine for seasonal differences in habitat use. In order to examine habitat selection, we compared habitat variables of used blocks with available blocks by maximum-likelihood logistic regression (Manly et al. 1993). A total of 1,999 random blocks were selected within the TSA to evaluate availability using the Arc View extension MOVEMENT (Hooge & Eichenlaub 1998). The proportions of random coordinates selected within 1 km and within 1–5 km of roads were chosen to be comparable to the proportions in the telemetry dataset (80 and 20%, respectively). All values were categorised for logistic regression analysis (Table 2). Snow depth was used only for winter analysis. Elevation was excluded from logistic regression analysis because it should correlate with snow depth or vegetation, and habitat selection of sika deer should be affected by snow depth, vegetation, or roads, not directly by elevation. We examined multi-collinearity among independent variables by a Variance Inflation factor (VIF). There was no significant multi-collinearity (all VIF < 3) according to Hayakawa (1986; if VIF > 10, there is a significant multi-collinearity). For each season, the model was developed using signif-icant (P < 0.05) variables by a backward stepwise elimination with the likelihood-ratio criterion. The model was ascertained by comparing the model containing all variables to the ‘no selection model’ (Hosmer & Le-meshow 1989). The logistic regression coefficients represented the contribution of each variable to explain resource selection relative to random availability. The association and power of each variable was ascertained by each coefficient' s sign (negative or positive) and strength (-2 log likelihood-ratio: -2LLR).
For each season, we calculated RSF using the equation from Manly et al. (1993):
where w(x) is the RSF, Xi is habitat variables from 1 to p, ßi is coefficient estimates for each habitat variables Xi, and ß1x1 + ... +ßpxp represents the available resource units. The RSF values represented the relative probability of use of each unique variable combination by sika deer. The RSF values were scaled as probability by dividing each RSF value by the largest RSF value for each season. Calculations were made with the GIS to produce a map containing RSF values within each 1-km block within the TSA, and then throughout the ESA. The extrapolation of TSA coefficients to the ESA assumed that the same factors controlling RSF levels in the TSA applied as well to the ESA. We considered the area of less than 1% RSF as the boundary of suitable and unsuitable habitat for sika deer in the TSA, because 1% RSF has an extremely low probability of use. We used the computer software SPSS (SPSS Inc. Chicago, Illinois, USA) for data analyses.Table 3.
Data on summer and winter habitat variables for radio-collared female sika deer in eastern Hokkaido, Japan, during 1997–1999.

Table 4.
Logistic regression coefficients for the final model for summer and winter with only significant variables (P < 0.05) entered.

Results
The summer and winter telemetry locations were contained in 240 and 54 blocks, respectively (Fig. 2). A comparison of habitat variables indicated no significant difference in coniferous forest, agricultural land, Japanese larch forest, alpine vegetation and bamboo grass variety used by sika deer between summer and winter (Table 3). Snow depth of blocks used during summer were significantly greater than those used during winter. Elevation of blocks used during summer was more variable than the elevation of blocks used during winter, and the elevation range during summer (10–977 m) included the elevation range during winter (129–486 m). Most radiocollared sika deer (71%) moved between high elevation summer and low elevation winter ranges, whereas some (29%) moved between low elevation summer and similar or high elevation winter ranges. Coniferous and mixed forests together comprised the majority (about 60%) of all vegetation types used during both summer and winter. Use of mixed and deciduous forests was significantly greater during winter than during summer. Agricultural land comprised a small proportion of used blocks, but was somewhat larger during summer than during winter. Alpine vegetation was a minor component of blocks used by deer in either season. Distance to roads was significantly greater for blocks used during summer than during winter.
The logistic regression model was significantly different from the ‘no selection model’ for summer (-2LL = 1335.252, x2 = 189.961, df = 6, P < 0.001, r2 = 0.165), and winter (-2LL = 367.283, x2 = 132.196, df = 5, P < 0.001, r2 = 0.289). Variables, eliminated by backward stepwise regression from the final model, were Japanese larch forest and alpine vegetation for winter, and deciduous forest and bamboo grass variety for both seasons (Table 4). During summer, females were more often in coniferous, mixed or Japanese larch forests and less often in agricultural land or alpine vegetation, and further from roads. During winter, females were more often in coniferous or mixed forests and less often in agricultural land or areas with deep snow, and closer to roads.
Figure 3.
Resource selection probability functions (RSF) during summer (A) and winter (B) for the telemetry study area (TSA) and the extrapolated study area (ESA). Each 1-km block represents the relative probability of use by female sika deer. White areas represent RSF values of less than 1%.

Within the TSA, the area of less than 1 % RSF was extremely small during summer (0.4% of TSA), but relatively large during winter (80.0% of TSA). During winter, the remaining 20% area of greater than 1% RSF appeared around the Shiranuka Hills (Fig. 3). For the ESA, the area of less than 1% RSF was 1.2% and 81.9% for summer and winter, respectively. Although almost all of the ESA was over 1% RSF during summer, mountainous area with coniferous or mixed forest and without roads tended to have a higher RSF value. During winter, the area of more than 1% RSF appeared around the Akkeshi and Nemuro area in addition to the area around the Shiranuka Hills (see Fig. 3). Areas of more than 1% RSF rarely appeared in the north portion of the ESA.
Discussion
Our study shows that female sika deer select the areas with less snow during winter, but are not attracted to bamboo grass variety during summer or winter. However, Kaji et al. (2000) demonstrated that sika deer distribution on Hokkaido was limited by snow depth and bamboo grass availability with a great forage value. This apparent conflict is probably due to the difference of scale on habitat selection (Johnson 1980, Powell 1994). Bamboo grass variety may be an important factor limiting sika deer distribution on Hokkaido, but not in eastern Hokkaido, because a large part of TSA of eastern Hokkaido is dominated by the species with a high forage value (Sn and Ss). We demonstrate that snow depth is the predominant factor affecting winter habitat selection of sika deer within the area dominated by bamboo grass variety with a high forage value. This pattern is consistent with the current knowledge of impacts of snow on ungulate behaviour. Deep snow reduces access to forage for ungulates (Gilbert, Wallmo & Gill 1970) and increases the cost of locomotion (Parker, Robbins & Hanley 1984, Fancy & White 1987). Both these factors could contribute to increased mortality (Robinette, Julander, Gashwiler & Smith 1952, Klein & Olson 1960, Uno, Yokoyama & Takahashi 1998), and/or could affect subsequent maternal care of young (Langenau & Lerg 1976, White & Luick 1984).
Garrott et al. (1987) suggested that migration is obligatory for mule deer occupying high elevation moun tainous areas during summer, as snow accumulation precludes use of these areas during winter. Similarly, many studies have reported the movement between high elevation summer and low elevation winter ranges with less snow (Schoen & Kirchhoff 1985, Albon & Langvatn 1992). Previous studies on sika deer in the central and northern parts of the island of Honshu in Japan also reported this type of movement (Maruyama 1981, Ito & Takatsuki 1987). Nevertheless, some radio-collared sika deer inhabiting eastern flatlands moved between low elevation summer and similar or high elevation winter ranges in our study. This pattern of movement has not been reported before in sika deer or other ungulates. It demonstrates the importance of habitat variables in selecting winter ranges, within the constraints imposed by climatic variables such as snow depth.
During winter, the habitat selection of sika deer was positively associated with coniferous and mixed forests. These are likely important winter habitats for sika deer to access food resources in Hokkaido (Hokkaido Institute of Environmental Sciences 1994, Kaji 1995). There are many studies indicating the importance of mature coniferous forest to deer during winter. Old-growth forest has been shown to provide protection from deep snow, access to nutritional forage, and essential winter habitat for Sitka black-tailed deer O. hemionous sitkensis on Admiralty Island, southeast Alaska (Bloom 1978, Wallmo & Schoen 1980). Vegetation cover may provide rest sites, protection and/or refuge from severe climatic conditions, predators and human activities (Hamilton, Drysdale & Euler 1980, Takatsuki 1989, 1992b). Only coniferous and mixed forests provide cover during winter in this area. We suggest that higher elevation habitats (over 400 m) are not suitable for deer during winter because of deep snow, and eastern lower elevation habitats (below 200 m) are not suitable for deer because of the lack of coniferous and mixed forests. Therefore, during winter, sika deer selected middle elevation habitats (200–400 m) with both a relatively low snow depth and the presence of coniferous and mixed forests. It is therefore likely that our data confirm the importance of vegetation cover for northern ungulates. During summer, sika deer were widely distributed regardless of elevation within the TSA, although highest elevation habitats (over 1,000 m) dominated by alpine vegetation were not selected and had less than 1% RSF value. This is probably because deciduous and Japanese larch forests in eastern low elevation regions, in addition to coniferous and mixed forests, also play a role in providing cover for deer during summer. Our data also likely support the importance of vegetation cover for deer during summer as well as during winter. Deer were less often in agricultural lands during summer, and this is probably due to avoidance of hunters (Sakuragi, Igota, Uno, Kaji, Kaneko, Akamatsu & Maekawa 2002).
Considering the limitation of our radio-telemetry location from roads, we chose the random coordinates to evaluate availability within 1 km and within 1 -5 km of roads with the proportions in the telemetry data set (80 and 20%, respectively). Thus the limitation of our radio-telemetry location from roads could be cleared by this treatment because habitat selection analysis showed that the deer were further from roads during summer and closer to roads during winter. During summer, deer may be able to select habitats further from roads to conceal from human disturbance because a wide area is available for deer. Female deer may especially avoid human disturbance during calving (Unsworth, Kuck, Carton & Butterfield 1998). During winter, however, sika deer may be forced to select habitats closer to roads because the selection is strongly affected by snow depth and the presence of coniferous and mixed forests. Studies of the responses of white-tailed deer to vehicles (Behrend & Lubeck 1968) and to snowmobiles (Dorrance, Savage & Huff 1975) suggested that hunted populations were less apt to habituate to these disturbances than were unhunted populations. If sika deer were free to select habitats, they might be distributed further away from roads even during winter, at least during the hunting season (November to January), because during this period many hunters use roads in this area to hunt sika deer (Hokkaido Institute of Environmental Sciences 1995). Rost & Bailey (1979) indicated that mule deer and elk C. canadensis avoided roads, and it might be primarily a characteristic of hunted populations. On the other hand, deer may use roads during winter after the hunting season because roads are easier to travel for animals (Mladenoff et al. 1999).
Within the TSA, the abundance and distribution of suitable habitat for deer changes considerably between summer and winter, with winter habitat of sika deer being very limited. The origin of migration patterns may be understood on the basis of the difference in abundance and distribution of suitable habitat between summer and winter. Deer may move from summer to winter ranges if they cannot overwinter on their summer ranges. This leads to the conclusion that seasonal migration could have developed as a way to avoid undesirable conditions at a particular time of year (Vaughan, Ryan & Czaplewski 2000). This provides an explanation why deer move from summer to winter ranges, but the question why deer move from winter to summer ranges remains undetermined. Resident deer can and do oversummer on winter ranges without migration. Current evidence suggests that migration is generally selected for as a means of enhancing access to high-quality food and/or reducing the risk of predation (Fryxell & Sinclair 1988). Selection should favour individuals that migrate, if by migrating, their reproductive success is enhanced (Baker 1978). Because of environmental fluctuation and individual differences in costs of migration, several strategies related to migration can occur in the same species or in the same population (Fretwell 1972). To answer the question, we should consider other habitat variables and/or some other factors in relation to migration events, such as survivorship and birth rate (Caughley 1977), quality and availability of forage (Taylor & Taylor 1977, Fryxell & Sinclair 1988) and selection of favourable environment for breeding and parenting (Bowyer, Kie & Ballenberghe 1998). Coefficient of determination of logistic regression models was higher for winter (0.289) than for summer (0.165), although the low coefficients of determination of logistic regression models for summer and winter suggested that the majority of the variance in the deer distribution might have been left unexplained. This implies that suitable winter habitats could be relatively well predicted from spatial landscape features at this scale (i.e. 1-km block), but suitable summer habitats are less predictable based on only these spatial landscape features.
The extrapolation of coefficients for the TSA to the ESA reveals that additional other winter habitats are suitable for deer in the southern portion of the ESA. Deer use these areas during winter (Kaji et al. 2000), indicating that the prediction of this model agrees with the actual sika deer distribution. However, our prediction of suitable winter habitats may be underestimated, especially in the northern portion of the ESA. A possible reason may be the climatic difference between the northern and southern portions: snowfall is heavier in the northern portion than in the southern portion.
The extrapolation study gives us some suggestions for winter and summer habitat selection in relation to the expansion of the population. Although some of the migratory radio-collared deer have their summer ranges closer to Akkeshi or Nemuro than to the Shiranuka Hills, all of them move to the Shiranuka Hills regardless of the presence of closer winter habitats. This phenomenon may reflect the expansion of the sika deer population in eastern Hokkaido, because the eastern Hokkaido deer population originated from the Akan population (Nagata et al. 1998). Other winter habitat areas such as Akkeshi and Nemuro may be newly selected winter habitats in recently reoccupied ranges.
Although a large part of the ESA is suitable for deer during summer, the extent of the summer range of the radio-collared deer is limited and no deer had their summer range south and west of the Shiranuka Hills. A possible reason for this may be the presence of other sika deer populations or male sika deer in these areas. Summer habitat selection might be affected by various factors such as, for example, the connectivity of corridor habitat from the winter to the summer ranges and costs of migration. However, what drives sika deer movements from winter to summer habitats still needs to be investigated. In addition, fidelity to seasonal ranges and migration route or migration tradition may be an important factor affecting the present pattern of seasonal migration. Matrilineal family groups of female and fawn white-tailed deer migrate together (Nelson & Mech 1981) and show fidelity and often exclusive use of seasonal ranges (Tierson, Mattfeld, Sage & Behrend 1985, Aycrigg & Porter 1997). Migration tradition is maintained through social interactions (Van Deelen et al. 1998), and fawns and yearlings learn the migration traditions of their social groups and remain faithful to them throughout life (Nelson & Mech 1981, 1984, Nixon, Hansen, Brewer & Chelsivig 1991). These assumptions need to be examined in future studies of sika deer and other species.
Conclusions
Our study examined the summer and winter habitat selection of an expanding sika deer population at a landscape scale, and presents a landscape-scale evaluation of sika deer habitat for summer and winter in eastern Hokkaido, Japan. Movements of some individuals between low elevation summer and similar or high elevation winter ranges in our study were exceptions to the movement reported from previous studies (Maruyama 1981, Schoen & Kirchhoff 1985, Garrott et al. 1987, Ito & Takatsuki 1987, Albon & Langvatn 1992). Snow depth and coniferous and mixed forests affected the winter habitat selection of sika deer, and migration from summer to winter ranges was dependent on the abundance and distribution of suitable winter habitat on a landscape scale. However, we could not explain migration from winter to summer ranges from the population-landscape scale study above.
Our study supported the assumption that the eastern Hokkaido deer population might have originated from the Akan subpopulation. Although some authors have already demonstrated that a landscape-scale habitat evaluation of a target species is useful for management and conservation (Pereira & Itami 1991, Mace et al. 1996, Mlandenoff et al. 1999, Kaji et al. 2000), our study is the first in which this approach has been applied to an expanding population, and demonstrates the relationships among seasonal habitat selection, seasonal migration and the expansion of the population.
Acknowledgements -
we thank D.R. McCullough for providing suggestions and advice and for critically reviewing this manuscript. We are very grateful to S. Hamasaki, M. Kishimoto, Y. Kamiyama, the staff of the Wildlife Management Office, EnVision, and the Hokkaido Government for capturing deer and supporting our work. We also thank T. Suzuki, Buheaosier, S. Shimada, S. Kameyama and K. Waseda for helping us with the data analysis, and Y. Furukawa, H. Takahashi, J. Tanaka and H. Hamamoto for assisting with collection of location data. We thank the Akan Town Office, Shiranuka Town Office, Teshikaga Town Office, and H. Sugawara for supporting our work. We acknowledge K. Tamada, T. Saitoh and the students of the Laboratory of Boreal Forest Conservation, Hokkaido University for stimulating discussions. This study is supported by the Hokkaido Government and by a Grant-in-Aid for Scientific Research (14000353 to MS and 14340240 to KM) from the Ministry of Education, Culture, Sports, Science and Technology, and the Agency of Japan Society for the Promotion of Science.