-
5·12汶川8.0级地震发生在龙门山断裂带上[1]，致使在其范围内形成了大量的滚石、崩塌、滑坡、泥石流、堰塞湖等次生地质灾害[2]。仅北川羌族自治县境内就有滑坡、崩塌、泥石流等共计581处。北川县在汶川特大地震中受灾严重，其地表破坏强烈。多种次生地质灾害的发生，导致大量耕地、林地被淹没毁坏，山地生态环境极度退化，群落结构遭到严重破环，生物多样性急剧降低[3]。地震区域多处地表被大量松散碎石覆盖，受到暴雨冲刷极易加重水土流失，造成该区域植被恢复困难。因此，灾害迹地植被恢复工作需受到重视。群落结构是群落中植物与植物、植物与环境间相互关系的标志，从群落结构方面来研究物种多样性有一定生态学意义[4-5]。物种多样性指数能够反映群落种类组成、结构水平及植物群落的稳定性和复杂性，关系到群落成熟度、生产力、进化时间、捕食竞争与空间异质性。植物多样性的恢复是受损生态系统重建的重要内容，在植被恢复过程中其变化说明了植被的恢复进程[6-8]。植被多样性研究有利于预测植物群落的演替方向及群落稳定程度[9]。近年来，关于震后灾区植被恢复的研究较多[10-11]，大多偏重于立地条件、植被自然恢复的研究[12-14]，而地震引发的不同次生地质灾害对植物生物多样性破坏程度和植被恢复效果的研究较少[15-16]。本研究对5·12汶川地震引发的崩塌、滑坡、泥石流、落石、堰塞湖等5种次生灾害迹地的植被恢复情况进行调查，了解地震诱发的次生灾害迹地植被自然恢复群落组成、结构特征及生物多样性特征，对进一步开展灾区植被恢复具有理论指导意义。
HTML
-
北川羌族自治县位于四川省西北部（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以下以砾石土为主。
-
试验地选取5·12地震重灾区北川地区的崩塌、滑坡、泥石流、落石、堰塞湖等5种次生地质灾害点，在尽量靠近重灾区的附近选择未受损林地桤木Alnus cremastogyne-樟树Cinnamomum camphora混交林进行对比研究。崩塌迹地样地设置在坡脚，地表为小块碎石及少量大粒碎石；泥石流迹地样地设置在泥石流形成区与堆积区连接处，地表为松散碎石；滑坡迹地样地设置在滑坡体上，地表覆盖砂土及碎石；落石迹地样地设置在坡脚偏上的大块落石滚过处，土层较薄；堰塞湖淹没迹地样地设置在坡中，地表被砂石碎屑覆盖。桤木-樟木混交林样地设置在坡中，植被保存完好，乔木层在10 m以上，冠幅大于2 m，林龄大于10 a。样地基本情况见表 1。
样地类型 纬度 经度 坡度/(°) 坡向 海拔/m 总覆盖度/6 崩塌自然恢复地(B1) 31°51"37.7"N 104°30"18.0"E 36 南偏东 624 50 泥石流自然恢复地(B2) 31°52"17.1"N 104°31"53.1"E 24 东偏南 735 80 滑坡自然恢复地(B3) 31°49"25.1"N 104°27"13.2"E 28 东 690 68 落石自然恢复地(B4) 31°49"01.0"N 104°27"07.4"E 22 西偏北 736 75 堰塞湖自然恢复地(B5) 31°50"36.8"N 104°25"31.4"E 21 西偏南 739 83 桤木-樟树混交林(ck) 31°36"26.9"N 104°29"17.0"E 22 西偏北 585 89 Table 1. Basic characteristics of plots
-
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个。调查并记录每个样方内物种名称、数量、平均高度、盖度等指标。
-
物种的优势度由其重要值表示，重要值（VI）=（相对密度+相对频度+相对盖度）/3。本研究草本层使用相对高度代替相对密度，乔木层中的相对盖度主要为相对基盖度。
采用Simpson多样性指数、Shannon-Wiener多样性指数、Hurlbert均匀度指数和Margalef丰富度指数进行计算[17]。
Simpson多样性指数D：
Shannon-Wiener多样性指数H：
Hurlbert均匀度指数E：
Margalef丰富度指数R：
式（1）~式（4）中：Pi为物种i的数量百分比，即Pi=Ni/N，其中Ni为样方中第i物种的个体数，N为全部种的个体总数；S为样方面积群落中植物种数。
Jaccard物种相似性指数I：
式（5）中：c表示A和B 2个样地群落的共有植物种数；a表示样地A群落的植物种数；b表示样地B群落的植物种数。
1.1. 研究区概况
1.2. 试验样地布设
1.3. 植被调查
1.4. 数据处理
-
研究发现，不同次生灾害迹地植被恢复较好，植被恢复初期以草本为主，有少量的灌木或乔木种。受损迹地群落结构简单，其植物相对比较单一，主要是菊科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）。
样地类型 物种组成 生长型 科/个 属/个 种/个 草本/种 灌木/种 乔木/种 B1 10 20 20 13 6 1 B2 13 20 26 22 2 3 B3 14 25 26 17 9 0 B4 29 39 40 25 15 0 B5 23 28 29 19 7 3 ck 24 30 30 17 9 4 Table 2. Species composition and growth form of diiferent slash communities
-
重要值是衡量群落中植物是否为优势种的一项重要指标。对重要值的分析可找出群落中的主要优势树种[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位的分别为山麦冬、蜈蚣蕨、地果。
样地类型 乔木层 灌木层 草本层 物种 重要值 物种 重要值 物种 重要值 桤木 1.00 鸡骨柴 0.55 飞蓬 0.26 B1 红麸杨 0.13 千里光 0.20 水麻 0.13 小蓬草 0.10 桤木 0.82 水麻 0.77 千里光 0.18 B2 银杏 0.09 长叶水麻 0.23 五节芒 0.17 铁木 0.09 蜈蚣蕨 0.14 水麻 0.37 千里光 0.18 B3 木蓝 0.25 五节芒 0.16 铁木 0.11 蒙古蒿 0.15 马桑 0.22 荩草 0.22 B4 水麻 0.16 五节芒 0.11 披针叶荚莲 0.12 野棉花 0.10 盐肤木 0.42 水麻 0.30 葛藤 0.17 B5 桤木 0.34 马桑 0.22 白荀蒿 0.16 铁木 0.24 披针叶荚莲 0.19 薄荷 0.16 樟树 0.51 刺花椒 0.21 山麦冬 0.28 ck 桤木 0.38 豪猪刺 0.19 蜈蚣蕨 0.21 八角枫 0.07 野蔷薇 0.19 地果 0.13 说明：银杏Ginkgo biloba，盐肤木Rhus chinensis，铁木Ostrya japonica，鸡骨柴Elsholtzia fruticosa，红麸杨Rhus punjabensis，长叶水麻Debregeasia longifolia，木蓝Indigofera tinctoria，马桑Coriaria nepalensis，披针叶荚蒾Viburnum lancifolium，刺花椒Zanthoxylum acanthopodium，飞蓬Erigeron acer，小蓬草Conyza canadensis，蜈蚣蕨Pteris vittata，荩草Arthraxon hispidus，葛藤Argyreia seguinii，薄荷Mentha haplocalyx，白苞蒿Artemisia lactiflora，山麦冬Liriope spicata，地果Ficus tikoua Table 3. A sorting of species important value of different slash communities(Top 3)
-
植物种种间相关性是植物群落重要的数量和结构特征之一，是衡量2个物种相似性的一种尺度。由表 4可知：次生灾害迹地物种与未受损林地物种存在较大差异，Jaccard物种相似性指数均在0.1左右，其中落石区（B4）物种与未被破坏地物种较为接近，相似性系数为0.19；崩塌地（B1）植物种与未被破坏地只有1种共有物种，相似性系数为0.02。次生灾害迹地与未受损林地物种相似性系数大小依次为B4＞B5＞B2＞B3＞B1。
样地类型 B1 B2 B3 B4 B5 ck B1 1.00 B2 0.21 1.00 B3 0.24 0.30 1.00 B4 0.18 0.27 0.27 1.00 B5 0.20 0.28 0.28 0.28 1.00 ck 0.02 0.10 0.06 0.19 0.13 1.00 Table 4. Species similarity index between different slash communities
相比而言，不同次生灾害迹地之间物种更为相似，其物种相似性指数为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。由于未受损林地植被经过多年发育，植被群落达到一个相对稳定状态，而受损地区的林地是在原生演替初期，所处的发育阶段不同，因此群落物种相似性低。
-
地震加剧了生境的破碎化程度，并引发次生灾害，造成了大片的植被破坏，减弱了生态功能。植物多样性在一定程度上能够反映植被恢复程度，是衡量植被演替进程的一项重要指标。对于地震后不同次生灾害多样性指数的计算和分析能够反映植物群落的生境差异、群落的结构类型、演替阶段和稳定性程度[19]。研究样地的植物分乔木、灌木和草本3种生长型，其多样性指数见表 5。
样地类型 乔木层 灌木层 草本层 D H E R D H E R D H E R B1 0 0 0 0.56 1.23 0.68 1.89 0.78 1.93 0.75 2.46 B2 0.26 0.51 0.46 0.76 0.28 0.45 0.65 0.56 0.89 2.51 0.83 3.66 B3 0.73 1.64 0.74 1.98 0.87 2.29 0.81 2.88 B4 0.89 2.41 0.89 3.82 0.89 2.52 0.78 4.00 B5 0.65 1.07 0.98 0.72 0.81 1.65 0.84 1.84 0.85 2.31 0.78 3.22 ck 0.57 0.96 0.70 0.82 0.83 1.92 0.87 1.79 0.81 1.97 0.70 2.87 说明：D为Simpson指数，H为Shannon-Wiener指数，E为Hurlbert均匀度指数，R为Margalef丰富度指数 Table 5. Diversity indices of different slash communities
乔木层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群落多样性指数和丰富度指数并不高，原因可能是稳定的乔木和灌木群落影响到林下草本层植物的多样性。
2.1. 不同灾害迹地群落组成及结构特征
2.2. 不同灾害迹地植物群落优势种
2.3. 不同灾害迹地群落物种相似性
2.4. 不同灾害迹地群落物种多样性特征
-
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样地最大，但不同次生灾害迹地之间差别较小。
-
5·12地震引起的次生灾害迹地土壤表层堆积着大量的松散碎石，土壤养分缺乏，土壤肥力较差，保水能力弱，水土流失严重，不利于植被恢复，这与大多数研究结果一致[3, 13]。菊科、禾本科植物在困难立地早期植物群落里占显著优势，对受损后植被的恢复和延续有重要作用，是典型的先锋植物种[20]。据中国地震局公布的资料，仅2008年5-10月4级以上的余震，北川县就发生了22次[21]。在地震发生之后的3 a内，北川县因地表破碎、暴雨等影响，次生灾害时常发生[22]，又因为缺失植被自然恢复的优良条件，使得植被自然恢复进程缓慢。
在植物多样性指数方面，次生灾害迹地群落多样性指数并不都低于未被破坏的桤木-樟木混交林。在草本层除崩塌迹地外，其他受损林地多样性指数均大于未被破坏的桤木-樟木混交林，在灌木层只有落石迹地多样性指数大于未被破坏的桤木-樟木混交林，在乔木层只有堰塞湖迹地多样性指数大于未被破坏的桤木-樟木混交林。就植被演替进程来看，受损林地仍处在自然演替的初期即草本群落阶段，该阶段群落环境较差，以1~2年生植物为主，但随着时间的推移群落种类逐渐增多，多年生草本植物逐渐占优势，并发展为较为茂密的草本植物群落。演替总是朝着适应性强的方向演进，即按照草本—灌木—乔木的顺序变化着[23]。对于未受损林地而言，已经处于基本稳定的阶段，乔木的生长，冠层郁闭度的增加，使得林下光照减少，灌木和草本的生存受到一定的限制且趋于稳定[24]。
植被恢复是次生灾害迹地重建的一项重要工作，植被的自然恢复与演替又是一个长期的过程，所以对于研究区应进行连续长期的观测与调查，更加深入地研究植被自然恢复与演替的过程，为次生灾害迹地植被的恢复与重建提供科学依据。比如采用合适的自然恢复与人工恢复相结合的方法，加快次生灾害迹地的植被恢复工作；其次，应选取乡土植物种类，同时兼顾植物群落的多样性，尽量做到乔灌草搭配适宜，使受损立地的植被得到快速有效的恢复。
DownLoad: