Seasonal and annual fluctuations of deer populations estimated by a Bayesian state–space model

Deer overabundance is a contributing factor in the degradation of plant communities and ecosystems worldwide. The management and conservation of the deer-affected ecosystems requires us to urgently grasp deer population trends and to identify the factors that affect them. In this study, we developed a Bayesian state–space model to estimate the population dynamics of sika deer (Cervus nippon) in a cool-temperate forest in Japan, where wolves (Canis lupus hodophilax) are extinct. The model was based on field data collected from block count surveys, road count surveys by vehicles, mortality surveys during the winter, and nuisance control for 12 years (2007–2018). We clarified the seasonal and annual fluctuation of the deer population. We found two peaks of deer abundance (2007 and 2010) over 12 years. In 2011 the estimated deer abundance decreased drastically and has remained at a low level then. The deer population increased from spring to autumn and decreased from autumn to winter in most years. The seasonal fluctuation we detected could reflect the seasonal migration pattern of deer and the population recruitment through fawn births in early summer. In our model, snowfall accumulation, which can be a lethal factor for deer, may have slightly affected their mortality during the winter. Although we could not detect a direct effect of snow on population dynamics, snowfall decrease due to global warming may decelerate the winter migration of deer; subsequently, deer staying on-site may intensively forage evergreen perennial plants during the winter season. The nuisance control affected population dynamics. Even in wildlife protection areas and national parks where hunting is regulated, nuisance control could be effective in buffering the effect of deer browsing on forest ecosystems.