宗學政 田曉瑞
(中國林業科學研究院森林生態環境與保護研究所 國家林業和草原局森林保護學重點開放性實驗室 北京 100091)
火是大興安嶺森林生態系統中重要的生態影響因子,對該區域森林景觀格局及演替起著重要作用(劉志華等,2009)。1987年以后,大興安嶺地區加強了林火管理,逐步建立了完善的森林防火機構和撲火隊伍,森林防火投入不斷增加,林火管理能力顯著提高,火發生次數也明顯減少(Tianetal., 2020)。但嚴格的林火管理政策改變了該地區的火動態,林內可燃物不斷積累,未來高強度森林火災的概率增加(田曉瑞等,2005)。可燃物是森林燃燒的物質基礎(Bradstocketal., 2010),也是燃燒三要素(可燃物、氧氣和溫度)中最易調控的因子。因此,需開展可燃物的科學管理,應對未來的火災風險,預防森林大火發生。
不同類型可燃物的火行為有差異(Ageeetal., 2005),可燃物載量及其空間結構影響火行為(蔓延速度、火強度等)(Anderson, 1970)。草地和灌木林的地表多細小可燃物,易形成地表火,火強度低,對這類可燃物多通過計劃火燒降低地表可燃物載量(Bradstocketal., 2006)。中幼齡針葉樹的自然整枝快,林下存有大量枯落物,且樹冠呈金字塔形,地表火易沿可燃物梯蔓延至樹冠形成樹冠火(張思玉等,1998)。老齡針葉林樹冠稀疏,缺少可燃物梯,較難形成樹冠火。對于針葉林的可燃物處理主要以機械清理為主,包括清理林內可燃物梯及地表粗可燃物(趙鳳君等,2010)。闊葉林內無可燃物梯,地表可燃物積累也較少,以清理地表可燃物為主。大興安嶺地區的森林以落葉針葉林和興安落葉松(Larixgmelinii)白樺(Betulaplatyphylla)混交林為主,可采用計劃火燒和機械方式清理林內可燃物梯與地表可燃物,減少林火發生(Changetal., 2007; 趙彬清等,2018)。
可燃物處理可有效減少林火的頻率和強度(Chionoetal., 2017)。野外可燃物處理試驗和火行為預測模型是評估可燃物管理措施效果的有效手段。McGinnis等(2010)通過模擬美國內華達山區可燃物不同處理(未處理、砍伐火燒木、使用除草劑)下的火行為,發現清理火燒木短期內會增加地表死可燃物載量,但在長時間尺度上對地表火行為沒有影響。Moghaddas等(2010)利用FlamMap和FARSITE模擬了內華達一些區域可燃物處理前后的燃燒概率和潛在火行為,發現可燃物處理后的森林燃燒概率降低,過火面積、火焰高度及樹冠火發生比例均下降。Salis等(2016)模擬了地中海地區可燃物處理后的火行為變化,表明隨可燃物處理強度增強和處理面積增大,燃燒概率、過火面積、火焰高度等顯著降低。Barnett等(2016)研究了可燃物處理與林火發生的相關性,提出了合理的可燃物處理措施。
可燃物管理是降低林火危害程度的有效措施(Fangetal., 2015)。可燃物特征與火行為顯著相關(吳志偉等,2012)。景觀模型模擬表明,可燃物處理后林火發生及火行為都明顯降低。劉志華等(2009)利用LANDIS模擬了可燃物處理對大興安嶺呼中林業局森林燃燒性的影響,認為經常進行計劃火燒和機械清理可燃物可減少災難性火災。陳宏偉等(2011)利用LANDIS模擬了不同采伐方式對森林燃燒性影響,發現擇伐可有效降低森林燃燒性及火行為。Wu等(2013)認為基于火險進行可燃物優化處理能進一步發揮可燃物處理的效益。林分尺度上的火行為模擬可有效評估可燃物處理效果(Ottmaretal., 2012)。
基于可燃物特征預測潛在火行為,是開展可燃物管理的基礎(Hollingsworthetal., 2012)。可燃物特征分類系統(fuel characteristic classification system,FCCS)是在林分尺度上模擬火行為的系統。該系統將每個林分看作景觀中的均質單元,基于Rothermel模型計算火行為指標,定量描述可燃物的燃燒特征(Ottmaretal., 2007),包括火強度、蔓延速度及火焰高度等(Sandbergetal., 2007)。該模型包括多種可燃物床層類型,用戶可選擇與研究對象特征相近的可燃物床層,并通過修改特征參數獲得所需的可燃物床層(Riccardietal., 2007; Arroyoetal., 2008)。
本文以大興安嶺地區主要林型為對象,設立標準地調查獲取主要林型森林的可燃物分布特征及環境變量,利用FCCS模擬不同強度的可燃物處理對林分潛在火行為的影響,為該區域開展可燃物科學管理提供參考。
大興安嶺(50°54′―53°33′N,121°11′―127°02′E)位于黑龍江省西北部。研究區總面積為646.2萬hm2。地勢較平坦,平均海拔573 m,屬淺山丘陵地帶。該區氣候為溫帶大陸性季風氣候,年均氣溫-2 ℃,冬長夏短,無霜期約90~110天,年均降水量450~500 mm(徐化成, 1998)。針闊混交林多以興安落葉松-白樺為主,常綠針葉林較少。主要樹種為興安落葉松、山楊(Populusdavidiana)、白樺、黑樺(Betuladahurica)等為主,灌木有興安杜鵑(Rhododendrondauricum)、越桔(Vacciniumvitisidaea)等,草本植物以杜香(Ledumpalustre)和苔蘚為主(周道瑋等,2010)。
根據統計資料(表1),2000—2018年共發生1 039起,平均每年發生55起,年均過火面積21 351.3 hm2。2007年火災次數最多(100起),但森林過火面積為817 hm2。2003年過火面積最多(251 373.4 hm2)。2013年火災次數最少,只有14起,過火面積11 hm2。

