The Chinese boreal forest is an important forest resource in China. However, it has been suffering serious disturbances of forest fires, which were caused equally by natural disasters (e.g., lightning) and human activities. The literature on humancaused fires indicates that climate, topography, vegetation, and human infrastructure are significant factors that impact the occurrence and spread of humancaused fires. But the studies on humancaused fires in the boreal forest of northern China are limited and less comprehensive. This paper applied the spatial analysis tools in ArcGIS 10.0 and Logistic regression model to investigate the driving factors of humancaused fires. Our data included the geographic coordinates of humancaused fires, climate factors during year 1974-2009, topographic information, and forest map. The results indicated that distance to railway (x1) and average relative humidity (x2) significantly impacted the occurrence of humancaused fire in the study area. The logistic model for predicting the fire occurrence probability was formulated as P=1/[1+e-(3.026-0.00011x1-0.047x2)] with an accuracy rate of 80%. The above model was used to predict the monthly fire occurrence during the fire season of 2015 based on the HADCM2 future weather data. The prediction results showed that the high risk of humancaused fire occurrence concentrated in the months of April, May, June and August, while April and May had higher risk of fire occurrence than other months. According to the spatial distribution of possibility of fire occurrence, the high fire risk zones were mainly in the west and southwest of Tahe, where the major railways were located.
全 文 :我国北方针叶林人为火发生的预测模型∗
郭福涛1 苏漳文1 王光玉2 王 强3 孙 龙3∗∗ 杨婷婷1
( 1福建农林大学林学院, 福州 350002; 2加拿大不列颠哥伦比亚大学林学院, 温哥华 V6T 1Z4; 3 东北林业大学林学院, 哈尔
滨 150040)
摘 要 我国北方针叶林带是重要的森林资源储藏地,也是林火发生的重灾区,其自然火和
人为火所占比例相当. 气象因子、地形特征、植被条件、人为基础设施等因素对人为火发生具
有显著影响,国内目前应用空间分析技术对北方针叶林带人为火影响因子的研究还存在一定
不确定性. 本文基于 1974—2009年间人为火的空间地理坐标,结合研究地的气象因子、基础
地理信息及矢量化林相图,应用 ArcGIS 10.0中的空间分析工具和 SPSS 19.0的逻辑斯蒂回归
模型对影响人为火发生的主要驱动因子进行分析,并建立人为火发生的概率模型. 利用 HAD⁃
CM2模式下研究区域未来气象数据对塔河地区 2015 年人为火发生情况进行计算.结果表明:
距离铁路距离(x1)和平均相对湿度(x2)对研究区域人为火发生具有显著影响,并得到火险概
率模型 P= 1 / [1+e-(3.026-0.00011x1-0.047x2)] . 模型校验结果显示,模型的准确度可达到 80%.林火发
生预测结果表明,塔河地区 2015年 4—6月、8 月为人为火高发期,其中,4—5 月的林火发生
概率最高.从火险空间分布来看,高火险主要集中在塔河西部和西南部,铁路线路主要包含在
此区域.
关键词 北方针叶林; 人为火; 逻辑斯蒂回归; 空间分析
文章编号 1001-9332(2015)07-2099-08 中图分类号 S762.2 文献标识码 A
Prediction model of human⁃caused fire occurrence in the boreal forest of northern China.
GUO Fu⁃tao1, SU Zhang⁃wen1, WANG Guang⁃yu2, WANG Qiang3, SUN Long3, YANG Ting⁃ting1
( 1College of Forestry, Fujian Agriculture and Forestry University, Fuzhou 350002, China; 2Faculty
of Forestry, University of British Columbia, Vancouver V6T 1Z4, Canada; 3School of Forestry, Uni⁃
versity of Northeast Forestry, Harbin 150040, China) . ⁃Chin. J. Appl. Ecol., 2015, 26(7): 2099-
2106.
Abstract: The Chinese boreal forest is an important forest resource in China. However, it has been
suffering serious disturbances of forest fires, which were caused equally by natural disasters (e.g.,
lightning) and human activities. The literature on human⁃caused fires indicates that climate, topog⁃
raphy, vegetation, and human infrastructure are significant factors that impact the occurrence and
spread of human⁃caused fires. But the studies on human⁃caused fires in the boreal forest of northern
China are limited and less comprehensive. This paper applied the spatial analysis tools in ArcGIS
10.0 and Logistic regression model to investigate the driving factors of human⁃caused fires. Our data
included the geographic coordinates of human⁃caused fires, climate factors during year 1974-2009,
topographic information, and forest map. The results indicated that distance to railway ( x1) and
average relative humidity ( x2) significantly impacted the occurrence of human⁃caused fire in the
study area. The logistic model for predicting the fire occurrence probability was formulated as P= 1 /
[1+e-(3.026-0.00011x1-0.047x2)] with an accuracy rate of 80%. The above model was used to predict the
monthly fire occurrence during the fire season of 2015 based on the HADCM2 future weather data.
The prediction results showed that the high risk of human⁃caused fire occurrence concentrated in the
months of April, May, June and August, while April and May had higher risk of fire occurrence
than other months. According to the spatial distribution of possibility of fire occurrence, the high fire
risk zones were mainly in the west and southwest of Tahe, where the major railways were located.
Key words: boreal forest; human⁃caused fire; Logistic regression; spatial analysis.
∗福建农林大学林学院青年科学基金项目(6112C035K)、霍英东基金基础项目(131029)和黑龙江省博士后科研启动经费项目资助.
∗∗通讯作者. E⁃mail: sunlong365@ 126.com
2014⁃08⁃12收稿,2015⁃03⁃02接受.
应 用 生 态 学 报 2015年 7月 第 26卷 第 7期
Chinese Journal of Applied Ecology, Jul. 2015, 26(7): 2099-2106
林火(包括雷击火和人为火)是北方针叶林中
十分重要的生态干扰因子,影响森林更新和演
替[1], 还将直接导致森林资源损失,并大大增加森
林经营管理成本[2] .林火是一个全球性问题,每年有
大约 45000 次林火发生在欧洲,过火面积达到 26×
105 hm2[3];加拿大平均每年大约发生 10000 次林
火,过火面积大约 25×105 hm2[4];我国每年发生林火
大约 10000次,年均过火面积为 82×104 hm2[5] .在森
林火灾中,人为火(由于人类活动或人为基础设施
而引起的火灾)占很大比重[6-7] .位于北美阿拉斯加
的北方针叶林带,年均人为火发生次数约占到总林
火次数的 60%以上,过火面积随着人为扑救措施和
技术手段的进步,近年来呈递减趋势,从 38%
(1950—1959年)降到 5%(2000—2005 年),但火灾
次数却一直居高不下[8] .我国南方地区人为火比重
也非常大,如江西、福建省人为火多在 90% 以
上[9-10] .但由于南方全年湿度和降水量较大,人为火
发生后如及时进行扑救,不容易造成大面积森林损
失.我国北方针叶林带(大兴安岭地区)的植被和气
候特征与加拿大及美国北方针叶林带十分相似,林
分树种组成多为针叶林,可燃性较强,春秋季节区域
气候较为干旱,一旦发生森林火灾,扑救难度很大,
容易造成大面积森林损失.研究表明,我国北方针叶
林带人为火的比重约为 50%~60%[11] .因此,了解该
区域人为火发生的主要影响因子, 对人为火的预测
及管理具有重要意义.
近年来国内很多学者对人为火的时空分布格
局、火源组成、人为火与气象因子的关系进行了广泛
研究, 取得了大量研究成果[12-17] .国外研究显示,除
了气象因素外,林内可燃物类型及人为基础设施等
因素对人为火的发生也具有显著影响[18-22] .综合分
析气象、可燃物、地形、人为基础设施等因素与人为
火的关系,对理解我国北方针叶林带人为火发生特
性具有重要意义,也是对该地区林火进行预测的基
础.为此,本文利用空间分析手段,探讨了大兴安岭
塔河地区人为火发生的驱动因子,并利用所获得的
驱动因子和预测模型,计算了 2015年该地区防火期
内的火险分布情况.
1 研究地区与研究方法
1 1 研究区概况
塔河地区位于我国黑龙江大兴安岭林区(图
1),属于典型的北方针叶林带,具有植物多样性
低、乔木组成树种简单等特点,植被以兴安落叶松
图 1 研究区域示意图
Fig.1 Sketch map of the study area.
(Larix gmelinii)林和白桦(Betula platyphylla)林为
主.该区属寒温带季风气候,年均气温-2 ~ 4 ℃,年
降水量 350~ 500 mm,相对湿度 70% ~ 75%.年积雪
期达 5个月,林内积雪深可达 30 ~ 50 cm.历史上塔
河地区是火灾多发区,人为火源和雷击火源共同存
在,从次数来说,人为火约占火灾总数的 2 / 3,从面
积来说,雷击火的过火面积高于人为火.该地区雷击
火在空间尺度上多呈现聚集分布,且多集中在铁路
和居民区附近,时间尺度上多集中在 5—9 月,4—6
月是人为火高发期[23-24] .
1 2 数据获取
本文数据来源包括 3部分,分别为林火数据、矢
量化林相图和气象数据.林火数据来源于大兴安岭
地区塔河县森林防火办公室,包括 1974—2009 年塔
河地区林火发生情况数据(起火地理坐标、起火原
因、发生时间、过火面积等).
1 2 1基础地理和林相数据 研究中所需基础地理
地形图和矢量化林相图等数据来源于东北林业大学
森林经理教研室.1 ∶ 10万的基础地理图和矢量化林
相图数据库分别建于 2000 和 2003 年.其中,林相图
数据库精确到小班,并对林下植被类型和土壤条件
进行了很好的记录,为本研究提供了有力的数据支
持.本文参考国外相关研究[6,18,20,22],通过 ArcGIS 软
件的空间分析技术计算各火点和对照点到河流、街
区、铁路和机耕路等人为基础设施距离,以及火点和
对照点所对应的林型、坡度、坡向、地被物盖度和优
势树种等林分和地形特征.
1 2 2历史气象数据 火点和随机点对应的当日气
象数据(历史数据)来源于中国气象数据共享网络
(http: / / cdc.cma.gov. cn / ) .由于空间异质性的普遍
0012 应 用 生 态 学 报 26卷
存在,原则上在提取研究区内火点相应的每日气象
数据时,需要对该研究区域内所有气象站点的数据
进行插值.然而,由于本研究地只有一个基础气象站
(50246),无法进行区域插值.为了尽可能减少空间
异质性对每日气象数据的影响,本文利用邻近气象
站数据,对所用气象数据进行矫正.步骤如下:1)确
定了塔河附近所有国家级基础气象站的位置(5 个
气象站覆盖了整个大兴安岭地区)(图 2),发现如果
以塔河地区东西长度为半径画圆,作为气象站的有
效覆盖范围,则 50136(国家标准气象站号)、50349、
50353和 50442这 4个气象站对塔河地区林火影响
微弱.但若以塔河气象站(50246)为圆心,则火点和
随机点距离气象站点的距离有显著差异.2)为减少
偏差,我们以气象站 50136 和 50246 为圆心的两个
圆的切点为新站点(虚拟气象站点),其气象数据为
两个气象站的平均值.3)以新站点和 50246 为基础,
对所有火点和随机点进行最近距离计算,火点和随
机点的气象数据将由与其最近的气象站决定.
气象数据包括 1974—2009 年的每日极大风速
等 22个气象因子,为保证数据的完整性,本文对气
象数据进行预处理, 将由自然灾害或监测设施损坏
等原因而缺失的气象数据进行剔除.并将使用的气
象站原始气象数据与周边气象站数据进行对比,排
查异常点,保证数据的准确性.最终确定进入分析的
气象数据包括日最低气温(℃)、日最高气温(℃)、
24 h降水量(mm)、平均风速(m·s-1)、平均气温
(℃)、平均水汽压(hPa)、平均相对湿度(%)、平均
气温(℃)、日照时数(h)、最大风速(m·s-1)、最小
相对湿度(%)共 11个气象因子.
图 2 塔河周边区域国家级气象站分布
Fig.2 Distributions of national meteorological stations around
Tahe station.
1 2 3未来气象数据 为提高预测模型的准确性,
本文应用每日气象数据建立林火发生概率模型.然
而,对研究区域未来的火险形势进行预测,以月或年
为时间尺度则更具实际意义.因此,本文以每日气象
尺度建立的预测模型为基础,将每日尺度模型中涉
及的气象因子进行月尺度转化,即以月尺度气象数
据代替日尺度数据,并通过 ArcGIS 的栅格运算器,
计算出未来研究区域的每月林火发生概率.
月尺度未来气象数据来源于 HADCM2 模式中
国区域数据集.数据包括研究区域未来降水、温度和
相对湿度情况.HADCM2 模式是 Hadley 气候预测和
研究中心发展的海气耦合模型.其中,大气模式水平
分辨率为 2.50°×3.75°(经度×纬度) [25] .模式模拟的
开始时间为 1860 年,情景预测至 2099 年,共模拟
240年.在模拟过程中,1860—1989 年温室气体排放
采用历史实测资料;自 1990 年起采用 IS92a 排放情
景,即温室气体排放综合效果相当于 CO2浓度(以下
简称相当 CO2浓度)逐年递增 1%.HADCM2 中国区
域数据集建立在 HADCM2基础上,原始数据来源于
IPCC数据发布中心,数据为 GRIB 格式,先将其转
化为文本格式,提取中国区域的数据点,结合中国的
DEM,利用 ANUSPLIN软件包进行空间插值.数据集
成为 NETCDF格式,空间分辨率为 10 km,时间分辨
率为月.数据由南京大学国际地球系统科学研究所
编辑,通过地球系统科学数据共享平台(http: / / geo⁃
data.nju. edu. cn / Portal / index. jsp)提供.已有研究表
明,HADCM2模式较好地模拟了中国区域地面气温
和降水的分布特征[26-28] .
1 3 数据处理
1 3 1火点与对照点的创建及空间属性提取 应用
逻辑斯蒂回归模型对人为火发生的概率及驱动因子
进行分析.根据模型使用原理(二项逻辑斯蒂模型原
理介绍)需要随机创建一定数量(通常大于等于火
点数量)的对照点与火点一同进行模型拟合.对照点
的创建追寻时间和空间上的双向随机原则,对照点
创建后,按对照点的具体时间提取相应的气象因子,
同时根据对照点的空间地理坐标提取林型、坡度、到
铁路距离等地形和植被因子,并建立对照点及影响
因子数据库.对照点的创建数量或比例,国际上并没
有统一标准.有学者按照指定时间范围内平均火点
数量选取非火点[24],也有学者以非火点之间的最小
距离为标准在整个研究区域内选取随机点[25] .这些
方法各有利弊,但对照点的“随机性”是选择的基本
原则.此外,由于实际情况往往是研究区域的火点空
10127期 郭福涛等: 我国北方针叶林人为火发生的预测模型
间分布少于或大大少于对照点,所以在随机选取对
照点时,其数量通常不能低于火点数量,以反映林火
发生的客观规律.本文按照火点与随机点 1 ∶ 2 的原
则选取对照点,共统计得到样本 150个.
根据火点和对照点的时间提取相对应的气象因
子,并应用 ArcGIS 10.0对火点和对照点与基础地理
和林相图数据进行空间分析,计算火点及对照点到
河流面、街区、铁路、机耕路及等级公路的距离.并将
林相图中的字符型变量进行属性转换以便进行统计
分析,包括坡向、林型、经营措施和林下优势树种.此
外,坡度级、地被物盖度、龄级和郁闭度也是本文考
虑的因子,但由于其本身以数值型记录,因此不需要
进行属性转换.通过对时间和空间的信息提取,建立
人为火点及影响因子的数据库.
1 3 2建模应用数据的特征描述与分析 本文通过
定位火点、对照点、河流面、街区、铁路、机耕路、等外
公路及等级公路的空间分布(图 3),运用 ArcGIS 栅
格计算器计算火点和对照点到河流等其他设施的距
离.并将火点和对照点图层与林相图属性图层进行
叠加运算,合并各图层属性.同时根据火点和对照点
的时间选取对应的气象因子,建立数据属性数据库.
1 3 3模型拟合及优度检验 二项逻辑斯蒂模型
(Logistic)是目前被广泛应用于林火发生概率模型
研究的方法[21,29-31] .在 Logistic 回归模型中,对两类
判别问题进行取值.第一类为有林火发生,记为 y =
1,第二类为无林火发生,记为 y= 0.根据指标可以建
立 LogitP关于自变量(气象因子、相关林相属性及
火点到河流、铁路、公路等设施的距离)的 Logistic
回归方程:
LogitP= b0+b1x1+b2x2+…+bmxm (1)
图 3 火点、对照点、河流及其他人为基础设施空间分布
Fig. 3 Spatial distributions of fire points, random control
points, rivers and other infrastructures in the study area.
得到概率估计公式:
P= e
b0+b1x1+b2x2+…+bmxm
1+eb0+b1x1+b2x2+…+bmxm
本文根据估计概率进行判别归类,并计算模型
的预测准确率.以往在运用逻辑斯蒂分析林火发生
概率的研究中,通常以 0.5 为临界值,认为高于 0.5
就有林火发生,低于 0.5 则无林火发生.但此种临界
值多为主观选择,缺少统计支持.因此,有学者运用
医学领域用来判断肿瘤是否发生的约登指数,来判
断林火发生的最佳临界值[31-33] .约登指数 =敏感性
(sensitivity)+特异性( specificity ) -1,其中,敏感性
为把实际为真值的判断为真值的概率,特异性为把
实际为假值的判断为假值的概率.一般来说,以约登
指数最大,即(敏感性+特异性-1)达到最大所对应
的值为最佳诊断界值.本文将数据随机分为 70%和
30%两部分,分别用来建模和验证.
此外,本文应用 ROC( receiver operating charac⁃
teristic) [34]检验方法对逻辑斯蒂模型进行拟合优度
分析.ROC曲线用于二分类判别效果的分析与评价,
一般自变量为连续变量,因变量为二分类变量.基本
原理是:通过判断点( cutoff point / cutoff value)的移
动,获得多对敏感性和误判率(1-特异性),以敏感
性为纵轴,以误判率为横轴,连接各点绘制曲线,然
后计算曲线下的面积(area under the curve,AUC),
面积取值范围为 0~1,越接近 1 说明判断价值越高.
本研究应用 SPSS 10.0软件进行建模和模型评估.
1 3 4研究区域未来火险分布 根据建模结果, 结
合 HADCM2 模式下中国 2015 年的月平均气象数
据,应用 ArcGIS 的图层运算功能计算出研究区域
2015年的林火概率分布图.
2 结果与分析
2 1 模型拟合结果及校验
2 1 1模型拟合结果 本文对火点和对照点进行二
元逻辑斯蒂回归拟合.模型的全变量拟合(全模型)
结果显示,模型整体拟合结果 R2为 0.523,模型参数
拟合结果显示距铁路距离、日平均相对湿度、日平均
气温和日最高气温与人为火发生之间均表现出显著
相关性(P<0.05),而其他自变量因子在全模型中与
人为火发生之间并未表现出显著相关(表 1).此外,
本文应用向后逐步回归法,将不显著变量逐一从全
模型中剔除,最后得到最优模型(剔除不显著自变
量后的模型).最优模型整体拟合度指标 R2 为
0 481,模型中距离铁路距离和平均相对湿度仍然表
2012 应 用 生 态 学 报 26卷
表 1 逻辑斯蒂全模型的参数拟合
Table 1 Parameters estimation of Logistic regression full model
自变量因子
Independent variable
估值系数
Estimation
coefficient
标准误
SE
Wald卡方值
Wald chi⁃
square value
P
常量 Constant -5.512 81.895 0.005 0.946
地理及人为 到河流面距离 Distance to river 0.003×10-3 0.009×10-3 0.001 0.973
基础设施因素 到街区距离 Distance to residential area -0.008×10-3 0.011×10-3 0.530 0.467
Geographical and human 到铁路距离 Distance to railway -0.016×10-3 0.005×10-3 8.612 0.003
infrastructure factors 到机耕路距离 Distance to Jigeng Road -0.008×10-3 0.011×10-3 0.533 0.465
到等级公路距离 Distance to Dengji Road 0.006×10-3 0.011×10-3 0.331 0.565
坡度 Slope -0.362 1.030 0.124 0.725
坡向 Aspect -0.182 0.226 0.650 0.420
可燃物因子 腐质层厚度 Depth of humus layer 3.227 1.780 3.285 0.070
Fuel factors 地被物盖度 Coverage of litter -0.057 0.030 3.591 0.058
林型 Forest type -0.771 0.486 2.510 0.113
经营措施 Management measure -0.603 1.008 0.358 0.550
优势树种 Dominant tree -0.231 0.547 0.179 0.673
龄级 Tree age -0.552 0.400 1.905 0.168
郁闭度 Canopy 5.599 3.255 2.959 0.085
气象因子 20:00—20:00降水量 20:00-20:00 precipitation 0.028 0.034 0.662 0.416
Weather factors 日平均气压 Daily average air pressure 0.015×10-2 0.008 0.032×10-2 0.986
日平均风速 Daily average wind speed 0.110 0.060 3.325 0.068
日平均气温 Daily average temperature -0.108 0.042 6.583 0.012
日平均相对湿度 Daily average relative humidity -0.035 0.051 0.467 0.046
日照时数 Sunshine hours 0.039 0.019 3.978 0.494
日最低气温 Daily minimum temperature 0.035 0.021 2.676 0.102
日最高气温 Daily maximum temperature 0.078 0.027 8.547 0.003
日最大风速 Daily maximum wind speed -0.037 0.036 1.048 0.306
日最小相对湿度 Daily minimum relative humidity -0.033 0.057 0.330 0.566
现出与人为火发生有显著相关性,此外这两个变量
的显著性水平较全模型平均明显提高(表 2).
建模数据 ( 70%)和校验数据 ( 30%)的模型
ROC曲线(图 4)下面积(AUC)值分别为 0.847 和
0 914,且显著性水平 P<0.0001(表 3),说明本文建
立的逻辑斯蒂模型具有统计学意义且拟合效果较
好.人为火发生概率(P)模型为:
P= 1 / [1+e-(3.026-0.00011x1-0.047x2)] (3)
式中:x1、x2分别为距离铁路距离和平均相对湿度.
2 1 2 模型预测准确度检验 本文根据约登指数计
表 2 逻辑斯蒂最优模型的参数拟合
Table 2 Parameters estimation of Logistic regression best
model
自变量因子
Independent variable
估值系数
Estimation
coefficient
标准误差
SE
Wald卡方值
Wald chi⁃
square value
P
常量 Constant 3.026 1.050 8.303 0.004
距铁路距离
Distance to railway
-0.0001 0 15.256 <0.0001
日平均相对湿度
Daily average relative
humidity
-0.047 0.017 7.517 0.006
算出最佳临界值为 0.278,以此为基础检验模型的预
测准确率.结果表明,模型建模数据(70%)部分的整
体预测准确率为 79.1%,校验数据(30%)部分的预
测准确率为 80.0%(表 4),说明模型具有较好的预
测能力.
2 2 研究区域林火发生概率分布
根据模型拟合结果,通过 HADCM2 模式下研究
区域 2015 年的月平均相对湿度和平均气温,应用
ArcGIS的图层运算功能计算出塔河地区 2015 年防
火期(3—11月)各月的林火概率分布.由图 5 可以
看出,塔河地区 2015年人为火发生概率较高的月份
为 4—6月和8月,其中,4、5月是林火高发期 .从空
表 3 AUC检验结果
Table 3 Result of AUC test
数据类型
Data type
面积
Area
标准误差
SE
P 95%置信区间
95% confidence interval
下限
Upper
上限
Lower
建模数据
Building data
0.847 0.044 <0.001 0.760 0.934
校验数据
Validation data
0.914 0.045 <0.001 0.825 1
30127期 郭福涛等: 我国北方针叶林人为火发生的预测模型
表 4 基于最佳临界值的林火建模数据(70%)和校验数据(30%)预测准确率
Table 4 Prediction accuracy of building data (70%) and validation data (30%) based on cutoff value
观测值
Observed value
预测值 Predicted value
建模数据 Building data (70%)
火点 Fire
0 1
预测准确率
Percentage correct
校验数据 Validation data (30%)
火点 Fire
0 1
预测准确率
Percentage correct
火点 Fire 0 59 18 76.6 28 7 80.0
1 4 24 85.7 2 8 80.0
总预测准确率 Overall percentage (%) 79.1 80.0
图 4 模型的 ROC拟合曲线
Fig.4 ROC curve of Logistic model.
a)70%拟合数据的计算结果 Result of building data (70%); b)30%
校验数据的计算结果 Result of validation data (30%).
间分布来看,高火险区将集中分布在南部和西南部,
也是铁路主要分布区.
3 讨 论
本文应用塔河地区历年人为火信息数据,结合
研究地气象、基础地理信息及矢量化林相图,应用
ArcGIS 10.0和 SPSS 19.0等工具,通过空间分析、图
层叠加、属性转换、逻辑斯蒂回归模型运算等一系列
手段,分析了对我国北方针叶林带人为火发生有显
著影响的驱动因子.结果显示,距离铁路距离和平均
相对湿度对人为火发生有显著影响.本文根据逻辑
斯蒂模型结果,运用 HADCM2模式下研究区域的未
来月平均气象因子,对大兴安岭塔河地区 2015 年逐
月人为火发生概率分布进行预测,结果显示,2015
年 6—9月为人为火高发期且主要集中在塔河南部
沿铁路沿线.
图 5 研究区域 2015年火险期人为火发生概率
Fig.5 Probability of human⁃caused fire during fire season of 2015 in the study area.
4012 应 用 生 态 学 报 26卷
本文主要研究结论与 Cardille 等[6]在北美针叶
林带所进行的研究结论相似,他们的研究指出,铁路
密度对研究区域人为火有显著影响.Maingi 等[19]在
美国肯塔西州的研究结果表明,到公路距离和公路
密度与人为火发生显著相关,这也说明了人为基础
设施的空间分布对人为火的重要影响.此外,本文关
于平均相对湿度是人为火发生主要驱动因子的结论
也与其他相关研究一致[35] .
由于数据本身的局限性,我们只获得 2003 年的
塔河地区矢量林相图.本研究假设塔河地区在
1974—2009年间主要林型的空间分布保持一致.实
际上塔河地区每年(尤其 1974—1996 年)都会进行
计划采伐,采伐形式多为择伐和抚育间伐,主伐比例
相对较小.择伐和抚育间伐通常不会改变群落优势
树种,而主伐后区域补植多以采伐树种为主,因此本
研究的假设在一定程度上具有合理性.但更多连续
时间梯度上的塔河地区矢量林相图,可以更准确地
反映林型的动态变化,进而提高预测模型的精度,这
也需要研究团队进一步深入研究.此外,由于塔河地
区属于我国北方林区,受林地规划因素影响,林区的
道路、铁路发展相对比较平缓,因此本研究假设林区
铁路、道路在空间分布上保持相对稳定.
本研究主要选取气象因子及一些人为基础设施
为研究对象,与国外同类研究相比[6,19-21],本文没有
考虑当地人口密度、GDP 等经济发展因素,主要由
于我国公开的 GDP 或人口密度等统计数据的最小
尺度单位为县,而本研究区域位于塔河县内, 因此
无法体现出空间异质性, 进而无法同其他因子进行
分析.目前关于林火与影响因子的空间分析研究,主
要有网格取点和随机取点两种方式.网格取点指对
研究区域进行网格化,规定网格大小.每个网格代表
一个火点或对照点(非火点),如果网格面积大(如
5 km×5 km),那么坡向、林型等因子将不具备唯一
性,没有办法考虑进模型拟合中,若网格面积小(如
1 km×1 km 或更小),则会出现对照点与火点比例严
重失衡现象,会导致模型拟合的偏差,从而无法判断
影响因子.因此,在火点和非火点比例失衡的情况
下,随机选取对照点是一个比较好的选择.
参考文献
[1] Podur J, Martell DL, Csillag F. Spatial patterns of light⁃
ning⁃caused forest fires in Ontario, 1976-1998. Ecologi⁃
cal Modelling, 2003, 164: 1-20
[2] Chas⁃Amill ML, Touza J, Prestemon JP. Spatial distri⁃
bution of human⁃caused forest fires in Galicia ( NW
Spain). Second International Conference on Modelling,
Monitoring and Management of Forest Fires. Greece,
2010: 247-258
[3] Lasaponara R, Santulli A, Telesca L. Time⁃clustering
analysis of forest⁃fire sequences in southern Italy.
Chaos, Solitons & Fractals, 2005, 24: 139-149
[4] BC Ministry of Forests and Range Wildfire Management
Branch. Climate Change and Fire Management Research
Strategy. London: BC Ministry of Forests and Range
Wildfire Management Branch, 2009
[5] Zhong MH, Fan WC, Liu TM, et al. Statistical analysis
on current status of China forest fire safety. Fire Safety
Journal, 2003, 38: 257-269
[6] Cardille JA, Ventura SJ, Turner MG. Environmental
and social factors influencing wildfires in the Upper Mid⁃
west, United States. Ecological Applications, 2001, 11:
111-127
[7] Prestemon JP, Butry DT. Time to burn: Modeling wild⁃
land arson as an autoregressive crime function. American
Journal of Agricultural Economics, 2005, 87: 756-770
[8] Todd SK, Jewkes HA. Wildland Fire in Alaska: A His⁃
tory of Organized Fire Suppression and Management in
the Last Frontier. Agricultural and Forestry Experiment
Station Bulletin No. 114. Fairbanks, AK: University of
Alaska Fairbanks, 2006
[9] Sun K⁃H (孙科辉), Xiao J⁃X (肖金香), Feng M⁃Y
(冯敏玉), et al. Analysis on forest fire and fire⁃danger
weather and integrated preventive techniques. Acta Agri⁃
culturae Universitatis Jiangxiensis (江西农业大学学
报), 2004, 2(3): 394-398 (in Chinese)
[10] Zhao H (赵 辉). Spatial and Temporal Distribution of
Forest Fires and Its Impact in Sanming of Fujian Pro⁃
vince from 2000 to 2009. Master Thesis. Beijing: Bei⁃
jing Forestry University, 2012 (in Chinese)
[11] Guo F⁃T (郭福涛). Forest Fire Distribution and Fore⁃
casting Modeling Basis on Spatial Analysis Modeling
Theory in Daxing’ an Mountains. PhD Thesis. Harbin:
Northeast Forestry University, 2010 (in Chinese)
[12] Jin S (金 森), Hu H⁃Q (胡海清). Study on forest
fire regime of Heilongjiang Province.Ⅰ. Forest fire spa⁃
tial and temporal dynamics and statistical distribution.
Scientia Silvae Sinicae (林业科学), 2002, 38 (1):
88-94 (in Chinese)
[13] Wang M⁃Y (王明玉), Shu L⁃F (舒立福), Tian X⁃R
(田晓瑞), et al. Spatial fluctuation of forest fires and
their response to global change. Fire Safety Science (火
灾科学), 2003, 12(3): 165-176 (in Chinese)
[14] Zhao F⁃J (赵凤君), Wang M⁃Y (王明玉), Shu L⁃F
(舒立福), et al. Progress in studies on influences of
climate change on forest fire regime. Advances in Climate
Change Research (气候变化研究进展), 2009, 5(1):
50-55 (in Chinese)
[15] Wei S⁃J (魏书精), Hu H⁃Q (胡海清), Sun L (孙
龙). The effects of climate change on the laws of forest
fire occurrence in China. Forest Fire Prevention (森林防
火), 2011(1): 30-34 (in Chinese)
[16] Sun Y⁃R (孙玉荣), Zhang G (张 贵). Characteris⁃
50127期 郭福涛等: 我国北方针叶林人为火发生的预测模型
tics of the spatial⁃temporal distribution for artificial fire
sources on forest fires in Hunan. Chinese Agricultural
Science Bulletin (中国农学通报), 2011, 27 ( 28):
103-108 (in Chinese)
[17] Zhang Z⁃X (张正祥), Zhang H⁃Y (张洪岩), Li D⁃X
(李冬雪), et al. Spatial distribution pattern of human⁃
caused fires in Hulunbeir grassland. Acta Ecologica Sini⁃
ca (生态学报), 2013, 33(7): 2023-2031 ( in Chi⁃
nese)
[18] Syphard AD, Radeloff VC, Keeley JE, et al. Human in⁃
fluence on California fire regimes. Ecological Applica⁃
tions, 2007, 17: 1388-1402
[19] Maingi JK, Henry MC. Factors influencing wildfire
occurrence and distribution in eastern Kentucky, USA.
International Journal of Wildland Fire, 2007, 16: 23-
33
[20] Romero⁃Calcerrada R, Novillo CJ, Millington JDA, et
al. GIS analysis of spatial patterns of human⁃caused
wildfire ignition risk in the SW of Madrid ( central
Spain). Landscape Ecology, 2008, 23: 341-354
[21] Martínez J, Vega⁃Garcia C, Chuvieco E. Human⁃caused
wildfire risk rating for prevention planning in Spain.
Journal of Environmental Management, 2009, 90:
1241-1252
[22] Padilla M, Vega⁃Garcia C. On the comparative impor⁃
tance of fire danger rating indices and their interaction
with spatial and temporal fire occurrences in Spain. In⁃
ternational Journal of Wildland Fire, 2011, 20: 46-58
[23] Yu C⁃L (于成龙), Hu H⁃Q (胡海清), Wei R⁃H (魏
荣华). Dynamic analysis of meteorological conditions of
forest fire in Tahe forestry bureau of Daxing’an
Mountains. Journal of Northeast Forestry University (东
北林业大学学报), 2007, 35(8): 23- 27 ( in Chi⁃
nese)
[24] Guo F⁃T (郭福涛), Hu H⁃Q (胡海清), Ma Z⁃H (马
志海), et al. Applicability of different models in simu⁃
lating the relationships between forest fire occurrence
and weather factors in Daxing’ an Mountains. Chinese
Journal of Applied Ecology (应用生态学报), 2010, 21
(1): 159-164 (in Chinese)
[25] Johns TC, Carnell RE, Gregory JM, et al. The second
Hadley Centre coupled ocean⁃atmosphere GCM: Model
description, spinup and validation. Climate Dynamics,
1997, 13: 103-134
[26] Xu Y⁃L (许吟隆). Analyses on scenario simulations of
the 21st century climate change in China. Journal of
Nanjing Institute of Meteorology (南京气象学院学报),
2005, 28(3): 323-329 (in Chinese)
[27] Martell DL, Bevilacqua E, Stocks BJ. Modelling seasonal
variation in daily people⁃caused forest fire occurrence.
Canadian Journal of Forest Research, 1989, 19: 1555-
1563
[28] de Vasconcelos MJP, Sllva S, Tome M, et al. Spatial
prediction of fire ignition probabilities: Comparing logis⁃
tic regression and neural networks. Photogrammetric En⁃
gineering & Remote Sensing, 2001, 67: 73-81
[29] Garcia CV, Woodard PM, Titus SJ, et al. A Logit
model for predicting the daily occurrence of human
caused forest fires. International Journal of Wildland
Fire, 1995, 5: 101-111
[30] Andrews PL, Loftsgaarden DO, Bradshaw LS. Evalua⁃
tion of fire danger rating indexes using logistic regression
and percentile analysis. International Journal of Wild⁃
land Fire, 2003, 12: 213-226
[31] Chang Y, Zhu ZL, Bu RC, et al. Predicting fire occur⁃
rence patterns with logistic regression in Heilongjiang
Province, China. Landscape Ecology, 2013, 28: 1989-
2004
[32] del Hoyo VL, Isabel MPM, Vega FJM. Logistic regre⁃
ssion models for human⁃caused wildfire risk estimation:
Analyzing the effect of the spatial accuracy in fire occur⁃
rence data. European Journal of Forest Research, 2011,
130: 983-996
[33] Jiménez⁃Valverde A. Insights into the area under the re⁃
ceiver operating characteristic curve ( AUC) as a dis⁃
crimination measure in species distribution modelling.
Global Ecology and Biogeography, 2012, 21: 498-507
[34] Fielding AH, Bell JF. A review of methods for the as⁃
sessment of prediction errors in conservation presence /
absence models. Environmental Conservation, 1997, 24:
38-49
[35] Wotton BM, Martell DL, Logan KA. Climate change and
people⁃caused forest fire occurrence in Ontario. Climatic
Change, 2003, 60: 275-295
作者简介 郭福涛,男,1982 年生,博士. 主要从事林火生
态、林火预测预报研究. E⁃mail: guofutao@ 126.com
责任编辑 杨 弘
6012 应 用 生 态 学 报 26卷