森林生态系统具有涵养水源、调节气候、保持水土等生态功能, 在维护生态平衡中具有重要作用。草地生态系统是我国陆地上面积最大的生态系统, 在发展畜牧业、保护生物多样性、固持水土等方面占有重要地位^[1]。森林火灾、草原火灾对生态系统破坏力强, 是最普遍的干扰方式之一。2019年年初至同年8月, 巴西境内森林火灾发生次数超过7万起^[2], 多发生在雨林边缘, 湿度相对较低且人类活动较为频繁的地区。2019年9月, 澳大利亚山火肆虐^[3], 至2020年7月底, 澳大利亚丛林大火导致30亿动物死亡, 经济损失更加无法估量。2020年3月30日, 我国四川凉山西昌发生森林大火, 因空气湿度极低, 气温较高, 火势蔓延迅猛, 造成18名消防员和1名林场职工牺牲的惨剧。2018年数据显示, 全球草原面积达32.73亿公顷, 草原所提供的生产力不容小觑, 草原火灾将威胁牧民的生命财产安全以及牲畜生存, 影响畜牧业发展, 带来巨大的经济损失。

可燃物是火灾发生的物质基础和重要载体, 可燃物含水率的高低程度决定了火的发生和蔓延^[4]。可燃物与其所处的环境进行水汽交换是影响可燃物含水率变化的主要因素, 含水率的多少与其所处的自然环境的干湿程度和气温度高低相关^[5-6]。森林、草原火灾危害巨大, 探究可燃物含水率能够发掘林火/草原火点燃的难易程度和蔓延速率, 预测林火、草原火行为^[7], 防患于未然, 降低火灾发生的可能性, 减少经济损失, 为生产生活提供便利。当前国内外有关可燃物含水率的研究已基本成熟, Ertugrul^[8]等模拟了正常储存的卡拉布里亚松(Pinus brutia Ten)地表可燃物的水分含量与土耳其穆格拉省的天气条件(即温度、相对湿度和风速)的关系。Lee等^[9]比较了基于3种10 h可燃物含水率预测模型的性能。胡海清^[7]等分析了防火期和坡向对大兴安岭地区南瓮河混交林地表可燃物含水率预测模型精度的影响。于宏洲^[10]等应用时滞平衡含水率法及气象要素回归法, 分析了大兴安岭地区阳坡上部落叶松林可燃物含水率的动态变化。国内外可燃物含水率研究大多探究含水率的影响因素和预测评估模型, 亦有关于可燃物含水率研究的综述性文章, 但是, 这些综述一般是基于对已有研究的总结, 鲜有利用文献计量的方法理清研究领域研究进展、发展脉络的文章。

文献计量是一种常见的定量分析方法, 它基于数学、统计学原理, 以可获取文献的诸多特征为研究对象, 探究其分布结构、数量关系和变化规律等^[11], 通过知识图谱将文献计量可视化结果展示出来^[12], 并以此为依据分析有关研究的发展特征和规律, 是一种成熟的文献分析和信息挖掘方法^[13]。CiteSpace文献计量分析软件由美国德雷赛尔大学陈超美博士设计开发^[14-15], 用于文献信息的计量与可视化, 在诸多领域中被广泛使用^[16-17]。因此, 本文以CiteSpace 5.7.R2为辅助手段, 探究当前国际上关于森林和草原可燃物含水率的研究进展, 旨在理清当前相关研究的主要方向, 寻求新的突破, 从而促进森林和草原可燃物含水率研究的未来发展, 更好地预测火行为, 增进人类福祉。

1 数据来源与研究方法 1.1 数据来源

研究基于WEB OF SCIENCE数据库核心合集, 检索主题中同时包含“fuel moisture”和“forest”, 以及同时包含“fuel moisture”和“grassland”的文献, 输入检索式TS=((fuel moisture) AND (forest)), 记为#1, TS=((fuel moisture) AND (grassland)), 记为#2, 在WEB OF SCIENCE中组配两个检索式检索到的结果, 记为#1OR#2, 由于web of science核心合集所收录文献最早为1980年, 因此检索的时间范围是1980至2020年10月, 剔除无效、重复的文献, 共计检索到文章1298篇。

1.2 研究方法

可视化软件CiteSpace能够清晰准确地解读某一知识领域的发展趋势、研究前沿、研究热点等信息。目前比较常见的可视化有CiteSpace, Bibexcel, Sati, VOSviewer等, CiteSpace因其可视化效果好, 功能全面, 布局合理, 在众学科领域中被学者们广泛使用^[18-19]。以此1298篇有关森林/草原可燃物含水率研究的文章为研究对象, 导出文章题录信息, 再基于web of science内的统计数据结合Excel对文章数量进行统计, 基于Java环境下的CiteSpace 5.7.R2进行学科领域分析、研究力量分析、对文章关键词根据时间序列进行分析, 以及根据引文的聚类分析和突现分析。

