基于光能利用率模型估算大興安嶺北部火燒跡地生產(chǎn)力1)

胡海清 蘇志杰 孫 龍 魏書精

(東北林業(yè)大學(xué)，哈爾濱，150040)

森林作為陸地生態(tài)系統(tǒng)的主體［1］，占陸地面積的27%，以及全球80%以上的地上碳儲量和40%以上的全球土壤碳儲量［2］，對于維持全球碳平衡起著不可替代的作用［3］。森林火災(zāi)是森林生態(tài)系統(tǒng)的重要干擾因子，是導(dǎo)致植被和土壤碳儲量發(fā)生變化的重要原因［4］。植被凈第一性生產(chǎn)力(Net primary productivity，NP，P)是生態(tài)系統(tǒng)中植被碳固定與碳釋放之差，是大氣中的碳進(jìn)入生物群落的主要途徑?；馃E地植被NP，P的變化直接影響該區(qū)域碳平衡和碳循環(huán)，因此，了解其變化規(guī)律對理解火干擾對森林碳平衡的影響有重要意義，也是目前國外對火燒跡地植被恢復(fù)的研究熱點。我國林火研究起步較晚，但發(fā)展較快，許多學(xué)者在森林可燃物、林火規(guī)律、林火碳釋放、林火對森林生態(tài)的影響方面進(jìn)行了卓有成效的研究［5-7］。隨著遙感的應(yīng)用，以火干擾后NP，P的變化研究森林恢復(fù)已成為新的熱點，我國在這方面還幾近空白。本文以我國林火發(fā)生典型地區(qū)大興安嶺北部為例，對火燒跡地森林分布和NP，P進(jìn)行定量分析，旨在為火干擾后森林恢復(fù)的研究提供參考數(shù)據(jù)。

1 研究區(qū)概況

研究區(qū)域位于大興安嶺林區(qū)北部(52°20' ～53°30'N，121°40' ～124°20'E)，是1987年“5·6”大火主要過火區(qū)域，面積16 799 km2(圖1)。

圖1 研究區(qū)域位置

該區(qū)海拔300 ～1 400 m，坡度一般小于15°，為寒溫帶季風(fēng)區(qū)，最低溫度-52.3 ℃，最高氣溫39.0℃。年降水量350 ～500 mm，降水集中于7—8月份。相對濕度70% ～75%，積雪期達(dá)5 個月，林內(nèi)雪深30 ～50 cm。區(qū)內(nèi)廣泛分布著大片連續(xù)多年凍土和島狀多年凍土，土壤以棕色針葉林土和暗棕壤為主，植被類型簡單，優(yōu)勢樹種是興安落葉松(Larix gmelinii)、樟子松(Pinus sylvestris)和白樺(Betula platyphylla)。興安落葉松是該地區(qū)最穩(wěn)定的氣候演替頂級群落，樟子松為干旱瘠薄生境上的土壤演替頂級群落，白樺往往以先鋒樹種在森林火燒跡地和采伐跡地形成天然更新次生林［8］。該區(qū)為我國森林火災(zāi)高發(fā)區(qū)，年均森林過火面積居全國之首，是我國森林火災(zāi)危害最嚴(yán)重的地區(qū)。

2 研究方法與數(shù)據(jù)來源

2.1 遙感數(shù)據(jù)

本研究所利用的ETM +影像由美國NASA 陸地衛(wèi)星計劃Landsat7 采集，來源于中國科學(xué)院對地觀測與數(shù)字地球科學(xué)中心，圖像空間分辨率為30 m×30 m，時間為1999年8月30日。Landsat7 攜帶的傳感器在成像時受到太陽位置、大氣氣溶膠、傳感器自身特性、地形起伏等因素的影響，使遙感圖像產(chǎn)生了幾何畸變和輻射失真，從而降低了圖像的質(zhì)量，影響了圖像應(yīng)用和價值，因此對數(shù)據(jù)進(jìn)行了幾何校正、大氣校正、去云處理以及傳感器訂正等預(yù)處理。

2.2 氣象數(shù)據(jù)和植被分類數(shù)據(jù)

