神农架林区华山松大小蠹危害风险分析洪承昊1,2，柏冰洋2，易家喜2，孙宜平2，谢志强3，向劲松2 (1.湖北省林业科学研究院，湖北武汉430079；2.神农架林区林业科学研究所，湖北神农架442499；
3.神农架森林病虫害防治检疫站，湖北神农架442499)
摘要：利用景观生态学的研究方法，以神农架林区华山松大小蠹Dendroctonus armandi为研究对象，从华山松Pinus armandii景观空间格局分析林区华山松大小蠹发生情况，并对林区各林场乡镇的风险程度进行分析。结果表明：华山松大小蠹发生与寄主华山松的面积分布占比、斑块密度呈正相关，相关系数R值分别为0.504和0.379；与华山松斑块个数、景观形状指数不相关。林区华山松面积占比最大的地区是温水林场，斑块密度最大的地区是红坪林场，因此温水林场和红坪林场华山松大小蠹危害风险最大，危害风险最低的是宋洛乡和红花朵林场。风险预测结果显示林区中部和西部的华山松斑块受侵害的风险最大。 关键词：华山松大小蠹；景观格局；风险分析
中图分类号：S763.38文献标志码:A文章编号：1671-0886(2021)02-0017-06
DOI：10.19688/jki.issn1671-0886.20200057
Risk analysis of Dendroctonus armandi in Shennongjia Forestry District/HONG Chenghao1,2,et al.
武汉木兰草原一日游攻略
(1.Hubei Academy of Forestry,Wuhan430079,China； 2.Shennongjia Institute of Forestry, Shennongjia442499,China)
Abstract:The occurrence of Dendroctonus armandi and the risk level of each forest farm and town in Shennongjia Forestry District were analyzed from the landscape spatial pattern of Pinus armandii using research methods of landscape ecology.The results showed that the situation of D.armandi had positive correlation with the percentage of landscape type(PLAND)and patch density(PD)of P.armandii，with Pearson correlation coefficient R of0.504and0.379,and had no correlation with number of patches(NP)and landscape shape index(LSI).The region with the biggest PLAND was Wenshui Forest Farm，and the region with the biggest PD was Hongping Forest Farm，thus the two regions had the highest risk，and the risk in Songluo Town and Honghuaduo Forest Farm were the lowest.The risk forecast results showed that P.armandii in the central and western part of Shennongjia Forestry District had the greatest possibility of being invaded.
Keywords:Dendroctonus armandi；landscape pattern；risk analysis
华山松大小蠹Dendroctonus armandi Tsai and Li，又名大凝脂小蠹，隶属鞘翅目Coleoptera小蠹科Scolytidae大小蠹属Dendroctonus，为中国特有种，主要分布在陕西、四川、甘肃、河南和湖北等地
区，是侵入华山松Pinus armandii的先锋虫种，主要危害30年以上的华山松健康树［1-2］，也是我国重要的毁灭性森林害虫之一［3］。华山松又名五针松，是松科中的著名常绿乔木之一，原产于中国，因集中产于陕西的华山而得名。华山松是我国温带及亚热带地区的重要山地针叶林树种⑷，广泛分布于我国北方的宁夏、山西、陕西、河南、甘肃等省市区，南方主要分布在湖北西部、四川、贵州、云南直至西藏东南部山
收稿日期：2020-09-11；修回日期：2020-12-16；网络首发：2020-12-17
基金项目：国家重点研发计划“人工林重大灾害的成灾机理和调控机制”(2017YFD0600100)
第一作者：洪承昊(1982—)，男，湖北武汉人，博士，从事森林保护及景观生态调控工作,E-mail：haohao730@163 通信作者：向劲松，高级工程师，主要从事森林资源培育相关研究,E-mail：1911338058@qq 。
区［5］。世界范围内，在缅甸、日本也有引种分布。华山松是我国秦巴山区的主要造林树种，也是我国特有的乡土树种，在森林生态系统中发挥着涵养水源、净化空气等重要生态功能，是该区域森林植物落演替的重要指标性树种，在我国经济、社会、环境建设和地理分布等方面具有举足轻重的地位。
我国最早于1932年在秦岭巴山林区首次发现华山松遭华山松大小蠹侵害［6］。近年来，由于华山松大小蠹的区域性危害，给我国湖北、陕西、甘肃、四川、云南和贵州等秦巴山区华山松的生长生存造成了前所未有的威胁，特别是在重要水源涵养区和风景名胜地的神农架出现了暴发危害，截至2018年成灾面积约0.67万hm2(神农架森林病虫害防治检疫站数据)，造成的影响和破坏更加严重，仅2012—2018年，华山松大小蠹导致湖北神农架地区直接经济损失1000多万元，间接生态效益损失上亿元。
由于华山松大小蠹的危害极为隐蔽，幼虫在韧皮部下进行取食，不易察觉，加之世代重叠［7］，造成防治不及时，防治方法相对落后，因此，研究和评价华山松大小蠹的扩散风险则显得尤为重要。景观斑块的分布和空间格局是影响有害生物流行扩散的直接因素，其次是人为活动、经济因素等［8-9］o近年来，有不少学者从华山松大小蠹的生理生化、行为学和引诱剂研制等方面做了许多研究［10-12］，然而从景观斑块层面研究和控制华山松大小蠹的报道却较少［13］。景观结构和功能的关系是研究森林保护学的重要方向［14-16］o研究表明，有害生物灾害发生程度与森林景观结构关系密切［17-18］。
神农架林区华山松分布广且垂直落差大，因此在神农架林区研究华山松大小蠹具有很好的代表性。笔者通过研究神农架林区的华山松空间分布，首先从景观层面分析华山松林斑的空间格局差异，对华山松大小蠹年均发生情况进行拟合分析，其次对神农架林区各林场乡镇进行风险评估和分析，为控制和预防华山松大小蠹的蔓延扩散提供理论和技术支持。
1研究地概况
神农架林区是我国唯一以“林区”命名的行政地区，湖北省直辖管理。地处湖北省西北部，属秦巴山区。地跨东经109°56'~110°58‘，北纬31°15'~ 31°75'，全区总面积3253平方公里，林地占区域面积85%以上，生物资源和矿产资源丰富。林区平均海拔1800m，最高点神农顶海拔约3106m，最低点为林区下谷坪乡约398m，最大海拔相对落差为2707.4m，因此造就林区独特的气候环境条件。
2数据来源和研究方法
2.1数据来源数据本底资料来源于湖北省林业调查规划院2019年二类调查数据，通过Arcmap［19］建立属性数据库，采取定点踏查核实，数据精确。将原始矢量图像进行坐标转换，将景观类型重新定义划分，分类后将数据资料生成分辨率为10m的栅格图像，进行景观指数运算。华山松大小蠹发生数据来源于神农架森林病虫害防治检疫站,2012—2017年神农架森林病虫害防治检疫站调查数据显示，木鱼林场、红坪林场、温水林场和红花朵林场华山松大小蠹危害最为严重，新华林场、宋洛乡、新华镇和徐家庄林场没有华山松大小蠹发生，林区其余地区无寄主或有极小量面积分布。
2.2分析方法采用一、二级结合分类体系重新划分土地利用类型，参照《土地利用现状分类》(GB/T 21010—2017)，并考虑研究对象危害特点，将华山松和针阔叶混交林划分为一类，其他非华山松大小蠹寄主划分为一类。神农架林区华山松大部分为20世纪70年代人工造林，均已达成熟林，各林场寄主
分布具体情况见表1。采用Fragstats景观格局指标系统描述斑块类型景观水平特征［20］，选取斑块个数(NP)、斑块密度(PD)和景观形状指数(LSI)分析华山松斑块景观格局差异。
斑块数(NP)是指类型水平或者景观水平上总的斑块数目，是标志景观破碎化程度的重要指标之一；其计算公式为NP=N,N为景观中的斑块总数。
斑块密度(PD)是单位面积上的斑块数，是描述景观破碎化的重要指标之一,PD值越大，斑块破碎化程度越大；其计算公式为PD=N/A X10000,N为景观中的斑块总数,A为景观总面积。
景观形状指数(LSI)是结合景观面积对景观总边缘长度或边缘密度的标准化度量，也可看作是对斑块聚集或离散程度的度量；其计算公式为LSI= E/minE,E指景观的边缘总长度,minE是E的最小可能值。
通过景观格局指标分析各林场乡镇的华山松大小蠹风险等级，并根据分析结果给各林场乡镇赋予
一定风险评估值，从景观层面分析华山松大小蠹在 神农架林区的风险危害分布情况。
表1各林场乡镇华山松分布情况
Tab. 1 The distribution of Pinus armandii in each forest farm and town
地区斑块总面积平均斑块面积林龄结构
Site
CA/hm 2Mean patch area/hm 2
Forest age structure
温水林场 Wenshui Forest Farm 2 430. 9230. 77
成熟林木鱼林场 Muyu Forest Farm 3 338. 1121 . 96成熟林红坪林场 Hongping Forest Farm 437. 628. 26
成熟林红花朵林场 Honghuaduo Forest Farm
535. 0315.29成熟林新华镇 Xinhua Town
778. 63
13.20成熟林徐家庄林场 Xujiazhuang Forest Farm
103. 03 1.34
成熟林新华林场 Xinhua Forest Farm 301 . 21
15.06成熟林宋洛乡 Songluo Town
368. 31
金茂大厦图片
5. 34
成熟林
3 结果与分析
3.1景观类型水平差异寄主华山松斑块面积百
分占比最大的是温水林场，达到13. 14%，最小的是
宋洛乡，仅为1. 46% ；斑块数量最多的是木鱼林场,
有152个，最少的是红花朵林场，有35个（表2）。 平均斑块面积最大的是温水林场，达到30.77 hm 2；
表2 各林场乡镇华山松景观格局指数及华山松大小蠹年均发生面积
最小的是徐家庄林场，仅有1.34 hm 2（表1）。斑块
密度最大的是红坪林场，为54个/ hm 2，斑块密度最 小的是红花朵林场，仅为13个/ hm 2；形状复杂程度 最高的是木鱼林场，景观形状指数为20. 98 ，形状复
杂程度最低的是红花朵林场，景观形状指数仅为
9. 10（表 2）。
Tab. 2 The landscape pattern of P. armandii and annual occurrence area of D. armandi in each forest farm and town
地区Site
斑块面积占比PLAND/%斑块数量/个
NP 斑块密度/ （个/hm 2）
PD 景观形状指数
LSI 年均发生面积/hm 2Annual occurrence area
温水林场 Wenshui Forest Farm 13.14794314. 93206. 67木鱼林场 Muyu Forest Farm 9. 22152
42
20. 98107. 80
红坪林场 Hongping Forest Farm 4. 45
535411.93
238. 87
红花朵林场 Honghuaduo Forest Farm 2. 0035
139.10138.87新华镇 Xinhua Town
6. 435949
14. 27
河南自驾游大全景点排名0徐家庄林场 Xujiazhuang Forest Farm 4. 17
77南充
3112.320新华林场 Xinhua Forest Farm 2. 352016
6. 82
0宋洛乡 Songluo Town
1 . 4669
2714. 00
3.2景观指数分析横向比较华山松大小蠹发生
区和未发生区，在斑块面积百分比均值方面，发生区
是7.20%，显著高于未发生区的3.60% （表2）；在
斑块个数、斑块密度和景观形状指数方面也都高于
未发生区。Pearson 相关性分析表明，斑块面积百分 比和斑块密度与华山松大小蠹的年均发生面积有弱
正相关，相关系数分别为0.504和0. 379，斑块个数
和景观形状指数与华山松大小蠹的年均发生面积无
相关性（ 表 3） 。
表3景观格局指数与华山松大小蠹年均发生面积的相关性Tab. 3 The correlation of P. armandii landscape pattern with
annual occurrence area of D. armandi
景观格局指数Landscape indices
相关性 Correlation
R （ Pearson 相关性）Sig （双侧）
斑块面积占比 PLAND 0. 5040. 138斑块数量 NP
0.169
0. 641斑块密度 PD 0. 3790. 280景观形状指数 LSI
0.193
0. 594
面积百分比和斑块密度与华山松大小蠹的危害
程度呈正相关，因此可建立相关性模拟。模拟结果
显示，面积百分比和华山松大小蠹发生面积之间呈
正相关，相关方程为Y = 24.61 +11.46X ；斑块密度
与华山松大小蠹发生面积之间呈正相关，相关方程 为Y = 2.83+2.41 X （图1）。斑块个数和景观形状 指数与华山松大小蠹的危害程度之间无相关性，因
此不做线性模拟分析。
-50050052 2
江西两日游适合去哪里
1
1
7=24.61+11.46^
0 2 4 6 8 10 12 14面积百分比/%
Percentage oflandscpae(PLAND)
o w o w o
5 o 501522 1 1ZUJZQoualm 。。。
o q H
z
日
畠
M
册叙
、F=2.83+2.41X
0 10 20 30 40 50 60
斑块密度/(个/hm?)Patch density(PD)
图1基于面积百分比和斑块密度的相关性模拟
Correlation simulation based on area percentage and patch density
坪林场、徐家庄林场、新华林场、红花朵林场和宋洛
乡，且各林场乡镇间存在显著差异（sig < 0. 05）
（图2）。斑块密度分析结果显示，华山松大小蠹危害 风险最高的是红坪林场，高于其他各林场乡镇，往后
依次是新华镇、温水林场、木鱼林场、徐家庄林场、宋
洛乡、新华林场和红花朵林场，且各林场乡镇间存在 显著差异（sig <0. 05）（图2）。
Fig.1
3.3基于景观指数的华山松大小蠹危害风险评估 华山松斑块的面积百分比和斑块密度与华山松大 小蠹发生面积呈正相关。因此，分析各个乡镇华山松
面积和斑块密度可预测神农架林区华山松大小蠹的
发生情况。温水林场华山松面积百分比最高，占比为 13. 14%，显著高于其他林场或乡镇，受到华山松大小
蠹危害的风险等级最高；其次是木鱼林场、新华镇、红
(Q N V u d v E d o s p u z c o O M E g o J a d
各乡镇林场
Forest farm and town
(Z
日悭图张靂
图2各林场乡镇华山松面积百分比和斑块密度横向比较
Compared the area percentage and patch density in each forest farm and town
Fig.2
综合华山松面积百分比和斑块密度指数两项景 观指标，对各个林场乡镇华山松大小蠹危害风险进 行赋值。基于面积百分比的赋值从小到大依次为:
宋洛乡，1；红花朵林场,2 ；新华林场,3 ；徐家庄林场,
4 ；红坪林场,
5 ；新华镇,
6 ；木鱼林场,
7 ；温水林场,8。
基于斑块密度的赋值从小到大依次为： 红花朵林场, 1；新华林场,2；宋洛乡,3；徐家庄林场,4；木鱼林场,
5； 温水林场,6； 新华镇,7； 红坪林场,8。 赋值后, 进 行归一化处理（表4）
。
表4各林场乡镇风险指数归一化处理值Tab.4Normalize value in each forest farm and town
地区Site
风险系数Risk factor 斑块面积占比
PLAND
斑块密度
PD
温水林场Wenshui Forest Farm10.71
木鱼林场Muyu Forest Farm0.860.57
新华镇Xinhua Town0.710.86
红坪林场Hongping Forest Farm0.571
徐家庄林场Xujiazhuang Forest Farm0.430.43
新华林场Xinhua Forest Farm0.290.14
红花朵林场Honghuaduo Forest Farm0.140
宋洛乡Songluo Town00.29
牡丹江百姓网
对各林场乡镇的华山松大小蠹危害风险指数进行赋值后，采取克里金插值（kriging interpoloction）分析方法，得出林区华山松大小蠹危害风险预测图（图3）。
直观分析图3可知，基于华山松面积百分比的华山松大小蠹危害风险预测显示风险区集中在林区西部地区和中部地区（图3A），该区域也是林区华山松主要分布区。基于华山松斑块密度的华山松大小蠹危害风险预测结果与上述一致，风险高的区域主要集中在林区中部地区和西部地区（图3B）。由于各林场乡镇在寄主面积占比和斑块密度两个指标的大小有差异，因此分析结果有一定差异。基于寄主面积分布占比的分析结果显示风险最高的在最西部，其次是中部地区；基于斑块密度的分析结果显示风险最咼的是中部偏北地区，其次是西部和中部偏南地区。结合两图可知，林区风险高的区域集中在西部和中部地区，其他位置风险较低。
N N
A.基于华山松面积百分比Based on area percentage of P.armandii；
B.基于华山松斑块密度Based on patch density of P.armandii
图3神农架林区华小松大小蠹风险分析
Risk analysis of Dendroctonus armandi in Shennongjia Forestry District
Fig.3
4讨论
笔者针对神农架林区华山松大小蠹寄主植株华山松进行了景观格局层面的分析研究，从华山松斑块的空间格局分析了与华山松大小蠹发生程度之间的关系。结果显示，华山松在林场乡镇的面积分布和斑块密度与华山松大小蠹发生程度之间有正相关性。寄主分布面积越广泛，会提供给华山松大小蠹更多的资源，发生的程度也会更大。斑块的密度越高，单位面积内的斑块数量更大，华山松大小蠹在扬飞期的传播扩散成功率更高，定殖能力和概率也会更高，进一步传播扩散的风险也更大［21］，因此风险等级也更高。
通过景观格局控制森林病虫害是森林保护学的一个重要手段，本研究可以为综合控制华山松大小蠹的危害和扩散提供一定基础［22］。在华山松斑块密度较大的区域，可以通过改造寄主斑块或减小斑块数量，切断邻近度高的斑块联系，形成抑制华山松大小蠹扩散的屏障廊道，从而达到控制华山松大小蠹危害的目的。
笔者仅从华山松的景观层面进行了分析，景观中的其他异质性斑块，人为活动程度和公路道路的建设程度等也是影响生物扩散和蔓延的重要因子之一，实际发生过程中也会对昆虫的行为产生较大的影响，如何定量分析景观生态学和森林有害生物之间的关系还需要进行更为深入的研究。
志谢：感谢湖北神农架森林生态系统国家定位观测研究站给予经费支持并协助开展工作
。