表1 2000—2018年大興安嶺地區火災發生次數及森林過火面積Tab.1 Fire number and burned areas in Daxing’anling during 2000-2018
選擇寒溫帶典型林分設置樣地。大興安嶺是我國重要的天然林區,天然林主要分布在大興安嶺中部和北部; 但也分布著大量人工林,主要在南部的加格達奇林業局。呼中林業局位于大興安嶺中部,主要森林為天然林,屬典型的寒溫帶針葉林。興安落葉松天然林(51°55′N,123°34′E)、白樺林(52°39′N,123°32′E)、興安落葉松-白樺混交林(52°20′N,123°35′E)樣地設在呼中林業局,而興安落葉松人工林(50°17′N,124°7′E)和樟子松(Pinussylvetrisvar.mongolica)人工林(50°19′N,124°6′E)樣地設在加格達奇林業局。
天然林內的灌木和草本生長密集,地表凋落物多,腐殖質層較厚; 落葉松人工林下無灌木,草本較少,針葉堆積較厚,地表散布前期清理的少量枯枝; 樟子松人工林清理較少,林下多灌木及針葉凋落物(圖1)。

圖1 5種林型植被狀況Fig.1 Vegetation of the five forest types
在距離道路或林緣20 m以上的林地內,每種林型設置3個樣地(20 m×20 m),樣地間隔10 m以上,立地環境相近,共設置15塊樣地(表2)。在每塊樣地內調查林分的樹種組成和郁閉度(包括活喬木及雖然死亡但仍有葉片殘余的死亡木)。對樣地內的喬木每木檢尺,分別測定活立木及死木的胸徑、樹高、可燃物梯高度(最低和最高)。郁閉度為樹冠垂直投影到樣地對角線上的長度與2條對角線長度和的比值(梁瀛等,2017)。死木分為葉片完整、無葉片但枝干完整、無側枝但主枝干完整、主枝干不完整4類(Riccardietal.,2007)。
在每個樣地內設置4塊10 m×10 m小樣地,其中1塊為對照,其他3塊做低、中和高強度的可燃物處理。低強度處理:割除林內易燃及枯死灌草,打掉林內3 m以下全部枯枝并移出林地,將地表大于10 h時滯的可燃物清理出林地; 中強度處理:清理林內3 m以下枯枝,連同林內倒木全部移出林地,割除林內可能影響樹冠火的灌木及小喬木,并清理地表未分解的枯落物等可燃物; 高強度處理:清理林內3 m以下枯枝并移出林地,割除小樣地內全部的灌草,移除林內枯死木及腐爛樹樁,將地表可燃物清理至腐殖質層。
對每塊10 m×10 m樣地進行可燃物處理前后的調查。地表凋落物包括枯枝、葉片、球果及枯死草本; 腐殖質是地表凋落物與土壤層之間的物質(Riccardietal., 2007)。倒死木質可燃物是距地面2 m以下的木質可燃物,分為1 h時滯(0~0.063 5 cm)、10 h時滯(0.063 5~2.54 cm)、100 h時滯(2.54~7.62 cm)、1 000 h時滯(7.62~22.86 cm)、10 000 h時滯(22.86~50.8 cm)及>10 000 h時滯(>50.8 cm)6個等級(Keaneetal., 2012)。
灌木分成主要和次要灌木,草本分成主要和次要草本(Prichardetal., 2011)。灌木和草本的高度、活植株的比例和可燃物及腐殖質厚度調查采用樣方法。將每塊10 m×10 m的樣方劃分為4塊5 m×5 m的樣方,并在每塊5 m×5 m樣方中,均勻設置5塊1 m×1 m樣方,在每個小樣方內分類調查,求得各類可燃物的特征平均值。各種類型可燃物蓋度的調查采用樣線法(Parresoletal., 2012),沿10 m×10 m樣方的2條對角線分別按10 cm間隔,記錄每段出現的灌木、草本、地表可燃物、腐殖質及各等級木質可燃物,各類型可燃物在2條對角線的頻度平均值作為該類可燃物的蓋度。
收集每個10 m×10 m樣地內的灌木、草本、地表可燃物及腐殖質,分別稱重計算載量。在每塊樣地內采集喬木冠層葉片、灌木、草本及1~100 h時滯的木質可燃物,每種樣品量在100 g以上,采樣后馬上稱取鮮質量后再帶回實驗室測定含水率。室內含水率測定條件是103 ℃烘干至恒質量。
FCCS系統包括全球458個不同的可燃物床層數據(Riccardietal., 2007)。可燃物床層參數包括喬木、灌木、草本、倒木、地衣苔蘚和土壤等6個層次的信息,每一層又包括2類以上的信息。通過樣地調查,獲取了FCCS系統的所需參數,在FCCS系統中建立大興安嶺地區典型林分的可燃物床層,預測不同天氣情景下的火行為。
可燃物處理只改變了林內的可燃物載量,而環境條件保持一致。根據野外調查數據設置可燃物床層的參數。在FCCS系統中分別選擇落葉闊葉林、落葉針葉林、針闊混交林及常綠針葉林為原始模板,基于調查結果調整相關參數,準確表達白樺純林、落葉松林、落葉松白樺混交林和樟子松林的林分狀況。
環境變量輸入包括樹冠、灌木和草本的活葉片、1~100 h時滯可燃物的含水率、林地坡度(表3)和風速。這些環境變量可根據每種林型的調查結果設定,風速采用該區域歷史氣象數據的均值,模擬各林分在火險期內平均(一般)天氣條件下的火行為。輸出結果包括地表火的蔓延速度(m·s-1)、火強度(kW·m-2)和火焰高度(m)。
FCCS還可以模擬BehavePlus模型中的16種標準可燃物含水率情景。D和L分別為林內死可燃物和活可燃物含水率,D1、D2、D3和D4分別表示死可燃物含水率極低、低、中和高; L1、L2、L3和L4分別表示草本植物全部干枯、2/3草本植物干枯、1/3草本植物干枯和草本沒有干枯(Prichardetal., 2011)。D2L2是系統默認的一種干旱情景,表示死可燃物含水率低且2/3草本植物干枯狀態。模型只能基于該情景模擬潛在火行為,包括潛在地表火和樹冠火行為指數(表4)(這些指數用0~9等級值表示),數值越大表示火行為越強(Prichardetal., 2011)。
不同林型林冠層下可燃物載量之間存在顯著差異(F檢驗,P<0.01)(圖2)。其中興安落葉松白樺混交林的灌木、草本、地表凋落物和倒死木質可燃物載量分別為2.6、0.7、2.3和4.3 t·hm-2。白樺林內平均灌木、草本、地表凋落物和倒死木質可燃物載量分別為0.9、1.1、2.0和4.5 t·hm-2。興安落葉松天然林內平均灌木、草本、地表凋落物和倒死木質可燃物載量分別為0.6、1.6、0.4和0.5 t·hm-2。興安落葉松人工林內平均灌木、草本、地表凋落物和倒死木質可燃物載量分別為0、0.1、5.1和3.9 t·hm-2。樟子松人工林內平均灌木、草本、地表凋落物和倒死木質可燃物載量分別為0.3、0.8、3.9和0.2 t·hm-2。

