5·12地震后北川次生灾害迹地植被的自然恢复与更新

崔艳红，史常青，孙丽文，彭贤锋，张 艳，赵廷宁

（1.北京林业大学 水土保持学院，北京100083；2.乌鲁木齐市高新区林业园林管理局，新疆 乌鲁木齐830001；3.广东省水利水电科学研究院，广东 广州 510610）

5·12汶川8.0级地震发生在龙门山断裂带上［1］，致使在其范围内形成了大量的滚石、崩塌、滑坡、泥石流、堰塞湖等次生地质灾害［2］。仅北川羌族自治县境内就有滑坡、崩塌、泥石流等共计581处。北川县在汶川特大地震中受灾严重，其地表破坏强烈。多种次生地质灾害的发生，导致大量耕地、林地被淹没毁坏，山地生态环境极度退化，群落结构遭到严重破环，生物多样性急剧降低［3］。地震区域多处地表被大量松散碎石覆盖，受到暴雨冲刷极易加重水土流失，造成该区域植被恢复困难。因此，灾害迹地植被恢复工作需受到重视。群落结构是群落中植物与植物、植物与环境间相互关系的标志，从群落结构方面来研究物种多样性有一定生态学意义［4-5］。物种多样性指数能够反映群落种类组成、结构水平及植物群落的稳定性和复杂性，关系到群落成熟度、生产力、进化时间、捕食竞争与空间异质性。植物多样性的恢复是受损生态系统重建的重要内容，在植被恢复过程中其变化说明了植被的恢复进程［6-8］。植被多样性研究有利于预测植物群落的演替方向及群落稳定程度［9］。近年来，关于震后灾区植被恢复的研究较多［10-11］，大多偏重于立地条件、植被自然恢复的研究［12-14］，而地震引发的不同次生地质灾害对植物生物多样性破坏程度和植被恢复效果的研究较少［15-16］。本研究对5·12汶川地震引发的崩塌、滑坡、泥石流、落石、堰塞湖等5种次生灾害迹地的植被恢复情况进行调查，了解地震诱发的次生灾害迹地植被自然恢复群落组成、结构特征及生物多样性特征，对进一步开展灾区植被恢复具有理论指导意义。

1 研究地区与方法

1.1 研究区概况

北川羌族自治县位于四川省西北部（31°14′～32°14′N， 103°44′～104°42′E）， 隶属北亚热带山地湿润季风气候区，全年气候温和，雨水充沛，四季分明，大陆性季风气候明显。研究区年平均气温为15.6℃，多年平均降水量为1 400.0 mm，年平均无霜期为276.0 d，常年日照时数为924.3 h。土壤和植被随海拔呈带状分布，海拔1 800.0 m以下以砾石土为主。

1.2 试验样地布设

试验地选取5·12地震重灾区北川地区的崩塌、滑坡、泥石流、落石、堰塞湖等5种次生地质灾害点，在尽量靠近重灾区的附近选择未受损林地桤木Alnus cremastogyne-樟树Cinnamomum camphora混交林进行对比研究。崩塌迹地样地设置在坡脚，地表为小块碎石及少量大粒碎石；泥石流迹地样地设置在泥石流形成区与堆积区连接处，地表为松散碎石；滑坡迹地样地设置在滑坡体上，地表覆盖砂土及碎石；落石迹地样地设置在坡脚偏上的大块落石滚过处，土层较薄；堰塞湖淹没迹地样地设置在坡中，地表被砂石碎屑覆盖。桤木-樟木混交林样地设置在坡中，植被保存完好，乔木层在10 m以上，冠幅大于2 m，林龄大于10 a。样地基本情况见表1。

1.3 植被调查

2015年4月26日至5月20日，用全球定位系统（GPS）测得每块样地的三维坐标，用罗盘测得坡度和坡向。在样地中顺坡设置5 m×30 m的样带，灾害自然恢复地样带设置在土壤堆积物上，未受破坏林地样带设在坡体的中间位置，共选样地6个。在设置好的样带里布置6个5 m×5 m的样方，进行乔木调查。在5 m×5 m样方的对角线位置取2个2 m×2 m的小样方进行灌木调查，同时随机取2个1 m×1 m的小样方进行草本调查，共设大小样方180个。调查并记录每个样方内物种名称、数量、平均高度、盖度等指标。

