生态学报  2020, Vol. 40 Issue (18): 6659-6671

文章信息

程志斌, 刘定震, 白加德, 钟震宇, 林润生, 田东晓, 王丽斌, 刘田
CHENG Zhibin, LIU Dingzhen, BAI Jiade, ZHONG Zhenyu, LIN Runsheng, TIAN Dongxiao, WANG Libin, LIU Tian
麋鹿鹿角脱落、群主更替、产仔的年节律及其环境影响因子
Circannual rhythms in antler casting, harem master alternation, calves birth of Pére David's deer (Elaphurus davidianus) and their environmental influencing factors
生态学报. 2020, 40(18): 6659-6671
Acta Ecologica Sinica. 2020, 40(18): 6659-6671
http://dx.doi.org/10.5846/stxb201912062653

文章历史

收稿日期: 2019-12-06
网络出版日期: 2020-07-12
麋鹿鹿角脱落、群主更替、产仔的年节律及其环境影响因子
程志斌1 , 刘定震2 , 白加德1 , 钟震宇1 , 林润生3 , 田东晓3 , 王丽斌1 , 刘田1     
1. 北京麋鹿生态实验中心, 北京 100076;
2. 北京师范大学生命科学学院, 教育部生物多样性与生态工程重点实验室, 北京 100875;
3. 北京市气象信息中心, 北京 100089
摘要: 于2014-2016年三个冬季(12月-翌年2月)收集了北京南海子麋鹿苑半散放麋鹿自然脱落的角,并观察和记录了2015-2017年发情期(5-9月)群主更替和2016-2018年产仔期(3-7月)麋鹿幼仔出生情况,结合2014-2018年年平均气温、季平均气温、月平均气温、年降雨量、雨季开始时间、种群密度等环境因子,对鹿角脱落、群主更替、产仔等繁殖特征的年节律及其环境影响因子进行了研究。结果表明:(1)麋鹿鹿角脱落、群主更替、产仔的年节律均存在年际差异。(2)鹿角脱落时间从12月开始,1月下旬或2月上旬结束。(3)发情期为5月下旬或6月上旬开始,9月上旬结束;2015-2017年发情期时间有延长的趋势。(4)产仔期为3月中旬或4月中旬开始,5月下旬或7月下旬结束。(5)麋鹿鹿角脱落、群主更替、产仔的年节律存在明显的同步关系,其中鹿角脱落开始时间、鹿角脱落高峰期、鹿角脱落结束时间、第一次发情期开始时间、群主更替高峰期、产仔期开始时间、产仔高峰期、产仔期结束时间与前一年度比较均出现同步提前的现象。(6)鹿角脱落年节律存在随着12月平均气温升高而提前的现象;产仔期开始时间和产仔高峰期存在随着前一年9月平均气温的升高而提前的现象。(7)鹿角脱落年节律表现出随着年降雨量的增多而提前的现象;第一次发情期开始时间、群主更替高峰期的年节律表现出随着前一年度年降雨量的增多而提前的现象。(8)麋鹿鹿角脱落、群主更替、产仔的年节律均不存在随着种群密度升高或降低而提前或延迟的现象。麋鹿繁殖特征的年节律是一个复杂的过程,受气候、营养、种群密度、纬度等环境因子的影响。
关键词: 麋鹿    年节律    繁殖    鹿角    群主    产仔    环境因子    
Circannual rhythms in antler casting, harem master alternation, calves birth of Pére David's deer (Elaphurus davidianus) and their environmental influencing factors
CHENG Zhibin1 , LIU Dingzhen2 , BAI Jiade1 , ZHONG Zhenyu1 , LIN Runsheng3 , TIAN Dongxiao3 , WANG Libin1 , LIU Tian1     
1. Beijing Milu Ecological Research Center, Beijing 100076, China;
2. Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China;
3. Beijing Meteorological Information Center, Beijing 100089, China
Abstract: In order to study the circannual rhythm of reproductive characteristics in Pére David's deer (Elaphurus davidianus) and its adaptability to the environment, the casted antler were collected during three annual winters (from December to February) of 2014 to 2016 in Beijing Milu Ecological Research Center (BMERC), Beijing, China. Then harem masters were observed and recorded during the rut time (from May to September) of 2015 to 2017. Offspring date of birth during the offspring period (from March to July) from 2016 to 2018 were also recorded. Combined with the environmental factors from 2014 to 2018 such as annual mean temperature, seasonal mean temperature, monthly mean temperature, annual precipitation, start and end of rainy season, and population density, we studied circannual rhythms in antler casting, harem master alternation, calves birth of Pére David's deer and their environmental influencing factors. The results showed that: (1) there were annual differences in the circannual rhythms antler casting, harem master alternation, calves birth of Pére David's deer. (2) Antler casting date started in December, ended in late January or early February. (3) There was a second estrus period in Pére David's deer population. The estrus period started in late May or early June, and ended in early September. The estrus duration from 2015 to 2017 were 91d, 102d, and 110d respectively, which had a tendency to prolong. (4) The offspring period started in mid-March or mid-April, and ended in late May or late July. (5) There was synchronization relationship among circannual rhythms in antler casting, harem master alternation, calves birth. Antler casting start date, median antler casting date, antler casting end date, rut start date, rut median date, offspring start date, median calving date and offspring end date were synchronization advanced than last year respectively. (6) The circannual rhythms in harem master alternation did not show the phenomenon of advance or delay with the increase or decrease of the annual mean temperature, season mean temperature or monthly mean temperature, while circannual rhythms in antler casting increased with the mean temperature in December, and offspring start date, median calving date increased with the mean temperature in September of the previous year. (7) The circannual rhythm in calves birth did not show a corresponding synchronous change with the precipitation factors. While circannual rhythms in antler casting showed advancement with the increase of annual precipitation. Rut start date and rut median date appeared to advance with the increase of annual precipitation in the previous year. (8) There was no synchronization advanced or delayed of circannual rhythms in antler casting, harem master alternation, calves birth as the population when the density increased or decreased. The circannual rhythm of the reproductive characteristics of Pére David's deer was a complex process, affected by the environmental factors such as climate, nutrition, population density, and latitude.
Key Words: Elaphurus davidianus    circannual rhythms    reproductive    antler    harem    offspring    environmental influencing factor    

