第一部分 研究背景及现状
全球气候变暖背景下,特大森林草原火灾等极端气候事件频发,给全球生态环境保护和百姓生命财产安全带来了严峻挑战。大面积的野火[自然的(natural)或非计划(unprescribed)的火[1,2]]会演变为森林草原火灾,破坏生态环境,威胁当地百姓生命财产安全[3]。森林草原火灾是一种突发性强、破坏性大、处置救助较为困难的自然灾害。森林草原防火和救援工作是中国防灾减灾工作的重要组成部分,是国家公共应急体系建设的重要内容,事关社会发展稳定的大局。因此,高精度地定量预估野火发生的潜在可能性(即野火风险预警[4]),将森林草原火灾消灭在萌芽状态,显得尤为迫切和重要。
针对野火风险预警,Pyne等提出了著名的野火环境三角(fire environment triangle)模型,指出森林野火的爆发与气象、地形和可燃物(fuel)三个方面密切相关[5]。目前,国内外相关研究[6-9]和业务化运行的森林火险系统,如中国森林火险气象预报系统(QX/T 77-2007)、加拿大森林林火天气指数系统(Canadian Fire Weather Index System)、美国国家火险等级系统(National Fire Danger Rating System)、澳大利亚森林火险等级系统(Australian McArthur Forest Fire Danger Rating System)等,虽然充分考虑了温度、相对湿度、降雨、风速以及坡度等气象和地形信息,但在可燃物信息方面仅考虑了部分死可燃物信息或绿度因子(greenness factor)[10-13]。其主要原因是气象和地形信息相对容易获取,而大范围、近实时的可燃物信息仍难以获取,导致森林野火风险预警仍存在理论缺陷。
鉴于此,作者团队十余年来一直致力于野火风险预警监测方面的研究,通过多学科交叉研究,构建了基于可燃物遥感定量反演和时空大数据挖掘技术的野火风险预警理论与方法体系,并研发了森林草原火灾预警监测系统,服务于我国森林草原火灾风险预警监测与防控。
第1章 森林草原可燃物含水率遥感反演方法研究进展
经典的野火环境三角模型中,气候和地形相关数据的获取比较容易,因此,准确地描述可燃物的空间分布特征和时间变化规律便成了基于野火环境三角模型评估和预警野火风险的关键。可燃物具体是指野火发生时所有可以被燃烧物质的总和。对于自然植被区域而言,占比*大的可燃物便是地表覆盖的植被(包括植被正常生长发育所产生的枯枝落叶等凋落物)。根据其组成部分的生物学特征,植被可燃物可具体划分为4部分:冠层可燃物、树枝可燃物、枝干可燃物以及地表枯枝落叶等死可燃物。表征植被冠层可燃物特征的参数包括冠层可燃物含水率(canopy fuel moisture content,CFMC)、冠层可燃物载量以及可燃物类别等。其中,植被冠层可燃物载量是指单位地表面积中所有植被冠层叶片的干物质重量,是冠层火发生和燃烧的物质基础,与植被自身特点和生长状态密切相关。可燃物类别是指一个可识别的可燃物要素的组合,主要与植被物种属性相关。以上两种可燃物信息的时空分布均比较稳定,而植被冠层可燃物含水率定义为冠层叶片含水量(canopy water content,CWC)和冠层叶片干物质重量(dry matter content,DMC)的百分比值[式(1-1)],比较容易受到气候因素(如空气温度和湿度、风速、降水等)的影响,在一个很小的空间或时间尺度上都容易发生较为剧烈的起伏变化,导致植被冠层可燃物含水率的时空监测结果存在较大的不确定性[14-16],从而直接影响*终野火风险的评估与预警精度。因此,针对野火风险评估与预警的重点就是研究如何高精度地描述植被冠层可燃物含水率在时空尺度上的变化分布特点。
传统的植被冠层可燃物含水率监测方案以人工地面采样为主,具体的测量流程如下:首先选定地表植被空间异质性较低的采样区域;然后在选定的地面采样区域内通过破坏性采样方式采集所有植被冠层叶片并现场称重(记为植被湿重),将样本带回实验室利用烘箱等设备烘干后再称重(记为植被干重);*后根据植被冠层可燃物含水率计算[式(1-1)]植被叶片含水量和植被叶片干重的百分比值,即为植被冠层可燃物含水率(CFMC)。
(1-1)
这种基于地面调查的方法可以获取到比较准确的植被冠层可燃物含水率参数数据,但相应地需要大量的地面实测操作以及复杂的后续处理流程,实施过程不但要消耗大量的人力物力,且无法实现及时、宏观和动态的植被冠层可燃物含水率监测[17]。因此,这种方法只在早期的野火防治工作中使用,目前除在采集数量较少的地面验证数据过程中使用外,已基本不用。
遥感技术的兴起与发展使得传统的地面调查方法转换到了空间尺度上的定量描述。相对于传统的地面测量方法而言,遥感技术可以重复性地对地观测,从而获取到大面积、高空间分辨率的多时相对地观测信息,因此越来越多的研究人员尝试利用遥感技术进行植被冠层可燃物含水率参数大范围、多时相的高精度监测。按照传感器的探测波段,遥感技术可具体划分为可见光-近红外遥感、红外遥感以及微波遥感3类。由于植被冠层可燃物含水率参数在近红外波段和短波红外波段具有一定的光谱灵敏性[14,17,18],所以目前国内外已经进行了大量基于可见光-近红外遥感和红外遥感的植被冠层可燃物含水率监测研究。但是可见光-近红外遥感和红外遥感光谱信号的穿透能力较弱,极易受云雨雾等天气条件和日照条件的影响[19],无法精确监测植被冠层可燃物含水率数据的时空连续变化特点。与此相比,微波信号具有比较强的穿透能力,且不受光照条件的限制,可以全天时全天候地对地观测。同时,地表植被和土壤层对于微波信号的吸收和散射作用主要由介电常数决定,而介电常数又主要受水分含量信息的影响,所以微波信号对于植被水分信息比较敏感[19-22]。因此,微波遥感在植被冠层可燃物含水率监测方面同样具有极高的潜在应用价值。
1.1 基于光学-红外遥感的植被冠层可燃物含水率反演
基于可见光-近红外遥感或者红外遥感的植被冠层可燃物含水率反演方法主要可以划分为两大类:传统的基于经验统计模型的定量估算方法和基于植被辐射传输模型的定量反演方法。
传统的经验统计方法利用多元线性回归等经验统计模型直接对植被冠层可燃物含水率进行遥感估算,目的是寻找植被冠层可燃物含水率参数与多光谱反射率数据或相关典型植被指数[如叶面积指数(leaf area index,LAI)、归一化差值植被指数(normalized difference vegetation index,NDVI)、归一化红外指数(normalized difference infrared inelex,NDII)以及增强型植被指数(enhanced vegetation index,EVI)等]之间的线性(或非线性)关系。20世纪80~90年代,Chladil和Paltridge等[23,24]首次基于NOAA-AVHRR传感器获取的多光谱反射率数据发现了归一化差值植被指数(NDVI)与植被冠层可燃物含水率(CFMC)参数的正线性相关关系。2012年,Caccamo等[25]基于中等分辨率成像光谱仪数据(moderate resolution imaging spectroradiometer,MODIS),对比了4种光谱指数与植被冠层可燃物含水率的相关性差异,发现对于澳大利亚东南部3种容易发生火灾的地表植被而言,基于归一化红外指数(NDII)所建立的经验统计模型的植被冠层可燃物含水率遥感估算精度*好。2014年,Casas等[17]基于机载可见光/红外成像光谱仪AVIRIS和MODIS多光谱数据计算多种典型植被指数,采用经验回归的方式定量估算了包含植被冠层可燃物含水率在内的4种典型植被参数。传统的经验统计方法虽然可以估算出比较准确的植被冠层可燃物含水率,但基于地面实测数据建立的经验统计模型容易受地表植被类型、下垫面土壤质地特点以及传感器性能差异等因素的影响,具有较差的普适性[16],一般只适用于基于特定传感器数据和特定研究区域的植被冠层可燃物含水率遥感估算。
基于植被辐射传输模型的植被冠层可燃物含水率反演方法由于具有较好的普适性而受到野火风险评估专家的推崇。此类方法通过同时反演PROSPECT[26]叶片光学模型中叶片等效水分厚度(equivalent water thickness,EWT)及叶片干物质重量(DMC)两个参数,通过相除运算[式(1-2)]得到植被冠层可燃物含水率(CFMC)[18,27]。
(1-2)
目前常用的定量反演植被冠层可燃物含水率的植被辐射传输模型主要包括PROSAILH模型[28](SAILH模型[29]+PROSEPCT模型[26])、PROACRM模型[30](ACRM模型[30]+PROSEPCT模型[26])以及PROGEOSAILH模型(SAILH模型[29]+Jasinski模型[31,32]+PROSPECT模型[26])3种。近年来,以澳大利亚国立大学野火风险研究专家Marta Yebra为代表的众多研究人员都在尝试基于植被辐射传输模型遥感定量反演植被冠层可燃物含水率,并取得了一定的研究成果[16,18,32-34]。2004年,K?tz等[35]分别利用FLIGHT和PROGEOSAILH植被辐射传输模型对比反演了森林冠层可燃物含水率数据。2009年,Yebra等[32]利用PROSPECT叶片光学模型耦合SAILH冠层辐射传输模型,同时结合长期观测所获得的丰富的先验知识,成功定量反演了地中海区域草地及灌木的冠层可燃物含水率。2013年,Jurdan等[16]利用PROGEOSAILH植被辐射传输模型反演了地中海区域森林植被的冠层可燃物含水率。虽然基于植被辐射传输模型的植被冠层可燃物含水率数据取得了比较好的参数反演效果,且模型构建过程中考虑了光学-红外信号的辐射传输机制,使得此类方法具有较高的普适性。但是依然存在4个方面的主要问题而限制其在大区域尺度下的高精度工程化应用:①干物质重量的光谱弱敏感特性使得基于植被辐射传输模型的植被冠层可燃物含水率参数遥感反演总体精度不高,特别是对于垂直结构更为复杂和空间异质性更强的森林覆盖类型的植被区域而言[36];②植被辐射传输模型对于地表植被的结构假设比较理想化[37],导致模型只适用于单一物种和均匀分布结构的植被区域,对于更为复杂的混合分层结构植被(如森林植被),模型难以准确地描述植被冠层的光谱反射特征,从而影响目标参数的定量反演效果;③以上所有的植被辐射传输模型均需要较多的输入参数来具体描述植被的属性特征,而其中某些参数的参数化过程又比较困难[36],这种不精确的模型参数化过程对于植被冠层可燃物含水率参数的高精度反演也有一定的限制作用;④植被辐射传输模型未知参数的个数一般多于观测多光谱反射率波段数,并且不同输入参数的差异性组合有可能对应相同或相似的多光谱反射率曲线,导致基于植被辐射传输模型的植被关键参数遥感反演方案本质上是病态多解的[38]。
1.2 基于微波遥感的植被冠层可燃物含水率反演
基于微波遥感的植被冠层可燃物含水率参数定量反演同样可以分为两个类别:传统基于经验统计模型的定量估算方法和基于雷达辐射传输理论的半经验或物理模型的反演方法。
传统基于经验统计模型的植被冠层可燃物含水率遥感定量估算方法主要是统计分析雷达后向散射系数以及相关衍生参数(比值参数、雷达植被指数和极化特征分解参数等)与植被冠层可燃物含水率的线性(或非线性)相关性,并以此来构建相应的估算模型。此类方法同样只适用于特定的且范围比较小的研究区域,普适性较差。相比基于可见光-红外遥感的植被冠层可燃物含水率反演方法,由于微波信号较强的穿透能力以及对于水分信息的高敏感性[20,22],基于微波遥感的植被冠层可燃物含水
展开