Winter Habitat Use of Sichuan Sika Deer Cervus nippon sichuanicus in the Tiebu Sika Deer Nature Reserve

Cheng Zhao1,2, Qin Shi1*, Gen Yang3, Youyuan Tan2, Songwen Tan2, Jiahong Li2 and Marwan M.A. Rashed2

1Yibin Key Laboratory of Zoological Diversity and Ecological Conservation, Yibin University, Yibin 644007, P. R. China.

2Faculty of Agriculture, Forestry and Food Engineering, Yibin University, Yibin 644007, P. R. China.

3Sichuan Ruoergai Forestry Bureau, Aba 624500, P. R. China.

ABSTRACT

Sichuan sika deer Cervus nippon sichuanicus is a threatened subspecies of sika deer. Winter is the critical season affecting the survival of Sichuan sika deer. To get strategic information about the relationship between Sichuan sika deer and the environment in the winter, we conducted our study on habitat use of Sichuan sika deer in the Tiebu Nature Reserve from January to February in 2019 and 2020. The results showed that Sichuan sika deer preferred shrub vegetation and southern slope. Compared to control plots, lower elevation, longer distance to forest edge, higher herb cover and concealing condition appeared in used plots. The results of the logistic regression models showed that herb cover and elevation were the key factors. It was indicated that diet, safety and temperature were the main factors affecting the habitat use of Sichuan sika deer in winter.

Article Information

Received 16 July 2021

Revised 18 September 2021

Accepted 01 October 2021

Available online 03 March 2022

(early access)

Published 07 November 2022

Authors’ Contribution

CZ is the main contributor of this research. QS is the person in charge of this project. YT, ST, JL and YG are contributors in the field work. MMAR revised the manuscript.

Key words

Cervus nippon sichuanicus, Habitat use, Tiebu sika deer nature reserve, Winter

DOI: https://dx.doi.org/10.17582/journal.pjz/20210716080715

* Corresponding author: [email protected]

0030-9923/2023/0001-333 $ 9.00/0

Copyright 2023 by the authors. Licensee Zoological Society of Pakistan.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

INTRODUCTION

Sichuan sika deer Cervus nippon sichuanicus is a subspecies of sika deer Cervus nippon, which belongs to the Cervinae, Cervidae, Artiodactyla (Corbet and Hill, 1991). Sichuan sika deer is distributed in a narrow range along the eastern edge of the Qinghai-Tibet Plateau. This deer is threatened by habitat loss and human disturbance (poaching, grazing, road building) (Zhao et al., 2014a). The wild population is about 1000 individuals (Qi et al., 2014). Accordingly, it is classified as a category i protected wild animal species in China, and listed as critically Endangered in the Red List of China’s Vertebrates (Jiang et al., 2016).

Habitat is crucial for wildlife survival, and understanding the habitat use of target species is useful for designing an animal protection plan (Khadka and James, 2016). Spatial temporal variation affects wildlife habitat use (Dupke et al., 2017). It has been shown that the key factors for habitat use by wild animals varied over different seasons (Chen et al., 2012; Rahman et al., 2017; Dargan et al., 2019), as the environmental pressure altered seasonally. Guo (2002) pointed that due to food availability, winter was the critical season affecting survival of Sichuan sika deer. However, till now, little research was focused on the winter habitat use by Sichuan sika deer, which severely restricts the development of Sichuan sika deer conservation.

To understand the relationship between the winter environmental pressure and Sichuan sika deer habitat use and to assist managers to design more effective protection plans, we conducted the research on the winter habitat use of Sichuan sika deer in Tiebu Sika Deer Nature Reserve. Our goals were to determine which factors had a greater influence on the habitat use of Sichuan sika deer in winter, to detect the habitat characteristics which Sichuan sika deer preferred in winter, and to improve protection plans for Sichuan sika deer conservation.

MATERIALS AND METHODS

Study area