地球上的光、温度、湿度等生态因子具有明显的周期性, 因而在生物的进化和发展中形成了昼夜节律、月节律、年节律等生物节律[1-3]。动物的年节律如繁殖、脱角、换羽(毛)、迁徙、冬眠等行为是一种复杂的生物学现象, 它是动物对光照、温度、湿度、水等非生物条件, 食物、种内关系和种间关系等生物条件的适应, 每个物种均有各自的生物节律特点, 来适应外界复杂而变动的外界环境条件[4-8]。因此, 通过动物的生物节律来研究动物与环境之间的适应性关系始终是生态学家和物候学家的研究热点[9-12]。繁殖活动是哺乳动物生命中最为重要的组成部分, 已有大量关于鸟类[13-15]和兽类[16-18]繁殖年节律会着温度、降水等环境变化而产生适应性变化的报道, 有蹄类作为哺乳动物中的重要类群, 其繁殖的年节律受光周期、气候、营养、纬度、种群密度等环境因子的影响[19-20]。然而结合具体的平均气温和降水量来探讨气候对哺乳动物繁殖特征年节律影响方面的报道较少。

各物种采取的繁殖策略是多种多样的, 包括从高度季节性到完全无季节性的繁殖模式[21]。生活在高纬度地区的鹿科动物属于前者, 雄鹿每年的繁殖周期包含发情期、身体恢复期、鹿角脱落期和长茸骨化期四个阶段, 雌鹿每年的繁殖周期包含发情期、怀孕期、产仔期和哺乳期四个阶段[22]。首先, 鹿角是鹿亚科(Cervinae)和空齿鹿亚科(Odocoileinae)动物特有的性二型特征[23], 是繁殖期雄性争夺交配权和雌性选择雄性的重要指标[24-25], 它是哺乳动物唯一可再生的器官[23], 高纬度地区鹿科动物的鹿角每年周期性脱换[23, 26]。其次, 占群繁殖, 是一雄多雌的哺乳动物雄性繁殖期的重要行为特征之一, 如川金丝猴(Rhinopithecus roxellana)[27]、象海豹(Mirounga Leonina)[28]、藏羚羊(Pantholops hodgsonii)[29]、普氏野马(Equus przewalskii)[30]等。第三, 雌性鹿类动物的每次的发情持续时间较短, 通常只有几个小时[31], 不易于观察, 因此更为直观的产仔周期成为研究雌鹿的繁殖年节律特征的首选[17, 32]

尽管关于鹿类动物年节律特征的研究有许多, 但主要为单一或者两个繁殖特征, 如马鹿[17, 33-34](Cervus elaphus)、狍[35](Capreolus pygargus)、驯鹿[36](Rangifer tarandus)等广布物种的角生长、发情或者产仔的年节律, 以及白尾鹿(Odocoileus virginianus)血液中激素的年节律变化[37-38]。而同一群体的多个繁殖特征的变化差异、相互之间的关系及其环境影响因子等方面的研究甚少。麋鹿(Elaphurus davidianus)为我国特有物种, 国家一级保护动物, IUCN红皮书列为野外灭绝, 属于典型的季节性繁殖动物, 雄性具角, 具有自身独特的脱角周期[23]。麋鹿的婚配制度为典型的一雄多雌的后宫制[39-40], 雄性麋鹿采用序列占群的繁殖策略[40], 通过顶角争斗获得群主地位和交配权[41-43], 目前尚未见到麋鹿年节律研究的相关报道。本文通过研究麋鹿鹿角脱落、群主更替、产仔等繁殖特征的年节律及其环境影响因子, 为深入了解哺乳动物在进化过程中对环境的适应性, 以及为保护麋鹿这一珍稀物种提供科学依据。

1 材料与方法 1.1 研究地概况

研究地为北京南海子麋鹿苑(39°50′N, 116°30′E), 位于北京城南10 km, 占地60 hm2, 麋鹿半散放活动的保护区面积约40 hm2, 为古代元明清皇家苑囿南海子的一部分, 是麋鹿模式种产地、我国最后灭绝地和我国第一处重引入地。苑内优势植物有牛筋草(Eleusine indica), 大画眉草(Eragrostis cilianensis), 马唐(Digitaria sanguinalis)和狗尾草(Setaria viridis)等; 年均气温为13.1℃, 1月平均气温为-3.4℃, 7月平均气温为26.4℃; 年降雨量约为600 mm[44]。苑内无河流经过, 为人工湖, 水域靠降雨和抽取地下水补给。

北京南海子麋鹿种群是麋鹿在国内灭绝后建立的第一个半散放迁地种群。1985年和1987年北京南海子麋鹿苑从英国乌邦寺共引入38只麋鹿。回归35年来, 北京南海子麋鹿苑分别向湖北石首、江西鄱阳湖、浙江慈溪、河北滦河上游等地输送536只麋鹿, 建立了39个迁地保护种群, 在长江中下游建立了麋鹿野生种群。北京南海子麋鹿种群维持在130—210只之间(雄雌性比约1 : 1), 采取半散放的饲养方式, 由于夏秋两季地面草本植物生物量有限, 一年四季均有人工补饲。

1.2 数据采集及处理 1.2.1 鹿角脱落数据采集及处理

2014—2016年, 每一年度掉角期12月—翌年2月, 每天上午和下午各一次在麋鹿生活区域巡视, 收集麋鹿自然脱落的角, 且在角上标记脱落日期。选择鹿角脱落开始时间、鹿角脱落高峰期、鹿角脱落结束时间作为鹿角脱落年节律因子。

脱角数量达到半数为脱角中位值日期[45], 由于脱角中位值日期处于脱角集中时间段, 为了便于讨论, 将其定义为鹿角脱落高峰期; 由于2岁龄鹿角(称笔杆角, 无分叉)在脱落后掉于草丛中难以收集, 因此本文将2叉(3岁龄)及2叉以上的鹿角全部脱落完毕的日期定义为鹿角脱落结束时间。

鹿角回收率(%)=(鹿角回收数/2岁以上雄性数×2)×100%

1.2.2 群主更替数据采集及处理

2015—2017年, 每年发情期5—9月, 一周观察4—7 d, 每天观察0.5—2.5 h, 使用双筒望远镜(SWAROVSKI 8×42WB)、相机(Canon550)记录群主的特征, 通过鹿角形状、耳标、耳缺、项圈、体型等区分个体和年龄。群主为占有和控制发情雌鹿群(后宫)的等级序位高的雄鹿[41]。群主有两种方式控制发情雌鹿群:一是直接控制发情雌鹿群, 二是控制发情场[40]。由于北京南海子麋鹿为半散放种群, 发情场地有限, 在水边泥滩的高大乔木下, 占有发情场地是群主最喜欢的方式。本研究中群主为占有固定发情场的群主, 而占有临时雌性群(指离开固定发情场外出觅食形成的临时雌性群, 它们通常1 h内会被群主驱赶回固定发情场或者被群主控制)的临时群主不计算在内。