2 结果与分析 2.1 文献基本特征分析 2.1.1 文献数量变化规律分析

根据检索到的1298篇有关森林/草原可燃物含水率文章的出版年份绘制出发文数量变化图(图 1)。根据图 1的发文数量, 可以将研究的发展历程分为3个阶段: 1980—2000年为缓慢发展阶段, 年发文数量少, 增长缓慢, 增长率浮动较大, 20世纪末, 相关研究大多处于起步阶段, 技术手段落后, 国际合作交流较少, 导致文献数量较少。2001—2010年为快速发展阶段, 这一时期增长率趋于平缓, 除个别年份小幅下降外, 大部分年份处于增长态势, 21世纪以来网络通讯发展完善, 信息数据全球共享, 全球卫星遥感事业取得长足发展, 这些都为研究提供了极大的支持。2011至今为波动增长阶段, 这一时期发文数量整体上波动较大, 2016年和2019年为负增长率, 由于统计时间在2020年中叶, 因此2020文献数量不具有代表性。2010年7月俄罗斯多地发生森林大火, 同年9月美国科罗多州发生森林大火, 2013年美国亚利桑那州森林大火, 因此推测特大火灾的发生是导致次年研究数量大量增加的原因。

2.1.2 学科领域

在WEB OF SCIENCE数据库核心合集中, 每篇文章会涉及一个或多个学科类别的学科标签加以显示, 得出学科领域的分类图(图 2), 图中带有玫粉色圈层的节点代表了该学科突现性。由图 2可知, 在森林/草原可燃物含水率的研究中, Forestry(林学), Environmental science and ecology(环境科学与生态学), Environmental science(环境科学), Ecology(生态学), Remote sensing(遥感), Engineering(工程学), Meteorology and atmospheric science(气象和大气科学), Energy and Fuels(能源与燃料), Geosciences multidisciplinary(地理科学多学科), Geology(地质学), Agriculture(农学)等学科是其涉及得的主要研究领域。其中林学所占的数量最多, 说明对于森林可燃物含水率的研究居多, 这是由于森林火灾强大的破坏性和频发性使这类研究数量居多。在分析过程中得知, 环境科学这一学科的中心性最高, 达到0.43, 说明研究森林/草原可燃物含水率这一研究方向主要致力于环境研究, 为提高人类赖以生存的地球环境而服务。这一方向是一个综合了林学、生态学、环境科学、地理科学、遥感等学科的多学科交叉领域, 与地理科学和气象大气科学的学科交叉, 体现对含水率的研究不只是考虑小范围的区域生态而是更具有空间维度, 与能源燃料和工程学的交叉体现了研究向实际应用的迈进。

2.1.3 研究力量分析

利用CiteSpace 5.7.R2, 对文献来源国家(Country)进行分析, 将发文数量大于30的国家标签加以显示, 共13个国家, 得出国际合作关系图(图 3), 得出13个国家科研实力信息表(表 1)。

结合表 1和图 3可知, (1)在有关森林/草原可燃物含水率相关研究中, 美国的研究实力在世界各国中居于领先地位, 相关文献中由美国作者发表或合作发表的数量达到496次, 文献首次发表年份最早, 远超居第二位的澳大利亚180次, h-指数和被引次数也远超各国, 说明美国科研实力最高, 其中心性也最高, 代表美国与世界各国有着较为紧密的合作关系。(2)我国相关研究起步较晚, 国际文献最早出现在2004年, 按照文献数量, 我国位居世界第五, 但我国的h-指数和被引次数均低于文献数量少于我国的两个国家, 这一现象表明我国在这一领域的科研实力有待加强, 在提高研究数量的同时更要保证研究质量。(3)位居前3位的国家都是国土面积较大且森林面积广阔的发达国家, 此外, 在前13位的国家中, 只有中国和巴西两个发展中国家, 这与两个国家拥有大面积森林和草原的自然地理特征有关, 由此说明, 森林/草原可燃物含水率研究是各国专家学者为解决实际问题而产生, 有较强的地理相关性。

2.2 研究进展分析 2.2.1 关键词分析

