春季東北地區(qū)森林火險(xiǎn)氣象指數(shù)及其極值重現(xiàn)期特征*

劉吉敏，黃 泓，王學(xué)忠

(1.解放軍92493部隊(duì)，遼寧 葫蘆島 125000；2.國(guó)防科技大學(xué) 氣象海洋學(xué)院，江蘇 南京 211101)

春季東北地區(qū)森林火險(xiǎn)氣象指數(shù)及其極值重現(xiàn)期特征*

劉吉敏1，黃 泓2，王學(xué)忠2

(1.解放軍92493部隊(duì)，遼寧 葫蘆島 125000；2.國(guó)防科技大學(xué) 氣象海洋學(xué)院，江蘇 南京 211101)

以1961-2010年間東北地區(qū)逐日氣象觀(guān)測(cè)數(shù)據(jù)為基礎(chǔ)，計(jì)算春季逐日森林火險(xiǎn)氣象指數(shù)(FFDI)并用廣義帕累托分布(GPD)模型對(duì)FFDI極值進(jìn)行了擬合，結(jié)果表明：東北地區(qū)春季FFDI高值區(qū)分布于吉林西部、遼寧西部和內(nèi)蒙東四盟市南部區(qū)域，高值區(qū)基本位于山脈背風(fēng)坡一側(cè)。春季各站FFDI值的趨勢(shì)變化大致以大興安嶺山脊為界，山脊以東站點(diǎn)的FFDI趨勢(shì)系數(shù)大部分為負(fù)值，其中部分站點(diǎn)趨勢(shì)系數(shù)通過(guò)了95%置信水平檢驗(yàn)，山脊以西站點(diǎn)趨勢(shì)系數(shù)為正值，但趨勢(shì)均不顯著。GPD模型的概率分布能很好的擬合FFDI極值，96.5%站點(diǎn)的FFDI極值經(jīng)驗(yàn)分布和理論分布相一致。20年一遇和50年一遇的FFDI極值空間分布大體一致，在50年一遇極值高值區(qū)中FFDI值超過(guò)40，低值區(qū)的FFDI值普遍低于25。

東北地區(qū)；森林火險(xiǎn)；氣象指數(shù)；極值；重現(xiàn)期

森林火災(zāi)是一種突發(fā)性強(qiáng)、破壞性大、處置困難的自然災(zāi)害，嚴(yán)重的森林火災(zāi)給當(dāng)?shù)貛?lái)重大的經(jīng)濟(jì)損失、人員傷亡和環(huán)境污染。Mouillot等[1]統(tǒng)計(jì)表明森林火災(zāi)占每年全球燃燒面積的14%，火災(zāi)面積約50～90 hm2。Dempsey對(duì)安大略省森林火災(zāi)的環(huán)境效應(yīng)研究顯示森林火災(zāi)產(chǎn)生的煙和排放物會(huì)增加空氣中的PM2.5、O3、碳?xì)浠衔锖偷趸餄舛?，給風(fēng)向下游地區(qū)帶來(lái)嚴(yán)重的空氣污染[2]。森林火災(zāi)發(fā)生原因眾多[3]，如閃電等自然原因，或高壓電線(xiàn)接觸植被意外觸發(fā)，亦可由人為縱火引起。無(wú)論何種原因引發(fā)火災(zāi)，森林火災(zāi)的蔓延、傳播和撲滅都需要一定的有利氣象條件，為了對(duì)森林火險(xiǎn)等級(jí)進(jìn)行評(píng)估，國(guó)內(nèi)外研究機(jī)構(gòu)開(kāi)發(fā)了一系列利用氣象因子預(yù)測(cè)火險(xiǎn)指數(shù)的預(yù)報(bào)系統(tǒng)，如加拿大火險(xiǎn)等級(jí)預(yù)報(bào)系統(tǒng)[4]、澳大利亞McArthur火險(xiǎn)尺[5-6]、美國(guó)國(guó)家火險(xiǎn)等級(jí)預(yù)報(bào)系統(tǒng)[7]等，中國(guó)國(guó)家氣象局基于修正的布龍-戴維斯方法開(kāi)發(fā)了業(yè)務(wù)化運(yùn)行的森林火險(xiǎn)氣象指數(shù)[8]。