将第一只群主开始占有和控制雌鹿群, 且占有固定发情场的日期定义为第一次发情期开始时间(即群体而言雄性发情开始时间)[17]; 第一次发情期间群主更替次数达到半数的日期为群主更替中位值日期[46], 由于群主更替中位值日期处于群主更替集中时间段, 为了便于讨论, 将其定义为群主更替高峰期; 对于整个种群而言, 将当年第一次发情期的群主占群更替结束的日期定义为第一次发情期结束时间(通常群主全身为夏毛, 未开始更换冬毛); 将当年第二次发情期的群主开始占群的日期定义为第二次发情期开始时间(通常群主夏毛已经脱落完毕, 正在长冬毛); 将当年第二次发情期的群主占群更替结束的日期定义为第二次发情期结束时间(即群体而言雄性发情结束时间)。

选择第一次发情期开始时间、群主更替高峰期、第一次发情期结束时间、第二次发情期开始时间、第二次发情期结束时间作为群主更替年节律因子。

1.2.3 产仔数据采集及处理

2016—2018年, 每一年产仔期3—7月, 每天使用望远镜观察麋鹿群, 尤其是离群即将分娩或者正在分娩的母鹿, 上午和下午各观察鹿群及其活动场所一次, 做好新生仔鹿出生记录。选择产仔期开始时间、产仔高峰期、产仔期结束时间作为产仔的年节律因子。

产仔数量达到半数的日期为产仔中位值日期[47], 由于产仔中位值日期处于产仔集中时间段, 为了便于讨论, 将其定义为产仔高峰期。

1.2.4 环境数据采集及处理

2014年1月—2018年7月气象数据由北京市气象信息中心提供, 站点为大兴站(站号54594, 国家一般气象站), 距离研究地点10.3 km; 2018年7—12月的气象数据取自中国地面国际交换站气候资料日值数据集, 站点为北京站(站号54511, 国家基本气象站), 距离研究地点3.7 km。选择气温(3—9月、12月的月平均气温, 春、夏、秋、冬季4个季节的季平均气温, 年平均气温)、降水(春季第一场雨、春季第一场透雨、雨季开始时间、雨季结束时间、年降雨量)作为气候因子。

春季第一场透雨参照姚佩珍和张强[48]北京春季第一场透雨指标的方法, 由于2014年、2015年、2017年、2018年春季降雨量均超过70 mm, 降水充沛, 第一场透雨日的标准按照降雨量超过或者约等于当旬轻旱的平均标准值统计。雨季开始和结束时间参照于晓澄等[49]华北雨季开始早晚的方法。

由于麋鹿半散放区域的面积不变, 以年末种群数量作为种群密度因子。

1.3 作图

按照实际日期制作每一年度麋鹿鹿角脱落、群主更替、产仔和日降水量情况变化图。日降水量存在的缺失值在图中不显示。

2 结果 2.1 麋鹿鹿角脱落、群主更替、产仔的年节律 2.1.1 鹿角脱落的年节律

表 1可知, 2014—2016年麋鹿角的回收率为81.37%—90.70%, 说明可根据实际所收集到的麋鹿角脱落日期统计鹿角脱落高峰期。图 1表明, 整个种群而言, 麋鹿角脱落的年节律存在差异, 2014—2016年的鹿角脱落开始时间分别为12月16日、12月7日、12月3日, 鹿角脱落高峰期分别为12月31日、12月26日、12月21日, 鹿角脱落结束时间分别为2月7日、2月4日、1月17日。说明北京南海子麋鹿种群鹿角脱落时间开始于12月, 结束于1月下旬或2月上旬; 2015年、2016年的鹿角脱落开始时间、鹿角脱落高峰期、鹿角脱落结束时间与前一年比较均出现同步提前的现象。

表 1 北京南海子麋鹿苑2014—2018年麋鹿年末存栏、脱角、群主更替次数、产仔情况 Table 1 Population at the end of the year, number of collected antler, alternate time of harem master, number of calves from 2014 to 2018 in Beijing Milu Ecological Research Center (BMERC)
年份Year 年末种群数/只
Population at the end of the year
2岁以上雄性数/只
Number of male above 2 years old
鹿角回收数/支
Number of collected antler
鹿角回收率/%
Ratio of collected antler
群主更替次数/次
Alternate time of harem master
幼仔出生数/只
Number of calves
2014 150 51 83 81.37 -
2015 157 ≈48 85 88.54 17 -
2016 147 ≈43 78 90.70 18 31
2017 181 ≈52 - - 31 42
2018 172 52 - - 48

图 1 北京南海子麋鹿苑2014—2016年度麋鹿角脱落情况 Fig. 1 Antler casting date from 2014 to 2016 in Beijing Milu Ecological Research Center (BMERC) A是鹿角脱落高峰期
2.1.2 群主更替的年节律

图 2表明, 雄性麋鹿存在第二次发情期的现象; 2017年出现两个固定发情场所, 其中第二个固定发情场6月7日出现, 6月28日消失。整个种群而言, 群主更替存在年节律的差异, 2015—2017年的第一次发情期开始时间分别为6月9日、5月30日、5月17日, 群主更替高峰期分别为7月22日、7月4日、6月28日, 第一次发情期结束时间分别为9月1日、8月12日、8月15日, 第二次发情期开始的时间分别为9月1日、8月15日、8月17日, 第二次发情期结束时间分别为9月8日、9月9日、9月4日; 2015—2017年发情期时间分别为91d、102d、110d。说明北京南海子麋鹿苑麋鹿发情期开始于5月下旬或6月上旬, 结束于9月上旬; 2016年和2017年的第一次发情期开始时间、群主更替高峰期与前一年比较均出现同步提前的现象; 而第一次发情期结束时间、第二次发情期开始时间和第二次发情期结束时间不存在同步提前的现象; 从2015年至2017年, 整个鹿群发情的时间有延长的趋势。

图 2 北京南海子麋鹿苑2015—2017年麋鹿群主更替情况 Fig. 2 Alternate date of harem master from 2015 to 2017 in BMERC
2.1.3 产仔的年节律

表 1图 3可知, 整个种群而言, 麋鹿产仔年节律存在差异, 产仔期开始时间2016—2018年分别为4月12日、3月28日、3月16日; 产仔高峰期分别为4月23日、4月17日、4月3日, 产仔期结束时间分别为7月25日、7月15日、5月25日; 产仔期分别长达104 d、109 d、70 d。2016年70.97%的幼仔集中在高峰期前后的13 d内出生, 2017年64.29%的幼仔集中在高峰期前后的21 d内出生, 2018年56.25%的幼仔集中在高峰期前后的19 d内出生。说明北京南海子麋鹿苑麋鹿产仔期开始于3月中旬至4月中旬, 结束于5月下旬至7月下旬; 2017年和2018年产仔期开始时间、产仔高峰期、产仔期结束时间与前一年比较均出现同步提前的现象; 麋鹿产仔同期化变弱。