本研究所采用的氣象數(shù)據(jù)為1999年大興安嶺地區(qū)氣候資料日值數(shù)據(jù)集和輻射日值數(shù)據(jù)集，來源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cbc. cma.gov.cn/)。利用GIS 對氣象數(shù)據(jù)進(jìn)行了柵格化，并從空間上與遙感數(shù)據(jù)相匹配。

植被分類數(shù)據(jù)來源于黑龍江省二級土地覆蓋分類圖，由中國科學(xué)院編制。在此基礎(chǔ)上，本研究用較成熟的最大似然分類器將土地覆蓋類型分為針葉林、闊葉林、針闊混交林和其他4 種類型。其他類型是非林地，主要為沼澤、荒草地等。經(jīng)檢驗，分類精度達(dá)到70%以上。

2.3 火災(zāi)等級圖

1987年火災(zāi)等級圖為火后實地調(diào)查并結(jié)合遙感資料形成的手繪圖，火災(zāi)等級劃分為3 類:輕度火燒，燒死木蓄積量＜林分蓄積量30%;中度火燒，燒死木蓄積量占林分蓄積量31% ～60%;中度火燒，燒死木蓄積量＞林分蓄積量60%。在Arcgis 9.2 中利用已校正的ETM+影像對火災(zāi)等級圖進(jìn)行了配準(zhǔn)和數(shù)字化，裁剪得到研究區(qū)火災(zāi)等級分布圖(圖2)。

2.4 光能利用率模型

自20 世紀(jì)70年代初，Leith 首次提出估算全球NP，P值［9］。NP，P研究經(jīng)歷了站點實測、統(tǒng)計回歸及模型估算3 個研究階段。通過參數(shù)的選取可以利用模型估測區(qū)域或全球尺度的NP，P，目前已成為NP，P研究的主流方向。光能利用率模型以遙感數(shù)據(jù)、植被分類數(shù)據(jù)和氣象數(shù)據(jù)為基礎(chǔ)，利用植被吸收的光合有效輻射(Absorbed photo-synthetically active radiation，AP，A，R)和光能利用率(ε)來估算植被NP，P［10］。光能利用率模型基于植被的生理生態(tài)過程，在全球和區(qū)域NP，P研究及全球碳循環(huán)、碳平衡研究中被廣泛應(yīng)用［11］。國內(nèi)研究學(xué)者將模型應(yīng)用在我國時，對光能利用率模型進(jìn)行了適當(dāng)改進(jìn)，模擬結(jié)果得到了較好的驗證［12-14］。光能利用率模型的參數(shù)少，遙感數(shù)據(jù)主要由光合有效輻射和光能利用率這兩個變量驅(qū)動，避免了由于參數(shù)缺乏而人為簡化或者估計所帶來的誤差。模型算法為:

式中:x 表示空間位置;t 表示時間;AP，A，R(x，t)表示t月像元x 所吸收的光合有效輻射;ε(x，t)表示t月像元x 實際光能利用率。

圖2 數(shù)字化火災(zāi)等級圖

AP，A，R(x，t)可以根據(jù)遙感數(shù)據(jù)中植被對紅外和近紅外波段的反射特征來實現(xiàn)。光合有效輻射(photosynthetically active radiation，PA，R)是植物光合作用的驅(qū)動力，它與植物生長直接相關(guān)。植被吸收的光合有效輻射由太陽總輻射和植物特性決定，如式(2):

式中:SO，L(x，t)表示t月在像元x 處的太陽總輻射;FP，A，R(x，t)為植被層對入射光合有效輻射的吸收比例;常數(shù)0.5 表示植被所能利用的太陽總輻射的比例。

FP，A，R(x，t)與歸一化植被指數(shù)(ND，V，I)、比值植被指數(shù)(VI)有很好的線性關(guān)系。研究發(fā)現(xiàn):FP，A，R(x，t)可以利用遙感資料提取的歸一化植被指數(shù)(ND，V，I)和比值植被指數(shù)(VI)分別來估算［15］。但利用ND，V，I估 算 的FP，A，R(x，t)偏 高，而VI估 算 的FP，A，R(x，t)偏低，由兩者估算的AP，A，R(x，t)平均值效果最好。