東北地區(qū)森林資源豐富，在大興安嶺、小興安嶺、長(zhǎng)白山地區(qū)分布有我國(guó)最大的天然林區(qū)，每年春秋季節(jié)防火任務(wù)繁重，研究容易引發(fā)火災(zāi)的天氣背景對(duì)防火救災(zāi)有重要的現(xiàn)實(shí)意義。牛若蕓等[8]對(duì)國(guó)際通用的5個(gè)森林火險(xiǎn)指數(shù)進(jìn)行比較后發(fā)現(xiàn)FFDI指數(shù)在我國(guó)東北地區(qū)實(shí)用性最好。在實(shí)際應(yīng)用中FFDI值越大表明火險(xiǎn)等級(jí)越高，說(shuō)明在當(dāng)時(shí)氣象條件下火災(zāi)極易發(fā)生，且火災(zāi)發(fā)生后難以控制，對(duì)于FFDI極大值需要引起足夠的重視。在極值研究方法中，廣義帕雷托分布(Generalized Pareto Distribution，GPD)近三十年來(lái)在氣象學(xué)極值研究中得到了成功應(yīng)用[9]，不同于以往極值研究中選取年最大值A(chǔ)M(Annual Maximum)構(gòu)成時(shí)間序列，GPD擬合中通過(guò)設(shè)定門(mén)限值，建立超過(guò)門(mén)限值數(shù)據(jù)的極值概率分布，即POT(Peaks over Threshold)抽樣[9]。POT抽樣充分利用了有限的觀(guān)測(cè)信息，增加了極值的樣本量，使模型的參數(shù)估計(jì)更加合理、重現(xiàn)期精度更高。本文利用東北地區(qū)氣象站點(diǎn)觀(guān)測(cè)數(shù)據(jù)計(jì)算了春季FFDI值，分析東北地區(qū)春季FFDI的空間分布和變化趨勢(shì)，并用GPD模型對(duì)FFDI極值進(jìn)行擬合，研究了給定重現(xiàn)期下FFDI極值的分布特點(diǎn)。

1 資料和方法

1.1 資料說(shuō)明

本文中東北地區(qū)包括黑龍江省、吉林省、遼寧省和內(nèi)蒙古東四盟市(呼倫貝爾市、通遼市、赤峰市、興安盟)。東北春季防火期研究時(shí)段取為3月1日至5月31日。所用資料為中國(guó)氣象局國(guó)家氣象中心資料室提供的中國(guó)743測(cè)站1961-2010年逐日資料，資料中包括了溫度、相對(duì)濕度和降水等要素。選取東北地區(qū)氣象記錄時(shí)間長(zhǎng)、缺測(cè)率低的88個(gè)站點(diǎn)，各站點(diǎn)在1951-2010年間觀(guān)測(cè)要素缺測(cè)率均低于1%，東北地區(qū)地形和88個(gè)氣象站點(diǎn)分布如圖1所示。

圖1 東北地區(qū)氣象站點(diǎn)分布(圓點(diǎn))和地形(陰影區(qū))。

1.2 FFDI計(jì)算

森林火險(xiǎn)氣象指數(shù)(Forest Fire Danger Index，簡(jiǎn)稱(chēng)FFDI)計(jì)算中需要風(fēng)速、氣溫、相對(duì)濕度和干旱因子。Noble等[10]通過(guò)對(duì)McArthur火險(xiǎn)尺[5-6]的擬合得到FFDI計(jì)算公式：

FFDI=2D0.987exp(-0.45-0.0345H+0.0338T+0.023V)。

(1)

式中：D∈[0,10]為干旱因子；H為相對(duì)濕度(%)；T為攝氏溫度(℃)；V為距地表10 m高度處平均風(fēng)速(km/h)。實(shí)際計(jì)算中T取為日最高溫度，H應(yīng)取日最高溫度對(duì)應(yīng)的相對(duì)濕度，但資料中僅有日平均相對(duì)濕度，故將H取為日平均相對(duì)濕度。FFDI的計(jì)算主要依賴(lài)于干旱因子D，求取干旱因子D過(guò)程中涉及到Keetch-Byram drought index(KBDI)的復(fù)雜計(jì)算，KBDI是用降水和蒸散估計(jì)值的差異來(lái)表示森林土壤水分的虧缺。

第t日的KBDI可定義為：

KBDIt=KBDIt-1+ETt-1-NRt-1。

(2)

式中：ET為日蒸散量(mm)；NR為日降水量減去植被冠層攔截水量后的凈降水量(mm)。蒸散量ET的估計(jì)值為：

(3)

(4)