图 3 北京南海子麋鹿苑2016—2018年度麋鹿产仔情况 Fig. 3 Offspring date of birth from 2016 to 2018 in BMERC A是产仔高峰期
2.2 环境因子

表 1可知, 北京南海子麋鹿苑2014—2018年麋鹿年末种群数量分别为150只、157只、147只、181只、172只。

2014—2018年的月平均气温见图 4; 年平均气温分别13.88℃、13.62℃、13.36℃、13.65℃、13.34℃; 春季平均气温分别为16.42℃、15.54℃、15.40℃、16.27℃、15.04℃, 夏季平均气温分别为26.30℃、26.17℃、26.37℃、26.08℃、27.51℃, 秋季平均气温分别为13.56℃、12.77℃、12.89℃、12.57℃、13.36℃。2014—2017年冬季(12月、翌年1—2月)平均气温分别为-0.13℃、-1.36℃、-0.11℃、-2.00 ℃。

图 4 北京南海子麋鹿苑2014—2018年月平均气温 Fig. 4 Monthly mean temperature from 2014 to 2018 in BMERC

2016年春季3—5月共7场降雨, 降雨量总和为30.7 mm, 属于干旱, 7场降雨均没有达到春季第一场透雨的标准。2014—2018年的年降雨量分别416.8 mm、505.8 mm、722.7 mm、682.9 mm、537.1 mm。由图 5可知, 2014—2018年的第一场春雨分别为3月28日、3月31日、4月2日、3月20日、3月17日; 2014年、2015年、2017年、2018年的春季第一场透雨分别为4月17日、3月31日、3月23日、4月5日; 2014—2018年的雨季开始时间分别为7月2日、7月17日、7月19日、7月4日、7月8日; 2014—2018年的雨季结束时间分别8月13日、8月31日、7月28日、8月27日、8月13日。

图 5 北京南海子麋鹿苑2014—2018年日降水量 Fig. 5 Precipitation from 2014 to 2018 in BMERC A是春季第一场雨; B是春季第一场透雨; C是雨季开始时间; D是雨季结束时间
2.3 鹿角脱落、群主更替、产仔的年节律之间的同步关系

图 13可知, 当2014—2016年麋鹿鹿角脱落开始时间、鹿角脱落高峰期和鹿角脱落结束时间逐年提前时, 2015—2017年第一次发情期开始时间、群主更替高峰期也随之逐年提前; 2016—2018年的产仔期开始时间、产仔高峰期和产仔期结束时间也紧接着逐年提前。2015—2017年的第一次发情期开始时间至翌年产仔期开始时间间隔分别为308 d、302 d、303 d, 平均间隔为304.3 d; 2015—2017年的群主更替高峰期至翌年产仔高峰期间隔分别为276 d、287 d、279 d, 平均间隔为280.7 d。说明2014—2016年的麋鹿鹿角脱落年节律因子(鹿角脱落开始时间、鹿角脱落高峰期和鹿角脱落结束期)、2015—2017年群主更替年节律因子(第一次发情期开始时间和群主更替高峰期)、2016—2018年的产仔年节律因子(产仔期开始时间、产仔高峰期和产仔期结束时间)存在明显的同步关系。

2.4 鹿角脱落、群主更替、产仔的年节律与环境因子之间的关系

图 15可知, 麋鹿鹿角脱落的年节律与秋季平均气温、冬季平均气温、年平均气温不存在随着平均气温升高或降低而提前或者延迟的现象; 2015—2017年群主更替的年节律与2014—2016年的年平均气温和秋季平均气温, 与本年度夏季平均气温不存在随着平均气温升高或降低而提前或者延迟的现象; 2016—2018年产仔的年节律与2015—2017年秋季平均气温、冬季平均气温和年平均气温, 与本年度春季平均气温不存随着平均气温升高或降低而提前或者延迟的现象。产仔期开始时间与3月平均气温, 产仔期开始时间和产仔高峰期与4月平均气温, 第一次发情期开始时间与5月平均气温, 第一次发情期开始时间和群主更替高峰期与6月平均气温, 群主更替高峰期与7月平均气温, 第一次发情期结束时间和第二次发情期开始时间与8月平均气温, 第二次发情期结束时间与9月平均气温均不存在随着平均气温升高或降低而提前或者延迟的现象。但鹿角脱落开始时间、鹿角脱落高峰期和鹿角脱落结束时间存在随着12月平均气温升高而提前的现象; 产仔期开始时间和产仔高峰期存在随着前一年9月平均气温的升高而提前的现象。说明除群主更替年节律外, 产仔年节律和鹿角脱落年节律表现出受气温因子影响的现象。

鹿角脱落开始时间、鹿角脱落高峰期和鹿角脱落结束时间均与年降雨量存在随着年降雨量的增加而提前的现象。2015—2017年第一次发情期开始时间、群主更替高峰期的年节律存在随着2014—2016年年降雨量的增加而提前的现象。2016—2018年产仔期开始时间和产仔高峰期不存在随着2015—2017年年降雨量的增加或减少而提前或延迟的现象。产仔期开始时间、产仔高峰期和第一次发情期开始时间均不存在随着春季第一场雨、春季第一场透雨的提前或者延迟而提前或者延迟的现象。群主更替高峰期不存在随着雨季开始的提前或者延迟而提前或者延迟的现象; 第一次发情期结束时间、第二次发情期开始时间和第二次发情期结束时间均不存在随着雨季开始时间、雨季结束时间的提前或者延迟而提前或者延迟的现象。2015年的群主更替高峰期出现在雨季期间, 2016年和2017年的群主更替高峰期出现在雨季开始之前。说明, 产仔年节律没有表现出与降水因子的变化而出现相应的同步变化, 鹿角脱落年节律表现出随着年降雨量的增多而提前的现象; 第一次发情期开始时间、群主更替高峰期的年节律表现出随着前一年度年降雨量的增多而提前的现象。

图 13表 1可知, 麋鹿鹿角脱落年的节律、群主更替的年节律、产仔的年节律均不存在随着种群密度升高或降低而提前或者延迟的现象。

3 讨论 3.1 麋鹿鹿角脱落、群主更替、产仔的年节律

