低温对云南盈江北部高黎贡白眉长臂猿鸣叫行为的影响.pdf

四川动物 Sichuan Journal of Zoology 2023，42（4）：400-406低温对云南盈江北部高黎贡白眉长臂猿鸣叫行为的影响朱永亮1#，周佳1#，李金发2，杨国韬2，张利祥3，李彭3，崔亮伟1*，管振华1*（1.西南林业大学云南省高校极小种群野生动物保育重点实验室，昆明650224；2.德宏州生态环境局盈江分局，云南盈江679300；3.云南铜壁关省级自然保护区管护局，云南芒市678400）摘要：2017年12月2018年11月，首次对云南省德宏傣族景颇族自治州盈江县北部7个高黎贡白眉长臂猿Hoolock tianxing家庭和2只独猿开展了连续1年的鸣声监测。共记录352次鸣叫，高黎贡白眉长臂猿平均鸣叫起始时间在日出后117.97 min67.1 min（n=352），鸣叫持续时间平均为21.52 min12.8 min（n=352），独猿和家庭的日鸣叫频次分别为1.1次/天和0.424次/天。广义线性混合模型分析结果显示，日最低温度对高黎贡白眉长臂猿的鸣叫行为有显著影响，长臂猿在低温月份会推迟鸣叫起始时间，但日最低温度对鸣叫持续时间和日鸣叫频次的影响不显著；监测结果表明占据领域的独猿的鸣叫频次高于其他家庭。本研究反映了在生境破碎化栖息地中，低温对高黎贡白眉长臂猿鸣叫行为有显著影响，建议未来开展持续性种群监测和栖息地恢复工作。关键词：高黎贡白眉长臂猿；鸣叫行为；日最低温度；配对状态；盈江县中图分类号：Q958.1 文献标志码：A 文章编号：1000-7083（2023）04-0400-07Impact of Low Temperature on the Singing Behavior of Hoolock tianxingin Northern Yingjiang County,Yunnan ProvinceZHU Yongliang1#，ZHOU Jia1#，LI Jinfa2，YANG Guotao2，ZHANG Lixiang3，LI Peng3，CUI Liangwei1*，GUAN Zhenhua1*（1.Key Laboratory for Conserving Wildlife with Small Population in Yunnan，Southwest Forestry University，Kunming 650224，China；2.Yingjiang Branch of Dehong Prefecture Ecological Environment Bureau，Yingjiang，Yunnan Province 679300，China；3.Yunnan Tongbiguan Provincial Nature Reserve Management and Protection Bureau，Mangshi，Yunnan Province 678400，China）Abstract:From December 2017 to November 2018，a one-year continuous singing monitoring was conducted on 7 Hoolock tianxing family groups and 2 solitary individuals occupying territories in the northern part of Yingjiang County，Yunnan.The aim of the study was to investigate the singing behavior of H.tianxing and its relationship with environmental factors in their habitat.A total of 352 singing behaviors were recorded and analyzed，revealing that the average start time of the singing behavior was 117.97 min67.1 min（n=352）after sunrise，with an average duration of 21.52 min12.8 min（n=352）.The daily singing frequency of family groups was 0.424 times on average，while that of solitary individuals was 1.1 times on average.Generalized linear mixed model analysis was performed to investigate the impact of environmental factors on the singing behavior.The results showed that the daily minimum temperature had a significant impact on the start time of DOI：10.11984/j.issn.1000-7083.20220286收稿日期：2022-08-13 接受日期：2023-04-13基金项目：云南省科技厅基础研究面上项目（202101AT070165）；云南省高层次人才培养支持计划（YNWR-CYJS-2018-052）作者简介：朱永亮（1999），男，硕士研究生，研究方向：野生动物生态与保护，E-mail：；周佳（2000），硕士研究生，研究方向：动物生态学，E-mail：#同等贡献作者*通信作者Corresponding author，E-mail：cuilweastern-；400朱永亮等：低温对云南盈江北部高黎贡白眉长臂猿鸣叫行为的影响the singing behavior，with delayed onset during colder months.This suggests that H.tianxing may adopt an energy conservation strategy during cold periods.However，the effect of daily minimum temperature on singing duration and daily singing frequency was not significant.In addition，the study showed that the singing frequency of solitary individuals occupying territories was higher than that of other family groups，reflecting the attracting function of the singing behavior.Based on our study，we have observed the impact of low temperature and habitat fragmentation on the singing behavior of H.tianxing.Our findings suggest that in order to effectively protect the H.tianxing，it is recommended to implement continuous population monitoring and habitat restoration initiatives.Keywords:Hoolock tianxing；singing behavior；daily minimum temperature；pairing status；Yingjiang County长臂猿科 Hylobatidae 属小型猿类，栖息于东南亚的热带雨林和亚热带常绿阔叶林中，大部分长臂猿为一夫一妻小家庭群体，占据固定领域生活（Geissmann，2002）。鸣叫是长臂猿科动物的典型行为，除穆氏长臂猿Hylobates muelleri和克氏长臂猿H.klossii外，其他长臂猿的配对个体都可发出雌雄配合的二重唱（Geissmann，2002）。长臂猿的鸣叫行为主要发生在日出后的清晨，受温度、光照、捕猎风险等影响。东黑冠长臂猿Nomascus nasutus在雾天的鸣叫起始时间推迟，雌性激动鸣叫次数减少（费汉榄等，2010）；大风天气会抑制婆罗洲白须长臂猿 H.albibarbis（Cheyne，2008）和东黑冠长臂猿（费汉榄等，2010）的鸣叫行为；降雨带来的环境温度下降会推迟长臂猿的鸣叫起始时间（合趾猿Symphalangus syndactylus：Chivers，1974；克氏长臂猿：Whitten，1982；穆氏长臂猿：Clink et al.，2020）；捕猎风险会抑制高黎贡白眉长臂猿H.tianxing的鸣叫行为（Yin et al.，2016）。长臂猿的鸣叫行为具有多种功能（Cowlishaw，1992）：向周边家庭宣告对领域的占领，以及宣告配 对 个 体 之 间 的 稳 定 性 关 系 等（Fan et al.，2009；Ma et al.，2020，2022）；对于未配对长臂猿，鸣叫有吸引配偶的功能（Henwood&Fabrick，1979），是联系潜在配偶建立家群的重要途径。当长臂猿到达亚成年后，通常会离开出生群成为游离独猿（Brockelman et al.，1998；Hu et al.，2018），独猿通过2种方式寻找配偶建立家群：一种是取代原住的成年同性个体（范朋飞等，2010），另一种是离开家庭，取代其他家庭的成年同性个体（Brockelmann et al.，1998；Hu et al.，2018），或者和新配偶占据新的领域（Bryant et al.，2016）。当配对关系不稳定时，独猿可以取代原住个体（范朋飞等，2010），已有证据表明雌性独猿可以通过鸣叫识别出雄性个体的配对情况，选择是否与其建立家庭（Gittins，1984；Raemaekers et al.，1984；Mitani，1988；Mitani&Marler，1989；Reichard&Neumann，2007）。此外，在有新的独猿个体时，独猿可通过鸣叫联系并建立家群（Bryant et al.，2016）。未配对的游离独猿发出更高频次、更长时间的鸣叫有利于其寻找配偶或者被其他异性独猿发现。在对白掌长臂猿H.lar的研究中，游离独猿和未配对亚成年个体都比成年个体发出的鸣叫频次更高（Raemaekers et al.，1984）。我国分布有长臂猿 3 属 7 种（魏辅文等，2021），均为国家一级重点保护野生动物，目前仅分布于云南、广西和海南的斑块化生境中。高黎贡白眉长臂猿曾广泛分布于怒江以西，受偷猎和栖息地质量退化等威胁，现仅分布于云南省高黎贡山国家级自然保护区（腾冲大塘和保山隆阳区）、保山市腾冲县猴桥镇和德宏傣族景颇族自治州盈江县 3 个片区（Fan et al.，2011；Chan et al.，2017；Fan et al.，2017；Zhang et al.，2020）；2017年的种群数量与分布调查表明，国内高黎贡白眉长臂猿个体数量不超过150只（Zhang et al.，2020），被列入世界上最濒危的25种灵长类之一（Schwitzer et al.，2019；Fan et al.，2020）。关于高黎贡白眉长臂猿的401四川动物 2023 年 第 42 卷 第 4 期Sichuan Journal of Zoology Vol.42 No.4 2023鸣叫行为研究仅限于腾冲大塘、保山赧亢和板厂的种群，高黎贡白眉长臂猿全年的鸣叫频次较低（0.25次/天和0.50次/天，张岛等，2011；0.20次/天，0.074 次/天，0.313 次/天和 0.097 次/天，Yin et al.，2016）。张岛等（2011）初步报道了腾冲大塘的高黎贡白眉长臂猿鸣叫行为的月变化，发现降雨有可能导致鸣叫起始时间推迟或者不发生鸣叫。Yin等（2016）对比赧亢和板厂高黎贡白眉长臂猿的鸣叫行为发现，鸣叫频次在寒冷季节会降低，且捕猎和低家庭密度会抑制鸣叫行为。分布在盈江县北部的高黎贡白眉长臂猿种群数量占全国的一半以上（1820个家庭/独猿），该区域也是高黎贡白眉长臂猿重要的种群分布地（Zhang et al.，2020），大部分长臂猿家庭分布于保护区外。傈僳族对长臂猿的原始崇拜维持了这一种群的生存，目前存在最大威胁是栖息地破碎化和人类活动干扰（Fan，2017；Zhang et al.，2020）。本研究基于保护区外的 9 个家庭/独猿连续 1年的鸣声监测结果，报道了盈江县北部高黎贡白眉长臂猿鸣叫的时间特征，分析了其鸣叫行为与环境温度的关系，对比了独猿和家庭鸣叫行为之间的差异，填补了盈江县高黎贡白眉长臂猿鸣叫行为的研究空白，为其监测保护工作提供科学依据。1研究方法1.1研究地点与研究对象研究地点位于云南省德宏傣族景颇族自治州盈江县北部的苏典乡和支那乡（97459812E，25022521N，海拔1 6952 243 m），区域内高黎贡白眉长臂猿栖息地森林破碎化严重，种群之间被村寨和道路隔离。2017年4月620日采用鸣声定位法开展了2周的长臂猿种群数量与分布前期调查，明确了长臂猿种群的具体分布位点。基于调查结果，分别在苏典乡的腊马河和梨树2个寨子，支那乡的香柏、中山坝和白岩3个寨子附近选择8个听点对周边明确领域位置的9个长臂猿家庭/独猿进行监测。1.2调查监测方法监测时间为2017年12月19日2018年11月29日。每个听点2名监测队员，每月中下旬连续监测10 d。监测队员每天在天亮后30 min内到达听点，并持续监听至12 00。调查区域的日出时间使用R中的suncalc package（Thieurmel&Elmarhraoui，2019）获得，监听过程中记录鸣叫信息，包括：听点GPS坐标、天气（多云、阴、雨、雾、微风、大风等）、使用两步路（户外助手）APP结合卫星影像图估计鸣叫地点的方位角（沿听点的指北方向线起依顺时针方向至鸣叫方向线间的水平夹角）和听点与鸣叫位置的直线距离、鸣叫起始时间和鸣叫结束时间。鸣叫间隔定义为10 min以上，即鸣叫结束后10 min之后再次鸣叫，计为2次鸣叫。1.3温度数据数据来源于开放数据集国家生态科学数据中心19792018年中国近地表气温数据集（Fang et al.，2022），类型为栅格数据，空间分辨率为0.1。通过经纬度提取2017年12月1日2018年12月1日监测家庭/独猿所在位置的每日最低温度、每日最高温度和每日平均温度。1.4数据分析为确保结果的准确性，选择一年内监测到鸣叫次数 20 次以上的 9 个家庭/独猿作为研究对象（表1），按月份计算每个家庭的日鸣叫频次，即每个月的鸣叫总次数/监测总天数，使用 Kruskal-Wallis非参数检验鸣叫频次的群间差异，并使用多重比较检验分析不同家庭/独猿间的鸣叫频次差异，显著性水平设置为=0.05。采用广义线性混合模型（GLMM；Bolker et al.，2009）分析日最低温度和配对状态对鸣叫起始时间、鸣叫持续时间的影响，模型中设置鸣叫起始时间、持续时间为Gamma分布，家庭编号设置为随机变量。计算不同月份中日最低温度的平均值作为月平均最低温度，采用 GLMM分析月平均最低温度对鸣叫频次的影响，设置鸣叫频次为Gamma分布，家庭编号设置为随机变量。通过 R中的 lme4 402朱永亮等：低温对云南盈江北部高黎贡白眉长臂猿鸣叫行为的影响package 进行 GLMM 分析（Bates et al.，2015），所有统计检验均在R 4.2.0中完成（R Core Team，2022）。2研究结果2.1鸣叫的时间特征共监测170 d，平均每月监测14.2 d（1022 d）。监测过程中共获得352次鸣叫记录显示，平均鸣叫起始时间在日出后 117.97 min67.1 min（n=352），鸣叫持续时间平均为21.52 min12.8 min（n=352）。独猿的日鸣叫频次为1.1次/天0.545次/天（n=105），家庭的为0.424次/天0.355次/天（n=603）。2.2长臂猿鸣叫起始时间和温度的关系GLMM结果显示，日最低温度显著影响长臂猿的鸣叫起始时间（Z=7.138，P0.001），随日最低温度上升，鸣叫起始时间提前，但日最低温度并不显著影响鸣叫持续时间（Z=0.418，P=0.676）（图1）。月 平 均 最 低 温 度 对 鸣 叫 频 次 的 影 响 不 显 著（Z=0.487，P=0.627）；配对状态显著影响了鸣叫频次，独猿鸣叫频次更高（Z=3.081，P=0.002），但对鸣叫起始时间（Z=0.040，P=0.968）和鸣叫持续时间（Z=0.114，P=0.909）的影响不显著。2.3家庭/独猿间鸣叫频次差异鸣叫频次在不同家庭/独猿间存在极显著差异（2=28.968，df=8，P0.001），独猿（腊马河-1和腊马河-2）的鸣叫频次较高，且显著高于梨树-1、白岩-1和白岩-2（P0.05）（图2）。3讨论本研究显示温度对长臂猿鸣叫行为有显著影响。以往在高黎贡山的研究表明，高黎贡白眉长臂猿鸣叫行为主要发生在日出后 2 h 左右（张岛等，2011；Yin et al.，2016；Chan et al.，2017）。盈江县高黎贡白眉长臂猿的平均鸣叫起始时间在日出后的117.97 min，与上述研究结果相似。高纬度栖息地的温度季节性变化明显，在冬季，长臂猿面临低温和食物资源匮乏的生存挑战，采取能量保守的行为策略（Fan et al.，2013），在清晨延迟离开过夜树的时间（Fei et al.，2019），因此推迟了鸣叫起始时间。相比之下，低纬度地区的长臂猿鸣叫起始时间更早，南黄颊长臂猿Nomascus gabriellae的平均鸣叫起始时间为 06 14（Rawson，2004）；95%戴 帽 长 臂 猿 Hylobates pileatus 的 在 09 00 以 前（Brockelman&Srikosamatara，1993）。表1云南盈江北部高黎贡白眉长臂猿监测对象编号和配对状态Table 1ID and pairing status of Hoolock tianxingin northern Yingjiang County，Yunnan村寨Village腊马河梨树香柏中山坝白岩编号ID腊马河-1腊马河-2腊马河-3梨树-1香柏-1香柏-2中山坝-2白岩-1白岩-2配对状态Pairing status雌性独猿雄性独猿家庭家庭家庭家庭家庭家庭家庭图1日最低温度对云南盈江北部高黎贡白眉长臂猿鸣叫起始时间的影响Fig.1Impact of daily mimimum temperature on the start time of singing behavior of Hoolock tianxingin northern Yingjiang County，Yunnan403四川动物 2023 年 第 42 卷 第 4 期Sichuan Journal of Zoology Vol.42 No.4 2023Yin 等（2016）通过对比不同种群密度长臂猿群体的鸣叫行为，发现低种群密度会降低鸣叫频次。本次监测的长臂猿家庭被道路、河流、村寨、农业种植等分隔为3个小种群，所有家庭均分布于低种群密度的栖息地中。虽然腊马河种群的种群密度最高，但家庭的鸣叫频次并没有显著高于其他家庭。除了高纬度的低温影响外，捕猎风险也是抑制鸣叫行为的重要因素，Yin等（2016）推断长臂猿在高捕猎风险环境中为了自我保护，鸣叫频次低于 0.1 次/天。这一结果明显低于本研究（0.424次/天）。虽然本次监测的对象大部分分布在保护区外，但因在傈僳族传统生物多样性保护文化中占据特殊地位，面临捕猎的风险极低（Zhang et al.，2020），故鸣叫频次更高。长 臂 猿 的 鸣 叫 行 为 是 一 种 领 域 行 为（Cowlishaw，1992），独猿的鸣叫频次和个体的配对状态有关，多数情况下，独猿的鸣叫频次极低（范朋飞，个人通信），但当出现取代和配对机会时，会增加鸣叫频次（Raemaekers et al.，1984；范朋飞等，2010；Bryant et al.，2016）。本次监测的腊马河种群是盈江县北部最大的种群，包括3个家庭和2只独猿，雌性独猿（腊马河-1）和雄性独猿（腊马河-2）均占据固定领域，由于存在潜在异性配偶，独猿显著提高鸣叫频次以达到吸引异性配偶建立家庭的目的。虽然2只独猿均能听到对岸的鸣叫声，但由于栖息地被公路和河流阻隔，无法建立家庭。本次监测一方面增加了对种群的直接保护；另一方面，监测数据的积累有助于分析物种的生存现状。在破碎化生境中，除了独猿外，隔离小种群的鸣叫频次显著下降，影响了种群间的交流，也潜在影响到长臂猿的基因交流。盈江县高黎贡白眉长臂猿目前面临的最大威胁就是栖息地的破碎化和质量下降，开展长期种群监测（管振华等，2017；Liu et al.，2022）、通过栖息地恢复和廊道建设加强不同种群的连接，将是未来高黎贡白眉长臂猿保护工作的重点内容。致谢：感谢德宏州生态环境局盈江分局对本项目的支持；感谢云南铜壁关省级自然保护区管护局、盈江县林业和草原局的全力配合和帮助。参考文献:范朋飞，黄蓓，管振华，等.2010.西黑冠长臂猿雄性取代前后鸣叫行为的变化 J.兽类学报，30（2）：139-143.费汉榄，范朋飞，向左甫，等.2010.东黑冠长臂猿鸣叫特征及气象因子对鸣叫的影响 J.兽类学报，30（4）：377-383.管振华，阎璐，黄蓓.2017.中国长臂猿科动物种群监测现状分析 J.四川动物，36（2）：232-238.魏辅文，杨奇森，吴毅，等.2021.中国兽类名录（2021版）J.兽类学报，41（5）：487-501.张岛，袁胜东，崔亮伟，等.2011.云南高黎贡山大塘东部白眉长臂猿鸣叫行为 J.四川动物，30（6）：856-860.Bates D，Mchler M，Bolker B，et al.2015.Fitting linear mixed-effects models using lme4J.Journal of Statistical Software，67（1）：1-48.Bolker BM，Brooks ME，Clark CJ，et al.2009.Generalized 图2云南盈江北部高黎贡白眉长臂猿鸣叫行为的多重检验比较Fig.2Dunn test of singing behavior amongHoolock tianxing in northern Yingjiang County，Yunnan深色柱代表独猿，浅色柱代表家庭，不同字母代表组间差异显著（P0.05）The dark bars represent solitary individuals，the light bars represent families，different letters indicate there is a significant difference between two groups（P0.05）404朱永亮等：低温对云南盈江北部高黎贡白眉长臂猿鸣叫行为的影响linear mixed models：a practical guide for ecology and evolutionJ.Trends in Ecology&Evolution，24（3）：127-135.Brockelman WY，Reichard U，Treesucon U，et al.1998.Dispersal，pair formation and social structure in gibbons（Hylobates lar）J.Behavioral Ecology and Sociobiology，42（5）：329-339.Brockelman WY，Srikosamatara S.1993.Estimation of density of gibbon groups by use of loud songsJ.American Journal of Primatology，29（2）：93-108.Bryant JV，Brul A，Wong MHG，et al.2016.Detection of a new hainan gibbon（Nomascus hainanus）group using acoustic call playback J.International Journal of Primatology，37（4）：534-547.Chan BPL，Mak CF，Yang JH，et al.2017.Population，distribution，vocalization and conservation of the Gaoligong hoolock gibbon（Hoolock tianxing）in the Tengchong section of the Gaoligongshan National Nature Reserve，China J.Primate Conservation，31：107-113.Cheyne SM.2008.Effects of meteorology，astronomical variables，location and human disturbance on the singing apes：Hylobates albibarbis J.American Journal of Primatology，70（4）：386-392.Chivers DJ.1974.The siamang in Malaya.A field study of a primate in tropical rain forest J.Contributions to Primatology，4：1-335.Clink DJ，Ahmad HA，Klinck H.2020.Gibbons aren t singing in the rain：presence and amount of rainfall influences ape calling behavior in Sabah，Malaysia J/OL.Scientific Reports，10：1282 2022-06-10.https：/doi.org/10.1038/s41598-020-57976-x.Cowlishaw G.1992.Song function in gibbonsJ.Behaviour，121（1-2）：131-153.Fan PF，Ai HS，Fei HL，et al.2013.Seasonal variation of diet and time budget of eastern hoolock gibbons（Hoolock leuconedys）living in a northern montane forestJ.Primates，54（2）：137-146.Fan PF，He K，Chen X，et al.2017.Description of a new species of hoolock gibbon（Primates：Hylobatidae）based on integrative taxonomy J/OL.American Journal of Primatology，79（5）：e226312022-08-11.https：/doi.org/10.1002/ajp.22631.Fan PF，Turvey S，Bryant JV.2020.Hoolock tianxing EB/OL .（2019-07-10）2022-08-11 .https：/dx.doi.org/10.2305/IUCN.UK.2020-1.RLTS.T118355648A166597159.en.Fan PF，Xiao W，Huo S，et al.2009.Singing behavior and singing functions of black-crested gibbons（Nomascus concolor jingdongensis）at Mt.Wuliang，central Yunnan，ChinaJ.American Journal of Primatology，71（7）：539-547.Fan PF，Xiao W，Huo S，et al.2011.Distribution and conservation status of the vulnerable eastern hoolock gibbon Hoolock leuconedys in China J.Oryx，45（1）：129-134.Fan PF.2017.The past，present，and future of gibbons in China J.Biological Conservation，210b：29-39.Fang S，Mao K，Xia X，et al.2022.Dataset of daily near-surface air temperature in China from 1979 to 2018J.Earth System Science Data，14（3）：1413-1432.Fei HL，Thompson C，Fan PF.2019.Effects of cold weather on the sleeping behavior of skywalker hoolock gibbons（Hoolock tianxing）in seasonal montane forestJ/OL.American Journal of Primatology，81（9）：e23049 2022-08-11.https：/doi.org/10.1002/ajp.23049.Geissmann T.2002.Duet-splitting and the evolution of gibbon songs J.Biological Reviews，77（1）：57-76.Gittins SP.1984.The vocal repertoire and song of the agile gibbonM/Preuschoft H.The lesser apes：evolutionary and behavioural ecology.Edinburgh：Edinburgh University Press：354-375.Henwood K，Fabrick A.1979.A quantitative analysis of the dawn chorus：temporal selection for communicatory optimization J.The American Naturalist，114（2）：260-274.Hu NQ，Guan ZH，Huang B，et al.2018.Dispersal and female philopatry in a long-term，stable，polygynous gibbon population：evidence from 16 years field observation and geneticsJ/OL.American Journal of Primatology，80（9）：e22922 2022-08-11 .https：/doi.org/10.1002/ajp.22922.Liu GQ，Lu XL，Liu Z，et al.2022.The critically endan405四川动物 2023 年 第 42 卷 第 4 期Sichuan Journal of Zoology Vol.42 No.4 2023gered Hainan gibbon（Nomascus hainanus）population increases but not at the maximum possible rateJ.International Journal of Primatology，43：932-945.Ma HG，Ma CY，Fan PF.2022.Adult male-female social bond advertising：the primary function of coordinated singing intensity in a small ape J/OL.American Journal of Primatology，84（1）：e23351 2022-08-11.https：/doi.org/10.1002/ajp.23351.Ma HG，Ma CY，Fei HL，et al.2020.Cao vit gibbons（Nomascus nasutus）sing at higher elevation but not in peripheral areas of their home range in a karst forest J.International Journal of Primatology，41（5）：701-713.Mitani JC，Marler P.1989.A phonological analysis of male gibbon singing behavior J.Behaviour，109（1-2）：20-45.Mitani JC.1988.Male gibbon（Hylobates agilis）singing behavior：natural history，song variations and functionJ.Ethology，79（3）：177-194.R Core Team.2022.R：a language and environment for statistical computingEB/OL.（2022-06-23）2022-08-11.https：/www.R-project.org/.Raemaekers JJ，Raemaekers PM，Haimoff EH.1984.Loud calls of the gibbon（Hylobates lar）：repertoire，organisation and context J.Behaviour，91（1-3）：146-189.Rawson B.2004.Vocalisation patterns in the yellow-cheeked crested gibbon（Nomascus gabriellae）M/Ha T.Conservation of primates in Vietnam.Endangered Primate Rescue Center，Cuc Phuong National Park：Frankfurt Zoological Society，Vietnam Primate Conservation：130-136.Reichard UH，Neumann C.2007.Long-distance solo calls of white-handed gibbon