1. 引言
鸣声(vocalization)由鸟类特有的发声器官鸣管(syrinx)发出 [1]，是鸟类重要的通讯手段，包含丰富的生物学信息，在求偶、保卫领域以及个体间交流中发挥着重要的作用 [2] [3]。Lwing Koch于1889年即开始对鸟类鸣声进行录制 [1]，随后鸟声受到了广泛的重视与研究 [4]。鸣声作为一类生物学特征，可以作为系统分类学的依据，探索物种的地理变异 [5] 、起源与进化历程 [6] [7] [8]。由于探索鸣声的地理变异，需要较大地理尺度的数据来源，因此大多数鸣声地理变异的研究都是针对雀形目鸟类开展，而针对其他鸟类的研究相对缺乏 [9] [10] [11] [12]。
啄木鸟(Picidae)分布于除大洋洲和南极洲以外的各个大陆 [13] [14]，是重要的森林健康指示物种和森林生态系统的伞护物种 [15] [16]，也是森林中鸟类多样性的重要指标 [17] [18]。啄木鸟属(Dendrocopos)中大斑啄木鸟(D. major)，白背啄木鸟(D. leucotos)和叙利亚啄木鸟(D. syriacus) 3种鸟类广泛分布于古北界(Palaearctic) [19]，种群数量大，栖息地生境类型多样，为鸣声地理变异的研究提供充分的数据来源。
随着科学技术的进步和科学研究的需要，广大研究者和鸟类爱好者记录了大量的鸣声数据并分类汇总成数据库，形成了鸟类鸣声大数据，为开展大地理尺度鸟声的比较研究奠定了基础 [20]。本文利用鸟类鸣声数据库，分析3种啄木鸟的鸣声差异，有助于探索鸟类鸣声地理变异的一般规律。
2. 材料与方法
2.1. 数据来源
我们通过xeno-canto.org网站鸣声数据库下载啄木鸟鸣声音频文件。xeno-canto由Bob和Willem-Pier于2005年主办，致力于分享世界各地鸟类各种鸣声数据，鸣声音频文件由研究学者或鸟类爱好者上传，鸣声数据按鸣声类型、录制地点、录制时间、录制者、海拔高度等进行了系统的整理 [21]。其中啄木鸟属3种啄木鸟的鸣声音频文件从1972年积累至今，鸣声数据充足，来源于古北界各个地区。
啄木鸟属的分类系统依据世界自然保护联盟濒危物种红色名录IUCN Red List [19]。为了比较属内的鸣声地理差异，同时收集黑啄木鸟属(Dryocopus)黑啄木鸟(D. martius)的鸣声数据作为系统聚类分析的外群。
下载鸣声文件类型为鸣叫(call)，音频文件类型为MP3格式。根据鸣声数据的数量和来源，结合研究目标，将鸣声文件按亚洲、东欧、西欧、北欧、南欧、中欧六个自然地理区域分类，其中亚洲地区鸣声文件主要包括中国、日本、伊朗、土耳其、哈萨克斯坦、巴基斯坦六个国家。亚洲与欧洲的地理分界线为乌拉尔山脉(Ural Mountains)、大高加索山脉(Caucasus Mountains)、乌拉尔河(Ural River)、博斯普鲁斯海峡(Bosporus Strait)、达达尼尔海峡(Dardanelles Strait) [22]。中欧地区和南欧地区的地理分隔为阿尔卑斯山脉(Alps)，西欧和北欧的地理分隔为大西洋(Atlantic Ocean)、挪威海(Bosporus strait)和巴伦支海(Barents Sea)，气候均为温带阔叶林气候。东欧与南欧距离相对最远，其中相隔波德平原(Cental European Plains)、黑海(Black Sea)、亚里德亚海(Adriatic Sea)、喀尔巴阡山脉(Carpathian Mountains) [22]。
2.2. 鸣声数据特征提取
2.2.1. 鸣声音频处理
使用Audacity软件将啄木鸟的原始MP3音频文件中的空白录音或与鸣声内容不相关的片段剪切删除，并进行降噪处理 [23]。
2.2.2. 鸣声参数提取
使用声谱分析软件Wavesurfer提取以下6种声学参数 [24] (图1)：
音节时长Syllable length (单位：s)，一个音节持续的时长；
音节间隔时长Syllable interval length (单位：s)，两个相邻音节之间持续的空白时长；
(a) 
(b) 
(c) 
(d)注：(a)：原始波形；(b)：语图；(c)：声强图；(d)：功率频谱图。
Figure 1. Example of vocalization map parameters of the great spotted woodpecker
图1. 大斑啄木鸟鸣声语图参数图例
最低频率Lowest frequency (单位：Hz)，功率频谱图中幅值突然开始下降处对应的频率；
主频率Principal frequency (单位：Hz)，功率频谱图中幅值最大的峰频率；
最高频率Highest frequency (单位：Hz)，功率频谱图中幅值突然开始上升处对应的频率；
频率范围Frequency range (单位：Hz)，最高频率与最低频率之差。
每组鸣声数据分别独立提取三次，取平均值，以控制操作误差。
2.3. 鸣声数据分析
2.3.1. 啄木鸟属鸣声特征差异分析
由于所有鸣声数据均不服从正态分布，采用非参数检验Kruskal-Wallis检验方法，对不同地区啄木鸟鸣声的音节时长、音节间隔时长、主频率、最高频率、最低频率以及频率范围共6个声学参数数据进行统计检验，分析古北界啄木鸟属3种啄木鸟的鸣声地理差异。
2.3.2. 啄木鸟属鸣声相似性分析
以6个声学参数为变量矩阵，计算不同地区鸣声间的欧几里得距离(Euclidean distance)，进行系统聚类(hierarchical cluster method)，分析啄木鸟属3种啄木鸟鸣声地理差异的相似性。为了分析啄木鸟属内鸣声是否具有相似性，另外选择黑啄木鸟鸣声作为外群，合并4种啄木鸟鸣声数据进行系统聚类分析。
全部数据分析在R4.0.5环境下完成 [25]。聚类分析使用dist函数和hclust函数。统计检验显著性水平为α = 0.05。
3. 结果
共收集得到啄木鸟属3种啄木鸟鸣声数据1050条，黑啄木鸟鸣声数据303条(表1)，从每个地区采取随机抽样的方法，选取20条鸣声样本进行后续分析。
Table 1. Effective song data about woodpeckers in different regions of palaearctic boundary
表1. 古北界不同地区啄木鸟的有效鸣声数据
注：EE: Eastern Europe; WE: Western Europe; CE: Central Europe; SE: Southern Europe; NE: Northern Europe; AS: Asia; NR：无该地区的鸣声音频数据。
3.1. 古北界啄木鸟属鸣声特征比较
Kruskal-Wallis检验结果显示：古北界大斑啄木鸟、白背啄木鸟、叙利亚啄木鸟的鸣声特征均存在着显著地理差异(图2~图4)。亚洲地区大斑啄木鸟鸣声相比于其他地区持续时间最长，鸣声间隔时间最长，即鸣声最为缓慢，且最高频率在古北界6个地区中最高(图2)。北欧地区大斑啄木鸟间隔时间最短, 即鸣声最为急促，最低频率相比于其他地区最小(图2)。南欧地区频率范围明显窄于其他地区(图2)。西欧地区大斑啄木鸟的主频率存在较多异常值，且主频率相较于其他地区主频率高，振幅高，即鸣声相比于其他地区较洪亮(图2)。
注：EE：东欧Eastern Europe；WE：西欧Western Europe；CE：中欧Central Europe；SE：南欧Southern Europe；NE：北欧Northern Europe；AS：亚洲Asia。统计检验使用Kruskal-Wallis方法。
Figure 2. Comparison of vocalization parameters of the great spotted woodpeckers in different Palaearctic regions
图2. 古北界不同地区大斑啄木鸟的鸣声参数比较
注：EE：东欧Eastern Europe；CE：中欧Central Europe；NE：北欧Northern Europe；AS：亚洲Asia。统计检验使用Kruskal-Wallis检验方法。
Figure 3. Comparison of vocalization parameters of the white-baked woodpeckers in different palaearctic regions
图3. 古北界不同地区白背啄木鸟的鸣声参数比较
注：EE：东欧Eastern Europe；CE：中欧Central Europe；NE：北欧Northern Europe AS：亚洲Asia。统计检验使用Kruskal-Wallis检验方法。
Figure 4. Comparison of vocalization parameters of the Syrian woodpeckers in different palaearctic regions
图4. 古北界不同地区叙利亚啄木鸟的鸣声参数比较
白背啄木鸟、叙利亚啄木鸟与大斑啄木鸟鸣声地理差异特征相似。亚洲地区白背啄木鸟(图3)、叙利亚啄木鸟(图4)鸣声音节时长持续时间最长，最高频率相比于其他地区高。北欧地区白背啄木鸟、叙利亚啄木鸟鸣声间隔时间最短，即鸣声较其他地区急促，北欧地区白背啄木鸟最低频率最小。南欧地区叙利亚啄木鸟频率范围明显窄于其他地区。
3.2. 啄木鸟属鸣声地理差异及其相似性
为了比较古北界啄木鸟属不同地区鸣声的相似性，我们分别对古北界不同地区啄木鸟属3种啄木鸟的6种声学参数进行综合系统聚类分析。结果显示：按照啄木鸟鸣声特征的地理差异分别将大斑啄木鸟
(a) (b) (c)注：EE：东欧Eastern Europe；NE：北欧Northern Europe；WE：西欧Western Europe；AS：亚洲Asia；CE：中欧Central Europe；SE：南欧Southern Europe。
Figure 5. Similarity of vocalization of three species of woodpeckers in different regions
图5. 啄木鸟属三种啄木鸟在不同地区鸣声的相似性
聚为5类(图5(a))，白背啄木鸟聚为3类(图5(b))，叙利亚啄木鸟聚为2类(图5(c))。并发现西欧地区和北欧地区的大斑啄木鸟鸣声特征最相似，与亚洲地区和中欧地区的鸣声差异最大。亚洲地区和中欧地区的叙利亚啄木鸟鸣声特征相似，与东欧和南欧地区的鸣声存在着明显地理差异，3种啄木鸟均是鸣声特征相似的欧洲地区先聚类，后再与亚洲地区聚类(图5)。
引入黑啄木鸟鸣声作为外群，与啄木鸟属鸣声进行系统聚类分析，结果显示：啄木鸟属3种啄木鸟聚为一类，与黑啄木鸟明显区分，说明啄木鸟属内鸣声具有强烈的相似性。且在啄木鸟属一支中，大斑啄木鸟和叙利亚啄木鸟鸣声特征相似，与白背啄木鸟聚成两类。这两个属中，啄木鸟鸣声地理差异有相同特征，黑啄木鸟的东欧地区和南欧地区鸣声特征相近，大斑啄木鸟和叙利亚啄木鸟的东欧地区和南欧地区鸣声特征也相近(图6)。
注：BW：黑啄木鸟(Dryocopus martius)；GW：大斑啄木鸟(Dendrocopos major)；WW：白背啄木鸟(Dendrocopos leucotos)；SW：叙利亚啄木鸟(Dendrocopos syriacus)；EE：东欧Eastern Europe；WE：西欧Western Europe；CE：中欧Central Europe；NE：北欧Northern Europe；AS：亚洲Asia；SE：南欧Southern Europe。
Figure 6. Biogeographic pedigree of three species of palaearctic woodpecker genus based on vocalization
图6. 基于鸣声的古北界啄木鸟属三种啄木鸟的生物地理谱系
4. 讨论
4.1. 鸣声数据库的利用
原始收集鸣声的方法需要研究者去采样地对每一个鸟声进行录制并筛选，这样的方法采样范围小，效率低，需要大量的人力和经费支持 [26]。而进入大数据时代后，我们通过鸣声数据库收集世界各地的鸣声数据，探索大地理尺度上鸣声的生物地理变异，可以最大程度上节约数据收集的经费支出和时间精力，并获得可靠且充分的数据。本研究我们通过xeno-canto.org鸣声网站收集到古北界啄木鸟鸣声数据共1353条(表1)，并将获得的数据根据鸣声数据的数量和研究目标分为6个研究区域。虽然亚洲为古北界面积最大的一个洲，但由于鸟类分布和地形地势等原因使鸟类研究者和鸟类爱好者较少，导致我们在xeno-canto.org网站所收集的啄木鸟鸣声数据较少，因此将亚洲统归为一个研究区域。欧洲地区鸣声数据相对较多，并分布于各个国家，因此将欧洲划分为东欧、西欧、北欧、南欧、中欧5个自然地理区域。随着科学研究的需要与发展，鸟类研究者和鸟类爱好者可以通过加大对亚洲各部分鸟声录制以增加其数据量，为后续研究提供更为充分的数据支撑。
4.2. 啄木鸟属鸣声变异
鸟类鸣声会随着地理分布的不同而发生变异 [27] [28]。Lovell和Lein (2013)对桤木纹霸鹟鸣声地理变化的研究中采用单因素分析和多因素分析方法，结果显示桤木纹霸鹟鸣声在不同的地区存在着显著差异 [29]。我们的研究显示出与其相似的结果，啄木鸟属3种啄木鸟均存在显著的鸣声地理差异(图2~图4)。3种啄木鸟的鸣声地理差异有相似之处，如亚洲地区的3种啄木鸟鸣声音节时长持续时间均最长，在北欧地区的鸣声间隔时间均最短(图2~图4)。但其也存在着相异之处，如亚洲地区白背啄木鸟的鸣声最低频率最高，而叙利亚啄木鸟的鸣声最低频率最低(图2~图4)；在中欧地区和东欧地区的白背啄木鸟鸣声特征相似，而大斑啄木鸟和叙利亚啄木鸟在这两地区的鸣声特征并不相似(图5)。这个结果看似矛盾实则不然。出现这样结果的可能原因是，尽管3种啄木鸟均为啄木鸟属，但啄木鸟不同种间存在遗传结构的差异，使与鸣声相关的结构或器官及其生理机制具有差异，导致鸣声存有差异。在鸟声研究的早期有相关研究证明，鸟类鸣声结构既有其稳定遗传的部分，又存在变异的特征 [3] [28]。如李金林(2008)对云南柳莺鸣声特点和地理差异的研究发现，在7个研究地区中，云南柳莺的最高频率﹑最低频率和持续时间没有显著差异；而音节的间隔时间和音素数存在着显著差异 [11]。
因为很少有在较大的地理尺度上研究鸟类鸣声的多样性，所以这种鸣声特征的变异和稳定性之间的平衡是一个有趣的进化谜题。而在这种尺度上，鸣声可能是一个非常重要的信息识别信号 [20]。在鸟类进化过程中，由于自然选择、栖息环境和社会适应等影响，鸟类在形态特征、种群的遗传结构等方面会发生相应的变化，而鸣声结构的变异往往与形态特征变化方向相一致，阻碍种群间的信息传递与基因交流，长时间的地理隔离与生殖隔离可能会导致物种分化，并对物种形成产生一定影响 [30]。影响鸣声变异的因素有多种，需要我们继续对鸣声地理变异进行研究与探讨。
4.3. 啄木鸟鸣声变异的影响因素
鸟类地理分布格局的形成与维持是其经历长期演化与生境适应的结果 [31]。鸟类为了在自然选择的压力下获得最大的生存适合度，其在遗传结构、行为方式与生理特征等各个方面进行不断的进化以适应所处的环境，从而在长期的演化过程中形成较好的适应策略和遗传结构 [32]。鸣声变异也是鸟类对其目前地理环境分布的适应。我们对啄木鸟属3种啄木鸟的鸣声地理差异及其相似性分析，几个欧洲地区的大斑啄木鸟和白背啄木鸟鸣声特征均较为相似，与亚洲地区的啄木鸟鸣声特征明显不同(图5)。产生这种结果的原因可能是亚欧交界地区多山脉，多海峡，如乌拉尔山、大高加索山、博斯普鲁斯海峡、达达尼尔海峡等，从而导致啄木鸟因地理隔离引发鸣声地理变异。乌拉尔山介于东欧平原和西西伯利亚平原之间平均海拔500 m~1200 m [22]。大高加索山脉位于黑海、亚速海和里海之间。地形以高山、高原为主。一般海拔在4 km以下，山体宽度约为200 km [22]。博斯普鲁斯海峡全场30.4 km，最宽处为3.6 km。达达尼尔海峡长约61 km，最宽处为6.4 km，与马尔马拉海和博斯普鲁斯海峡一起组成黑海海峡 [22]。山脉的隆起和海峡的阻隔导致啄木鸟栖息地片段化，从而形成基因流屏障，对物种的种群历史和遗传分化有重要的影响 [33]。大高加索山脉的高海拔和黑海海峡可能是导致啄木鸟属亚欧两洲鸣声地理变异的根本性原因。
西欧地区和北欧地区的大斑啄木鸟鸣声特征最为相似，可能因为两地区均是温带阔叶林气候，具有相似的气候环境，大斑啄木鸟为适应所生存的气候环境，鸣声结构发生了相似的适应性变化，使两地区的鸣声特征也具有相似性。南欧地区和东欧地区的大斑啄木鸟鸣声特征相似，叙利亚啄木鸟和黑啄木鸟也存在同样的鸣声地理差异相似性(图6)。可能的原因是东欧与南欧两地由波德平原、黑海、亚里德亚海、喀尔巴阡山脉相隔，在长期的地理隔离情况下，演化过程中两地区均保留了保守性较强的鸣声要素，作为该种啄木鸟的基础鸣声，并把这种鸣声结构稳定的遗传下来。该结果与2013年Lovell和Lein (2013)对桤木纹霸鹟鸣声地理变化的研究结果类似，这些区域的鸣声特征差异有遗传基础，一些距离较远地区的鸟类鸣声可能比邻近地区的鸟类鸣声更相似 [29]。
此外，将黑啄木鸟鸣声作为外群引入后，进行综合系统聚类分析时，表现出啄木鸟属内鸣声特征的一致性以及与黑啄木鸟的差异性。并且发现大斑啄木鸟和叙利亚啄木鸟鸣声特征更为相似，白背啄木鸟各地区鸣声单独聚为一类(图6)。产生此现象的原因可能与大斑啄木鸟和叙利亚啄木鸟亲缘关系更近，遗传结构相似，所以两种啄木鸟的形态特征与鸣声特征也相似。这与鸟类鸣声的地理变异与形态学特征的变异在鸟类系统学研究结果相一致 [34]。
纵观对于啄木鸟鸣声的生物地理变异的研究，仅有几篇是基于分子生物学的研究，基于鸣声的地理变异研究目前知之甚少 [12] [35]。因此关于啄木鸟地理变异研究有必要进一步加强探索。通过对啄木鸟鸣声地理变异一般规律的研究，并与形态学、分子生物学等综合分析，将更有利于对啄木鸟的生物地理变异进行探究与应用。
5. 结论
1) 本研究应用Kruskal-Wallis检验，成功地对不同地区啄木鸟属3种啄木鸟的鸣声地理差异进行分析，得出古北界大斑啄木鸟、白背啄木鸟、叙利亚啄木鸟各种种内的鸣声特征均存在显著的地理差异。同时采用系统聚类分析方法，比较分析古北界啄木鸟属不同地区鸣声的相似性，得出：亚洲地区的3种啄木鸟鸣声特征与欧洲地区的鸣声特征均存在着地理差异。再次证明，无论物种种内还是中间，啄木鸟鸣声均存在生物地理差异。
2) 通过引入黑啄木鸟鸣声作为外群，采用系统聚类分析方法与啄木鸟属鸣声特征比较，证明了啄木鸟属内鸣声具有较强的相似性。且大斑啄木鸟和叙利亚啄木鸟鸣声特征更为相似，白背啄木鸟各地区鸣声单独聚为一类，这与其遗传结构、地理分布和气候环境有着密不可分的关系。
致谢
感谢东北林业大学野生动物与自然保护地学院徐笑和李梦奇两位同学前期对鸣声音频文件收集工作中给予的帮助，特此致谢。
基金项目
中央高校基本科研业务费专项资金资助(2572021AW03)；国家自然科学基金面上项目(31970385)。
参考文献
NOTES
*第一作者。
#通讯作者。