The Effects of Landscape Components, Wildlife Behavior and Hunting Methods on Hunter Effort and Hunting Efficiency of Sika Deer

Population management of overabundant deer is essential in many countries. Hunting is one of the ways to manage such deer so it is important to understand how hunting affects deer and numbers caught. The number caught is the product of hunter effort and hunting efficiency. However, factors affecting effort and efficiency have rarely been properly investigated. I used a Bayesian state—space model to examine the effects of landscape components and deer behavior on the effort and efficiency of two hunting methods in hunting sika deer Cervus nippon. With shooting by gun, effort increased as the deer abundance of previous year increased and the percentages of wildlife protection areas and city areas decreased. With trapping, effort increased as the percentage of wildlife protection areas decreased. Efficiency of shooting decreased as the effort of the previous year increased. With trapping, efficiency was not affected by the effort of the previous year. Efficiency of trapping plateaued with the increase of deer abundance and the decrease of efficiency in relation to deer abundance was stronger with trapping than with shooting. In conclusion, effort and efficiency were affected by landscape components, deer behavior, and hunting methods. I recommend intensive shooting in the initial phase, and after deer abundance decreases and deer vigilance increases, trapping should be adopted for a sustainable hunting efficiency.

In recent years, the management of deer populations has become increasingly necessary (Côté et al. 2004, Takatsuki 2009). Population control by game hunting is one of the ways to manage overabundant deer (Decker andConnelly 1989, Simard et al. 2013). Successful management requires that the effects of hunting be understood in order to avoid unwanted population explosions or depletion of the targeted wildlife. A direct effect of hunting is the number of animals caught, and this number is determined by the product of hunter effort (hereafter, effort) and hunting efficiency (hereafter, efficiency). Therefore, to conduct effective population control of deer by hunting, factors affecting effort and efficiency should be clarified.
Many studies have examined the factors affecting effort and efficiency (Davidson and Fraser 1991, Harden et al. 2005, Ward and Myers 2005. Davidson and Fraser (1991) showed that deer hunter effort was high in places of high accessibility and good visibility. Harden et al. (2005) showed that deer hunters avoided urban areas. However, each of these studies has two major problems. First is the treatment of two type errors. In addition to demographic process error, effort and efficiency fluctuate greatly with each hunting event because of uncontrollable conditions like weather, the equipment for hunting, and the skill of hunters. A state-space model that can explicitly treat both stochastic and measurement errors (Calder et al. 2003) is probably the most useful way to analyze effort and efficiency, but statespace models have rarely been applied to effort and efficiency. The second problem is the difficulty in distinguishing between change in wildlife abundance and change in wildlife behavior caused by previous effort, although this problem is specific to efficiency. It is known that changes in efficiency correspond to changes in wildlife abundance (Bigelow et al. 2002, Haggarty andKing 2006). However, at the same time, wildlife escape from hunters (Kilgo et al. 1998, Martin andBaltzinger 2002), which also decreases efficiency. Therefore, it is necessary to examine the factors affecting efficiency by discriminating between changes in abundance and changes in wildlife behavior (Harley et al. 2001). If the behavior of the target wildlife changes in response to the previous year's hunting effort, then hunting efficiency should not be used as an index of wildlife abundance.
Deer are a major target of game hunting, and data on effort and efficiency with respect to deer hunting have been collected from locations all over the world (e.g. Solberg et al. 1999 in Norway;Ueno et al. 2010 in Japan). Models to estimate deer abundance have been developed because of the demand for management of overabundant deer (Iijima et al. 2013). Thus, deer are a good study subject to use for examining the factors affecting effort and efficiency. Effort with respect to deer hunting is expected to vary depending on landscape characteristics like the type of land use (e.g. forest, farmland, and urban area), slope, and road density. The type of land use will affect outlook for hunters that would affect the hunters' selection of hunting point and the slope and road density will affect the accessibility for hunters. Hunters will select the place where deer abundance is high because they expect to hunt easily in such places. Furthermore, if the percentage of wildlife protection area and urban area is large, hunters will avoid hunting in such regions. Efficiency with respect to deer hunting is expected to vary depending on the degree of slope.
In addition to landscape components, the method of hunting is also expected to affect effort and efficiency. In Japan, there are two major hunting methods: shooting and trapping. Shooting is the most common method, but it is expected that the use of guns may make deer extremely vigilant if the deer had previously been shot at and escaped. On the other hand, trapping is superior to shooting because a trapped deer will rarely escape. However, because hunters are legally required to check their traps every day in Japan, it is expected that hunting efficiency of trapping will be more rapidly saturate with the increase of deer abundance than the efficiency of shooting. In Japan, deer populations are managed both by nuisance control and by game hunting. Population control of deer in Japan is characterized by little participation of professional hunters, so both game hunting and nuisance control are conducted by amateur hunters, and nuisance control is conducted using the same methods of game hunting. Furthermore, Japanese hunters are aging and the number of active hunters is decreasing, which limits the number of hunted deer. For these reasons, examination of factors affecting effort and efficiency of game hunting is useful for population control of deer species.
To contribute to the effective management of overabundant deer populations, this study aimed to clarify the effects of landscape components, deer behavior, and hunting methods on variations in effort and efficiency of hunting deer. In this study, I analyzed effort and efficiency with respect to hunting sika deer Cervus nippon. I examined the hypotheses that 1) effort and efficiency of hunting are affected by landscape components, 2) efficiency of shooting, but not trapping, decreases with the increase of previous year's hunter effort, and 3) hunting efficiency of trapping will be more rapidly saturate with the increase of deer abundance than the efficiency of shooting.