关键词是每一篇学术论文的组成部分, 是对文章的高度概括和凝练^[20]。同时, 关键词又代表着不同研究领域的研究焦点^[21]。因此, 透过关键词往往能够理清某一研究领域的发展脉络与趋势, 进而整体把握该领域的研究现状, 有助于进一步展开深入研究。本研究提取了检索到的1298篇森林/草原可燃物含水率研究的关键词, 导入CiteSpace 5.7.R2中, 文献筛选范围从1990—2020年, 设置两年为一个时间切片, 将出现频次大于30的关键词标签加以显示, 对共现情况进行可视化, 并选择时区视图加以呈现如图 4。

如图 4, 图中节点和线条不同的颜色代表不同的时间切片, 颜色越深, 年代越久。每个圆形节点代表一个关键词, 节点越大代表其出现频次越多, 节点所在位置代表其最早出现的时间, 关键词节点年轮的颜色由中心向四周逐渐变浅, 某一层越厚, 说明这一关键词在该颜色所代表的时间切片中出现次数越多；节点间的连线代表两个关键词共同出现在同一篇文献中的事件, 线条颜色代表两个关键词第一次出现在同一篇文章中的年份, 线条越粗, 共现次数越多^[12]。

结合图 4按照时间脉络, 可以归纳出森林/草原可燃物含水率研究的大致发展特点。(1)20世纪90年代的关键词主要是forest(森林)、fire(火)、climate change(气候变化)、biomass(生物量)、vegetation(植被)等, 此时的研究大多是探究森林地区的可燃物含水率与气象条件如温度和湿度、气候变化、区域生物量等一些变量的内在联系, 或火灾现象对生态系统造成的影响。California(加利福尼亚)作为关键词出现, 说明当时开始出现研究以州/省域为研究区域。Fire regime(管理体制)一词在20世纪90年代末出现, 这表明可燃物含水率研究不再停留于理论层面, 开始面向实际应用。(2)21世纪开始出现model(模型)、modi(程序)、index(指数)、remote sensing(遥感)、fire severity(火烈度)、fire danger(火险)等词, 说明这一时期伴随着计算机的普及、全球卫星事业的发展, 研究工具得以改进, 人们对于可燃物含水率的估测以及火灾的量化评估开始稳步推进。(3)2008年左右出现了United States(美国), 标志着研究区域已经从州/省域扩大到国家层面。(4)2010年以后logging residue(采伐残渣)一词的出现, 体现了森林可燃物含水率的研究中考虑了更加多样的可燃物种类, 将人类活动影响加以考虑, 研究更加客观、精确、符合实地情况。

2.2.2 聚类分析

每一篇文献都会引用数篇同领域的, 或与文章内容相近的文献, 频繁的共同引用表明它们具有共同相关的研究主题, 因此探究文章的共被引情况, 共引分析可以根据内容的相似程度将相关参考文献分组聚类, 从而分析这一领域的研究热点^[22]。Citespace的共被引分析是基于检索到的文章及其引用文章进行分析, 得出文献共被引网络可视图之后, 对所呈现的点进行聚类, 从施引文献的标题、关键词或摘要中提取聚类命名术语。本文对有关森林/草原可燃物含水率的1298篇文章进行聚类, 两年为一个切片, 得出13个群集, 群集名称选取LLR(Log-likelihood Ratio)算法计算得到, 采用时间线视图(Timeline)加以显示, 每个节点代表一篇文献, 节点间连线代表其共同被引用的关系, 得出引文聚类情况网络图(图 5)。

聚类后得出Q值(Modularity聚类模块值)和S值(Silhouette)聚类平均轮廓值, 一般认为Q＞0.3时, 聚类结构显著, S＞0.5聚类是合理的, S＞0.7, 聚类是可信的。本次聚类中Q值为0.8489, 代表聚类结构显著；S值为0.9318, 说明聚类可信度高。群集标签从0—5依次为#0 radiative transfer modeling(辐射传输建模), #1 assessing crown fire potential(树冠火潜力评估), #2 shoot flammability(枝条可燃性), #3 western USA(美国西部), #4 fire science(火灾科学), #5 predicting temperate ecosystem sensitivity(生态系统敏感性预测), #6 inter-hemispheric synchrony(两半球同步), #7 fire risk assessment(火灾风险评估), #8 dead fuel(死可燃物), #10 wildfire consumption(野火消耗), #11 forest floor duff(森林半腐层), #12 beetle attack(甲虫攻击), #13 dominant forest cover(主要森林覆盖)。群集序号越小, 所包含的文献数量越多。由于群集标签是从影响力最强的文章题目中提取, 导致标签不能准确地概括群集中的所有文章内容, 因此本文选取了前5个群集对文章进行具体分析, 剩余群集中包含的文献由于数量少, 代表性不强, 因此不做具体分析。