光能利用率是在一定時期單位面積上生產(chǎn)的干物質(zhì)中所包含的化學(xué)潛能與同一時刻投射到該面積上的光合有效輻射能之比。在現(xiàn)實條件下，光能利用率ε 受溫度、水分及最大光能利用率εmax的影響:

式中:Tε1(x，t)和Tε2(x，t)表示溫度對光能轉(zhuǎn)化率的脅迫作用;Wε(x，t)為水分脅迫作用系數(shù)［10，16］。

Tε1和Tε2利用氣象資料計算，方法見文獻(xiàn)［12］。Wε(x，t)的計算涉及到土壤水分子模型，比較復(fù)雜，數(shù)據(jù)不易獲取。本文采用改進(jìn)的方法，利用實際蒸散量與區(qū)域潛在蒸散量的比值來反映水分對NP，P的影響，詳細(xì)方法見文獻(xiàn)［13］、［17］、［18］。月最大光能利用率εmax與植被類型有關(guān)［19］，本研究采用Zhu等［20］模擬的各植被類型的最大光能利用率。

3 結(jié)果與分析

3.1 不同強(qiáng)度過火區(qū)森林類型分布

不同火強(qiáng)度后林型分布如表1所示。針葉林在未過火區(qū)以及輕、中、重度過火區(qū)的比重分別為40%、36.8%、31.6%、20.8%，并隨火燒強(qiáng)度的增加而降低，表明針葉林的分布與火燒強(qiáng)度密切相關(guān)。輕度過火區(qū)針葉林的比重接近未過火區(qū)，一方面是該地區(qū)針葉林中主要為興安落葉松林，對輕度火燒有一定的抵抗性，因此輕度火燒對針葉林造成的損失很小。另一方面，輕度火燒后，利于興安落葉松種子的萌發(fā)，使得過火后針葉林恢復(fù)迅速。重度過火區(qū)針葉林的比重為20. 8%，相比未過火區(qū)針葉林40%的比重降幅較大。重度火燒造成了70%以上的樹木死亡，林下枯枝落葉層全部被燒毀，埋藏在土壤表層的種子大部分喪失萌發(fā)能力?；馃箅m然采取了人工恢復(fù)措施，但是由于立地條件的差異，人工恢復(fù)措施效果一般［21］，針葉林的生長恢復(fù)最差。中度過火區(qū)針葉林恢復(fù)介于輕度過火區(qū)和重度過火區(qū)之間。

表1 1999年8月不同火強(qiáng)度過火區(qū)林型分布 %

闊葉林在輕、中、重度火燒區(qū)的比重分別為7.9%、13.6%、8.8%，相對未過火區(qū)17.0%的比重均有降低?；旖涣值谋戎卦谖催^火區(qū)為11.4%，在輕、中、重度過火區(qū)的比重分別為13.0%、11.1%、17.4%，除中度過火區(qū)與未過火區(qū)比重接近外，其它都有增加?；馃髥棠颈淮罅繜溃馃E地上透光率增加，陽光充足，白樺、山楊等先鋒樹種迅速生長，形成闊葉林或與未燒死的針葉林形成針闊混交林。隨著林內(nèi)喬木均勻度的增加，土壤養(yǎng)分在短時間內(nèi)無法補(bǔ)充白樺、山楊等速生樹種的需求，在土壤養(yǎng)分多的地方喬木分布比較集中，從而又使喬木層的生長受到限制。十幾年林齡的闊葉林出現(xiàn)了大量的枯爛現(xiàn)象，使適合本地生長的針葉林有了更好的生長環(huán)境，這是混交林比重增加的原因。

過火區(qū)林地的NP，P相比未過火區(qū)明顯下降，與此對應(yīng)的是非林地的比重明顯增加。在火燒前處于沼澤邊緣和水分條件優(yōu)越的河谷地帶的林地，由于立地條件較差，林木生長不良，火燒后立地條件進(jìn)一步惡化，很難重新形成林地，導(dǎo)致林地對降水的截留量和蒸散量減少。而該地區(qū)土壤下層普遍存在永凍層，水分不能下移，使過火區(qū)地表積水增加，土壤含水量增大，土壤沼澤化增加［22］，而隨著火后地表徑流量的增加，逐漸發(fā)育成沼澤。在陽坡和陡坡地帶，火燒使土壤失去了植被的保護(hù)作用，造成了水土流失和土壤干旱化、貧瘠化。貧瘠的土壤更不利于植被的恢復(fù)，使森林發(fā)生了逆行演替。在火后12 a時，過火區(qū)非林地面積比未過火區(qū)提高了20%多，這也可能是造成闊葉林在過火區(qū)的比重下降、針闊混交林的比重較低的原因。