The study area is in the Tiebu Sika Deer Nature Reserve (abbreviation, TNR), China, which is located on the eastern edge of the Qinghai-Tibet Plateau. Elevation ranges from 2,450 m to 3,800 m above sea level. The annual precipitation is about 650 mm, with only ten percent of rainfall appearing from November to March. The mean annual temperature is 6.7oC and extreme low temperature recorded is -18.9oC. The TNR covers about 260 km2 (Yang et al., 2012), which is a mosaic of forest, bush, and meadow. Forest is distributed mainly in shady slopes, bush distributed mainly on sunny slopes, and meadow occurs mainly on gentle areas (Zhao et al., 2014b).

Sampling method and variable design

We established transects along altitudinal gradients in the TNR, ensuring the representative vegetation types could be sampled with a transect (Morrison et al., 1992). Sampling plots were centered on signs of sika deer, were treated as used plots, with a distance not less than 100 m between them along transects. Control plots were established along transects at every 100 m change in elevation.

After each sampling point location was established, the vegetation type was determined. Then, seventeen biological and physical variables (Table I) were measured following the method of Wei et al. (2000), which was used in former sika deer research (Zhao et al., 2014b). We sampled 132 used plots and 62 control plots from January to February 2019 and 2020 (Fig. 1).

Data analyses

For categorical variables, exact goodness-of-fit tests were used to detect the difference between used and control plots, then post-hoc tests were conducted. Since we could not confirm that every variable in our research was normally distributed, both independent sample T tests and Mann–Whitney U tests were used to find the difference between used and control plots.

 

Table I. Description and definition of variables in the habitat sampling plots.

Variable	Definition and description
Vegetation type	Five categories: coniferous forest, broad-leaved forest, bush, meadow, cultivated field
Slope aspect	Aspect of each 400m2 plot, defined by five categories: northern slope (315–45°), eastern slope (45–135°), southern slope (135–225°), western slope (225–315°), and no slope aspect (the plot in a flat land)
Slope (°)	The slope measured at the center of the 400 m2 plot
Elevation (m)	The elevation measured at the center of the 400 m2 plot
Concealing condition	Mean greatest distance looking eastward, southward, westward, and northward at a height of 1.0m at the center of the 400 m2 plot, four categories: >40 m, 20–40 m, 10–20 m, and <10 m (as the concealing condition from low to high)
Tree height (m)	Average height of trees in 400 m2 plot
Tree density (n)	Number of trees in 400 m2 plot
Tree size (cm)	Average diameter at breast height of the trees in 400 m2 plot
Shrub height (m)	Average height of shrubs in two 20 m2 rectangular subplots
Shrub density (n)	Average number of shrubs in two 20 m2 rectangular subplots
Shrub cover (%)	Average percent shrub cover in two 20 m2 rectangular subplots
Herb height (cm)	Average height of herbs in five 1.0 m2 subplots
Herb cover (%)	Average percent herbs cover in five 1.0 m2 subplots
Distance to forest edge (m)	Estimated straight-line distance from the sampling plot to the nearest forest edge, if the plot was in the forest, record as +x, otherwise, record as –x
Distance to water (m)	Straight-line distance from the sampling plot centre to the nearest water source
Distance to roads (m)	Straight-line distance from the sampling plot centre to the nearest road
Distance to habitation (m)	Straight-line distance from the sampling plot centre to the nearest residential area

 

To identify the key factors affecting sika deer habitat use, information-theoretic approach (Burnham and Anderson, 2002) was used. This investigation has conducted a correlation analysis to test independence between variables firstly. For those variables with a correlation coefficient (Pearson) above 0.5, only one variable remained in the subsequent analysis in order to reduce multicollinearity (Zhang et al., 2011; Torres et al., 2012). Then, GLM models with binomial distribution and logit-link were built. We calculated the Akaike information criterion (corrected for sample size, AICc) to evaluate model fit, using the differences in AICc values (ΔAICc) to assess a candidate set of models of all of the possible combinations of variables. For models with ΔAICc<2, we used model averaging to obtain estimates of each regression coefficient and confidence limits to identify the key factors (Creel et al., 2014).

RESULTS