在KBDI的計(jì)算中需要初值和初始化期，將每年2月1日的FFDI初值設(shè)為0，每年的2月份作為初始化期，經(jīng)初始化后可計(jì)算3-5月的逐日FFDI。由KBDI計(jì)算干旱因子 時(shí)參考Griffiths定義[11]，即D=max[min(D*,10),0]，式中D*公式為：

(5)

其中λ=max{ψ*,ψt,ψt-1,…,ψt-19}，ψt的算法為

(6)

式中：i=0,1,…，19，ψ*為Finkele等[12]提出的修正項(xiàng)，公式為：

(7)

2 FFDI氣候特征

2.1 FFDI和降水空間分布

東北地區(qū)各站春季FFDI日平均值分布如圖2a所示，從圖中可看出春季FFDI大值區(qū)主要分布在內(nèi)蒙東四盟市南部、遼寧西部和吉林西北部，這些大值區(qū)域?qū)?yīng)于東北平原和科爾沁草原，站點(diǎn)FFDI日均值超過(guò)4.4，F(xiàn)FDI日均值最大值為6.7。在遼寧中部、黑龍江西部和呼倫貝爾西部地區(qū)的FFDI日均值在3.2～4.4，而大興安嶺、小興安嶺和長(zhǎng)白山山脈海拔較高區(qū)域的FFDI日均值小于3.2。

圖2 1961-2010年春季平均值空間分布。

在春季各站日均降水分布(圖2b)中，日均降水量大值區(qū)集中在遼寧東部和吉林南部，日均降水量超過(guò)1mm，日均降水量最大值為1.7mm，黑龍江大部和呼倫貝爾北部的日均降水量在1.0～1.0mm，吉林西部和內(nèi)蒙東四盟市中南部的日均降水量低于0.6mm。因?yàn)樵诖杭久晒艊?guó)被高壓脊控制(圖略)，冷高壓攜帶冷空氣頻繁南下，致使東北地區(qū)春季降水量為四季中最小值，同時(shí)春季東北地區(qū)高空盛行西北氣流，內(nèi)蒙東四盟市南部、吉林西部位于大興安嶺背風(fēng)坡，遼寧西部也位于奴魯爾虎山背風(fēng)坡一側(cè)，背風(fēng)坡一側(cè)下沉氣流也導(dǎo)致這些區(qū)域春季干燥少雨，而在遼寧東部和吉林南部位于長(zhǎng)白山迎風(fēng)坡，降水相對(duì)較多。由FFDI和降水量分布比較可看出，F(xiàn)FDI大值區(qū)與降水量低值區(qū)能較好對(duì)應(yīng)。

2.2 FFDI變化趨勢(shì)

對(duì)各站50年的春季FFDI平均值進(jìn)行線(xiàn)性趨勢(shì)分析(圖3)，從圖中可看出東北地區(qū)春季FFDI變化趨勢(shì)分布大致以大興安嶺山脊為界，在山脊以東區(qū)域大部分站點(diǎn)的FFDI趨勢(shì)系數(shù)為負(fù)值，表明這些地區(qū)FFDI有減小趨勢(shì)，尤其在東北平原和遼寧西部地區(qū)站點(diǎn)的趨勢(shì)系數(shù)通過(guò)了95%置信水平的F檢驗(yàn)。在山脊西側(cè)大部分站點(diǎn)的FFDI趨勢(shì)為正，但各站點(diǎn)趨勢(shì)系數(shù)并不顯著。

圖3 1961-2010年春季FFDI趨勢(shì)系數(shù)。(注：紅色表示趨勢(shì)為正值，綠色為負(fù)值，填色圓點(diǎn)表示趨勢(shì)系數(shù)通過(guò)95%置信水平檢驗(yàn)。)

3 FFDI極值GDP擬合

廣義帕雷托分布(GPD)的分布函數(shù)為：

(8)

概率密度函數(shù)為：

(9)

式中：ζ為門(mén)限值；a稱(chēng)為尺度參數(shù)；κ為線(xiàn)型參數(shù)或形狀參數(shù)。GPD公式中參數(shù)的估計(jì)采用極大似然估計(jì)法(MLE)，具體算法參考文獻(xiàn)[13-14]。

參數(shù)確定后可得到GPD的極值分位數(shù)，該公式為：

(10)

式中：λ=n/m，n為大于門(mén)限值的極值數(shù)量，m為總樣本量；T為選定的重現(xiàn)期值。

3.1 FFDI極值空間分布