表3 一般天氣情景下地表火模擬環境Tab.3 Simulating environment for surface fires under average weather scenario %

表4 干旱情景下潛在地表火及樹冠火指標Tab.4 Indexes of the surface fire and potentialcrown fire under drought scenario
與未進行可燃物處理情景相比,低強度可燃物處理情景下,5種林型的灌木、草本和地表凋落物降低30%,倒死木質可燃物載量變化明顯。其中興安落葉松白樺混交林和白樺林內存有較多的大于10 h的倒死木質可燃物,低強度處理清除了這些可燃物,導致這2種林型內倒死木質可燃物載量分別降低88.4%和91.1%(圖2c)。隨著可燃物處理強度的增加,可燃物載量明顯降低(F檢驗,P=0.000)。在中強度可燃物處理情景下,各林型灌木、草本和地表凋落物載量降低60%。在高強度可燃物處理情景下,5種林型林冠下可燃物載量降低90%。其中灌木和草本可燃物全部清除,地表凋落物及1 h時滯可燃物零散分布在樣地內,無法支持火的持續燃燒。

圖2 不同可燃物處理強度的可燃物載量變化Fig.2 Change of fuel loading after the fuel treatments
一般天氣情景下,興安落葉松天然林平均地表火蔓延速度、火焰高度及火強度分別為2.1(2.1~2.2)m·min-1、1.3(1.2~1.4)m和1 450.9(536.1~2 365.8)kW·m-2。落葉松人工林相應分別為0.4(0.2~0.6)m·min-1、0.3(0~0.6)m和208.3(18.5~398.1)kW·m-2。白樺林、興安落葉松白樺混交林和樟子松人工林的地表火蔓延速度分別為0.9(0.7~1.1)、0.7(0.6~0.8)和0.8(0.7~0.9)m·min-1; 火焰高度分別為0.6、0.7和0.6 m; 火強度分別為391.5(213.3~569.7)、931.1(701~1 161.3)和761.2(604.3~918.2)kW·m-2。
隨著可燃物處理強度增加,5種林型的地表火蔓延速度、火強度和火焰高度都明顯降低(F檢驗,P<0.001)(圖3)。與可燃物未處理的對照相比,低強度處理情景的興安落葉松天然林和人工林的地表火行為降幅明顯:平均蔓延速度分別降低51.6%和42.8%,平均火焰高度分別降低33.6%和39.4%,平均火強度分別降低22.8%和34%。白樺林平均蔓延速度和火焰高度分別降低15.7%和12.3%。樟子松人工林的平均蔓延速度和火焰高度降幅最小,分別為3.9%和7.9%。
中強度可燃物處理后,興安落葉松天然林地表火平均蔓延速度、火焰高度及火強度分別降低79.1%、67.2%和54.6%,興安落葉松人工林相應分別降低83.3%、69.7%和54.5%,最快蔓延速度為0.1 m·min-1,最高火焰高度為0.1 m。興安落葉松白樺混交林、白樺林和樟子松人工林地表火平均蔓延速度分別降低29.4%、37.1%和19.7%,火焰高度分別降低33.3%、29.8%和38.1%,火強度分別降低50.7%、22.5%和54.5%,其中興安落葉松白樺混交林、白樺林和樟子松人工林最快蔓延速度分別為0.6、0.7和0.6 m·min-1,最高火焰高度分別為0.5、0.4和0.4 m。
高強度可燃物處理后,興安落葉松天然林和人工林地表火平均蔓延速度分別降低95.3%和97.6%,平均蔓延速度都低于0.1 m·min-1; 平均火焰高度分別降低93.1%和93.9%,都低于0.1 m。興安落葉松白樺混交林和白樺林地表火平均蔓延速度為0.1 m·min-1,分別降低85.7%和88.9%; 平均火焰高度分別降低92.6%和87.6%。樟子松人工林地表火平均蔓延速度、火焰高度及火強度分別降低77.6%、87.3%和93.4%。