3.2 不同強(qiáng)度過火區(qū)NP，P

表2所示為火干擾12 a 后不同強(qiáng)度過火區(qū)針葉林、闊葉林、針闊混交林的NP，P。輕度和中度過火區(qū)針葉林NP，P超過了未過火區(qū)針葉林。輕度和中度火燒只燒毀了灌草層、喬木層中的小徑木和幼林，保留了很多針葉林母樹，為過火區(qū)針葉林生產(chǎn)力提供了基礎(chǔ)。而火災(zāi)后過火區(qū)林窗加大，充足的陽光促進(jìn)了林下植被生長恢復(fù)。重度火燒使大部分針葉林母樹死亡，其自然恢復(fù)情況要比輕度和中度過火區(qū)差，但火后進(jìn)行了大量的人工恢復(fù)措施，使植被恢復(fù)速度大大增快。到1999年8月，重度過火區(qū)針葉林達(dá)到未過火區(qū)NP，P的99%，植被生長恢復(fù)情況良好。

表2 1999年8月不同強(qiáng)度過火區(qū)的NP，P g·m -2

闊葉林耐火性較差，但輕度火燒由于強(qiáng)度較低，只燒毀了些幼樹，大徑木都得以保留，火后林內(nèi)立木間空間較大，土壤養(yǎng)分供應(yīng)充足，生長情況比同時期未過火區(qū)闊葉林好。中度和重度過火區(qū)喬木層死亡很多，火后白樺和山楊等快速生長，并在演替層成為優(yōu)勢樹種。火干擾12 a 時，其生長情況與未過火區(qū)基本相同，比輕度過火區(qū)差。

輕度和中度火燒對針闊混交林的影響較小，雖然灌草層全部被燒毀，但喬木層針葉樹和闊葉樹大徑木均保留了下來，因此其林內(nèi)生長環(huán)境要好于未過火區(qū)，火后12 a 時輕度和中度過火區(qū)NP，P超過了未過火區(qū)，重度過區(qū)NP，P比未過火區(qū)低。重度過火區(qū)只保留下來一部分針葉樹，林內(nèi)空地上白樺和山楊快速萌生，鑲嵌在針葉樹間，形成了大量針闊混交林，并很快郁閉。由于闊葉樹受制于養(yǎng)分和空間，生長放緩，因此重度過火區(qū)針闊混交林NP，P比未過火區(qū)低，但其生長恢復(fù)也很快。

4 結(jié)論

通過火后森林分布和NP，P的研究發(fā)現(xiàn)，在1987年“5·6”大火12 a 后，與未過火區(qū)相比，輕、中、重度過火區(qū)林地的比重下降，且隨火強(qiáng)度的增加而降低?；鸶蓴_造成了部分森林的逆行演替，以重度過火區(qū)最嚴(yán)重。林地內(nèi)輕度和中度過火區(qū)NP，P超過了未過火區(qū)，重度過火區(qū)林地NP，P雖低于未過火區(qū)，但已接近未過火區(qū)NP，P的水平。

本文研究的是火后植被覆蓋和NP，P的變化狀況，未探討立地條件和林齡對火后NP，P的影響，這是以后需要進(jìn)一步研究的問題。

［1］ 魏書精，胡海清，孫龍，等.氣候變化背景下我國森林防火工作的形勢及對策［J］.森林防火，2011(2):1 -4.

［2］ 王邵軍，阮宏華.全球變化背景下森林生態(tài)系統(tǒng)碳循環(huán)及其管理［J］. 南京林業(yè)大學(xué)學(xué)報:自然科學(xué)版，2011，35(2):113 -116.

