贵州花溪大学城破碎化林地蝴蝶多样性与嵌套分布格局.pdf

四川动物 Sichuan Journal of Zoology 2023，42（5）：579-585贵州花溪大学城破碎化林地蝴蝶多样性与嵌套分布格局伍素，李义，杨开慧，李琳，曾维波，蔡明飞，刘宇，张凯，徐雨*（贵州师范大学生命科学学院，贵阳550025）摘要：为揭示城镇化进程中生境破碎化对蝴蝶多样性及分布格局的影响，本研究于20172022年每年的58月对贵州省花溪大学城26个残存自然林地中的蝴蝶多样性进行调查，共记录蝴蝶5科51属81种。不同斑块中的蝴蝶物种数为253种，平均每个斑块10.2种9.7种。线性回归分析表明：蝴蝶物种丰富度与斑块面积呈边缘显著正相关，符合岛屿生物地理学理论提出的“面积效应”；物种丰富度与距最大斑块距离呈显著负相关，表明存在隔离效应；进一步分析揭示不同斑块间的蝴蝶群落不存在嵌套格局或呈现反嵌套格局，具体形成机制有待进一步研究。基于研究结果，建议在城镇规划建设中应注重维持栖息地的完整性，对不同面积的林地斑块均加以保护。关键词：生境破碎化；蝴蝶；面积；隔离度；嵌套格局；城市化中图分类号：Q958.1 文献标志码：A 文章编号：1000-7083（2023）05-0579-07Butterfly Diversity and Nestedness on Fragmented Woodlots in the Huaxi University Town,GuizhouWU Su，LI Yi，YANG Kaihui，LI Lin，ZENG Weibo，CAI Mingfei，LIU Yu，ZHANG Kai，XU Yu*（School of Life Sciences，Guizhou Normal University，Guiyang 550025，China）Abstract:To examine the impact of habitat fragmentation on the diversity and distribution pattern of butterflies during the process of urbanization，the diversity of butterfly in 26 fragmented woodlots in the Huaxi University Town（Guizhou Province）was investigated.Through six-round surveys，81 species of butterflies belonging to 51 genera，5 families were recorded.The number of butterfly species in different woodlots ranged from 2 to 53，with an average of 10.29.7 species per woodlot.The result of linear regression analysis showed that the species richness of butterflies was positively correlated with the area of woodlots，and this was accordant to the“area effect”proposed by the island biogeography theory.On the other hand，the correlation of species richness with the distance to the largest woodlot was significantly negative，indicative of an isolation effect.Furthermore，butterfly assemblages among woodlots were not significantly nested or even anti-nested.Based on the results of the current study，it is recommended that habitat integrity should be maintained in urban planning and construction，and woodlots with different sizes should be protected.Keywords:habitat fragmentation；butterflies；area；isolation；neatedness；urbanization生境破碎化是指在人为活动或自然干扰作用下，大片连续分布的自然生境被其他非适宜生境DOI：10.11984/j.issn.1000-7083.20210270收稿日期：2021-08-02 接受日期：2023-05-16基金项目：国家自然科学基金委员会-贵州省人民政府喀斯特科学研究中心项目（U1812401）；国家自然科学基金项目（32270540）；贵州省科技计划项目（黔科合基础-ZK 2021 一般098）作者简介：伍素（1999），女，硕士研究生，主要从事昆虫群落多样性研究*通信作者Corresponding author，E-mail：579四川动物 2023 年 第 42 卷 第 5 期Sichuan Journal of Zoology Vol.42 No.5 2023分隔成许多面积较小斑块的过程（Fahrig，2003）。如水位上升导致部分陆地与大陆隔离而形成陆桥岛屿（Wang et al.，2010；Hu et al.，2011），森林砍伐后形成的片断化森林斑块（Hill et al.，2011），城镇化过程中道路、房屋修建对原生生境破坏和分割后形成的斑块（Hill et al.，2011；郑进凤等，2021）。目前，已有大量研究分析了生境破碎化对生物多样性的影响，认为生境破碎化是导致生物多样性减少的重要因素之一（武晶，刘志民，2014；吴倩倩等，2017；Miller-Rushing et al.，2019）。在众多研究中，基于岛屿生物地理学理论提出的面积效应、隔离效应是最基础的研究内容。面积效应、隔离效应从物种丰富度角度预测：随着生境被分割、缩小和隔离，其中的物种丰富度可能降低（Watling&Donnelly，2006）。在此基础上，一些学者开始从群落组成角度研究生境破碎化的影响，报道破碎化生境中的生物群落呈现嵌套格局，即小斑块中的物种是大斑块中物种集合的子集（Wang et al.，2010）。如浙江千岛湖陆桥岛屿中的大中型兽类、鸟类和两栖爬行类（张竟成等，2008；王熙等，2012；Chen et al.，2019），广西喀斯特破碎化森林斑块中的鸟类（Li et al.，2013），浙江杭州城市绿地中鸟类呈现明显的嵌套格局（Wang et al.，2013）。嵌套格局的产生可能与生境被分割、缩小和隔离过程中的物种选择性灭绝（Patterson&Atmar，1986）、生境嵌套（Blake，1991）、选择性迁入（Cook&Quinn，1995）等因素有关。但是，一些研究也指出，受竞争作用、环境过滤和生境异质性等因素的影响，破碎化生境中的一些动物群落不存在嵌套格局，甚至出现反嵌套格局（Chen et al.，2018；郑进凤等，2021；田璐嘉等，2022）。蝴蝶是一类体型较大、易观察和辨别的昆虫类群。它们虽然有一定的飞行能力，但幼虫以植物为食，对寄主植物较为专一，对生境恶化做出反应的速度比其寄主植物快（张立微，张红玉，2016）。在生境丧失和破碎化时，蝴蝶的生存和繁殖会受到较大的影响。因此，在生物多样性监测与评价中，蝴蝶常被作为主要的指示生物类群之一（Chen et al.，2021）。当前，一些学者以海洋性群岛（英国和爱尔兰大陆附近的岛屿、中国舟山群岛）为样地，分析了蝴蝶多样性的种面积关系和嵌套分布格局，揭示随着岛屿面积减小、隔离度增加，受选择性灭绝和选择性迁入等因素的影响，蝴蝶群落分布呈现嵌套格 局（Dennis et al.，2012；Xu et al.，2017）。但是，有关破碎化生境中的蝴蝶多样性与群落分布的研究案例还比较缺乏（MacDonald et al.，2018）。本研究以贵州省花溪大学城破碎化自然林地为研究样地，对蝴蝶多样性进行调查，研究斑块面积、隔离度与物种丰富度的关系，并分析群落是否存在嵌套格局，以此丰富生境破碎化研究理论，并为城市规划建设和蝴蝶多样性保护提供科学依据。1研究方法1.1研究区概况研究区位于贵州省贵安新区东南部，紧邻贵阳市花溪区（1063610640E，26222624N），面积约10 km2。研究区属于典型的喀斯特地貌特征，海拔 1 0001 200 m，属亚热带湿润温和型气候，年均温15，年均降水量1 130 mm。该区域在规划修建大学城以前，为农田环绕的自然林地。2009年大学城动工修建至投入使用，自然林地被人工设施、建筑物、道路等不透水面分割。残留林地斑块中的植被类型相对简单，主要由暖性松林、常绿落叶阔叶林、落叶阔叶灌丛和禾草灌草丛组成，受喀斯特地貌特征影响，不同斑块间的小生境尤其是植物种类存在较大差异（曹朝阳，2021）。本研究选择其中 26 个斑块作为研究样地（面积0.3290.4 hm2；表1）。1.2斑块地理参数的测量从Google Earth下载研究区的卫星地图，导入ArcGIS 10.2中勾绘出各个斑块的边界，采用“计算几何”的方式计算斑块面积，以“生成近邻表”的方580伍素等：贵州花溪大学城破碎化林地蝴蝶多样性与嵌套分布格局式计算各斑块到最大斑块的距离，以反映斑块的生境破碎化程度。1.3蝴蝶采集及鉴定参考Schoereder等（2004）的研究，在不同斑块中设置长度与面积成正比的固定样线。于20172022年的58月分别对每个斑块进行了6次蝴蝶调查与样品采集，每条样线每次调查1次。选择天气晴朗、蝴蝶活跃度较高的时间段（10 0017 00）采用扫网法捕捉蝴蝶。每条样线 2名调查人员，调查行进速度为10 m/min；用捕虫网捕捉样线两侧5 m范围内的蝴蝶，将采集到的个体置于硫酸纸三角袋中，记录采集日期、地点和采集人等相关信息（张雪梅等，2016；Xu et al.，2017）。物种鉴定参考 中国蝴蝶分类与鉴定（周尧，1998）、中国昆虫生态大图鉴（张巍巍，李元胜，2011）、贵州蝴蝶（陈正军，2016）和 中国蝴蝶图鉴（武春生，徐堉峰，2017）。对形态鉴定存疑的种类，利用DNA条形码进行分类（Hebert et al.，2003）。提取蝴蝶腿部总DNA，扩增线粒体CO基因的部分片段（658 bp）测定序列，后与 BOLD 数据库（http：/www.boldsystems.org/）进行比对和鉴定。1.4统计与分析（1）物种丰富度：分别统计26个斑块中的蝴蝶物种数。（2）计算 Berger-Parker优势度指数（d）（Berger&Parker，1970）：d=ni/N，式中，ni为物种i的个体数量，N 为全部物种的总个体数量。d0.05 为优势种，0.005d0.05为常见种，d0.005为少见种或偶见种。（3）建立线性回归模型：按照岛屿生物地理学理论提出的种面积函数公式，对物种丰富度、斑块面积、距最大斑块距离采取10为底数的对数转换（由于最大斑块到自身的距离为0，因此在对数转换前对各个距离值进行+1处理）。利用R.4.2.1构建26个斑块中蝴蝶物种丰富度与斑块面积、距最大林地斑块距离的线性回归模型，分析是否具有面积效应和隔离效应。使用visreg程序包绘制模型的偏残差图（Breheny&Burchett，2020）。（4）分析群落嵌套格局：构建斑块与蝴蝶物种多度矩阵，分析群落的嵌套格局（Almeida-Neto&Ulrich，2011）。使用NODF2.0中的WNODF方法计算物种（WNODFr）、斑块（WNODFc）和整体的嵌套值（Almeida-Neto&Ulrich，2011）。设置行列固定（rc）、总数固定（aa）和总物种丰富度固定（ss）等3种重采样零模型，依据软件随机产生1 000个矩表1贵州花溪大学城26个破碎化林地斑块的地理特征、样线数量和长度Table 1Characteristics of 26 fragmented woodlots and number and length of transect lines recordedin the Huaxi University Town，Guizhou林地斑块编号ID of the fragmented woodlots1234567891011121314151617181920212223242526面积Area/hm22.46.29.43.62.10.90.68.62.63.70.92.31.51.53.21.2290.425.56.519.034.21.41.11.411.50.3距最大林地斑块距离Distance tothe largestwoodlot/m962.7826.5506.61 290.9974.51 221.5688.3828.71 223.5865.3663.7438.81 662.9244.12 188.191.30504.01 666.92 114.12 390.62 561.12 772.32 553.256.6271.4样线数量Number oftransectlines11111111111111114111111111样线长度Length oftransectlines/m442.6915.5926.3336.3298.6172.1126.2617.1282.3392.8153.6293.2250.9254.1252.6159.17 073.71 302.8558.6892.21 064.4183.2204.5244.3585.482.5581四川动物 2023 年 第 42 卷 第 5 期Sichuan Journal of Zoology Vol.42 No.5 2023阵对嵌套度进行显著性检验。上述统计分析中，P0.05表示差异显著。2结果共采集蝴蝶个体 586只，属于 5科 51属 81种（https：/ Nymphalidae物种数最多（32种），其次是灰蝶科 Lycaenidae（18 种）、凤蝶科 Papilionidae（13种）、粉蝶科Pieridae（12种）和弄蝶科Hesperiidae（6 种）。优势种有 4 种，分别是东方菜粉蝶Pieris canidia、东亚豆粉蝶Colias poliographus、琉璃灰蝶Celastrina sp.和酢浆灰蝶Zizeeria maha。常见种44种，少见种或偶见种33种。不同林地斑块中的蝴蝶物种数为253，平均为 10.2种9.7种。线性回归分析表明，蝴蝶物种丰富度与林地斑块面积呈边缘显著正相关（=0.149，P=0.057）；蝴蝶物种丰富度与距最大林地斑块距离呈显著负相关（=0.201，P=0.007；图1）。3种重采样零模型的WNODF嵌套分析显示，与随机产生的 1 000 个分布矩阵相比，rc 模型和aa模型表明没有显著嵌套格局，ss模型表明存在反嵌套格局（表2）。3讨论本研究在贵州花溪大学城26个破碎化自然林地中记录到蝴蝶81种，数据分析表明蝴蝶物种丰富度与林地斑块的面积呈边缘显著正相关，符合岛屿生物地理学理论提出的“面积效应”（McArthur&Wilson，1967）。由于土地大面积开发，贵州省花溪大学城的原生植被正不断被房屋、道路等人工建设物分割。斑块面积在一定程度上反映了可供物种利用的空间分布范围与资源量；小斑块中的生境丰富度和资源量可能较低，难以满足一些物种的生存需求，因此蝴蝶物种丰富度较低（张雪梅等，2016；Xu et al.，2017；MacDonald et al.，2018）。蝴蝶物种丰富度与距最大斑块距离呈显著负相关，表明存在隔离效应。虽然本研究区内各林地斑块的隔离度较低，但蝴蝶的飞行或扩散图1贵州花溪大学城26个破碎化林地中蝴蝶物种丰富度与林地斑块面积和距最大斑块距离的线性回归偏残差Fig.1Partial residual for linear regressions of between the species richness and the woodlot area andthe distance to the largest woodlot in the 26 fragmented woodlots in the Huaxi University Town，Guizhou黑线表示拟合值，灰色区域表示95%置信区间The black lines are model fits，the grey bands represent 95%confidence intervals582伍素等：贵州花溪大学城破碎化林地蝴蝶多样性与嵌套分布格局能力较弱，且人流量大和公路行驶车辆多等一定程度上增强了隔离效应。不同林地斑块间蝴蝶群落未呈现出嵌套结构，甚至可能出现反嵌套格局。这一结果与一些海洋群岛中蝴蝶群落的研究结果不同（Dennis et al.，2012；张雪梅等，2016；Xu et al.，2017），也不同于一些破碎化城市生境中其他动物群落的研究结果（Gonzlez-Oreja et al.，2012；Wang et al.，2013；de la Hera，2019）。一般认为，嵌套格局的产生与选择性灭绝（Patterson&Atmar，1986）、生境嵌套（Blake，1991）和选择性迁入（Cook&Quinn，1995）等因素有关。选择性灭绝理论认为，物种分布与斑块面积密切相关，一些种群较小、最小需求面积较大的物种通常具有较高的灭绝风险，会率先从小斑块中消失。选择性迁入理论认为物种的分布与其扩散能力有关，一些扩散能力较弱的物种可能会有序从隔离度较大的斑块中消失。生境嵌套理论认为，物种的空间分布与生境结构密切相关，斑块间的生境呈嵌套分布将导致物种的嵌套分布。虽然本研究揭示蝴蝶物种丰富度随林地斑块面积减小、隔离度增加而降低，但喀斯特特殊的地质地貌造成的局部小生境和植物种类差异（曹朝阳，2021；郑进凤等，2021；朱芸等，2022），以及蝴蝶对寄主植物的较强选择性（Krmer et al.，2012），可能使物种出现在各自适应的特定生境中，从而削弱了嵌套格局，甚至可能导致反嵌套格局（Gutirrez-Cnovas et al.，2013；郑进凤等，2021）。综上所述，本研究以贵州省花溪大学城为例，揭示了生境破碎化影响蝴蝶物种丰富度，即斑块面积越小，蝴蝶物种丰富度越低。但是，不同林地斑块间的蝴蝶群落没有由此产生显著的嵌套格局，甚至可能呈现反嵌套格局，这可能与喀斯特生境异质性造成的斑块间物种组成差异较大等因素有关，具体机制有待进一步研究。基于本研究结果，建议在城市规划建设时应注重维持栖息地的完整性，对不同面积的林地斑块均加以保护。致谢：感谢申力欣、张秋秋、高丽桦、张健、殷枭、梁晓珊、罗叶、李洁、朱芸在野外调查和蝴蝶形态鉴定工作中付出的辛劳，感谢汪小芳、赵娟在DNA条形码技术中所付出的努力。参考文献:曹朝阳.2021.贵州花溪大学城破碎化自然林地中木本植物多样性及空间分布格局研究 D.贵阳：贵州师范大学.陈正军.2016.贵州蝴蝶 M.贵阳：贵阳科技出版社.田璐嘉，杨小波，李东海，等.2022.海口和三亚两城市破碎化林地中鸟类群落多样性与嵌套分布格局 J.生物多样性，30（6）：92-102.王熙，王彦平，丁平.2012.千岛湖两栖爬行类动物群落结构嵌套分析 J.动物学研究，33（5）：439-446.吴倩倩，梁宗锁，刘佳佳，等.2017.中国生境片段化对生物多样性影响研究进展 J.生态学杂志，36（9）：2605-2614.武春生，徐堉峰.2017.中国蝴蝶图鉴 M.福州：海峡书局出版社.武晶，刘志民.2014.生境破碎化对生物多样性的影响研究综述 J.生态学杂志，33（7）：1946-1952.张竟成，王彦平，蒋萍萍，等.2008.千岛湖雀形目鸟类群落嵌套结构分析 J.生物多样性，16（4）：321-331.表2贵州花溪大学城26个破碎化林地中蝴蝶物种多度矩阵的嵌套格局分析Table 2Nestedness analyses for species-by-sitesabundance matrix of butterflies recordedin 26 fragmented woodlotsin the Huaxi University Town，Guizhou嵌套度量化Nestedness metric零模型：rc零模型：aa零模型：ss整体斑块物种整体斑块物种整体斑块物种观测值WNODFobs18.0118.8417.9318.0118.8417.9318.0118.8417.93期望值WNODFexp16.6518.6416.4516.1721.2515.6737.6438.4437.57SD1.661.561.751.722.301.835.044.745.10Z0.820.130.841.071.041.233.904.143.86P0.2050.4490.2000.1430.1480.1080.0010.0010.001583四川动物 2023 年 第 42 卷 第 5 期Sichuan Journal of Zoology Vol.42 No.5 2023张立微，张红玉.2016.蝶类对生境的指示作用研究进展J.生物学杂志，33（3）：88-91.张巍巍，李元胜.2011.中国昆虫生态大图鉴 M.重庆：重庆大学出版社.张雪梅，韩徐芳，刘立伟，等.2016.舟山群岛蝶类群落嵌套分布格局及其影响因素 J.生物多样性，24（3）：321-331.郑进凤，唐蓉，贺霜，等.2021.贵州花溪大学城破碎化林地鸟类多样性与嵌套分布格局 J.生物多样性，29（5）：661-667.周尧.1998.中国蝴蝶分类与鉴定 M.郑州：河南科学技术出版社.朱芸，盛尚，郑进凤，等.2022.贵州花溪大学城破碎化林地中鸟类群落的小岛屿效应 J.动物学杂志，57（2）：205-212.Almeida-Neto M，Ulrich W.2011.A straightforward computational approach for measuring nestedness using quantitative matricesJ.Environmental Modelling&Software，26（2）：173-178.Berger WH，Parker FL.1970.Diversity of planktonic foraminifera in deep-sea sediments J.Science，168（3937）：1345-1347.Blake JG.1991.Nested subsets and the distribution of birds on isolated woodlotsJ.Conservation Biology，5（1）：58-66.Breheny P，Burchett W.2020.Package visreg version 2.7.0 EB/OL 2020-06-04 .https：/cran.r-project.org/web/packages/visreg/.Chen CW，Holyoak M，Si XF，et al.2018.Do seasonal species assemblages differ in their biogeography？Evidence from the spatial structure of bird communities on landbridge islandsJ.Journal of Biogeography，45（2）：473-483.Chen CW，Xu AC，Ding P，et al.2019.The small-island effect and nestedness in assemblages of medium-and large-bodied mammals on Chinese reservoir land-bridge islandsJ.Basic and Applied Ecology，38：47-57.Chen CW，Xu AC，Wang YP.2021.Area threshold and traitenvironment associations of butterfly assemblages in the Zhoushan Archipelago，China J.Journal of Biogeography，48（4）：785-797.Cook RR，Quinn JF.1995.The influence of colonization in nested species subsets J.Oecologia，102（4）：413-424.de la Hera I.2019.Seasonality affects avian species distribution but not diversity and nestedness patterns in the urban parks of VitoriaGasteiz（Spain）J.Animal Biodiversity and Conservation，42（2）：279-291.Dennis RLH，Hardy PB，Dapporto L.2012.Nestedness in island faunas：novel insights into island biogeography through butterfly community profiles of colonization ability and migration capacityJ.Journal of Biogeography，39（8）：1412-1426.Fahrig L.2003.Effects of habitat fragmentation on biodiversity J.Annual Review of Ecology，Evolution，and Systematics，34（1）：487-515.Gonzlez-Oreja J，Fuente-Daz-Ordaz AD，Hernndez-Santn L，et al.2012.Can human disturbance promote nestedness？Songbirds and noise in urban parks as a case study J.Landscape&Urban Planning，104（1）：9-18.Gutirrez-Cnovas C，Milln A，Velasco J，et al.2013.Contrasting effects of natural and anthropogenic stressors on beta diversity in river organisms J.Global Ecology&Biogeography，22（7）：796-805.Hebert PDN，Cywinska A，Ball SL，et al.2003.Biological identifications through DNA barcodesJ.Proceedings of the Royal Society B：Biological Sciences，270（1512）：313-321.Hill JK，Gray MA，Khen CV，et al.2011.Ecological impacts of tropical forest fragmentation：how consistent are patterns in species richness and nestedness？J.Philosophical Transactions of the Royal Society of London，Series B，Biological Sciences，366（1582）：3265-3276.Hu G，Feeley KJ，Wu JG，et al.2011.Determinants of plant species richness and patterns of nestedness in fragmented landscapes：evidence from land-bridge islandsJ.Landscape Ecology，26（10）：1405-1417.Krmer B，Poniatowski D，Fartmann T.2012.Effects of landscape and habitat quality on butterfly communities in pre-584伍素等：贵州花溪大学城破碎化林地蝴蝶多样性与嵌套分布格局alpine calcareous grasslandsJ.Biological Conservation，152：253-261.Li ZT，Lu Z，Shu XL，et al.2013.Nestedness of bird assemblages in the karst forest fragments of southwestern Guangxi，China J.Chinese Birds，4（2）：170-183.MacDonald ZG，Anderson ID，Acorn JH，et al.2018.Decoupling habitat fragmentation from habitat loss：butterfly species mobility obscures fragmentation effects in a naturally fragmented landscape of lake islands J .Oecologia，168（1）：11-27.McArthur R，Wilson E.1967.The theory of island biogeography M.Princeton，NJ：The Princeton University Press.Miller-Rushing AJ，Primack RB，Devictor V，et al.2019.How does habitat fragmentation affect biodiversity？A controversial question at the core of conservation biologyJ.Biological Conservation，232：271-273.Patterson B，Atmar W.1986.Nested subsets and the structure of insular mammalian faunas and archipelagos J.Biological Journal of the Linnean Society，28（1-2）：65-82.Schoereder JH，Galbiati C，Ribas CR，et al.2004.Should we use proportional sampling for speciesarea studies？J.Journal of Biogeography，31（8）：1219-1226.Wang YP，Bao YX，Yu MJ，et al.2010.Nestedness for different reasons：the distributions of birds，lizards and small mammals on islands of an inundated lake J.Diversity and Distributions，16（5）：862-873.Wang YP，Ding P，Chen SH，et al.2013.Nestedness of bird assemblages on urban woodlots：implications for conservation J.Landscape and Urban Planning，111：59-67.Watling JI，Donnelly MA.2006.Fragments as islands：a synthesis of faunal responses to habitat patchiness J.Conservation Biology，20（4）：1016-1025.Xu AC，Han XF，Zhang XM，et al.2017.Nestedness of butterfly assemblages in the Zhoushan Archipelago，China：area effects，life-history traits and conservation implicationsJ.Biodiversity and Conservation，26（6）：1375-1392.585