對(duì)門(mén)限值的選取是進(jìn)行GPD模擬的關(guān)鍵步驟，如果門(mén)限值選取過(guò)高，超門(mén)限樣本數(shù)就太少，估計(jì)的參數(shù)方差很大；若門(mén)限值選取過(guò)小，則超門(mén)限值樣本數(shù)過(guò)多，會(huì)導(dǎo)致估計(jì)值產(chǎn)生大的偏差。目前關(guān)于門(mén)限值的選取有多種，如超額限望法、試湊法、擬合優(yōu)度法和Hill圖解法等[15-17]。萬(wàn)仕全等[18]對(duì)南京過(guò)去100年極端日降水模擬表明，對(duì)于長(zhǎng)度不小于50年的氣候序列，第91%百分位點(diǎn)上的數(shù)據(jù)作為門(mén)限值基本能滿(mǎn)足GPD研究需要，本文采用該方法來(lái)確定各站FFDI的門(mén)限值。Bonsal等[19]提出一種百分位閾值的計(jì)算方法：

(11)

式中：n為站點(diǎn)序列{X}的總樣本量，m為按升序排列的序列{x1,x2,…,xn}中xn的序號(hào)(m=1,2,…,n)，而p定義為變量X的某個(gè)觀(guān)測(cè)值小于等于xm的概率。由此，選取p為91%時(shí)的xm作為各站的門(mén)限值，若序列中沒(méi)有p值恰好為91%的觀(guān)測(cè)值，則門(mén)限值由p為91%相鄰的觀(guān)測(cè)值線(xiàn)性?xún)?nèi)插得到。

圖4為超過(guò)各站門(mén)限值的FFDI極值平均值分布，從圖可看出FFDI極值高值區(qū)主要分布在吉林白城市到內(nèi)蒙古赤峰市之間的帶狀區(qū)域內(nèi)，遼寧西部的朝陽(yáng)市也有高值區(qū)存在。圖中也可發(fā)現(xiàn)在長(zhǎng)春、哈爾濱、大慶一帶FFDI極值平均值也較高，這可能與這一帶為人口密集區(qū)有關(guān)，研究表明[20-21]在人口稠密區(qū)空氣溫度比郊區(qū)偏高，可能導(dǎo)致這一地帶FFDI較高。

圖4 超過(guò)門(mén)限值的春季FFDI指數(shù)平均值空間分布。

3.2 GPD擬合檢驗(yàn)

為了對(duì)GPD模型來(lái)擬合FFDI極值的合理性進(jìn)行診斷，首先查看2個(gè)站點(diǎn)的GPD擬合診斷圖。圖4為齊齊哈爾站FFDI極值的GPD擬合效果，可見(jiàn)圖4a中(PP圖)的點(diǎn)幾乎都落在了直線(xiàn)上，在圖4b中(QQ圖)中除了高分位數(shù)有個(gè)別點(diǎn)誤差較大外，其余點(diǎn)均能落在直線(xiàn)上，表明GPD模型能很好的擬合樣本數(shù)據(jù)，圖4c中重現(xiàn)期水平圖上可看出重現(xiàn)期都在95%置信區(qū)間內(nèi)，也說(shuō)明了模型的合理性。圖5為長(zhǎng)春站的模型診斷圖，圖中結(jié)果說(shuō)明GPD模型對(duì)長(zhǎng)春站也同樣適用。

圖4 GPD方法對(duì)齊齊哈爾站FFDI極值的擬合診斷圖。

圖5 GPD方法對(duì)長(zhǎng)春站FFDI極值的擬合診斷圖。

圖6 給定重現(xiàn)期下FFDI極值的空間分布。

Kolmogorov-Smirnov正態(tài)性檢驗(yàn)(K-S檢驗(yàn))用來(lái)比較一個(gè)經(jīng)驗(yàn)分布F(x)和理想分布G(x)是否相同。假定F(x)和G(x)中樣本數(shù)分別為m和n，定義：

(12)

表1 FFDI極端值GPD分布的K-S檢驗(yàn)

用K-S檢驗(yàn)法對(duì)85個(gè)站點(diǎn)FFDI極值的GPD擬合進(jìn)行檢驗(yàn)(表1)，可在86個(gè)站點(diǎn)中只有3個(gè)站點(diǎn)(約3.5%)的經(jīng)驗(yàn)分布和理論分布在90%置信水平下不一致，在95%置信水平下不一致的僅有2個(gè)，表明用GPD分布來(lái)描述各站FFDI極值是切實(shí)可行的。