［3］ IPCC. Land-use，land-use change and Forestry:A special report of the IPCC［M］. Cambridge:Cambridge University Press，2000.

［4］ 胡海清，魏書精，金森，等.森林火災(zāi)碳排放計量模型研究進(jìn)展［J］.應(yīng)用生態(tài)學(xué)報，2012，23(5):1423 -1434.

［5］ 胡海清，魏書精，孫龍.1965 -2010年大興安嶺森林火災(zāi)碳排放的估算研究［J］.植物生態(tài)學(xué)報，2012，36(7):629 -644.

［6］ 段建田，魏世強(qiáng)，于萍萍，等.林火跡地森林生態(tài)系統(tǒng)恢復(fù)研究［J］.安徽農(nóng)業(yè)科學(xué)，2007，35(1):182 -184.

［7］ 王緒高，李秀珍，孔繁花，等.大興安嶺北坡火燒跡地自然與人工干預(yù)下的植被恢復(fù)模式初探［J］.生態(tài)學(xué)雜志，2003，22(5):30 -34.

［8］ 梁鳳仙，顧鐘偉.遙感技術(shù)在大興安嶺火燒跡地森林凍土環(huán)境變化調(diào)查中的應(yīng)用［J］.冰川凍土，1993，15(1):27 -33.

［9］ Leith H，Whittaker R H. Modeling the primary productivity of the world［M］. Primary Productivity of the Biosphere，New York:Springer Verlag，1975:237 -263.

［10］ Potter C S，Randerson J T，F(xiàn)ield C B，et al. Terrestrial ecosystem production:a process model based on globe satellite and surface data［J］. Global Biogeochemical Cycles，1993，7(4):811 -841.

［11］ Ichii K，Matsui Y，Yamaguchi Y，et al. Comparison of global net primary production trends obtained from satellite-based normalized difference vegetation index and carbon cycle model［J］.Global Biogeochemical Cycles，2001，15(2):351 -363.

［12］ 樸世龍，方精云，郭慶華.利用CASA 模型估算我國植被凈第一性生產(chǎn)力［J］.植物生態(tài)學(xué)報，2001，25(5):603 -608.

［13］ 朱文泉，潘耀忠，張錦水.中國陸地植被凈初級生產(chǎn)力遙感估算［J］.植物生態(tài)學(xué)報，2007，31(3):413 -424.

［14］ 高志強(qiáng)，劉紀(jì)遠(yuǎn). 中國植被凈生產(chǎn)力的比較研究［J］. 科學(xué)通報，2008，53(3):317 -326.

［15］ Los S O，Justice C O，Tucker C J. A global by NDVI data set for climate studies derived from the GIMMS continental NDVI data［J］. International Journal of Remote Sensing，1994，15(17):3493 -3518.

［16］ Field C B，Randerson J T，Malmstrom C M. Global net primary production:combining ecology and remote sensing［J］. Remote Sensing of Environment，1995，51(1):74 -88.

［17］ 周廣勝，張新時.全球變化的中國氣候—植被分類研究［J］.植物學(xué)報，1996，38(1):8 -17.

［18］ 周廣勝，張新時.全球變化的中國植被的凈第一性生產(chǎn)力研究［J］.植物生態(tài)學(xué)報，1996，20(1):11 -19.

［19］ Ruimy A，Saugier B，Dedieu G. Methodology for the estimation of terrestrial net primary production from remotely sensed data［J］. Journal of Geophysical Research，1994，99(D3):5263 -5283.

［20］ Zhu Wenquan，Pan Yaoxhong，He Hao，et al. Simulation of maximum light use efficiency for some typical vegetation types in China［J］. Chinese Science Bulletin，51(4):457 -463.

［21］ 解伏菊，肖篤寧，李秀珍，等.基于NDVI 的不同火燒強(qiáng)度下大興安嶺林火跡地森林景觀恢復(fù)［J］. 生態(tài)學(xué)雜志，2005，24(4):368 -372.

［22］ 解伏菊，肖篤寧，李秀珍，等. 大興安嶺西北坡火燒跡地濕地格局變化及其影響因子研究［J］.水科學(xué)進(jìn)展，2007，18(2):176 -181.