群集#0所包含的文章主题相似, 主要是以遥感技术为依托, 评估叶片含水量, 从而对森林火灾进行监测^[23]。Emilio^[24]基于AVHRR数据的多时相分析, 提出了一种提取地中海地区草地和灌木物种燃料含水率的方法。Kötz^[25]通过对针叶林上的成像光谱仪数据的采集, 获取植被冠层结构和叶面含水量的空间分布信息, 用于评估森林火灾对树木的影响。这一群集的文献平均发表年份在2002年, 由于20世纪90年代国际遥感技术发展成熟, 影像分辨率得到较大提升, 使遥感影像在各个领域中广泛使用, 导致这一时期通过遥感影像评估植被冠层含水率的方法模型的相关研究大量涌向。

群集#1中的文章以评估形式为主^[26], 探究森林火灾的发生与不同变量之间的联系, 其中对于气候变化与火灾间的关系研究较为常见。如Flannigan^[27]分析了加拿大火灾天气指数(FWI)与加拿大生态区域历史燃烧面积之间关系, 并结合加拿大和哈德利中心GCMs的产品, 预测未来的燃烧面积。Cochrane^[28]探究了亚马逊地区火灾与土地利用、气候变化间的协同作用。Dymond^[29]提出了一种基于卫星火灾探测(hotspot)数据的火险分级系统校准的方法, 探究西印度尼西亚和马来西亚森林火灾问题与加拿大森林火灾天气指标体系的两个组成部分之间的关系。这部分文章多以辩证思维展开探究, 对已有的方法展开评估或验证其准确程度。

群集#2包含的文章以植被或其凋落物的可燃性为主要研究对象^[30]。研究可燃物含水率的目的是为了探究可燃物燃烧的可能性, 从而对燃烧行为加以利用或规避, 因此可以说, 植被或凋落物的可燃性是森林/草原可燃物含水率的研究的核心部分。专家学者以不同的视角、方法分析不同植被的可燃性。Essaghi^[31]以可燃物含水率的两个相关参数——点火时间和火焰高度作为可燃性的指数, 对摩洛哥森林燃料的可燃性进行评估。Alam^[32]分析了新西兰的43中本地多年生植物的叶片和茎部的可燃性数据, 解耦了茎部可燃性和叶片可燃性间的联系。Cawson^[33]分析了潮湿气候条件下桉树林可燃性的驱动因素, 结合了从21位火灾专家技术中获取的启发, 建立了湿润桉树林区可燃性的概念模型。

群集#3软件生成的标签为western USA, 主要是由于这一群集的文章中大部分以美国西部为研究区, 通过上文分析, 美国在这一方向的研究中一直处于遥遥领先的地位, 因此美国学者也往往会选择本土作为自己的研究区。美国西部地形高大的山地、高原为主, 科迪勒拉山系高度大、宽度广, 包括海岸山脉、内华达山脉和落基山脉, 森林覆盖面积大, 物种多样, 生态环境复杂, 是森林火灾的频发地区。几十年来美国西部野火的频繁出现, 人们逐渐发现火灾与气候之间的联系, Abatzoglou^[34]探究了美国西部八个地理区域协调中心的地区气候与大型火灾面积的关系。Gabrielle^[35]探究了科学地管理野火对内华达山脉植被变化的影响、植被变化对当地水文的影响以及植被变化对森林火灾燃烧特征和生态系统弹性的影响。

群集#4所包含的文章大多以火灾科学为主题, 或最终目的致力于火灾的防控。这些文章发表时间较晚, 平均年份为2016年, 是在可燃物含水率研究发展到一定时期必然兴起的产物, 经过了早期的对可燃物的可燃性研究, 以及对燃烧过后的生态恢复研究, 消防任务显得尤为重要。因此, Collins^[36]使用陆地卫星图像生成准确的火灾严重程度地图, 这对于记录火灾状态、为火灾后管理反应以及在研究应用设定优先级具有重要意义。Chuvieco^[37]总结了卫星遥感在消防科学与火灾管理上的贡献。全球气候变暖、人类活动增强的大背景下, 森林/草原消防工作的重要性日益凸显, 消防工作需要科学化系统化的管理才能将不必要的火灾防患于未然, 一旦火灾出现, 可以以最迅速地扑灭, 将损失降到最低, 因此有关消防科学的研究应运而生。

综上, 森林/草原可燃物含水率相关研究在研究方法上以遥感反演为主流研究方法, 气候变化是最常见的研究视角, 植被或枯落物的可燃性是研究的核心内容, 美国西部是研究最多的区域, 消防科学是最主要的应用领域。