3.3 給定重現(xiàn)期極值的分布

圖6為給定重現(xiàn)期下FFDI極值的空間分布，在20年一遇和50年一遇重現(xiàn)期下FFDI極值分布基本相似。圖中可見(jiàn)FFDI極值大值中心分布不均勻，在遼寧鐵嶺，吉林白城、長(zhǎng)春，黑龍江嫩江有大值中心存在，此外內(nèi)蒙東四盟市西部地區(qū)FFDI極值也較高。在50年一遇重現(xiàn)期下，大值中心的FFDI極值都超過(guò)40，其中黑龍江嫩江的FFDI極值高達(dá)57。圖中也可看出，在大興安嶺和長(zhǎng)白山區(qū)域的FFDI極值普遍較低，在50年一遇重現(xiàn)期下這些區(qū)域的FFDI極值低于25。

4 結(jié)論與討論

以東北地區(qū)85個(gè)站點(diǎn)的春季逐日氣象觀(guān)測(cè)數(shù)據(jù)為基礎(chǔ)，研究了春季各站FFDI的空間分布和變化趨勢(shì)，并用GPD擬合研究了FFDI極值分布和重現(xiàn)期特征，結(jié)論如下：

(1)東北地區(qū)春季FFDI高值區(qū)分布于吉林西部、遼寧西部和內(nèi)蒙東四盟市南部區(qū)域，這些區(qū)域基本位于山脈背風(fēng)坡一側(cè)。各站50年春季FFDI值的趨勢(shì)變化大致以大興安嶺山脊為界，山脊以東站點(diǎn)的FFDI趨勢(shì)系數(shù)大部分為負(fù)值，其中部分站點(diǎn)趨勢(shì)系數(shù)通過(guò)了95%置信水平檢驗(yàn)，山脊以西站點(diǎn)趨勢(shì)系數(shù)為正值，但趨勢(shì)不顯著。

(2)利用GPD模型的概率分布對(duì)各站FFDI極值的擬合效果較好，96.5%站點(diǎn)的FFDI極值分布和理論分布相一致。

(3)20年一遇和50年一遇的FFDI極值空間分布大體一致，在50年一遇中FFDI極值高值區(qū)FFDI值超過(guò)40，而大興安嶺和長(zhǎng)白山區(qū)域的FFDI極值相對(duì)較低，這些區(qū)域的FFDI值普遍低于25。

[1] Mouillot F, Field CB.Fire history and the global carbon budget: a 1°×1° fire history reconstruction for the 20th century [J].Glob Chang Biol., 2005, 11:398-420.

[2] Dempsey F.Forest fire effects on air quality in Ontario: Evaluation of several recent examples [J].Bull.Amer.Meteor.Soc., 2013, 94:1059-1064.

[3] Beale J and W Jones.Preventing and reducing bushfire arson in Australia: A review of what is known [J].Fire Technol., 2011, 47:507-518

[4] Stocks B J, B D Lawson, M E Alexander, et al.The Canadian forest fire danger rating system: An overview [J].For.Chron., 1989, 65:450-457.

[5] McArthur A G.Weather and grassland fire behavior [M].Forestry and Timber Bureau Leaflet, 1966, 100:23.

[6] McArthur, A.G.Fire behavior in eucalypt forests [M].Forestry and Timber Bureau Leaflet, 1967, 107:36.

[7] Deeming J E, R E Burgan and J D Cohen.The National Fire-Danger Rating System [M].USDA Forest Service General Tech.Rep.INTUS, 1978, 39-49.

[8] 牛若蕓,翟盤(pán)茂,佘萬(wàn)明.森林火險(xiǎn)氣象指數(shù)的應(yīng)用研究[J].應(yīng)用氣象學(xué)報(bào),2007,18(4):479-489.

[9] 丁裕國(guó),李佳耘,江志紅.極值統(tǒng)計(jì)理論的進(jìn)展及其在氣候變化研究中的應(yīng)用[J].氣候變化研究進(jìn)展,2011, 7(4):248-252.

[10] Noble I R, G A V Bary and A M Gill.McArthur’s fire danger meters expressed as equations [J].Aust.J.Ecol., 1980, 5:201-203.

[11] Griffiths D.Improved formula for the drought factor in McArthur’s forest fire danger meter [J].Aust.For., 1999, 62:202-206.