高纬度地区的有蹄类, 均具有各自的繁殖节律特征来适应高度季节性的生境。本研究中, 北京南海子麋鹿种群鹿角脱落时间从12月开始, 1月下旬或2月上旬结束, 这与极地物种驯鹿、欧亚大陆高纬度地区的广布物种狍的鹿角脱落时间相一致, 而马鹿、梅花鹿(Cervus nippon)、黇鹿(Dama dama)等大部分鹿科动物为春夏之际开始脱角[50]。有蹄类喜欢选择在春夏之际气候和食物条件最适的时候集中产仔, 研究显示麋鹿产仔期为3—7月, 集中于4—5月, 这与温带地区和寒带地区的其他鹿类相一致[50]。本研究表明, 麋鹿发情期为5月下旬或6月上旬开始, 9月上旬结束, 麋鹿存在二次发情期现象; 这与李春旺等[42]研究结果相同。本研究发现群主更替高峰期与翌年产仔高峰期间隔平均为280.7 d, 而人工圈养的麋鹿妊娠期平均为280 d[51], 因此可以根据前一年的群主更替高峰期推测第二年的产仔高峰期。有研究表明, 黑尾鹿(Odocoileus hemionus)发情高峰期可以预测翌年产仔时间的情况[52]; 可通过马鹿鹿角脱落年节律来预测紧接着的雄性发情年节律[17]

高纬度地区鹿类动物的鹿角脱落、发情、产仔等繁殖周期具有同步性[53]。雄鹿为了提高等级序位, 增加自身繁殖后代的贡献率, 可通过使角提前骨化而先进入发情占群期[46]。南非一些反刍动物在反捕食策略上表现为产仔行为和交配行为同步[54]。本研究发现麋鹿鹿角脱落、群主更替、产仔的年节律存在明显的同步关系, 其中鹿角脱落开始时间、鹿角脱落高峰期、鹿角脱落结束时间、第一次发情期开始时间、群主更替高峰期、产仔期开始时间、产仔高峰期、产仔期结束时间与前一年度比较均出现同步提前的现象。这与马鹿繁殖特征的年节律相似[33, 55], Clements等[46]对野生马鹿长达28年的研究表明, 马鹿群体鹿角脱落开始时间、鹿角脱落结束时间、雄性发情开始时间、雄性发情高峰期、雄性发情结束时间、雌性发情开始时间和产仔期开始时间等七个繁殖特征之间均呈显著正相关关系, 并呈线性逐年同步提前现象[17]

3.2 环境影响因子 3.2.1 气候因子

光周期具有恒定性和规律性, 随着地理纬度和季节而很有规律地变动, 它是引起动物繁殖活动年节律变化的直接原因[3, 19, 56], 通过神经——垂体——性腺轴控制动物繁殖期[42, 57]。根据动物繁殖与光照长短的关系, 可将动物分为长日照动物和短日照动物, 羚羊类和鹿类只有在白昼逐步缩短的秋冬之际才开始性腺活动, 它们属于短日照动物[3], 然而麋鹿是鹿类中唯一的夏季发情的短日照发情动物[42, 50, 53], 这与本研究的结果相一致, 麋鹿发情期为5月下旬—9月上旬, 发情高峰期集中在光照时长最长的6—7月。通过适应光周期的变化, 温带地区的鹿类动物产生繁殖年节律, 它们可以优化调整产仔时间, 使其与最适的环境气候因子(如食物)相同步[58]。本研究表明, 麋鹿的鹿角脱落、群主更替、产仔等繁殖特征存在明显的年节律, 恰好说明它们受光周期的影响, 形成周期性变化的年节律。有研究表明, 将南半球的普度鹿(Pudu pudu)迁移至北半球, 它的繁殖期会很快与当地的光周期同步[58]

光在决定动物繁殖年节律中起着重要作用, 但并不是说它是唯一的信号, 温度、湿度、降水等气候原因对动物繁殖年节律也有重要的影响[3]。非洲的有蹄类如斑马(Equus quagga)的产仔年节律同环境的气温、降水量的变动有关[59]。地中海马鹿雌性发情开始时间在更加干旱和炎热的春季会推迟; 在更加凉爽和降雨量多的秋季会提前[60]。我们的研究中, 麋鹿的鹿角脱落、群主更替、产仔等繁殖特征连续3个年度的年节律不是固定不变的, 而是存在年际差异, 说明它们受到其他因子的影响。

温度可以影响动物的生活史特点, 如使狍的产仔年节律提前[35]。气候变暖伴随着有效积温的逐年升高, 使得苏格兰Rum岛植物的生长量增加, 给当地马鹿带来更多的食物, 导致鹿角脱落开始时间、鹿角脱落结束时间、产仔期开始时间、雄性发情开始时间、雄性发情结束时间呈逐年提前的趋势[17, 46]。而且马鹿鹿角开始脱落的时间与前一年秋季温度有关, 因为秋冬季发情结束之后的食物多少与秋季温度密切相关[61], 而且前一年9月和当年4月平均气温升高会使马鹿鹿角脱落年节律提前[62]。我们的研究结果也显示产仔期开始时间和产仔高峰期存在随着前一年9月平均气温的升高而提前的现象, 鹿角脱落年节律存在随着每年12月平均气温升高而提前的现象。

降水是陆生动物的重要环境因子, 降水会改变哺乳动物繁殖的年节律, 本研究表明, 麋鹿第一次发情期开始时间、群主更替高峰期的年节律表现出随着前一年度年降雨量的增多而提前的现象。雌性驯鹿在4月降水量大的年份, 表现出产仔期开始时间推迟[18]; 但是本研究中, 产仔年节律未表现出随着降水因子的变化而发生相应的同步变化。有研究认为, 当冬季降雪量增多时, 第二年马鹿鹿角脱落开始时间会推迟[62-63], 本研究中麋鹿鹿角脱落年节律表现出随着年降雨量的增多而提前的现象。

3.2.2 营养和种群密度因子