Study site
The study site was the whole of Yamanashi Prefecture, central Japan. In Yamanashi Prefecture, the population density of sika deer has increased in recent years (Iijima and Ueno 2016), resulting in extensive debarking of trees (Nagaike andHayashi 2003, Iijima andNagaike 2015b) and browsing of understory vegetation in forests (Iijima and Nagaike 2015a) and herbaceous grasses in subalpine grasslands (Nagaike 2012). The hunting season in Yamanashi Prefecture is from November to March, and registered hunters hunt sika deer by shooting or trapping (snaring the legs). Furthermore, nuisance control of sika deer is conducted from April to November. The numbers of hunted deer by game hunting in 2005 was much higher than that of culled deer by nuisance control, but the number of hunted deer in 2010 was similar to that of culled deer (Fig. 1). Thus, game hunting also has important role in population control. Yamanashi Prefecture is characterized by complex topography. The elevation of Yamanashi Prefecture ranges from 36 m to the peak of Mt. Fuji at 3376 m.

Data
A rectangular grid (each cell was 5.5  4.6 km  25.3 km 2 ) was established to cover the entire prefecture (a total of 216 grid cells). The size and shape of the cells was set by the Ministry of the Environment of Japan and is uniform all over Japan. All data were obtained at the scale of a cell.
The hunter effort (the number of hunting days multiplied by the number of shooters or the number of days that traps were set) and the number of caught deer from 2005 to 2010 in each cell were obtained from hunters' reports. Hunters (both shooters and trappers) voluntarily submit their efforts and the number of deer caught to the local government office at the end of the hunting season. The submission of hunters' reports is voluntary, but  90% of registered hunters submit their reports in Yamanashi Prefecture because public officers repeatedly and persistently request the submission of these reports through the Yamanashi branch of the Japan Hunting Association, and the Japan Hunting Association has a strong influence on hunters in Japan. Therefore, these reports can be regarded as unbiased. In contrast, there is no system to submit hunters' report in nuisance control in Yamanashi Prefecture. Then, I could not analyze effort and efficiency of nuisance control. However, as already stated, nuisance control is also conducted by amateur hunters who go to hunting in the same region and the number of hunted deer by game hunting is similar or larger than that of culled deer by nuisance control. Then, in this study, I concentrate on the game hunting. Deer abundance from 2005 to 2010 was obtained from Iijima et al. (2013), who estimated deer abundance in each cell (mean and standard deviation) from three types of monitoring data and took into consideration measurement error and spatial autocorrelation of deer abundance. Please refer to Iijima et al. (2013) for further details of the methods of estimation. In this study, I treat deer abundance in each cell as a random variable with a normal distribution with the estimated mean and standard deviation.
For each of the cells, the areas (to account for cells located at the boundary of Yamanashi Prefecture, km 2 ) and the percentages of the cell occupied by evergreen forests, deciduous forests, wildlife protection areas, artificial grassland (e.g. farmland and golf course) areas, city areas, the mean slope of the cell, and the forest road density (forest road length within a cell/area of the cell) were calculated by the QGIS software package (< http://qgis.osgeo.org/en/site/ >, accessed 19 January 2016). The areas of the landscape components were obtained from the Natural Environment Information geographic information system (< www.biodic.go.jp/ trialSystem/top.html >, accessed 19 January 2016) provided by the Biodiversity Center of Japan. A digital elevation model with 10-m mesh was used to calculate mean slope; this was obtained from the Geospatial Information Authority of Japan (< http://fgd.gsi.go.jp/download/ >, accessed 19 January 2016). Forest road data were obtained from the Yamanashi prefectural government.

Statistical analysis
I adopted a state-space model (Calder et al. 2003) to examine the factors affecting effort and efficiency. The state-space model is composed of a process model, which describes latent processes like temporal changes in efficiency, and an observation model, which describes the relationship between latent processes and observed data. Thus, stochastic and observation errors could be incorporated explicitly. I analyzed effort and efficiency by hunting method (i.e. shooting and trapping).

Process model
The process model for hunter effort was as below: where Effort t,c is latent (true) effort in year t in cell c; b AB is the coefficient of D t-1,c ; D t-1,c is deer abundance in year t-1 in cell c; s 1 is an variance parameter; a 1 is an intercept term; b ER is the coefficient of the percentage of evergreen forest within cell c (ER c ); b DR is the coefficient of the percentage of deciduous forest within cell c (DR c ); b WR is the coefficient of the percentage of wildlife protection area within cell c (WR c ); b SL1 is the coefficient of mean slope within cell c (SL c ); b AR is the coefficient of area within cell c (AR c ); b GR is the coefficient of artificial grassland area within cell c (GR c ); b CR is the coefficient of city area within cell c (CR c ); b RD1 is the coefficient of forest road density within cell c (RD c ); and e 1 is the random effect of cell c. Prior distributions of b AB , a 1 , b ER , b DR , b WR , b SL1 , b AR , b GR , b CR , and b RD were a vague Gaussian distribution with mean 0 and variance 1000. Prior distribution of e 1 was a Gaussian distribution with mean 0 and variance σ 2 2 . Prior distributions of s 1 and s 2 were a uniform distribution from 0 to 100 (Gelman 2006 where Efficiency t,c is latent (true) efficiency in year t in cell c; b EF is the coefficient of effort in the previous year; s 3 is an variance parameter; a 2 is an intercept term; b SL2 is the coefficient of SL c ; and e 2 is the random effect of cell c. A negative value of b EF indicates that effort in the previous year decreased efficiency. In previous studies, efficiency has been determined as a mixture of target wildlife abundance, wildlife behavior, and landscape components, as noted in the Introduction. However, in this study, deer abundance was already known. Therefore, deer abundance was not included in the process model for efficiency and was treated in the observation model. Prior distributions of b EF , a 2 and b SL2 were a vague Gaussian distribution with mean 0 and variance 1000. Prior distribution of e 2 was a Gaussian distribution with mean 0 and variance σ 4 2 . Prior distributions of s 3 and s 4 were a uniform distribution from 0 to 100 (Gelman 2006).
where E t,c is the observed effort (the number of hunting days multiplied by the number of shooters or the number of days that traps were set) in year t in cell c; C t,c is the number of deer caught in year t in cell c; and b D is the strength of the deer abundance dependence of hunting efficiency (shape parameter; Tsuboi and Endou 2008). Prior distribution of b D was a uniform distribution from 0 to 2. If b D is smaller than 1, the number of deer caught will not increase linearly with an increase in deer abundance. As stated above, I treat deer abundance of each cell as a random variable with normal distribution with the estimated mean ( µ t c , ) and variance ( σ t c

Parameter estimation
Parameter estimation of the state-space model was conducted using a Bayesian framework. Posterior samples of parameters were obtained by the Markov chain Monte Carlo (MCMC) method (Calder et al. 2003). I ran three parallel MCMC chains and retained 200 000 iterations after an initial burn-in of 50 000 iterations. I thinned the sampled values to 0.5% (i.e. obtained 1000 samples as posterior distributions

Discussion
The factors affecting effort differed between shooting and trapping. Effort for shooting was high under high deer abundance, low percentages of wildlife protection area, and low percentages of city area (Table 1). The increase of effort by high deer abundance and the suppression of effort by wildlife protection areas and city area (Table 1) has been reported in other areas (Davidson andFraser 1991, Harden et al. 2005). Contrary to my prediction, effort of shooting was not affected by forest, artificial grassland, slope and road density (Table 1). It is common in Yamanashi Prefecture to use dogs to flush out deer for shooting, so it is not necessary for hunters to walk long distances to find deer. The general usage of dogs in shooting might also have contributed to why I did not detect an effect of forest, artificial grassland, slope and road density on effort. Increased percentage of wildlife protection areas decreased effort for trapping as well (Table 1). However, city area did not affect effort for trapping (Table 1). Trap hunters would like to set their traps near their home because hunters must patrol set traps every day in Japan. Hunters live in both city and rural areas. Then, there was no relationship between the city area and effort for trapping. The factors affecting efficiency also differed between shooting and trapping. Effort of the previous year decreased efficiency of shooting of the current year (Table 2, Fig. 2). Because deer abundance was incorporated into my model (Eq. 6) and considered estimation error (Eq. 7), the decrease in hunting efficiency by the previous year's high hunter effort for each chain). Convergence of MCMC sampling was judged by the criterion that R  was smaller than 1.1 (Gelman et al. 2004). To conduct MCMC sampling, I used JAGS (Plummer 2003) in the "rjags" package in R (< www. r-project.org >). I also used the package 'snow' to conduct multi-thread MCMC sampling. I concluded that a parameter was significant if the 95% credible intervals of the parameter did not overlap 0.

Results
Effort t,c for shooting significantly increased with the increase in D t,c and AR c and the decrease in WR c and CR c (Table 1). We failed to detect any significant difference in Effort t,c for shooting with the differences in ER c , DR c , SL c , GR c and RD c (Table 1). Effort t,c for trapping significantly increased with the increase in AR c and the decrease in WR c but was not affected by other factors (Table 1).
Efficiency t,c for shooting in a certain year significantly decreased following a high hunter effort in the previous year (Table 2, Fig. 2). During the study period (six years), Efficiency t,c for shooting of the cell where the sum of effort during six years was maximum among cells decreased more than 40%. Mean (95% credible interval) of b D for shooting was 0.79 (0.69-0.89; Table 2, Fig. 3). Efficiency t,c for trapping was not affected by any factors including hunter effort in the previous year (Table 2, Fig. 2). Mean (95% credible interval) of b D for trapping was 0.46 (0.22-0.74; Table 2, Fig. 3). Table 1. Posterior summary of estimated parameters for latent effort. b AB is the coefficient of D t,c ; D t,c is deer abundance in year t in cell c; b ER is the coefficient of the percentage of evergreen forest within cell c (ER c ); b DR is the coefficient of the percentage of deciduous forest within cell c (DR c ); b WR is the coefficient of the percentage of wildlife protection area within cell c (WR c ); b SL1 is the coefficient of mean slope within cell c (SL c ); b AR is the coefficient of area within cell c (AR c ); b GR is the coefficient of artificial grassland area within cell c (GR c ); b CR is the coefficient of city area within cell c (CR c ); and b RD1 is the coefficient of forest road density within cell c (RD c ). A parameter was considered significant if the 95% credible intervals of the parameter did not overlap 0. if too many traps are set. Therefore, trapping is not suitable for taking large numbers of deer.
The decrease in efficiency of shooting by the previous year's high effort also indicates the difficulty in using catch per unit effort (CPUE) as a deer abundance index. In previous studies about deer population dynamics, CPUE was sometimes used as the index of deer abundance (Noss et al. 2005, Uno et al. 2006. However, efficiency decreased by the increase in the previous year's effort even though the change of deer abundance was considered (Fig. 2). Thus, if the managers of deer abundance only monitor CPUE as the index of deer abundance, they can overestimate the effect of hunting which will lead to the failure of population control.
These characteristics of each hunting method clarified by my analysis indicate the most effective method of managing populations of overabundant deer. Intensive hunting by gun is recommended in the initial phase because the decrease of efficiency of shooting in relation to deer abundance was not strong (Table 2, Fig. 3). After deer abundance decreases and deer vigilance increases (Table 2, Fig. 2), trapping should be adopted for a sustainable hunting efficiency (Table 2, indicates a change in the behavior of deer. To the best of my knowledge, this is the first study to show directly a change in deer behavior as a result of human behavior (previous year's effort) although such a possibility has been suggested from ethological studies or the change in vegetation after hunting (Martin and Baltzinger 2002). In Yamanashi Prefecture, shooters use dogs to flush out deer, and it is known that the percentage of hunting failure per unit effort (where hunters encounter their quarry, but the quarry escapes) increases when using dogs (Godwin et al. 2013). Thus, the increase of effort increased the number of high-vigilance deer scared by dogs, resulting in lower efficiency the following year. In contrast to shooting, efficiency for trapping was not affected by the previous year's hunting effort (Table 2, Fig. 2). Basically, trapped deer cannot escape and so a change in deer behavior as a result of the previous year's effort would rarely occur. From this point of view, trapping is superior to shooting. However, the decrease of efficiency in relation to deer abundance (b D ) was stronger with trapping than with shooting (Table 2, Fig. 3). The number of traps that can be set is limited by law in Japan, and patrolling of set traps is difficult  Fig. 2). However, hunting efforts of shooting and trapping were low around wildlife protection area (Table 1). Then, control of nuisance deer around wildlife protection areas should be actively introduced. The combination of game hunting and nuisance control is an effective way of managing overabundant deer (Kaji et al. 2010). In conclusion, effort and efficiency were affected by landscape components, deer behavior, and hunting methods. For best results, it may be necessary to use of a variety of methods to control overabundant wildlife.