[12] Finkele K, G A Mills, G Beard, et al.National gridded drought factors and comparison of two soil moisture deficit formulations used in prediction of Forest Fire Danger Index in Australia[J].Aust.Meteor.Mag., 2006, 55:183-197.

[13] Martins E S and Stedinger J R.Generalized maximum likelihood extreme value quantile estimators for hydrologic data [J].Water Resources Research, 2000, 36 (3):737-744.

[14] Martins E S and Stedinger J R.Generalized maximum likelihood Pareto-Poisson estimators for partial duration series [J].Water Resources Research, 2001,37 (10): 2551-2557.

[15] 司波,余錦華,丁裕國(guó).四川盆地短歷時(shí)強(qiáng)降水極值分布的研究[J].氣象科學(xué),2012,32(4):403-410.

[16] 江志紅,丁裕國(guó),朱蓮芳,等.利用廣義帕雷托分布擬合中國(guó)東部日極端降水的試驗(yàn)[J].高原氣象,2009, 28(3):573-580.

[17] 蔡敏,丁裕國(guó),江志紅.我國(guó)東部極端降水時(shí)空分布及其概率特征[J] .高原氣象, 2007, 26(2):309-317

[18] 萬(wàn)仕全,周?chē)?guó)華,潘柱,等.南京過(guò)去100年極端日降水量模擬研究[J].氣象學(xué)報(bào),2010,68(6):790-799.

[19] Bonsal B R, ZHANG X, Vincent L A et al.Characteristics of Daily and Extreme Temperatures over Canada [J].2001, 14:1959-1976.

[20] 侯學(xué)英.中國(guó)城市化水平與城市氣溫變化的灰色關(guān)聯(lián)分析[J].生態(tài)經(jīng)濟(jì),2009(1):59-62.

[21] 邵全琴,孫朝陽(yáng),劉紀(jì)遠(yuǎn),等.中國(guó)城市擴(kuò)展對(duì)氣溫觀(guān)測(cè)的影響及其高估程度[J].地理學(xué)報(bào),2009,64(11):1292-1302.

Forest Fire Danger Index of Northeast China Region in Summer and its Characteristics of Extreme Value

LIU Jimin1, HUANG Hong2and WANG Xuezhong2

(1.No.92493ofPLA,Huludao125000,China;2.SchoolofMeteorologyandOceanography,NationalUniversityofDefenceTechnology,Nanjing211101,China)

Based on daily observational data in northeast China region in the period 1961-2010, the forest fire danger index (FFDI) were calculated and Generalized Pareto Distribution (GPD) was also used to fit the extreme value of FFDI, the main results show as following: the high value area of FFDI were located in western Jilin province, western Liaoning province and Chifeng, Tongliao, Hinggan league in Inner Mongolia, and all these areas are basically located in the leeward slope. The distribution of trend coefficient of FFDI were divided into two parts, the stations located to the east of the ridge of the Greater Khingan presented a negative trend and the trend were significant for some stations, the trend were positive for stations west of the ridge of Greater Khingan, but none of them was significant. The probability distribution of GPD model can be well fitted to the extreme value of FFDI, and the empirical distributions of 96.5% station were consistent with the theoretical distribution. The spatial distributions were substantial agreement for the 20-year and 50-year return values of extreme FFDI, the high value of 50-year return value were greater than 40 and low value were generally less than 25.

northeast China region; forest fire risk; meteorological index; extreme value; return period

劉吉敏，黃泓，王學(xué)忠.春季東北地區(qū)森林火險(xiǎn)氣象指數(shù)及其極值重現(xiàn)期特征[J].災(zāi)害學(xué)，2018，33(2)：76-80.[LIU Jimin, HUANG Hong and WANG Xuezhong.Forest Fire Danger Index of Northeast China Region in summer and Its Characteristics of Extreme Value[J].Journal of Catastrophology，2018，33(2)：76-80.

10.3969/j.issn.1000-811X.2018.01.015.]

2017-07-24

2017-09-18

國(guó)家自然科學(xué)基金(41375049；40905021；41275099)；中國(guó)博士后科學(xué)基金(2011M500894)；國(guó)家公益性行業(yè)(氣象)科技專(zhuān)項(xiàng)(GYHY201206005)

劉吉敏(1987-)，男，甘肅武威人，助理工程師，碩士，從事短期氣候預(yù)測(cè).E-mail:atmosljm@163.com

P49；X43

A

1000-811X(2018)01-0076-05

10.3969/j.issn.1000-811X.2018.01.015