表1 样地基本信息Table 1 Basic characteristics of plots

1.4 数据处理

物种的优势度由其重要值表示，重要值（VI）=（相对密度+相对频度+相对盖度）/3。本研究草本层使用相对高度代替相对密度，乔木层中的相对盖度主要为相对基盖度。

采用Simpson多样性指数、Shannon-Wiener多样性指数、Hurlbert均匀度指数和Margalef丰富度指数进行计算［17］。

Simpson多样性指数D：

Shannon-Wiener多样性指数H：

Hurlbert均匀度指数E：

式（1）～式（4）中：Pi为物种i的数量百分比，即Pi=Ni/N，其中Ni为样方中第i物种的个体数，N为全部种的个体总数；S为样方面积群落中植物种数。

Jaccard物种相似性指数I：

式（5）中：c表示A和B 2个样地群落的共有植物种数；a表示样地A群落的植物种数；b表示样地B群落的植物种数。

Margalef丰富度指数R：

2 结果与分析

2.1 不同灾害迹地群落组成及结构特征

研究发现，不同次生灾害迹地植被恢复较好，植被恢复初期以草本为主，有少量的灌木或乔木种。受损迹地群落结构简单，其植物相对比较单一，主要是菊科Compositae，禾本科Gramineae植物。未受损林地植物种类多样化，群落结构相对复杂，乔木有樟树，桤木，构树Broussonetia papyrifera，八角枫Alangium chinense，灌木有豪猪刺Berberis julianae，野蔷薇Rosa multiflora，小叶鼠李Rhamnus parvifolia，悬钩子Rubus corchorifolius和胡颓子Elaeagnus pungens等。

在次生灾害迹地中，落石迹地（B4）植物种最多，共有植物29科39属40种，草本层和灌木层最为发达，草本和灌木种类均大于ck群落，但缺少乔木层；而崩塌地（B1）则最少，共有植物10科20属20种，其中乔木为桤木；堰塞湖淹没地（B5）植物群落层次明显；泥石流迹地（B2）草本层最为发达，主要有五节芒Miscanthus floridulus，蒙古蒿Artemisia mongolica，千里光Senecio scandens和野棉花Anemone vitifolia等，灌木仅有水麻Debregeasia orientalis；B3和B4群落组成只有草本和灌木，缺少乔木层植物（表 2）。

表2 不同迹地群落物种组成和生长型Table 2 Species composition and growth form of different slash communities

2.2 不同灾害迹地植物群落优势种

重要值是衡量群落中植物是否为优势种的一项重要指标。对重要值的分析可找出群落中的主要优势树种［18］。从表3可见：①对乔木层来说，桤木-樟树混交林（ck）以桤木和樟树为优势种，2个物种的重要值相差不大，其次是八角枫。而次生灾害迹地，崩塌迹地 （B1）上仅有桤木其重要值为1；泥石流迹地（B2）有3种乔木，其重要值从大到小依次为桤木＞银杏＞铁木；滑坡体（B3）和落石迹地（B4）没有乔木；堰塞湖淹没地（B5）优势种为盐肤木、桤木、铁木。由此可以看出，桤木为该地区主要乔木优势种。②对灌木层来说，水麻为次生灾害迹地主要灌木优势种。B1样地灌木物种重要值位于前3位的依次为鸡骨柴、红麸杨、水麻；B2样地灌木只有2种水麻，其中水麻重要值为0.77，长叶水麻重要值为0.23；B3群落灌木物种重要值位于前3位的分别为水麻、木蓝、铁木（幼树）；B4和B5样地灌木优势种一样，都为马桑、水麻、披针叶荚蒾；未受损林地ck灌木优势种为刺花椒、豪猪刺、野蔷薇。③对草本层来说，次生灾害迹地上优势植物主要有千里光、五节芒，其次是飞蓬、蒙古蒿、荩草、葛藤、白苞蒿等；未受损林地（ck）重要值位于前3位的分别为山麦冬、蜈蚣蕨、地果。