营养(或食物)也是动物繁殖年节律变化的重要原因[59]。植物提前生长或者增加生物量, 可为食草动物在关键的繁殖期提供更多的食物[64-65]。增加或者提前的春季食物可使雌鹿妊娠期和雄鹿鹿角脱落时间提前; 植物生长时间延长, 可提高哺乳期和长茸期食物的质量, 让它们提前进入秋季的发情期[17]。动物能够根据产仔的营养需要, 调整产仔时间, 将产仔时间与环境气候最适时间(如食物供应高峰期)同步[35, 52, 63], 提高幼崽存活率[32]。本研究中, 2015—2017年麋鹿鹿群发情时间有延长的趋势, 2016—2018年麋鹿产仔同期化变弱, 分析原因, 可能与北京南海子麋鹿种群为半散放状态, 食物营养充分, 食物压力较小有关。圈养动物与野生动物的产仔年节律不同, 受气候和食物等因子的影响, 圈养的马麝(Moschus chrysogaster)产仔同期化比野生种群更弱[66]。白尾鹿的动物园种群产仔期开始时间为2月, 而野生种群4月才开始产仔[25]。南美草原鹿(Ozotoceros bezoarticus), 阿根廷干旱草原生境中的种群全年产仔, 产仔80%集中在10月至翌年4月出生; 而同纬度相邻的乌拉圭半干旱大陆性气候生境中的种群为季节性产仔, 产仔同期化更加集中, 高峰期推迟至夏季末或秋季初食物充沛的时期; 而热带环境下的巴西种群, 由于明显的旱季和雨季, 食物变化更为显著, 产仔同期化最显著, 产仔高峰期为8—9月, 与当地的雨季完全同步[32]。但是, 有学者认为, 这一同步现象不是立即发生的[67-68], 有时会存在时间误差, 对于狍而言它的胚胎滞育策略就显得尤为重要[35]

野外条件下食物和生存空间有限, 种群密度大, 意味着个体之间的竞争力大。地中海马鹿雌性发情开始时间在种群密度大的年份会推迟; 在种群密度小的年份会提前[60]。苏格兰Rum岛马鹿种群也存在高密度种群会推迟角脱落开始时间和角脱落结束时间的现象[62]。本研究表明, 麋鹿鹿角脱落、群主更替、产仔的年节律均不存在随着种群密度升高或降低而提前或延迟的现象, 也许和半散放状态下, 每天人工补饲, 不存在食物压力有关。

3.2.3 纬度

随着纬度的变化鹿类动物繁殖年节律会随之变化[52]。江苏大丰麋鹿脱角开始时间和结束时间早于北京南海子种群[23]。本研究中, 北京南海子麋鹿群主更替高峰期为6—7月, 而低纬度地区的江苏大丰麋鹿种群雄性发情高峰为提前为5—6月[42]。这与其他鹿类动物相似, 墨西哥的白尾鹿发情开始时间为7月, 而美国和加拿大地区为10月[25], 高纬度和低纬度地区马鹿发情和产仔年节律研究的结果也表现出这一特征[33]。然而, 本研究发现麋鹿鹿角脱落开始时间、发情开始时间、产仔期开始时间逐年提前, 而江苏大丰麋鹿种群2004—2010年鹿角脱落开始时间、发情开始时间、产仔期开始时间出现了延后的现象[69], 与本研究结果相反, 也许与江苏大丰为滨海湿地生境, 而北京南海子为内陆湿地生境有关。

种群繁殖特征年节律是一个复杂的规律, 受气候、营养、种群密度等诸多因子影响, 甚至种群内个体间的相互作用及各饲养区的管理模式也会对它造成影响[17, 60, 66]。本研究的数据为3个年度3种繁殖特征的年节律, 时间相对较短, 对麋鹿繁殖特征与环境因子之间的关系还需要更长时间的观察研究。