圖3 不同強度可燃物處理后的地表火行為Fig.3 Surface fire behavior with different fuel treatments不同字母代表同一林型不同處理強度間在0.05水平上差異顯著。Different letters indicate significant difference at 0.05 level among the different fuel treatment intensity in the same forest.
干旱(D2L2)情景下,5種林型之間潛在地表火行為(SFBP)指數差異明顯(F檢驗,P<0.01)(圖4A)。興安落葉松天然林平均地表火行為指數為6.3(3.5~9.0),顯著高于其他4種林分。興安落葉松人工林平均地表火行為指數最小,為3.3(1.9~4.8)。興安落葉松白樺混交林和樟子松人工林平均地表火行為指數均為4,但混交林平均地表火強度(RP)指數(5.6)高于樟子松人工林(5.2),蔓延速度(SP)指數(3.9)低于樟子松林(4.1)。白樺林平均地表火行為指數為4.3(2.9~5.8)。
白樺林內無樹冠火發生可能,而其他4種林型之間樹冠火行為(CFBP)指數差異明顯(F檢驗,P=0.00)(圖4B)。平均樹冠火行為指數最大出現在興安落葉松天然林(5),其中平均樹冠火發生指數(IC)和蔓延速度指數(RC)分別為4.5(3.4~5.6)和4.5(2.5~6.5); 最小值出現在興安落葉松人工林(3.3),平均樹冠火發生指數和蔓延速度指數分別為2.2(1.4~3.1)和2.5(0.9~4.1)。興安落葉松白樺混交林和樟子松人工林平均樹冠火行為指數均為4,而平均樹冠火發生指數和蔓延速度指數分別為4.1(2.7~5.5)、2.4(2.0~2.7)和3.1(2.9~3.2)、3.0(2.2~3.8)。4種林型樹冠火蔓延指數(Tc)差異不明顯(F檢驗,P=1.00)。
可燃物處理可有效降低干旱情景下各林型的潛在地表火和樹冠火行為指數,并隨著可燃物處理強度的增加而明顯(F檢驗,P<0.01)(圖5)。在低強度可燃物處理情景下,白樺林平均地表火行為指數降低7.7%。興安落葉松白樺混交林平均潛在地表火指數無明顯降低,但平均地表火蔓延速度指數、火焰長度指數和火強度指數分別降低5.1%、6.9%和12%。興安落葉松天然林和人工林平均地表火行為指數分別降低26.3%和20%。樟子松人工林潛在地表火行為指數無明顯變化,但平均地表火蔓延速度指數、火焰長度指數和火強度指數分別降低0.8%、4.9%和9%。興安落葉松天然林和人工林、樟子松人工林潛在樹冠火行為指數分別降低6.7%、10%和8.3%; 雖然混交林潛在樹冠火行為指數未明顯降低,但樹冠火發生指數和蔓延速度指數分別降低4.3%和17%。
中強度可燃物處理后,興安落葉松白樺混交林、白樺林、興安落葉松天然林和人工林地表火行為指數分別降低16.7%、23.1%、52.6%和60%。而樟子松人工林地表火行為指數仍未明顯變化,但平均地表火蔓延速度指數、火焰長度指數和火強度指數分別降低31.6%、8.9%和20.8%。興安落葉松白樺混交林、落葉松天然林、樟子松人工林樹冠火行為指數分別降低16.7%、26.7%和16.7%; 興安落葉松人工林樹冠火行為指數降低10%,雖與低強度處理情景下無顯著差異,但樹冠火發生指數和蔓延速度指數降幅增加,分別比未處理情景下降低41.8%和6.4%。
高強度可燃物處理后,興安落葉松白樺混交林、白樺林、興安落葉松天然林和人工林、樟子松人工林地表火行為指數分別降低75%、61.5%、75.8%、76%和50%。其中地表火蔓延速度指數、火焰長度指數和火強度指數降幅均超過60%。興安落葉松白樺混交林、落葉松天然林和人工林、樟子松人工林樹冠火行為指數分別降低41.7%、40%、30%和25%,其中樹冠火發生指數降幅超過60%,樹冠火蔓延指數降幅超過30%。
但可燃物處理并未降低興安落葉松白樺混交林、興安落葉松天然林和人工林、樟子松人工林4種林分的樹冠火蔓延指數(F檢驗,P=1.00)。
FCCS模型以世界各地區典型林分的可燃物床層信息作為模板,相關參數的調整主要基于區域內可燃物特征信息。該系統采用Rothermel公式計算火行為(Rothermel, 1972)。當前,FCCS已廣泛應用于火險評估、火行為預測及研究可燃物消耗等方面(Pettinarietal., 2016),在我國云南、廣西等地也得到了較好的使用(黃小榮等,2014; 馬振宇等,2020)。這些研究結果表明,在FCCS中選擇與研究地點生態區相近、植被類型相同的可燃物床層,利用標準地調查數據進行相關參數的調整,可以準確地描述研究區主要林型的可燃物特征及分布,對于林分火行為的預測貼合實際,在實際林火管理工作中具有一定指導意義。