表3 不同迹地群落物种重要值排序（前3位）Table 3 A sorting of species important value of different slash communities（Top 3）

2.3 不同灾害迹地群落物种相似性

植物种种间相关性是植物群落重要的数量和结构特征之一，是衡量2个物种相似性的一种尺度。由表4可知：次生灾害迹地物种与未受损林地物种存在较大差异，Jaccard物种相似性指数均在0.1左右，其中落石区（B4）物种与未被破坏地物种较为接近，相似性系数为0.19；崩塌地（B1）植物种与未被破坏地只有1种共有物种，相似性系数为0.02。次生灾害迹地与未受损林地物种相似性系数大小依次为B4＞B5＞B2＞B3＞B1。

相比而言，不同次生灾害迹地之间物种更为相似，其物种相似性指数为0.18～0.30。其中：以B2和B3样地间物种相似性系数最大，为0.30；B1与B2，B3，B4，B5物种相似性系数相对较小；而B2，B3，B4和B5样地之间物种相似性系数在0.28左右，相差不大，说明其物种较为相似。不同次生灾害迹地之间主要共有物种有五节芒，千里光，飞蓬，苦荬菜Ixeris polycephala，白苞蒿，荩草，野青茅Deyeuxia arundinacea，蕨Pteridium aquilinum和凤尾蕨Pteris cretica。由于未受损林地植被经过多年发育，植被群落达到一个相对稳定状态，而受损地区的林地是在原生演替初期，所处的发育阶段不同，因此群落物种相似性低。

表4 不同迹地群落物种相似性系数Table 4 Species similarity index between different slash communities

2.4 不同灾害迹地群落物种多样性特征

地震加剧了生境的破碎化程度，并引发次生灾害，造成了大片的植被破坏，减弱了生态功能。植物多样性在一定程度上能够反映植被恢复程度，是衡量植被演替进程的一项重要指标。对于地震后不同次生灾害多样性指数的计算和分析能够反映植物群落的生境差异、群落的结构类型、演替阶段和稳定性程度［19］。研究样地的植物分乔木、灌木和草本3种生长型，其多样性指数见表5。

乔木层Simpson指数，Shannon-Wiener指数和Hurlbert均匀度指数在B5样地最大，分别是0.65，1.07和0.98，而在其余受损地较小。由于崩塌、滑坡、泥石流、落石灾害对原本乔木破坏较大，多数乔木由于次生地质灾害的破坏倒塌死亡，乔木在短期难以恢复。由于混交林中乔木种类较多，Margalef丰富度指数却是未受损林地（ck）的最大，为0.82。

灌木层Simpson指数，Shannon-Wiener指数，Hurlbert均匀度指数和Margalef丰富度指数均在B4样地最大，其次是ck。因为B2样地仅有2种灌木，而B2样地多样性指数最小。次生灾害迹地群落Simpson指数，Shannon-Wiener指数小于ck，说明次生灾害迹地植被处于恢复初期，物种多样性必定小于长期稳定的未受损林地。然而均匀度与未受损林地差距不大，甚至B5均匀度大于ck均匀度。

草本层Simpson指数，Shannon-Wiener指数和Margalef丰富度指数均在B4样地最大，为0.89，2.52和4.00，Hurlbert均匀度指数B2样地（0.83）最大。林地受次生地质灾害破坏后，最先恢复的是草本层，因此，各样地草本多样性指数总体上差异不大。次生灾害迹地草本层多样性指数高于其他群落，而ck群落多样性指数和丰富度指数并不高，原因可能是稳定的乔木和灌木群落影响到林下草本层植物的多样性。

表5 不同迹地群落多样性指数Table 5 Diversity indices of different slash communities

3 结论与讨论

3.1 结论

5·12地震引起的次生灾害迹地植被破坏较为严重，目前仍处于恢复期，且植被恢复程度不同，其中草本层自然恢复与更新较快，其次为灌木层，乔木层恢复较慢。在演替早期，草本植物尤其是禾本科和菊科植物，可作为次生灾害迹地植被恢复的先锋植物。