致谢: 数据收集中得到北京麋鹿生态实验中心王建昌、包海元的帮助, 气象数据分析中得到国家气候中心中国气象局气候研究开放实验室赵俊虎的帮助, 特此致谢。
参考文献
[1]
Liu Q S, Li J Y, Wang D H. Ultradian rhythms and the nutritional importance of caecotrophy in captive Brandt's voles (Lasiopodomys brandtii). Journal of Comparative Physiology B, 2007, 177(4): 423-432. DOI:10.1007/s00360-006-0141-4
[2]
Helm B, Lincoln G A. Circannual rhythms anticipate the earth's annual periodicity//Kumar V, ed. Biological Timekeeping: Clocks, Rhythms and Behaviour. New Delhi: Springer, 2017: 545-569.
[3]
孙儒泳, 王德华, 牛翠娟, 刘定震, 张立. 动物生态学原理(第四版). 北京: 北京师范大学出版社, 2019: 61-67.
[4]
Chen J F, Zhong W Q, Wang D H. Seasonal changes in body mass, energy intake and thermogenesis in Maximowiczi's voles (Microtus maximowiczii) from the Inner Mongolian grassland. Journal of Comparative Physiology B, 2012, 182(2): 275-285. DOI:10.1007/s00360-011-0608-9
[5]
马瑞俊, 蒋志刚. 全球气候变化对野生动物的影响. 生态学报, 2005, 25(11): 3061-3066.
[6]
Goymann W, Helm B. Seasonality of life histories in tropical birds: circannual rhythms and Zeitgeber//Numata H, Helm B, eds. Annual, Lunar, and Tidal Clocks. Tokyo: Springer, 2014: 247-275.
[7]
Sokolov L V, Tsvey A L. Mechanisms controling the timing of spring migration in birds. Biology Bulletin, 2016, 43(9): 1148-1160. DOI:10.1134/S1062359016110145
[8]
Hazlerigg D, Blix A S, Stokkan K A. Waiting for the Sun:the circannual programme of reindeer is delayed by the recurrence of rhythmical melatonin secretion after the arctic night. Journal of Experimental Biology, 2017, 220(21): 3869-3872. DOI:10.1242/jeb.163741
[9]
王文君, 杨万勤, 谭波, 刘瑞龙, 吴福忠. 四川盆地亚热带常绿阔叶林不同物候期凋落物分解与土壤动物群落结构的关系. 生态学报, 2013, 33(18): 5737-5750.
[10]
李佳, 李言阔, 缪泸君, 谢光勇, 袁芳凯, 黄燕, 许鹏. 越冬地气候条件对鄱阳湖自然保护区白琵鹭种群数量的影响. 生态学报, 2014, 34(19): 5522-5529.
[11]
Loe L E, Bonenfant C, Mysterud A, Gaillard J M, Langvatn R, Klein F, Calenge C, Ergon T, Pettorelli N, Stenseth N C. Climate predictability and breeding phenology in red deer:timing and synchrony of rutting and calving in Norway and France. Journal of Animal Ecology, 2005, 74(4): 579-588. DOI:10.1111/j.1365-2656.2005.00987.x
[12]
Visser M E, Caro S P, Van Oers K, Schaper S V, Helm B. Phenology, seasonal timing and circannual rhythms:towards a unified framework. Philosophical Transactions of the Royal Society B:Biological Sciences, 2010, 365(1555): 3113-3127. DOI:10.1098/rstb.2010.0111
[13]
Price T, Kirkpatrick M, Arnold S J. Directional selection and the evolution of breeding date in birds. Science, 1988, 240(4853): 798-799. DOI:10.1126/science.3363360
[14]
Crick H Q P, Dudley C, Glue D E, Thomson D L. UK birds are laying eggs earlier. Nature, 1997, 388(6642): 526.
[15]
Dunn P O, Winkler D W. Climate change has affected the breeding date of tree swallows throughout North America. Proceedings of the Royal Society B:Biological Sciences, 1999, 266(1437): 2487-2490. DOI:10.1098/rspb.1999.0950
[16]
Lane J E, Kruuk L E B, Charmantier A, Murie J O, Dobson F S. Delayed phenology and reduced fitness associated with climate change in a wild hibernator. Nature, 2012, 489(7417): 554-557. DOI:10.1038/nature11335
[17]
Moyes K, Nussey D H, Clements M N, Guinness F E, Morris A, Morris S, Pemberton J M, Kruuk L E B, Clutton-Brock T H. Advancing breeding phenology in response to environmental change in a wild red deer population. Global Change Biology, 2011, 17(7): 2455-2469. DOI:10.1111/j.1365-2486.2010.02382.x
[18]
Paoli A. Breeding Phenology of a Semi-Domesticated Reindeer (Rangifer tarandus) Population in Response to Climatic Variability[D]. Montreal: Concordia University, 2019.
[19]
Curlewis J D, Loudon A S I, Coleman A P M. Oestrous cycles and the breeding season of the Père David's deer hind (Elaphurus davidianus). Reproduction, 1988, 82(1): 119-126. DOI:10.1530/jrf.0.0820119
[20]
Ungerfeld R, González-Sierra U T, Piaggio J. Reproduction in a semi-captive herd of pampas deer Ozotoceros bezoarticus. Wildlife Biology, 2008, 14(3): 350-358. DOI:10.2981/0909-6396(2008)14[350:RIASHO]2.0.CO;2
[21]
Asher G W. Reproductive cycles of deer. Animal Reproduction Science, 2011, 124(3/4): 170-175.
[22]
Nussey D H, Kruuk L E B, Morris A, Clements M N, Pemberton J M, Clutton-Brock T H. Inter-and intrasexual variation in aging patterns across reproductive traits in a wild red deer population. The American Naturalist, 2009, 174(3): 342-357. DOI:10.1086/603615
[23]
程志斌, 白加德, 钟震宇. 麋鹿鹿角生长周期及影响因子. 生态学报, 2016, 36(1): 59-68.
[24]
Pozo R A, Schindler S, Cubaynes S, Cusack J J, Coulson T, Malo A F. Modeling the impact of selective harvesting on red deer antlers. The Journal of Wildlife Management, 2016, 80(6): 978-989. DOI:10.1002/jwmg.21089
[25]
Ahuja-Aguirre C, Pérez-Gutiérrez O, López-deBuen L, Rojas-Maya S, Montiel-Palacios F. Early onset of the reproductive season in captive white-tailed deer (Odocoileus virginianus veraecrucis) in eastern Mexico. Journal of Zoo and Wildlife Medicine, 2017, 48(4): 1031-1038. DOI:10.1638/2016-0148.1
[26]
Sempéré A J, Boissin J. Relationship between antler development and plasma androgen concentrations in adult roe deer (Capreolus capreolus). Journal of Reproduction and Fertility, 1981, 62(1): 49-53. DOI:10.1530/jrf.0.0620049
[27]
Li B G, Zhao D P. Copulation behavior within one-male groups of wild Rhinopithecus roxellana in the Qinling Mountains of China. Primates, 2007, 48(3): 190-196. DOI:10.1007/s10329-006-0029-7
[28]
McMahon C R, Howe H, Van Den Hoff J, Alderman R, Brolsma H, Hindell M A. Satellites, the all-seeing eyes in the sky:counting elephant seals from space. PLoS One, 2014, 9(3): e92613. DOI:10.1371/journal.pone.0092613
[29]
Shi Y, Chen J R, Su J P, Zhang T Z, Wasser S K. Genetic resilience of a once endangered species, Tibetan antelope (Pantholops hodgsonii). bioRxiv, 2019. DOI:10.1101/628727
[30]
Burger D, Dolivo G, Wedekind C. Ejaculate characteristics depend on social environment in the horse (Equus caballus). PLoS One, 2015, 10(11): e0143185. DOI:10.1371/journal.pone.0143185
[31]
钟立成, 卢向东, 李艳冰, 张仁明. 半散放东北马鹿人工授精技术. 经济动物学报, 2002, 6(3): 19-22.
[32]
Ungerfeld R, González-Pensado S, Bielli A, Villagrán M, Olazabal D, Pérez W. Reproductive biology of the pampas deer (Ozotoceros bezoarticus):a review. Acta Veterinaria Scandinavica, 2008, 50(1): 16. DOI:10.1186/1751-0147-50-16
[33]
Bonenfant C, Loe L E, Mysterud A, Langvatn R, Stenseth N C, Gaillard J M, Klein F. Multiple causes of sexual segregation in European red deer:enlightenments from varying breeding phenology at high and low latitude. Proceedings of the Royal Society B:Biological Sciences, 2004, 271(1542): 883-892. DOI:10.1098/rspb.2003.2661
[34]
Pettorelli N, Mysterud A, Yoccoz N G, Langvatn R, Stenseth N C. Importance of climatological downscaling and plant phenology for red deer in heterogeneous landscapes. Proceedings of the Royal Society B:Biological Sciences, 2005, 272(1579): 2357-2364. DOI:10.1098/rspb.2005.3218
[35]
Plard F, Gaillard J M, Coulson T, Hewison A J M, Delorme D, Warnant C, Bonenfant C. Mismatch between birth date and vegetation phenology slows the demography of roe deer. PLoS Biology, 2014, 12(4): e1001828. DOI:10.1371/journal.pbio.1001828
[36]
Pettorelli N, Weladji R B, Holand Ø, Mysterud A, Breie H, Stenseth N C. The relative role of winter and spring conditions:linking climate and landscape-scale plant phenology to alpine reindeer body mass. Biology Letters, 2005, 1(1): 24-26. DOI:10.1098/rsbl.2004.0262
[37]
Bubenik G A, Bubenik A B, Schams D, Leatherland J F. Circadian and circannual rhythms of LH, FSH, testosterone (T), prolactin, cortisol, T3 and T4 in plasma of mature, male white-tailed deer. Comparative Biochemistry and Physiology Part A:Physiology, 1983, 76(1): 37-45. DOI:10.1016/0300-9629(83)90289-X
[38]
Bubenik G A, Smith P S. Circadian and circannual rhythms of melatonin in plasma of male white-tailed deer and the effect of oral administration of melatonin. Journal of Experimental Zoology, 1987, 241(1): 81-89. DOI:10.1002/jez.1402410110
[39]
Li C W, Jiang Z G, Zeng Y, Yan C E. Relationship between serum testosterone, dominance and mating success in Père David's deer stags. Ethology, 2004, 110(9): 681-691. DOI:10.1111/j.1439-0310.2004.01003.x
[40]
李春旺, 蒋志刚, 曾岩, 游章强. 不同密度下和不同发情期时段中雄性麋鹿的繁殖计策. 生物多样性, 2005, 13(5): 424-431.
[41]
Jiang Z G. Age-dependent rut strategy in Milu. Ethology, 1999, 84(S1): 168.
[42]
李春旺, 蒋志刚, 房继明, 姜国华, 丁玉华, 沈华, 徐安红. 麋鹿繁殖行为和粪样激素水平变化的关系. 兽类学报, 2000, 20(2): 88-100.
[43]
蒋志刚, 李春旺, 曾岩. 麋鹿的配偶制度、交配计策与有效种群. 生态学报, 2006, 26(7): 2255-2260.
[44]
李春旺, 蒋志刚, 曾岩, 阎彩娥, 张林源, 夏经世, 唐保田. 圈养雄性麋鹿血清睾酮和皮质醇含量在发情期的变化. 动物学研究, 2003, 24(1): 49-52.
[45]
Johansson A. Territory establishment and antler cycle in male roe deer. Ethology, 1996, 102(4): 549-559.
[46]
Clements M N, Clutton-Brock T H, Guinness F E, Pemberton J M, Kruuk L E B. Variances and covariances of phenological traits in a wild mammal population. Evolution, 2011, 65(3): 788-801. DOI:10.1111/j.1558-5646.2010.01161.x
[47]
Patterson B R, Mills K J, Middel K R, Benson J F, Obbard M E. Does predation influence the seasonal and diel timing of moose calving in central Ontario, Canada?. PLoS One, 2016, 11(4): e0150730. DOI:10.1371/journal.pone.0150730
[48]
姚佩珍, 张强. 北京地区春旱和春季第一场透雨指标的确定. 气象, 1999, 25(5): 49-53.
[49]
于晓澄, 赵俊虎, 杨柳, 支蓉, 封国林. 华北雨季开始早晚与大气环流和海表温度异常的关系. 大气科学, 2019, 43(1): 107-118.
[50]
Goss R J. Deer Antlers:Regeneration, Function and Evolution. New York:Academic Press, 1983, 8-46.
[51]
白加德. 麋鹿生物学研究. 北京: 北京科学技术出版社, 2014: 177-177.
[52]
Bowyer R T. Timing of parturition and lactation in southern mule deer. Journal of Mammalogy, 1991, 72(1): 138-145. DOI:10.2307/1381988
[53]
Bubenik G A. Seasonal versus non-seasonal reproduction in deer: From the arctic to the tropics//Bartoš L, Dušek A, Kotrba R, Bartošová-Víchová J, eds. Advances in Deer Biology: Deer in a Changing World. Prague: Research Institute of Animal Production, 2006: 16-17.
[54]
Skinner J D, Van Jaarsveld A S. Adaptive significance of restricted breeding in southern African ruminants. South African Journal of Science, 1987, 83: 657-663.
[55]
Lincoln G A. Biology of antlers. Journal of Zoology, 1992, 226(3): 517-528. DOI:10.1111/j.1469-7998.1992.tb07495.x
[56]
Lu W Q, Meng Q J, Tyler N J C, Stokkan K A, Loudon A S I. A circadian clock is not required in an arctic mammal. Current Biology, 2010, 20(6): 533-537. DOI:10.1016/j.cub.2010.01.042
[57]
Reiter R J. The pineal and its hormones in the control of reproduction in mammals. Endocrine Reviews, 1980, 1(2): 109-131.
[58]
Vidal F, Smith-Flueck J A M, Flueck W T, Bartoš L. Variation in reproduction of a temperate deer, the southern pudu (Pudu puda). Animal Production Science, 2012, 52(8): 735-740. DOI:10.1071/AN11364
[59]
Leuthold W, Leuthold B M. Temporal patterns of reproduction in ungulates of Tsavo East National Park, Kenya. African Journal of Ecology, 1975, 13(3/4): 159-169.
[60]
Peláez M, San Miguel A, Rodríguez-Vigal C, Perea R. Climate, female traits and population features as drivers of breeding timing in Mediterranean red deer populations. Integrative Zoology, 2017, 12(5): 396-408. DOI:10.1111/1749-4877.12252
[61]
Clutton-Brock T H, Guinness F E, Albon S D. Red Deer:Behaviour and Ecology of Two Sexes. Chicago:University of Chicago Press, 1982, 1-377.
[62]
Clements M N, Clutton-Brock T H, Albon S D, Pemberton J M, Kruuk L E B. Getting the timing right:antler growth phenology and sexual selection in a wild red deer population. Oecologia, 2010, 164(2): 357-368. DOI:10.1007/s00442-010-1656-7
[63]
Clutton-Brock T H, Albon S D, Guinness F E. Fitness costs of gestation and lactation in wild mammals. Nature, 1989, 337(6204): 260-262. DOI:10.1038/337260a0
[64]
Post E, Stenseth N C. Climatic variability, plant phenology, and northern ungulates. Ecology, 1999, 80(4): 1322-1339. DOI:10.1890/0012-9658(1999)080[1322:CVPPAN]2.0.CO;2
[65]
Mysterud A, Bonenfant C, Loe L E, Langvatn R, Yoccoz N G, Stenseth N C. The timing of male reproductive effort relative to female ovulation in a capital breeder. Journal of Animal Ecology, 2008, 77(3): 469-477.
[66]
孟秀祥, 杨奇森, 冯祚建, 施哲梅, 蒋应文. 马麝(Moschus sifanicus)的分娩时间格局. 四川动物, 2005, 24(1): 30-35.
[67]
Stenseth N C, Mysterud A. Climate, changing phenology, and other life history traits:nonlinearity and match-mismatch to the environment. Proceedings of the National Academy of Sciences of the United States of America, 2002, 99(21): 13379-13381. DOI:10.1073/pnas.212519399
[68]
Visser M E, Both C, Lambrechts M M. Global climate change leads to mistimed avian reproduction. Advances in Ecological Research, 2004, 35: 89-110. DOI:10.1016/S0065-2504(04)35005-1
[69]
侯立冰, 丁晶晶, 丁玉华, 任义军. 气候变化与麋鹿生存之间的关系. 北京林业大学学报, 2011, 33(S2): 11-15.