圖4 D2L2情景下不同林型潛在火行為Fig.4 Potential fire behavior of each forest type under D2L2 scenario不同小寫字母代表同一火行為指標不同林型間在0.05水平上差異顯著。Different small letters indicate significant difference at 0.05 level among the different forest types in the same fire behavior index.

圖5 不同強度可燃物處理與潛在火行為Fig.5 The potential fire behavior with different fuel treatments不同大寫字母代表同一林型不同處理強度間在0.05水平上差異顯著。Different capital letters indicate significant difference at 0.05 level among the different fuel treatment intensity in the same forest.
本文模擬的環境條件包括一般天氣條件(用戶定制的)和D2L2干旱條件。其中在用戶定制的環境情景下,由于林分結構、天氣及地形等條件都會對可燃物含水率造成影響(胡海清等,2017),無法進行不同林型之間的火行為比較,而各林型模擬結果也只反映特定天氣條件下的地表火行為。而在D2L2情景下(環境變量相同),可以比較各種林分在相同天氣情景下的潛在地表火和樹冠火行為。這也是各林型在高火險天氣條件下的潛在火行為。興安落葉松天然林地表火和樹冠火行為指數最高,興安落葉松人工林地表火和樹冠火行為指數最低,而白樺林無樹冠火發生可能。
可燃物處理是大興安嶺地區主要的林火管理措施,現有研究多基于景觀尺度和野外試驗來研究可燃物處理對火發生的影響(劉志華等,2009; 高仲亮等,2015),但這些研究并未考慮林內可燃物對火行為的影響。
一般天氣情景下的模擬結果表明,進行中強度可燃物處理時,各林型地表火蔓延速度降低至1 m·min-1以下,火焰高度低于1 m,這種低強度的地表火可以很快被初始攻擊所撲滅。干旱(D2L2)情景下模擬結果也表明,中強度可燃物處理后,各林型地表火行為降幅明顯,興安落葉松天然林和人工林、興安落葉松白樺混交林和樟子松人工林的樹冠火行為也明顯降低。本次試驗主要對林冠層下的可燃物進行處理,未進行林冠層的處理,因此,模擬結果中的樹冠火蔓延可能性未發生明顯變化。
雖然經高強度可燃物處理后,火不能持續燃燒及蔓延,但清光林冠下易燃可燃物改變了林分結構,影響其水源涵養和維持生物多樣性等功能(李紅振等,2014; 陳百靈等,2017),對生態環境有較大影響。因此,在塊狀處理可燃物或進行可燃物阻隔帶清理時,采用中強度的處理方式比較合適,不但可顯著降低火行為,而且基本不影響生態環境。對重點防火區,可通過間伐改變林冠層的可燃物連續性,降低樹冠火風險(Polletetal., 2002)。
相同天氣情景下,在大興安嶺地區的5種典型林分類型中,興安落葉松天然林的地表火和樹冠火行為指標最高,興安落葉松人工林地表火行為指標最低。樟子松人工林和興安落葉松白樺混交林的火行為指標相近。白樺林地表火行為指標較低,且無樹冠火發生可能。大興安嶺的可燃物管理應優先處理興安落葉松天然林。
可燃物的低強度處理對興安落葉松天然林的火行為效果明顯,對其他林分不明顯。中強度處理后,5種林型的地表火都易于直接撲救。高強度處理后,火不能在林內持續燃燒,但對生態環境影響大。建議在林火管理活動中采用中強度的可燃物處理,可以滿足林火管理需要。