从物种重要值来看，次生灾害迹地乔木层的优势种主要是桤木；灌木层优势种以水麻为主，其次是马桑和披针叶荚蒾；草本层优势种以五节芒、千里光等禾本科和菊科植物为主。未受损林地植物优势种更为多元化。

在群落物种相似性方面，次生灾害迹地物种与未受损林地物种存在较大差异，相似性指数较低，在0.1左右；而不同次生灾害迹地之间物种更为相似，其物种相似性指数为0.18～0.30，且其共有植物种有五节芒、千里光、飞蓬、苦荬菜、白苞蒿、荩草、野青茅、蕨和凤尾蕨。

在群落物种多样性方面，乔木层Simpson多样性指数，Shannon-Wiener多样性指数和Hurlbert均匀度指数在B5样地最大，分别是0.65，1.07和0.98；灌木层Simpson多样性指数，Shannon-Wiener多样性指数，Hurlbert均匀度指数和Margalef丰富度指数均在B4样地最大，分别为0.89，2.41，0.89和3.82，其次是ck；草本层也是B4样地最大，但不同次生灾害迹地之间差别较小。

3.2 讨论

5·12地震引起的次生灾害迹地土壤表层堆积着大量的松散碎石，土壤养分缺乏，土壤肥力较差，保水能力弱，水土流失严重，不利于植被恢复，这与大多数研究结果一致［3,13］。菊科、禾本科植物在困难立地早期植物群落里占显著优势，对受损后植被的恢复和延续有重要作用，是典型的先锋植物种［20］。据中国地震局公布的资料，仅2008年5-10月4级以上的余震，北川县就发生了22次［21］。在地震发生之后的3 a内，北川县因地表破碎、暴雨等影响，次生灾害时常发生［22］，又因为缺失植被自然恢复的优良条件，使得植被自然恢复进程缓慢。

在植物多样性指数方面，次生灾害迹地群落多样性指数并不都低于未被破坏的桤木-樟木混交林。在草本层除崩塌迹地外，其他受损林地多样性指数均大于未被破坏的桤木-樟木混交林，在灌木层只有落石迹地多样性指数大于未被破坏的桤木-樟木混交林，在乔木层只有堰塞湖迹地多样性指数大于未被破坏的桤木-樟木混交林。就植被演替进程来看，受损林地仍处在自然演替的初期即草本群落阶段，该阶段群落环境较差，以1～2年生植物为主，但随着时间的推移群落种类逐渐增多，多年生草本植物逐渐占优势，并发展为较为茂密的草本植物群落。演替总是朝着适应性强的方向演进，即按照草本—灌木—乔木的顺序变化着［23］。对于未受损林地而言，已经处于基本稳定的阶段，乔木的生长，冠层郁闭度的增加，使得林下光照减少，灌木和草本的生存受到一定的限制且趋于稳定［24］。

植被恢复是次生灾害迹地重建的一项重要工作，植被的自然恢复与演替又是一个长期的过程，所以对于研究区应进行连续长期的观测与调查，更加深入地研究植被自然恢复与演替的过程，为次生灾害迹地植被的恢复与重建提供科学依据。比如采用合适的自然恢复与人工恢复相结合的方法，加快次生灾害迹地的植被恢复工作；其次，应选取乡土植物种类，同时兼顾植物群落的多样性，尽量做到乔灌草搭配适宜，使受损立地的植被得到快速有效的恢复。

［1］ 王根龙,张军慧,刘红帅.汶川地震北川县城地质灾害调查与初步分析［J］.中国地质灾害与防治学报,2009,20（3）：47-51.WANG Genlong,ZHANG Junhui,LIU Hongshuai.Investigation and preliminary analysis of geologic disasters in Beichuan County induced by Wenchuan Earthquake ［J］.Chin J Geol Hazard Control,2009,20（3）：47-51.

［2］ 崔鹏,韦方强,陈晓清，等.汶川地震次生山地灾害及其减灾对策［J］.中国科学院院刊,2008,23（4）：317-323.CUI Peng,WEI Fangqiang,CHEN Xiaoqing,et al.Geo-hazards in Wenchuan Earthquake area and countermeasures for disaster reduction ［J］.Bull Chin Acad Sci,2008,23（4）：317-323.

