武廣源,田玲玲,張運(yùn)林,劉 璐
(1. 東北林業(yè)大學(xué) 林學(xué)院,黑龍江 哈爾濱 150040;2. 貴州師范學(xué)院 a. 生物科學(xué)學(xué)院;b. 貴州高等學(xué)校林火生態(tài)與管理重點(diǎn)實(shí)驗(yàn)室,貴州 貴陽 550018;3.國家林業(yè)和草原局 林草調(diào)查規(guī)劃院,北京 100013)
對(duì)森林火險(xiǎn)預(yù)測(cè)具有重要意義的參數(shù)為凋落物含水率,其影響了凋落物被點(diǎn)燃的難易程度和點(diǎn)燃后可能表現(xiàn)的火行為特征[1-3]。得到最準(zhǔn)確含水率的方法為烘干法,該方法需要烘干凋落物,不能得到實(shí)時(shí)含水率值,在實(shí)際中無法應(yīng)用[4-5]。
目前國內(nèi)外使用最廣的森林凋落物含水率預(yù)測(cè)模型為氣象要素回歸模型和直接估計(jì)模型,本研究也是基于這兩種模型分析步長對(duì)預(yù)測(cè)精度的影響。其中氣象要素回歸模型是通過分析氣象要素對(duì)凋落物含水率動(dòng)態(tài)變化的影響,進(jìn)而建立兩者之間的關(guān)聯(lián)關(guān)系,但該方法受地形、凋落物類型和氣候區(qū)等的影響較大[6-7];直接估計(jì)法是以能夠表示物理機(jī)理的水分?jǐn)U散模型為主體方程,其中通過實(shí)驗(yàn)獲取模型參數(shù),該方法既具有統(tǒng)計(jì)學(xué)方法的簡潔性,又具有物理法的適用性[8-9]。雖然這兩種方法應(yīng)用最廣,但也存在一定的誤差,有些情況下的誤差甚至能導(dǎo)致火險(xiǎn)等級(jí)與實(shí)際情況相差1 ~2 級(jí)[10]。
對(duì)于特定類型的凋落物,其含水率預(yù)測(cè)模型產(chǎn)生誤差的原因主要是由于對(duì)凋落物含水率日變化過程描述過于簡單造成的。一天中空氣溫度和濕度變化過程相對(duì)穩(wěn)定,而森林凋落物含水率動(dòng)態(tài)變化受空氣溫濕度的驅(qū)動(dòng),因此凋落物含水率動(dòng)態(tài)變化也表現(xiàn)出穩(wěn)定的日變化過程。而實(shí)際中凋落物含水率在一天中表現(xiàn)更為復(fù)雜的動(dòng)態(tài)變化主要還受到風(fēng)速、降雨和其他氣象要素的影響。此前,由于受到氣象要素觀測(cè)技術(shù)的限制,只能獲取每天固定時(shí)間的氣象值,無法得到能夠反映氣象要素日變化的基本數(shù)據(jù),凋落物含水率也只能根據(jù)某一時(shí)刻的氣象要素值進(jìn)行預(yù)測(cè),技術(shù)上無法通過氣象要素預(yù)測(cè)凋落物含水率日變化情況。現(xiàn)有的森林火險(xiǎn)預(yù)報(bào)系統(tǒng)沿襲了過去的方法,凋落物含水率依舊是以天為步長進(jìn)行計(jì)算。然而,在火險(xiǎn)預(yù)測(cè)中卻出現(xiàn)了一定的誤差,這是因?yàn)闆]有考慮到凋落物每天含水率的變化過程。
隨著氣象要素觀測(cè)技術(shù)的快速發(fā)展,能夠以任意步長監(jiān)測(cè),林火研究學(xué)者也認(rèn)識(shí)到凋落物含水率對(duì)火發(fā)生的可能性和火行為模型的影響十分重要,特別是林火行為的預(yù)測(cè)。準(zhǔn)確地預(yù)測(cè)需要比天更短的時(shí)間尺度得到的含水率數(shù)據(jù),建立一種更短步長的含水率預(yù)測(cè)模型十分迫切。國內(nèi)外學(xué)者也進(jìn)行了廣泛研究,都集中在以小時(shí)為步長進(jìn)行分析,且認(rèn)為以時(shí)為步長能夠顯著提高含水率預(yù)測(cè)精度,并能夠滿足森林火險(xiǎn)預(yù)報(bào)需要[11-17]。
國內(nèi)外學(xué)者以時(shí)為步長建立預(yù)測(cè)模型,但這些研究大多數(shù)都是集中在分析不同凋落物含水率預(yù)測(cè)方法在進(jìn)行以時(shí)為步長預(yù)測(cè)時(shí)的適用性,預(yù)測(cè)步長都是1 h,這些研究都需要以1 h 或者更短的時(shí)間步長進(jìn)行野外監(jiān)測(cè)凋落物含水率,將會(huì)耗費(fèi)大量人力物力[18-21]。若能得到合適的步長進(jìn)行預(yù)測(cè),既能減少人力物力的消耗,也能達(dá)到火險(xiǎn)預(yù)測(cè)的需求,對(duì)于實(shí)際應(yīng)用具有重要意義。因此,本研究旨在分析凋落物含水率的多大步長時(shí)能夠較好地表現(xiàn)出日變化過程,并建立不同的預(yù)測(cè)模型,通過分析不同步長的預(yù)測(cè)效果,得到能夠滿足預(yù)報(bào)要求的最大步長。對(duì)于理解森林凋落物含水率日變化和預(yù)測(cè)精度研究具有現(xiàn)實(shí)意義。
研究區(qū)地處貴州省鳳岡縣天合山(107°43′30″~107°43′38″E, 27°44′53″~27°45′01″N),屬大婁山山脈,森林覆蓋率高,平均海拔720 m,年平均氣溫16 ℃,年均濕度63%,雨水多集中在6—8 月,氣候?qū)賮啛釒駶櫺约撅L(fēng)氣候。研究區(qū)喬木主要包括馬尾松Pinusmassoniana、白櫟Quercusfabri、香樟Cinnamomumcamphora、青岡Cyclobalanopsis glauca、杉木Cunninghamialanceolata等,灌木有海桐Pittosporumtobira、女貞Ligustrumlucidum和龜甲冬青Ilexchinensis等。
在研究區(qū)內(nèi)選擇具有代表性的白櫟林和馬尾松林,分別設(shè)置25.82 m×25.82 m 的標(biāo)準(zhǔn)地,白櫟林和馬尾松林標(biāo)準(zhǔn)地的詳細(xì)信息下表所示。貴州省防火期為當(dāng)年10 月至次年5 月,由于凋落物在第2 年經(jīng)過風(fēng)化后的水分動(dòng)態(tài)變化對(duì)氣象要素的響應(yīng)不同,而研究區(qū)在2—4 月森林火災(zāi)更為嚴(yán)重,因此本研究于3 月份進(jìn)行凋落物含水率野外監(jiān)測(cè)試驗(yàn)。
表1 樣地基本信息Table 1 The basic information of plot
研究表明,破壞性采樣得到的含水率最準(zhǔn)確,因此本研究選擇破壞性采樣監(jiān)測(cè)凋落物含水率。在白櫟和馬尾松樣地內(nèi),分別隨機(jī)設(shè)3 個(gè)樣點(diǎn),每隔1 h 選擇破壞性采樣方法采集凋落物,現(xiàn)場(chǎng)稱重,記錄為Wh(濕質(zhì)量)。將其置于信封并帶回實(shí)驗(yàn)室,選擇鼓風(fēng)式干燥箱烘干至沒有水分為止,記錄此時(shí)質(zhì)量為Wd(干質(zhì)量),根據(jù)公式(1)計(jì)算含水率。若遇降雨天氣,則需要用紙將其表面自由水擦拭,并在信封外套封口袋,防止相互之間洇濕,影響結(jié)果。此次凋落物含水率監(jiān)測(cè)試驗(yàn)歷時(shí)7 d,各樣點(diǎn)數(shù)據(jù)為7 d×24 h =168 組,取3 個(gè)樣點(diǎn)的算數(shù)平均值,即該時(shí)刻該樣地的凋落物含水率。
式中:M為凋落物含水率(%);Wh和Wd分別為凋落物濕質(zhì)量和干質(zhì)量(g)。
HOBO氣象站架設(shè)在白櫟和馬尾松樣地之間,與采樣時(shí)間同步,以10 min 為間隔,對(duì)樣地內(nèi)的氣象要素進(jìn)行同步監(jiān)測(cè),主要包括空氣溫度、相對(duì)濕度、風(fēng)速和降水量。氣象要素的監(jiān)測(cè)時(shí)間要提前于含水率的監(jiān)測(cè)時(shí)間。
研究證明,只有在森林凋落物含水率低于35%時(shí),才有可能被引燃發(fā)生林火并蔓延為森林火災(zāi),因此本研究中所有數(shù)據(jù)分析均為含水率低于35%的。對(duì)于白櫟凋落物數(shù)據(jù),從3 月6 日8:00至3 月12 日10:00,共147 組數(shù)據(jù);對(duì)于馬尾松凋落物,從3 月6 日10:00 至3 月12 日9:00,共144 組數(shù)據(jù)。
1.4.1 數(shù)據(jù)基本統(tǒng)計(jì)
統(tǒng)計(jì)監(jiān)測(cè)期內(nèi)白櫟和馬尾松凋落物含水率基本情況,包括凋落物含水率最小值、最大值、平均值和標(biāo)準(zhǔn)差等指標(biāo)。
1.4.2 凋落物含水率日變化最大步長分析
分別繪制不同監(jiān)測(cè)步長時(shí)凋落物含水率的日變化圖,以小時(shí)為橫坐標(biāo),不同步長對(duì)應(yīng)的凋落物含水率值為縱坐標(biāo),確定能夠反映凋落物日動(dòng)態(tài)變化的最大步長。
為保證不同時(shí)間步長建立的含水率預(yù)測(cè)模型具有可比性,數(shù)據(jù)總量應(yīng)該保持一致。本研究共連續(xù)監(jiān)測(cè)野外含水率7 d,白櫟和馬尾松凋落物含水率低于35%的數(shù)據(jù)共有147 和144 組。通過研究得到能夠反映凋落物日變化的最大步長為n小時(shí),則白櫟和馬尾松以n小時(shí)為預(yù)測(cè)步長時(shí)分別共有147/n和144/n組數(shù)據(jù),因此之后的研究中,其余步長時(shí)也隨機(jī)選擇連續(xù)的147/n和144/n組數(shù)據(jù)進(jìn)行研究。
1.4.3 相關(guān)性分析
對(duì)森林凋落物含水率動(dòng)態(tài)變化進(jìn)行氣象要素影響分析。選擇斯皮爾斯相關(guān)性分析得到凋落物采集前x(x<10)小時(shí)的平均空氣溫度、相對(duì)濕度和風(fēng)速,以及前x小時(shí)的累計(jì)降雨量對(duì)不同步長時(shí)的凋落物含水率動(dòng)態(tài)變化的影響。
1.4.4 森林凋落物含水率預(yù)測(cè)模型介紹
1)直接估計(jì)法
Catchople 在2001 年提出了直接估計(jì)法,是以水分?jǐn)U散模型為主體方程,并以實(shí)際的凋落物含水率值和氣象數(shù)據(jù)估計(jì)模型參數(shù)值。采用直接估計(jì)模型時(shí)的重要參數(shù)為平衡含水率,本研究選擇被廣泛使用的Nelson 模型和Simard 模型,計(jì)算平衡含水率,并建立含水率預(yù)測(cè)模型。直接估計(jì)法的水分?jǐn)U散模型是基于Byram 提出的水分微分?jǐn)U散方程得到的,公式詳見參考文獻(xiàn)[22-24]。
2)氣象要素回歸法
以不同步長的凋落物含水率數(shù)據(jù)為因變量,每小時(shí)和前x小時(shí)氣象要素值為驅(qū)動(dòng)因子,選擇向后逐步回歸法得到不同步長的凋落物含水率氣象要素回歸模型。
1.4.5 確定最適預(yù)測(cè)步長
針對(duì)不同步長凋落物含水率的Nelson 模型、Simard 模型和氣象要素回歸模型,選擇k 折交叉驗(yàn)證法檢驗(yàn)?zāi)P途?,并得到模型的平均絕對(duì)誤差和平均相對(duì)誤差。1∶1 圖以實(shí)測(cè)值為橫軸,預(yù)測(cè)值為縱軸,分析凋落物含水率在不同區(qū)間時(shí)預(yù)測(cè)效果。分別不同模型,以不同步長凋落物含水率預(yù)測(cè)誤差(MAE 和MRE)為因變量,進(jìn)行差異性分析,研究不同步長時(shí)預(yù)測(cè)模型的誤差之間是否有顯著差異。得到最適預(yù)測(cè)步長。
表2 給出監(jiān)測(cè)期內(nèi)白櫟和馬尾松凋落物含水率基本情況??梢钥闯?,研究期內(nèi)白櫟凋落物含水率變化范圍為16.99%~30.90%,標(biāo)準(zhǔn)差為2.16%。馬尾松凋落物含水率最小值為11.11%,最大值為23.73%。整個(gè)研究期間白櫟和馬尾松凋落物的均值分別為21.91%和16.09%。
表2 凋落物含水率情況統(tǒng)計(jì)Table 2 Information of litter moisture content
圖1 給出研究期內(nèi)白櫟含水率不同步長時(shí)的日變化情況,可以看出,對(duì)于白櫟凋落物,當(dāng)預(yù)測(cè)步長超過4 h 時(shí),在一天中凋落物含水率就沒有表現(xiàn)出明顯的日變化,因此對(duì)于白櫟凋落物選擇研究的最大步長為4 h。
圖1 白櫟凋落物含水率日變化Fig. 1 Diurnal change of litter moisture content of Q. fabri
對(duì)于馬尾松凋落物,同樣也是當(dāng)?shù)蚵湮锖时O(jiān)測(cè)步長超過4 h 時(shí),其含水率值沒有表現(xiàn)出明顯的日變化。因此本研究中馬尾松凋落物含水率日變化的最大步長也為4 h(圖2)。
圖2 馬尾松凋落物含水率日變化Fig. 2 Diurnal change of litter moisture content of P. massoniana
圖3 給出不同步長的白櫟和馬尾松凋落物含水率動(dòng)態(tài)變化與氣象要素的相關(guān)系數(shù)值??梢钥闯觯徽摰蚵湮锉O(jiān)測(cè)步長如何改變,隨著采樣時(shí)間延長,空氣溫度和相對(duì)濕度與凋落物含水率的相關(guān)性先增加后下降。隨著空氣溫度的增加,白櫟凋落物含水率呈下降趨勢(shì),隨著相對(duì)濕度的增加含水率也呈現(xiàn)增加趨勢(shì),風(fēng)速對(duì)其幾乎沒有影響。馬尾松凋落物含水率隨空氣溫度的升高而升高,相對(duì)濕度對(duì)其的影響并不顯著,隨著風(fēng)速的增加,馬尾松凋落物含水率基本都在增加。
圖3 白櫟和馬尾松凋落物含水率與氣象要素相關(guān)系數(shù)Fig. 3 The correlation coefficient of moisture content of litter and meteorological elements
圖4 不同步長時(shí)白櫟凋落物預(yù)測(cè)模型1∶1 圖Fig. 4 Comparison of measured and predicted value of the prediction model of litter moisture content of Q. fabri
2.4.1 預(yù)測(cè)模型
表3 給出不同步長時(shí)白櫟和馬尾松凋落物的氣象要素回歸模型和預(yù)測(cè)誤差。由此可見,對(duì)于白櫟凋落物而言,氣象要素回歸模型在預(yù)測(cè)步長為3 h 時(shí)MAE 最低,只有0.96%,MRE 等于5.94%,而在預(yù)測(cè)步長為1 h 時(shí)誤差最大。對(duì)于馬尾松凋落物,同樣也是預(yù)測(cè)步長為3 h 時(shí),預(yù)測(cè)精度最好。
表3 不同步長時(shí)氣象要素回歸模型Table 3 Meteorological element regression model of litter moisture content under different prediction step
表4 給出不同步長時(shí)白櫟和馬尾松的直接估計(jì)法預(yù)測(cè)模型參數(shù)和預(yù)測(cè)誤差??梢钥闯?,對(duì)于白櫟凋落物,隨著步長增加,預(yù)測(cè)誤差增大,步長為4 h 時(shí)誤差最大。步長為1 h 時(shí),預(yù)測(cè)精度最高,模型MAE 僅為0.86%,MRE 為3.98%;對(duì)于馬尾松凋落物,預(yù)測(cè)步長為2 h 時(shí),預(yù)測(cè)精度最高,模型MAE 僅為0.38%,MRE 為3.17%。
表4 不同步長時(shí)直接估計(jì)法預(yù)測(cè)模型Table 4 Direct estimation prediction model of litter moisture content under different prediction step
2.4.2 1∶1 圖
對(duì)于白櫟凋落物,不論是直接估計(jì)模型還是氣象要素回歸模型,都是含水率值較低時(shí),預(yù)測(cè)值偏高,當(dāng)含水率實(shí)測(cè)值較高時(shí)會(huì)低估。對(duì)于直接估計(jì)法,以1 h 為步長時(shí)實(shí)測(cè)值和預(yù)測(cè)值擬合線最接近1∶1 線,并與其他步長有明顯區(qū)別,而2、3 和4 h 步長時(shí)的預(yù)測(cè)效果都接近。對(duì)于氣象要素回歸模型,3 h 步長時(shí)預(yù)測(cè)效果相對(duì)較好,但不同步長的預(yù)測(cè)效果都差別不大。
對(duì)于直接估計(jì)法,不同步長的馬尾松凋落物含水率預(yù)測(cè)模型都基本與1∶1 線重復(fù),預(yù)測(cè)效果都較好;對(duì)于氣象要素回歸法,所有步長時(shí)都是當(dāng)含水率較低時(shí)高估,而當(dāng)含水率較高時(shí)則低估,當(dāng)預(yù)測(cè)步長為2 h 時(shí)預(yù)測(cè)最接近1∶1 線(圖5)。
圖5 不同步長時(shí)馬尾松凋落物預(yù)測(cè)模型1∶1 圖Fig. 5 Comparison of measured and predicted value of the prediction model of litter moisture content of P. massoniana
2.4.3 不同步長時(shí)預(yù)測(cè)精度差異性分析
圖6 給出白櫟和馬尾松凋落物含水率氣象要素回歸模型不同步長時(shí)MAE 和MRE 的差異性分析結(jié)果??梢钥闯?,不論是白櫟還是馬尾松,氣象要素回歸模型的MAE 和MRE 在不同步長間都沒有顯著差異。
圖6 氣象要素回歸模型差異性檢驗(yàn)結(jié)果Fig. 6 Results of difference test of meteorological regression model
圖7 給出白櫟和馬尾松凋落物含水率直接估計(jì)模型不同步長時(shí)MAE 和MRE 的差異性分析結(jié)果??梢钥闯?,對(duì)于白櫟凋落物,直接估計(jì)模型的MAE 和MRE 在1 h 為步長時(shí)的顯著低于其他步長,其他步長間沒有顯著差異;對(duì)于馬尾松凋落物,直接估計(jì)模型的MAE 和MRE 在不同步長間卻沒有顯著差異。
圖7 直接估計(jì)模型差異性檢驗(yàn)結(jié)果Fig. 7 Results of difference test of direct estimation model
不論步長如何改變,白櫟和馬尾松凋落物含水率動(dòng)態(tài)變化對(duì)氣象要素的響應(yīng)并不相同,這主要是由于凋落物類型不同導(dǎo)致。凋落物類型不同,其床層結(jié)構(gòu)、理化性質(zhì)等都不同,對(duì)水分變化的響應(yīng)情況也不相同[25-26]。隨著采樣時(shí)間增加,不論步長如何改變,白櫟和馬尾松凋落物含水率動(dòng)態(tài)變化與氣象要素的相關(guān)性都呈先增加后下降的趨勢(shì)[27]。
采用直接估計(jì)法得到白櫟和馬尾松凋落物含水率不同步長時(shí)的預(yù)測(cè)模型,參數(shù)α的變化范圍分別為0.249~0.257和0.029~0.093,Catchpole等[12]研究得到的變化范圍為0.26~0.37,Slijepcevic等[3]研究得到值范圍為0.28 ~0.41;Sun 等[14]以興安落葉松林為研究對(duì)象,得到值為0.087 ~0.594。研究結(jié)果不同可能與凋落物類型、試驗(yàn)條件及試驗(yàn)設(shè)定時(shí)間段等有關(guān)系。該值不同可能是由于凋落物類型、監(jiān)測(cè)步長及監(jiān)測(cè)時(shí)間段等有關(guān)系。模型中β值反映的是凋落物含水率動(dòng)態(tài)變化對(duì)外界的響應(yīng)情況,其絕對(duì)值越大,表示含水率動(dòng)態(tài)變化更容易受到外界環(huán)境的影響[23]??梢钥吹?,不論是白櫟還是馬尾松凋落物,不同步長時(shí)的變化不大,說明該值也主要與凋落物自身?xiàng)l件有關(guān)系。馬尾松凋落物的值高于白櫟,說明馬尾松床層的持水能力要弱于白櫟床層,其含水率動(dòng)態(tài)變化更容易受到環(huán)境的影響,在實(shí)際應(yīng)用中需要更加注意馬尾松含水率的變化。
對(duì)于氣象要素回歸法而言,白櫟凋落物在不同步長時(shí)的含水率預(yù)測(cè)模型的MAE 變化幅度在0.96%~1.65%之間,MRE 變化幅度在5.94%~7.11%之間,兩者在不同步長之間無顯著差異。不同步長時(shí)馬尾松凋落物含水率預(yù)測(cè)模型的MAE 和MRE 變化范圍分別為1.01%~2.96%和10.23%~15.34%,不同步長之間都沒有顯著差異。研究結(jié)果精度都低于其他研究,主要是由于本研究步長在4 h 以內(nèi)所致,模型誤差均能滿足火險(xiǎn)預(yù)報(bào),但無法滿足火行為預(yù)報(bào)精度要求。
對(duì)于直接估計(jì)法而言,白櫟凋落物在不同步長時(shí)的含水率預(yù)測(cè)模型MAE 變化范圍在0.86%~1.47%之間,MRE 在3.98%~6.23%之間,1 h 步長時(shí)的預(yù)測(cè)精度要顯著低于其他步長,其余步長之間并沒有顯著差異。不同步長時(shí)馬尾松凋落物含水率預(yù)測(cè)模型的MAE 和MRE 變化范圍分別為0.38%~0.60%和3.17%~4.00%,不同步長之間沒有顯著差異。結(jié)果與Sun 等[14]以興安落葉松為研究對(duì)象得到的結(jié)果相似(0.41%~1.30%),與Catchpole 等[12]的研究結(jié)果也相似(0.8%~1.9%)。因此,對(duì)于白櫟凋落物,當(dāng)用于火險(xiǎn)預(yù)報(bào)時(shí),4 h 步長以內(nèi)可以滿足,而對(duì)于火行為預(yù)報(bào)則是1 h 步長更適合;對(duì)于馬尾松凋落物類型,2 h 步長完全能夠滿足火險(xiǎn)預(yù)報(bào)和火行為預(yù)報(bào)要求。
研究有一定的局限性,如針對(duì)特定凋落物類型,并未考慮林分特征的影響,不同林分特征、床層密實(shí)度、坡度和坡向等對(duì)凋落物含水率有顯著影響;此外,直接估計(jì)法中的平衡含水率預(yù)測(cè)模型選擇有限,且僅在春季防火期進(jìn)行研究等。在今后研究中,適當(dāng)增加研究范圍,對(duì)提高凋落物含水率預(yù)測(cè)模型的準(zhǔn)確性和實(shí)際應(yīng)用具有重要意義,如在各種林分和地形條件下,綜合考慮凋落物含水率日變化規(guī)律,擴(kuò)大研究時(shí)間尺度等。
研究表明,白櫟和馬尾松凋落物含水率動(dòng)態(tài)變化以4 h 以內(nèi)步長時(shí),能夠反映日變化過程。隨著采樣時(shí)間延長,空氣溫度和相對(duì)濕度與凋落物含水率的相關(guān)性先增加后下降,風(fēng)速對(duì)白櫟凋落物含水率的變化影響很小,與空氣溫度呈負(fù)相關(guān),與相對(duì)濕度呈正相關(guān)。馬尾松凋落物含水率隨空氣溫度的升高而升高,相對(duì)濕度對(duì)其的影響并不顯著,隨著風(fēng)速的增加,馬尾松凋落物含水率也增加。對(duì)于氣象要素回歸法,不論是白櫟還是馬尾松,以4 h 之內(nèi)為步長進(jìn)行含水率預(yù)測(cè)都能滿足預(yù)報(bào)要求;對(duì)于直接估計(jì)法,以1 h 為步長預(yù)測(cè)白櫟凋落物含水率效果更好,而馬尾松凋落物則2 h為步長進(jìn)行預(yù)測(cè)即可。