For categorical variables, both vegetation type and slope aspect showed significant differences between habitat use plots and control plots (Table II). The outcome of post-hoc tests indicated that sika deer preferred bush and south slope, avoiding broad-leaved forest (Table II). For continuous variables, elevation, distance to forest edge, shrub height, shrub cover, shrub density, and herb cover showed significant differences between habitat use plots and control plots (Table III).

The top ten models contained 7 variables, herb cover, shrub cover, elevation, shrub height, concealed condition, distance to forest edge, distance to water respectively, of which, herb cover and elevation appeared in all the ten models (Table IV). The models averaged coefficients showed that, Sichuan sika deer preferred habitat with higher herb cover and lower elevation (Table V).

 

Table II. Categorical variables in used plots and control plots of Cervus nippon sichuanicus in Tiebu Nature Reserve during winter.

Categories	Frequency (%)		P (post-hoc tests)
	Control	Used
Vegetation type
Coniferous forest	19(0.30)	36(0.27)	0.450
Broad-leaved forest	11(0.17)	7(0.05)	<0.001
Bush	18(0.29)	66(0.50)	<0.001
Meadow	12(0.19)	22(0.16)	0.509
Cultivated field	2(0.03)	1(0.008)	0.136
Slope aspect
North	17(0.27)	26(0.19)	0.051
East	15(0.24)	30(0.22)	0.761
South	11(0.17)	41(0.31)	<0.001
West	15(0.24)	28(0.21)	0.477
No	4(0.06)	7(0.05)	0.724
P<0.001

 

Table III. Continuous variables in used plots and control plots of Cervus nippon sichuanicus in Tiebu Nature Reserve during winter.

Variable	Control plots	Used plots	F	U
	Mean±SD	Mean±SD
Elevation	2833.94±181.02	2758.94±161.57	0.004	0.007
Slope	28.18±16.36	24.25±13.61	0.104	0.198
Distance to forest edge	-9.05±221.26	-104.67±249.19	0.011	0.010
Distance to water	368.82±287.55	304.05±258.55	0.118	0.115
Concealing condition	2.16±1.08	1.77±0.93	0.012	0.021
Tree height	11.52±12.76	8.39±12.12	0.101	0.074
Diameter at breast height	14.90±16.90	10.90±16.42	0.119	0.068
Tree density	6.06±8.25	4.23±7.47	0.126	0.058
Shrub height	1.61±1.03	2.17±1.00	0.000	0.001
Shrub cover	0.10±0.29	0.27±0.44	0.002	0.008
Shrub density	4.98±6.04	6.20±5.55	0.167	0.006
Herb height	17.23±10.22	19.48±11.60	0.191	0.460
Herb cover	0.40±0.49	0.73±0.44	0.000	0.000
Distance to roads	214.19±338.44	188.68±245.17	0.552	0.654
Distance to habitations	929.55±736.27	868.94±630.25	0.555	0.928

 

Table IV. Summary of the top 10 logistic regression model sets predicting the habitat use of Cervus nippon sichuanicus.

Model construction	K	AICc	ΔAICc	Wi
Herb cover + Shrub cover + Elevation + Shrub height	5	208.04	0	0.10
Herb cover + Shrub cover + Elevation + Concealing condition	5	208.15	0.11	0.09
Herb cover + Shrub cover + Elevation	4	208.49	0.45	0.08
Herb cover + Shrub cover + Elevation + Shrub height + Concealing condition	6	208.51	0.47	0.08
Herb cover + Shrub cover + Elevation + Shrub height + Distance to forest edge	6	208.82	0.78	0.06
Herb cover + Shrub height + Elevation + Distance to forest edge + Concealing condition	6	209.08	1.04	0.06
Herb cover + Shrub cover + Elevation + Distance to forest edge + Concealing condition	6	209.17	1.13	0.05
Herb cover + Shrub cover + Elevation + Shrub height + Distance to forest edge + Concealing condition	7	209.18	1.14	0.05
Herb cover + Shrub cover + Elevation + Shrub height + Distance to water	6	209.26	1.22	0.05
Herb cover + Shrub cover + Elevation + Distance to water	5	209.34	1.30	0.05

 

