历史生物地理学在揭示地质历史变迁、解释生物多样性的演化历程方面起着重要的作用。生物演化以地质历史的变迁为客观条件，而地质历史的变迁又以生物分布规律及其演化进程为其佐证。随着板块的分裂、移动或连接，新的山脉随之隆起、海洋出现或消失、气候转变，生物物种分化、灭绝和扩散^[1]。我国地处亚洲大陆东南边缘，晚古生代以来泛大陆(Pangaea)的解体和重组，特别是印度板块的北移与青藏高原的急剧隆起，对我国现生生物的分布格局产生重大影响。Hsü等^[2]基于大陆内部的地缝合线将中国地壳划分为华北板块、塔里木板块、扬子板块、华夏板块、西伯利亚板块和冈瓦纳板块的一部分。受地壳运动的影响，至更新世，这六大板块与古西伯利亚和欧洲板块连为一体^[2-3]。受地壳运动的影响，中国的地形地势西高东低，以青藏高原最高，向东逐渐下降呈阶梯状分布。

中国复杂的地质演化历史(喜马拉雅造山运动和青藏高原隆起)、错综紊乱的地形(青藏高原，长江、黄河中下游平原等)及其历史演替过程和气候类型(从最北部的寒带至最南部的热带)的多样性不仅对现存生物的地理分布格局产生重大影响，而且也使得中国生物地理学的研究结果备受国际学术界的重视。张荣祖^[4]根据对现代陆栖脊椎动物的研究，将我国大陆动物区系划分为古北界和东洋界；古北界包括西北、东北和华北的大部分地区；东洋界包括以西南、华中区南部以及华南区为主的中国南部的省区。中国作为东亚地区的主体，是唯一横跨两大动物界的国家，其各生物类群的研究对于生物地理学的研究具有重要意义。近年来，依据两栖类、鸟类、陆生哺乳类和高等植物的分布进行生物地理学研究的报道越来越多^[5-7]。然而，基于昆虫类群进行分布格局的研究并不多见。

蝴蝶是自然生态系统的重要组成部分之一，与被子植物在长期的协同进化过程中形成了密切的关系。蝴蝶数量是衡量一个地区生物多样性丰富程度的重要指标。本文选用地理分布数据较为全面的弄蝶科Hesperiidae作为研究类群，采用特有性简约分析方法对其分布进行聚类分析，在此基础上对中国弄蝶的分布进行分区，推测弄蝶现有分布的成因，为中国的地质历史演化提供生物学方面的证据。

1.   材料与方法

1.1   研究区域

本文基于Meng等^[8]的方法，依据中国的地形、植被、弄蝶科昆虫的分布特点和已有的其他研究结果^[5]，对弄蝶的28个地理区域进行了稍微的调整，譬如：将四川分为川西高原(A15) 和四川盆地(A13)，将云南分为滇南山地(A08) 和云南高原(A27)，采纳冀晋山地(A07) 和秦巴山区(A06) 单独分区的观点，结果见图 1。

图  1  基于中国地形特点划分的28个弄蝶地理区

A01：淮北平原；A02：洞庭湖平原；A03：江汉平原；A04：浙闽丘陵北段；A05：西部荒漠；A06：秦巴山区；A07：河北平原；A08：滇南丘陵；A09：鄱阳湖平原；A10：内蒙古高原；A11：青藏高原东北部；A12：山西高原；A13：四川盆地；A14：长白山区；A15：川西高原；A16：贵州高原；A17：冀北辽西山地；A18：东北北部平原；A19：西北荒漠、高原；A20：山东丘陵；A21：西藏高原；A22：浙闽丘陵南段；A23：台湾地区；A24：广东沿海丘陵；A25：广西盆地；A26：海南地区；A27：西南山地；A28：陕晋高原。

Figure  1.  The 28 geographic areas of skippers in China based on topographical characteristics

下载: 全尺寸图片 幻灯片

1.2   方法

1.2.1   构建分布矩阵

查阅《中国动物志》及其他相关文献资料，结合华南农业大学鳞翅目昆虫研究室多年来的野外调查实物标本，统计中国已记录的弄蝶种类及分布，得到共计428种(包括亚种)的分布记录。将分布区有某物种存在编码为“1”，反之，编码为“0”，将不符合特有性简约分析原理的物种(分布区过于狭窄的物种)排除在外，最终用于分析的弄蝶为328种。以划定的28个区为内群，假设1个没有任何分布的A0为外群^[9-11]，构建1个29×328的矩阵。