［3］ 梁超,赵廷宁,史常青,等.基于NDVI的汶川大地震前后北川县次生地质灾害区植被破坏评估［J］.中国水土保持科学,2013,11（4）：86-92.LIANG Chao,ZHAO Tingning,SHI Changqing,et al.Vegetation damage assessment of Beichuan County after the earthquake in May 2008 based on NDVI ［J］.Sci Soil Water Conserv,2013,11（4）：86-92.

［4］ 张翠微.崇州市汶川地震灾后重建林草植被恢复工程建设监理实践与探讨［J］.四川林业科技,2013,34（5）：87-89.ZHANG Cuiwei.Practice and discussion on the construction of forest and grass vegetation restoration project after Wenchuan Earthquake in Chongzhou ［J］.J Sichuan For Sci Technol,2013,34（5）：87-89.

［5］ 常云妮,钟全林,程栋梁，等.闽西北地区不同林龄常绿阔叶混交林物种多样性比较［J］.生态环境学报,2013,22（6）：955-960.CHANG Yunni,ZHONG Quanlin,CHENG Dongliang,et al.Species diversity of ever-green broad-leaved mixed forest with different forest age in the north-west of Fujian Province ［J］.Ecol Environ Sci,2013,22（6）：955-960.

［6］ 彭少麟,方炜,任海，等.鼎湖山厚壳桂群落演替过程的组成和结构动态［J］.植物生态学报,1998,22（3）：245-249.PENG Shaolin,FANG Wei,REN Hai,et al.The dynamics on organization in the successional process of Dinghushan Cryptocarya community ［J］.Acta Phytoecol Sin,1998,22（3）：245-249.

［7］ 常学礼,邬建国.科尔沁沙地沙漠化过程中的物种多样性［J］.应用生态学报,1997,8（2）：151-156.CHANG Xueli,WU Jianguo.Species diversity during desertification on Kerqin Sandy Land ［J］.Chin J Appl Ecol,1997,8（2）：151-156.

［8］ 罗双,孙海龙,刘冲,等.四川道路边坡自然恢复的植被多样性研究［J］.水土保持研究,2011,18（6）：51-56.LUO Shuang,SUN Hailong,LIU Chong,et al.Study on vegetation diversity of road verge after natural restoration in Sichuan Province,China ［J］.Res Soil Water Conserv,2011,18（6）：51-56.

［9］ 朱晓勇,胡海波,鲁小珍，等.太湖西区公路两侧植物物种多样性的研究［J］.南京林业大学学报（自然科学版）,2006,30（3）：85-88.ZHU Xiaoyong,HU Haibo,LU Xiaozhen,et al.Study on the plant species diversity along the roadsides in Taihu West Zone ［J］.J Nanjing For Univ Nat Sci Ed,2006,30（3）：85-88.

［10］ ROTHAUS R M,REINHARDT E,NOLLER J.Regional considerations of coastline change,Tsunami damage and recovery along the southern coast of the Bay of Izmit（The Kocaeli（Turkey） Earthquake of 17 August 1999） ［J］.Nat Hazards,2004,31（1）：233-252.

［11］ CHOU W J,LIN W T,LIN C Y.Vegetation recovery patterns assessment at landslides caused by catastrophic earthquake：a case study in central Taiwan ［J］.Environ Monit Assess,2009,152（1/4）：245-257.

［12］ 田佳,田涛,赵廷宁,等.微立地因子植被恢复法在汶川地震植被重建中的应用［J］.中国水土保持科学,2008,6（5）：16-20.TIAN Jia,TIAN Tao,ZHAO Tingning,et al.Application of micro-site factors revegetation technology （MFRT） to revegetation of Wenchuan after earthquake ［J］.Sci Soil Water Conserv,2008,6（5）：16-20.

［13］ 刘锦春,何丙辉,徐小军,等.汶川草坡乡地震次生灾害迹地植物群落的恢复研究［J］.西南大学学报（自然科学版）,2013,35（4）：51-56.LIU Jinchun,HE Binghui,XU Xiaojun,et al.On vegetation restoration in Caopo of Wenchuan,a township destroyedby seismic secondary disaster ［J］.J Southwest Univ Nat Sci Ed,2013,35（4）：51-56.