Table V. Parameter estimates, standard errors, and 95% confidence limits (CL) of the model-averaged for estimating the habitat use of Cervus nippon sichuanicus.

Variable	Estimate	SE	Lower 95% CL	Upper 95% CL	P
Intercept	7.358	3.126	-0.813	7.065	0.019
Herb cover	2.758	0.687	-0.179	1.552	0.000
Shrub cover	2.130	1.279	-0.333	2.891	0.097
Shrub height	0.189	0.242	-0.063	0.547	0.436
Elevation	-3.355e-3	1.146e-3	-2.980e-4	2.590e-3	0.004
Concealing condition	0.146	0.2101	-0.055	0.475	0.490
Distance to forest edge	-3.865e-4	7.594e-4	-1.974e-4	1.716e-3	0.612
Distance to roads	-1.510e-4	4.470e-4	-1.162e-4	1.010e-3	0.736
Distance to habitations	-7.777e-6	6.635e-5	-1.725e-5	1.500e-4	0.907
Shrub density	-9.198e-4	8.924e-3	-2.320e-3	0.020	0.918
Slope	-1.659e-4	2.530e-3	-6.578e-4	5.718e-3	0.948
Vegetation type (Broad-leaved forest)	-0.100	0.390	-0.101	0.882	0.798
Vegetation type (Bush)	-0.111	0.422	-0.110	0.955	0.793
Vegetation type (Meadow)	-0.104	0.411	-0.107	0.930	0.801
Vegetation type (Cultivated field)	-0.2006	0.8096	-0.210496	1.830	0.805

 

Note: regard vegetation type (coniferous forest) as a fixed reference factor.

 

DISCUSSION

Diet and safety are the basic needs of wild animals (Rozhnov et al., 2015; Fattebert et al., 2019). Food availability varies over time, and winter is the diet shortage season for Sichuan sika deer (Guo, 2002). Previous research indicated that meadow was the main foraging patch of Sichuan sika deer habitat (Zhao et al., 2014b). Our research showed that herb cover was the key factor affecting the habitat use of Sichuan sika deer in winter (Table V), which indicated the importance of the diet need in winter. The safety needs caused Sichuan sika deer to utilize habitat with better concealing condition (Table III). The preference of shrub vegetation (Table II) was due to the fact that shrub could not only provide food but also acted as a shelter patch (Zhao et al., 2014b).

It has been demonstrated that human disturbance is an important factor affecting wildlife habitat use (Filla et al., 2017; Fattebert et al., 2019), while the influence of roads and habitations of Sichuan sika deer habitat use were not found in our results (Tables III and V). This may have been caused by the food availability stress, which forced the deer to utilize the habitat close to roads and habitations, and the protection of Sichuan sika deer by local Tibetans leads to the harmonious coexistence between man and deer.

Compared with summer, winter habitat use of Sichuan sika deer appeared different. In winter, Sichuan sika deer preferred habitat at lower elevation and south slope aspect, which lacked in summer (Zhao et al., 2020). These two factors may be highly related to temperature. Previous study demonstrated that temperature could affect habitat use behaviour in many kinds of large herbivores, such as alpine ibex Capra ibex (Aublet et al., 2009), black-tailed deer Odocoileus hemionus columbianus (Bowyer and Kie, 2009), mouflon Ovis gmelini musimon (Bourgoin et al., 2011), and moose Alces alces (van Beest et al., 2012). The lower elevation and south slope aspect may mainly be caused by warmer habitat selection of Sichuan sika deer. In addition, distance to water is an important factor affecting Sichuan sika deer habitat use in the summer (Zhao et al., 2020), which had little to contribute on habitat use of Sichuan sika deer in winter (Tables III, V). The snow scattered on the mountains served as a water source.

CONCLUSIONS

Our results clearly showed that the meadow was important as a food patch for Sichuan sika deer in winter, while water was not a restricting factor. Our results also indirectly demonstrated the influence of temperature on habitat use of Sichuan sika deer. Even though road and habitations had little influence on the habitat use of Sichuan sika deer in this research, considering the huge group of free-ranging livestock and intensive firewood cutting, we suggest the higher protection of the shrub and meadow vegetation especially in winter.