1.2.2   特有性简约分析

特有性简约分析(Parsimony analysis of endemicity, PAE)是生物地理学研究中普遍采用的方法，用于生物地理格局及其特有性分析。该方法对划定区域的大小和形状没有特殊的要求，且区域的大小和外形不会影响其分析结果^[12]。本研究根据28个区域(图 1)中的物种数，采用特有性简约分析方法对其分布格局进行分析。利用PAUP4.0^[13]构建简约树(Maximum parsimony, MP)，进行自展检验^[14]，并以重复1 000次抽样的自展检验值(Bootstrap proportion, BP)来评估其分支的稳定性。

2.   结果与分析

构建基于弄蝶分布特点的MP树(图 2)，步长=1 135，一致性指数(Consistency index，CI)=0.290 7，存留指数(Retention index，RI)=0.448 3。中国弄蝶分布区聚成12支：B₁[(A01+A03+A04+A09)+(A02+A13)]，B₂[(A08+A26+A22+A24+A25)+(A16)]，B₃(A23)，C₁(A11)，C₂(A05+A28)，C₃(A10+A19)，D₁(A20)，D₂(A07+A12)，D₃(A06)，D₄[(A14+A18)+(A17)]，E₁(A15+A27) 和E₂(A21)；而且每一支对应一种地形。其中B₁包括2个亚支：“A01+A03+A04+A09”和“A02+A13”，并得到一定的支持(BP=34)，它们分别与中国动物地理区划中的东部丘陵平原亚区和西部山地高原亚区相对应。B₂由“A08+A26+A22+A24+A25”和“A16”2个亚支组成，分别对应于滇南山地亚区、海南岛亚区、闽广沿海亚区和贵州高原南部地区，且得到一定的支持(BP=39)。D₂由A07和A12聚成，以冀晋山区为主。D₄对应于东北三省，包括长白山山区和小兴安岭山区，且其支持值(BP=50) 高于D₂的(BP=47)。C₂包括A05和A28，以黄土高原为主。C₃由A10和A19组成，包括中国西北荒漠草原。E₁分为A15和A27亚支，依次对应于川西高原和云南高原，且支持值(BP=43) 低于C₂(BP=52) 和C₃(BP=55)。B₃(A23)、C₁(A11)、D₁(A20)、D₃(A06) 和E₂(A21) 分别独立成支，依次对应于台湾地区、青海地区、黄淮海平原地区、秦巴山区和西藏高原地区。

图  2  特有性简约性分析所得的严格合意支序图

图中的数字为1 000次自展置信检验值；A01：淮北平原；A02：洞庭湖平原；A03：江汉平原；A04：浙闽丘陵北段；A05：西部荒漠；A06：秦巴山区；A07：河北平原；A08：滇南丘陵；A09：鄱阳湖平原；A10：内蒙古高原；A11：青藏高原东北部；A12：山西高原；A13：四川盆地；A14：长白山区；A15：川西高原；A16：贵州高原；A17：冀北辽西山地；A18：东北北部平原；A19：西北荒漠、高原；A20：山东丘陵；A21：西藏高原；A22：浙闽丘陵南段；A23：台湾地区；A24：广东沿海丘陵；A25：广西盆地；A26：海南地区；A27：西南山地；A28：陕晋高原；B₁：中部地区；B₂：南部地区；B₃：台湾地区；C₂：中北部地区; C₃：西部地区，包括内蒙古和新疆；D₂：北部地区；D₃：秦巴山区；D₄：东北地区；E₁：西南地区；E₂：西藏地区。

Figure  2.  The strict consensus cladogram based on parsimony analysis of endemicity

下载: 全尺寸图片 幻灯片

根据简约性分析结果和弄蝶物种的分布特点，将中国版图划分成10个主要区(图 3)：中部地区(B₁)；南部地区(B₂)；台湾地区(B₃)；中北部地区(C₂); 西部地区(C₃)包括内蒙古和新疆；北部地区(D₂)；秦巴山区(D₃)；东北地区(D₄)；西南地区(E₁)；西藏地区(E₂)。

图  3  基于特有性简约分析的中国弄蝶分布格局