2.3 研究前沿分析

为了了解当前森林/草原可燃物含水率的研究前沿, 本文利用CiteSpace 5.7.R2分析文献的突现性, 突现性是指某一篇文献在某一时期被引次数陡然增加, 被引次数越多, 说明在这一段时期里相关研究大量出现, 因此近期的突现性一定程度上可以代表研究前沿。本文筛选出突现性最强的25篇文献, 并选取其中突现阶段距今最近的7篇文献, 如表 2(每个条块代表一年, 较粗条块代表突现时期), 将其作为当前的研究前沿。

结合表 2以及相关文献总结出当前的研究前沿为以下3个方面:

1) 大尺度的全球火险预测研究。全球气候变化的大环境推进大尺度的全球火险预测研究, 大尺度全球火险预测研究反之将丰富全球气候变化对人类影响的研究, 对全球范围内人类共同应对气候变化提出新看法, 推动全球可持续发展战略部署。相比于遥感影像等地面数据, 气象数据往往更容易获取, 因此大尺度研究通常以大量的气象数据加以分析佐证, 结合气候模型进行分析^{[40, 43]}。

2) 可燃物可燃性的多维度分析。可燃性是可燃物含水率研究的核心内容, 早期研究大多集中在可燃物含水率预测模型的构建, 或探究其与气象因子的关系, 这些研究大多分析可燃物时空间变化, 进而实现火险的预测目的。进一步剖析可燃性的属性, 可将其划分为可燃性、热释放和火灾蔓延率3个维度^[42], 用于分析燃烧产生的环境效应以及从空间上探究火灾蔓延态势, 便于人们掌握不同尺度下可燃性的可变性, 从而为进一步研究提供基础。

3) 基于MODIS反演可燃物含水率。随着遥感卫星技术的成熟发展, 使得遥感影像愈加容易获取。MODIS数据可免费获取, 光谱范围广, 更新频率高, 在实时地球观测和应急处理方面(如森林和草原火灾监测和救灾)有较大的实用价值。因此利用MODIS数据反演地表可燃物含水率^{[39, 44]}在近期成为研究热点, 学者们的数据获取不再只局限于实地测量和利用气象数据进行模拟, 而是利用MODIS数据反演后利用实测数据进行验证。遥感反演有着比气象模拟更准确, 比实地测量覆盖面积更广、更加省时省力的优势。

3 结论

本文通过对森林/草原可燃物含水率相关文献进行计量分析, 分析研究数量的增长态势、所涉及的研究领域、主要研究力量、研究进展与前沿, 为未来学者对森林/草原可燃物含水率的深入研究提供借鉴。

从研究的基本特征上看, 相关研究在可以分成3个研究阶段, 缓慢发展阶段、快速发展阶段和波动增长阶段；这一方向综合了林学、生态学、环境科学、地理科学、遥感等学科具有一定的空间维度和实际应用性；美国的研究实力在世界各国中居于领先地位, 我国相关研究起步较晚, 且我国在未来的研究中要注重提高研究质量, 可燃物含水率研究数量较多的国家往往有较大的森林面积。

从研究进展上看, 伴随着科学技术的提高, 信息获取难度大大降低, 研究区面积逐渐扩大, 由州域扩大到全国, 进而扩大到全球范围；关于研究方法, 20世纪90年代相关研究主要以探究可燃物含水率与气象因子或其他环境条件之间的关系为主, 实际应用研究也于90年代末期开始出现。21世纪随着科技进步, 研究工具得以改进, 依赖于遥感影像、数学建模等方法大量出现, 对于可燃物含水率的估测的准确性得以提升, 火灾的量化评估开始稳步推进。整体上看, 森林/草原可燃物含水率相关研究在研究方法上以遥感反演为主流研究方法, 气候变化是最常见的研究视角, 植被或枯落物的可燃性是研究的核心内容, 美国西部是研究最多的区域, 消防科学是最主要的应用领域。

总结出当前的研究前沿主要有3个方面, 大尺度的全球火险预测研究、可燃物可燃性的多维度分析和基于MODIS反演可燃物含水率。与早期研究相比, 当前的前沿在研究尺度上更广阔, 在研究视角上更深入透彻, 在研究数据上的精度更高, 更准确。因此, 未来我国在森林/草原可燃物含水率研究中, 应注重对深度和广度的挖掘, 利用遥感空间信息, 增加对森林/草原可燃物含水率空间异质性的研究, 形成对整体区域上的火险预测, 注重研究的实用性。