ACKNOWLEDGMENTS

This study was supported by Key Fund Project of Sichuan Provincial Education Department (18ZA0543), Special Fund of Yibin University (2017RC04).

Statement of conflict of interest

The authors have declared no conflict of interest.

REFERENCES

Aublet, J.F., Festa-Bianchet, M., Bergero, D. and Bassano, B., 2009. Temperature constraints on foraging behaviour of male alpine ibex (Capra ibex) in summer, Oecologia, 159: 237-247. https://doi.org/10.1007/s00442-008-1198-4

Bourgoin, G., Garel, M., Blanchard, P., Dubray, D., Maillard, D., and Gaillard, J.M., 2011. Daily responses of mouflon (Ovis gmelini musimon x Ovis sp.) activity to summer climatic conditions. Can. J. Zool., 89: 765-773. https://doi.org/10.1139/z11-046

Bowyer, R.T. and Kie, J.G., 2009. Thermal landscapes and resource selection by blacktailed deer: implications for large herbivores. Calif. Fish Game, 95: 128-139.

Burnham, K., and Anderson, D., 2002. Model selection and multimodel inference: A practical information-theoretic approach. Springer, New York.

Chen, W., Wu, Q.G., Hu, J.C., Lu, X., and You, Z.Q., 2012. Seasonal habitat use of Chinese goral (Naemorhedus griseus) in a subtropical forest. Russ. J. Ecol., 43: 256-260. https://doi.org/10.1134/S106741361203006X

Corbet, G.B., and Hill, J.E., 1991. A world list of mammalian species. Oxford University, London.

Creel, S., Schuette, P., and Christianson, D., 2014. Effects of predation risk on group size, vigilance, and foraging behavior in an African ungulate community. Behav. Ecol., 25: 773-784. https://doi.org/10.1093/beheco/aru050

Dargan, J.R., Moriyama, M., Mella, V.S.A., Lunney, D., and Crowther, M.S., 2019. The challenge for koala conservation on private land: koala habitat use varies with season on a fragmented rural landscape. Anim. Conserv., 22: 543-555. https://doi.org/10.1111/acv.12487

Dupke, C., Bonenfant, C., Reineking B., Hable R., Zeppenfeld T., Ewald M., and Heurich M., 2017. Habitat selection by a large herbivore at multiple spatial and temporal scales is primarily governed by food resources. Ecography, 40: 1014–1027. https://doi.org/10.1111/ecog.02152

Fattebert, J., Morelle, K., Jurkiewicz, J., Ukalska, J., and Borkowski, J., 2019. Safety first: Seasonal and diel habitat selection patterns by red deer in a contrasted landscape, J. Zool., 308: 111-120. https://doi.org/10.1111/jzo.12657

Filla, M., Premier, J., Magg, N., Dupke, C., Khorozyan, I., Waltert, M., Bufka, L., and Heurich, M., 2017. Habitat selection by Eurasian lynx (Lynx lynx) is primarily driven by avoidance of human activity during day and prey availability during night, Ecol. Evol., 7: 6367-6381. https://doi.org/10.1002/ece3.3204

Guo, Y.S., 2002. Determination of sika deer’s food resources and loading capacity in Tiebu Nature Reserve Sichuan, China. Acta Theriol., 22: 257-262.

Jiang, Z.G., Jiang, J.P., Wang, Y.Z., Zhang, E., Zhang, Y.Y., Li, L.L., Xie, F., Cai, B., Cao, L., Zheng, G.M., Dong, L., Zhang, Z.W., Ding, P., Luo, Z.H., Ding, C.Q., Ma, Z.J., Tang, S.H., Cao, W.X., Li, C.W., Hu, H.J., Ma, Y., Wu, Y., Wang, Y.X., Zhou, K.Y., Liu, S.Y., Chen, Y.Y., Li, J.T., Feng, Z.J., Wang, Y., Wang, B., Li, C., Song, X.L., Cai, L., Zang, C.X., Zeng, Y., Meng, Z.B., Fang, H.X. and Ping, X.G., 2016. Red list of China’s vertebrates. Biol. Sci., 24: 500-551.