A01：淮北平原；A02：洞庭湖平原；A03：江汉平原；A04：浙闽丘陵北段；A05：西部荒漠；A06：秦巴山区；A07：河北平原；A08：滇南丘陵；A09：鄱阳湖平原；A10：内蒙古高原；A11：青藏高原东北部；A12：山西高原；A13：四川盆地；A14：长白山区；A15：川西高原；A16：贵州高原；A17：冀北辽西山地；A18：东北北部平原；A19：西北荒漠、高原；A20：山东丘陵；A21：西藏高原；A22：浙闽丘陵南段；A23：台湾地区；A24：广东沿海丘陵；A25：广西盆地；A26：海南地区；A27：西南山地；A28：陕晋高原；B₁：中部地区；B₂：南部地区；B₃：台湾地区；C₂：中北部地区; C₃：西部地区，包括内蒙古和新疆；D₂：北部地区；D₃：秦巴山区；D₄：东北地区；E₁：西南地区；E₂：西藏地区。

Figure  3.  Distribution patterns of skippers in China based on parsimony analysis of endemicity

下载: 全尺寸图片 幻灯片

中部地区(B₁)包括四川盆地、长江及其支流。气候属于亚热带暖湿类型，雨热同期。总体地貌特征属长江水系，河网密布，山地、丘陵、平原交错，地形复杂，弄蝶生物多样性较丰富(184种)，其中有1个特有种。南部地区(B₂)包括广东、广西和云南南部、福建东南沿海、海南、贵州，总的地势由西北向南、东南或西南倾斜，境内山地、丘陵、平原、河川纵横交错，地形复杂，气候类型众多，属南亚热带及热带，降水充沛，植被茂密，弄蝶生物多样性最为丰富，由于复杂的地质历史演化，该区特有种十分丰富，共计38种。台湾地区(B₃)因欧亚大陆板块、菲律宾海洋板块挤压而隆起，全岛地形东高西低，具有丰富的低地、低山及高山环境。高山环境为动物提供了合适的居所，有利于接纳不同水平地带的种类，故该区弄蝶生物物种相对丰富，共计61种。中北部地区(C₂)由黄土高原及河谷地带组成，受人为因素的影响，原有的天然森林和植被自然景观已被农田和零星分散的林地所替代，剩余的一点天然植被仅在黄土塬上的山地有分布，生物赖以生存的栖息环境的消失使得该区弄蝶种类少(54种)。西部地区(C₃)涵盖新疆和内蒙古，主要以草原和荒漠为主，气候属于半干燥和干燥型，东西部差异显著，西部受西来湿润气候的影响，形成草原。干旱的气候、荒漠和草原为主的植被条件使得弄蝶资源颇为贫乏(24种)，其中有1个特有种。北部地区(D₂)包括北京市、天津市、山西省、河北省及陕西东部的一部分，属于暖温带半湿润大陆性气候，四季分明，光照充足；冬季寒冷干燥且较长，夏季高温降水相对较多，春秋季较短，受人类活动的干扰较大，弄蝶种类和数量较少(42种)。秦巴山区(D₃)：秦岭、巴山是我国中西部的横断山系，自西向东，横亘甘肃南部、陕西南部、河南西部的广大地区，是世界动物地理区划东洋界与古北界的交汇处，属于亚热带和温带过渡地带，森林茂盛，弄蝶种类较多(97种)。东北地区(D₄)由大、小兴安岭，张广才岭，老爷岭，长白山以及松花江和辽河平原组成，属于温带湿润、半湿润大陆性季风气候。该区弄蝶少(34种)，无特有种。西南地区(E₁)包括四川西部、云南中北部，境内三江并流，南北走向的高山、峡谷纵横交错，气候复杂多变，干湿季明显。西南山地复杂的地形与有利的湿润条件使其生物多样性非常丰富，拥有大量的特有动植物种类，该区弄蝶种类十分丰富(180种)，占最终用于分析弄蝶种类的55%，含8个特有种。西藏地区(E₂)位于我国西南部，平均海拔4 000 m以上，属典型的大陆性气候，境内地理位置高低悬殊，生态环境极其复杂多样，由于地势高，各地平均气温都较低，且随着地形地势的变化气温差异较大，自东南至西北，从温、热带逐渐过渡到寒冷地区，本区记录了65种弄蝶。

3.   讨论与结论