［14］ 张翔,王庆安,王文国,等.汶川地震极重灾区植被恢复分区研究［J］.中国水土保持,2010,31（5）：48-51.ZHANG Xiang,WANG Qing’an,WANG Wenguo,et al.Study on vegetation recovery in the ponderosity disaster area of Wenchuan Earthquake ［J］.Soil Water Conserv Chin,2010,31（5）：48-51.

［15］ 林勇明,俞伟,刘奕,等.汶川震区典型区不同受损阶段群落多样性特征：以北川县为例［J］.北京林业大学学报,2012,34（6）：75-79.LIN Yongming,YU Wei,LIU Yi,et al.Community diversity features of different damaged stages in typical regions of Wenchuan Earthquake affected area：a case study in Beichuan,Sichuan Province of southwestern China ［J］.J Beijing For Univ,2012,34（6）：75-79.

［16］ 刘守江,张斌,杨清伟，等.汶川地震非规范滑坡体上植被的自然恢复能力研究：以彭州银厂沟谢家店子滑坡体为例［J］.山地学报,2010,28（3）：373-378.LIU Shoujiang,ZHANG Bin,YANG Qingwei,et al.Research on the natural recovery of vegetation on the non-normative landslide mass in Wenchuan Earthquake：take landslide mass in Xiejiadian of Yinchanggou in Pengzhou as an example ［J］.J Mount Sci,2010,28（3）：373-378.

［17］ 张金屯.数量生态学［M］.2版.北京：科学出版社,2011：372.

［18］ 郝建锋,王德艺,唐永彬，等.人为干扰对江油地区马尾松人工林群落结构和物种多样性的影响［J］.生态环境学报,2014,23（5）：729-735.HAO Jianfeng,WANG Deyi,TANG Yongbin,et al.Effects of human disturbance on species diversity of Pinus massoniana plantation in Jiangyou district,Sichuan Province ［J］.Ecol Environ Sci,2014,23（5）：729-735.

［19］ 秦伟,朱清科,刘中奇,等.黄土丘陵沟壑区退耕地植被自然演替系列及其植物物种多样性特征［J］.干旱区研究,2008,25（4）：507-513.QIN Wei,ZHU Qingke,LIU Zhongqi,et al.Study on natural series of vegetation and plant species diversity on returning land for farming to forests and grassplots in the hilly-gully regions of the Loess Plateau ［J］.Arid Zone Res,2008,25（4）：507-513.

［20］ 张广帅,邓浩俊,杜锟,等.汶川地震生态治理区土壤种子库及其与地上植被的关系［J］.中国生态农业学报,2015,23（1）：69-79.ZHANG Guangshuai,DENG Haojun,DU Kun,et al.Soil seed bank and its correlations with above ground vegetation in ecological restoration zones of Wenchuan Earthquake region ［J］.Chin J Eco-Agric,2015,23（1）：69-79.

［21］ 王青,李国蓉,梁斌,等.龙门山地震带5·12汶川地震余震空间分布特征［J］.地质科技情报,2009,28（2）：1-6.WANG Qing， LI Guorong,LIANG Bin,et al.Spatial distribution characters of aftershocks of 5·12 Wenchuan Earthquake in Longmenshan Seismic Belt［J］.Geol Sci Technol Inf,2009,28（2）：1-6.

［22］ 王娟.5·12地震北川县震后受损林地立地类型划分及其质量评价［D］.北京：北京林业大学,2012.WANG Juan.Site Type Division and Evaluation of Beichuan’s Damaged Forest after the 5·12 Earthquake in Wenchuan Country Town ［D］.Beijing： Beijing Forestry University,2012.

［23］ LI Liguang,HE Xingyuan,LI Xiuzhen,et al.Depth of edge influence of the agricultural-forest landscape boundary,Southwestern China ［J］.Ecol Res,2007,22（5）：774-783.

［24］ 王凯博.子午岭植被演替过程中物种多样性研究［D］.杨凌：西北农林科技大学,2008.WANG Kaibo.Study on the Species Diversity in Vegetation Succession in Ziwuling Area ［D］.Yangling：Northwest A＆F University,2008.