Khadka, K.K., and James, D.A., 2016. Habitat selection by endangered Himalayan musk deer (Moschus chrysogaster) and impacts of livestock grazing in Nepal Himalaya: Implications for conservation. J. Nat. Conserv., 31: 38-42. https://doi.org/10.1016/j.jnc.2016.03.002

Morrison, M.L., Marcot, B.G., and Mannan, R.W., 1992. Wildlife habitat relationships: Concepts and applications. University of Wisconsin, Madison.

Qi, W.H., Jiang, X.M., Yang, C.Z., and Guo, Y.S., 2014. Reproductive behavior of Sichuan sika deer, Acta Ecol. Sin., 34: 6548-6559. https://doi.org/10.5846/stxb201302210289

Rahman, D.A., Gonzalez, G., Haryono, M., Muhtarom, A., Firdaus, A.Y., and Aulagnier, S., 2017. Factors affecting seasonal habitat use, and predicted range of two tropical deer in Indonesian rainforest, Acta. Oecol., 82: 41-51. https://doi.org/10.1016/j.actao.2017.05.008

Rozhnov, V.V., Chistopolova, M.D., Lukarevskii, V.S., Hernandez-Blanco, J.A., Naidenko, S.V., and Sorokin, P.A., 2015. Home range structure and space use of a female Amur leopard, Panthera pardus orientalis (Carnivora, Felidae), Biol. Bull., 42: 821-830. https://doi.org/10.1134/S1062359015090095

Torres, R.T., Vorgós, E., Santos, J., Linnell, J.D.C., and Fonseca, C., 2012. Habitat use by sympatric red and roe deer in a Mediterranean ecosystem. Anim. Biol., 62: 351-366. https://doi.org/10.1163/157075612X631213

van Beest, F.M., van Moorter, B., and Milner, J.M., 2012. Temperature-mediated habitat use and selection by a heat-sensitive northern ungulate. Anim. Behav., 84: 723-735. https://doi.org/10.1016/j.anbehav.2012.06.032

Wei, F.W., Wang, Z.W., and Hu, J.C., 2000. Habitat use and separation between the giant panda and the red panda, J. Mammal., 80: 448–455. https://doi.org/10.1093/jmammal/81.2.448

Yang, C.Z., Xiao, Z., Guo, Y.S., Xiong, Y.Q., Zhang, X.Y., and Yue, B.S., 2012. Alarm signals of the Sichuan sika deer Cervus nippon sichuanicus, Zool. Sci., 29: 423-427. https://doi.org/10.2108/zsj.29.423

Zhang, Z.J., Swaisgood, R.R., Zhang, S.N., Nordstrom, L.A., Wang, H.J., Gu, X.D., Hu, J.C. and Wei, F.W., 2011. Old-growth forest is what giant pandas really need. Biol. Lett., 7: 403-406. https://doi.org/10.1098/rsbl.2010.1081

Zhao, C., Hu, J., Li, Y.H., Ran, J.H., Guan, J.K., Yang, C.Z., Xiong, Y.Q., Yue, B.S., 2014a. Habitat use of Sichuan Sika deer in forest, bush and meadows in the Tiebu Nature Reserve, Sichuan, China, Pakistan J. Zool., 46: 941-951.

Zhao, C., Li, Y.H., Li, D.Y., Guan, J.K., Xiong, Y.Q., and Hu, J., 2014b. Habitat suitability assessment of Sichuan sika deer in Tiebu Nature Reserve during periods of green and dry grass, Acta Ecol. Sin., 34: 135-140. https://doi.org/10.1016/j.chnaes.2014.03.001

Zhao, C., Li, Y.H., Shi, Q., Wang, B., Xu, K.N., Zhao, S.Q., Hu, J., and Xiong, Y.Q., 2020. Summer habitat selection of Cervus nippon sichuanicus in the Sichuan Tiebu sika deer Nature Reserve. Sichuan J. Zool., 39: 56-62.