Hsü等^[2]根据造山运动群岛模型将中国版图划分为六大地质省区：劳伦/华夏南部和西南边缘省区、西北省区、华中省区、华南省区、西藏省区和环太平洋省区。中国弄蝶的分布格局与这六大地质省区大致重叠。如：南部地区(B₂)与华南地质省区基本重合；中部地区(B₁)与华中地质省区大致相同；西部地区(C₃)(包括内蒙古和新疆)及华北(C₂和D₂)西南部区域与西北地质省区基本重叠；东北地区(D₄)与环太平洋北部省区基本一致；西藏地区(E₂)与西藏省区相吻合；西南地区(E₁)与劳伦/华夏南部和西南边缘省区基本重叠。

通过比对地质板块分布格局与弄蝶的分布格局可知：弄蝶主要分布在秦岭山脉的东部、藏东南地区、黄河和长江中下游以及南部地区。这些地区自更新世以来，伴随着青藏高原的隆起，经历了深刻的气候变化(东季风的形成，冰期和间冰期的频繁交替)以及河流系统的演变^[15-16]。最近的研究表明，在冰期期间，复杂的地形和山区可作为合适的潜在栖息地为温带物种提供相对稳定的微生态系统^[17]。中国中部、东部、南部于更新世冰川期间成为动植物的避难所，为生物提供适宜的栖息环境^[18-25]。避难所的存在可能对弄蝶的分布格局具有一定的作用，也可能是促进弄蝶种间和种内进一步分化的原因。

不同类群受同一地质历史的影响可能具有相同的分布格局。为检验这一假设，我们将本研究结果与其他类群进行了比较^{[4-5, 7]}，结果显示有5个区基本一致：中部地区(B₁)、南部地区(B₂)、西部地区(C₃)、北部地区(D₂)、东北地区(D₄)。同时，不同类群其分布格局亦有差异。张荣祖^[4]根据我国已有陆栖脊椎动物的研究，将云南南部山地作为滇南亚区处理，隶属于华南区。Xing等^[7]以动物区系内物种数为研究对象，对中国鸟类和哺乳类动物地理区划的19个亚区进行聚类分析，结果显示云南南部山地亚区与其他分支明显分开，独成一支，在此基础上，结合亚区的气候条件、植被类型和地质演化历程，将云南南部山地亚区提升为区，即滇南山区。本研究根据中国弄蝶的分布特点，利用特有性简约分析方法进行聚类，结果显示滇南山地亚区、闽广沿海亚区和海南亚区聚为一支，说明将滇南山地亚区放在华南区是非常合适的，故本研究支持张荣祖^[4]提出的视云南南部山地为1个亚区的观点。这种差异可能是因研究类群的不同所导致的，也可能与调查不全面有关。

在传统的中国动物地理区划中，台湾一直被视为1个亚区即台湾亚区，归属于华南区^[4]。在中国生物地理区划研究中，Xie等^[5]将台湾及其附近各小岛作为一个整体视为1个亚区。在对中国40 000种昆虫蜘蛛区系的多元相似性聚类分析研究中，申效诚等^[26]提出台湾单独为1个分布区。本研究结果显示台湾(A23) 单独聚为一支，显示了其物种分布格局的特殊性。台湾是我国东部最大的陆地性岛屿，因印度大陆北缘与亚洲大陆南缘挤压而隆起，其虽为岛屿性陆地，但有丰富的低地、低山及高山环境。小生境气候条件和高山环境为动物提供了合适的居所，因此也被赋予了“避难所”的称号^[27-29]。此外，台湾岛属于西太平洋岛弧系的一部分, 具备独特的地质和地球物理特征。因此，本研究认为台湾单独成为1个分布区比较合理，这一观点与申效诚等^[26]的研究结果一致。

西藏在我国动物地理区划中，存在较大分歧。陈宜瑜等^[30]从历史时空的角度对青藏高原地区鱼类区系的历史演化过程进行分析，认为高原地区在动物地理区划上应作为独立的区划单元，即为青藏区。王保海等^[31]根据西藏昆虫特有种和高山成分以及昆虫区系的垂直带谱，提出青藏高原区为单独的动物地理区，与古北区、东洋区等其他动物地理区并列。Xiang等^[6]依据中国及周边国家577种陆生哺乳动物的定量分析研究，提出青藏地区为1个亚区。本研究结果显示西藏(A21) 独立为一支，与其他分支明显分开，故本研究暂且将西藏区作为单独的1个区。由于对西藏地区的划分，仅考虑了物种分布上的特点，建议今后的研究充分考虑当地的自然气候条件、植被类型、不同海拔高度和地质历史的演变历程，结合已有的动植物区划体系多角度、多方位对西藏在动物地理区划中的划分作进一步的探讨。