吳金卓 孔琳琳 董希斌 沈微

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

低質林改造更新?lián)嵊到y(tǒng)養(yǎng)分的循環(huán)模擬1)

吳金卓 孔琳琳 董希斌 沈微

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

依據(jù)低質林改造更新?lián)嵊到y(tǒng)的養(yǎng)分循環(huán)過程，采用內部集散物流分析、物質流分析方法，建立低質林改造更新?lián)嵊到y(tǒng)養(yǎng)分元素流動數(shù)學模型；應用實地調查數(shù)據(jù)對模型有效性驗證后，運用模型對研究區(qū)內低質林更新階段(2008—2015年)的營養(yǎng)元素積累與循環(huán)情況進行模擬分析。結果表明：通過營養(yǎng)元素流動模型得到的土壤養(yǎng)分指標，預測值與實測值之間的相關系數(shù)，水解N為0.87、有效P為0.77、速效K為0.88，顯示了較高的模型精度。2008—2015年，更新系統(tǒng)通過降雨過程從系統(tǒng)外吸收的N、P、K，分別為320、28、120kg/hm2；同時隨徑流流出系統(tǒng)的N、P、K，分別為224.0、19.6、84.0kg/hm2。林分從土壤中吸收的N、P、K，分別為26.45、4.29、11.11kg/hm2；通過淋溶作用，歸還到土壤中的N、P、K，分別為0.36、0.13、0.56kg/hm2。另外，幼齡林撫育措施及逐年累積的凋落物，回流至土壤中的N、P、K，分別為4.17、0.63、1.31kg/hm2；土壤中累積的營養(yǎng)元素N、P、K，分別達到289.08、14.27、124.75kg/hm2；比改造當年均有較大幅度的增加，表明低質林改造對于林地內的養(yǎng)分循環(huán)具有一定的促進作用。

低質林；低質林改造；低質林更新；林分養(yǎng)分

養(yǎng)分循環(huán)是森林生態(tài)系統(tǒng)分析中的重要內容之一。養(yǎng)分循環(huán)是非常復雜的生物過程，不僅受到環(huán)境因素、樹種生物學特性的影響，而且各種營養(yǎng)元素的循環(huán)也具有明顯的差異[1-2]。養(yǎng)分元素的循環(huán)與平衡，直接影響生產(chǎn)力的高低和生態(tài)系統(tǒng)的穩(wěn)定和持續(xù)[3-4]。研究森林生態(tài)系統(tǒng)的養(yǎng)分循環(huán)，不僅有助于了解物質的循環(huán)與積累機制，而且對于指導生產(chǎn)實踐、改善和調節(jié)各種新環(huán)境限制因素，提高生態(tài)系統(tǒng)的生產(chǎn)力，具有十分重要的作用。自上個世紀國際生物學計劃(IBP)實施以來，世界范圍內已經(jīng)開展了大量的森林生態(tài)系統(tǒng)養(yǎng)分循環(huán)研究，并積累了大量的數(shù)據(jù)和資料[5-9]；但是，對于低質林改造林分的養(yǎng)分循環(huán)與積累方面的研究很少。

小興安嶺林區(qū)是我國重要的用材林生產(chǎn)基地之一，然而，由于多年的過度采伐和自然因素的干擾，這一地區(qū)已經(jīng)形成了大面積郁閉度較低、林相衰敗的低質林[10-12]。為了更好地掌握低質林改造更新階段營養(yǎng)元素的積累與循環(huán)流動情況，本研究采用內部集散物流分析(BIF)、物質流分析(MFA)方法，建立營養(yǎng)元素流動數(shù)學模型，并結合小興安嶺地區(qū)低質林改造實地調查數(shù)據(jù)進行比較分析，旨在為定量分析低質林改造更新?lián)嵊到y(tǒng)內營養(yǎng)元素流動的數(shù)量、結構、特征及為低質林改造的經(jīng)營管理和效果評估提供參考。

1 研究區(qū)概況

研究區(qū)為小興安嶺林區(qū)鐵力林業(yè)局馬永順林場500林班內低質林改造帶。該區(qū)北部地勢平緩，南部地勢稍有斜坡，平均坡度10°，海拔117～284 m，土壤為暗棕壤，土層平均厚度45 cm。2007年，對研究區(qū)內的低質林進行水平帶狀采伐改造，沿同一海拔高度分別設置6 m×100 m、8 m×100 m、10 m×100 m、15 m×100 m的改造帶，伐除帶內的全部非經(jīng)營目的喬木和灌木。2008年，在皆伐帶內分3段栽植紅皮云杉(Piceakoraiensis)、興安落葉松(Larixgmelinii)、紅松(Pinuskoraiensis)幼苗，栽植比例為紅皮云杉∶興安落葉松∶紅松=6∶2∶2。

文獻[13]研究結果表明，在帶狀8 m×100 m、10 m×100 m改造帶內造林成活率及生長率較高，造林更新效果較好，且改造帶的設置可有效提高林木攔截降雨能力，利于林區(qū)內水土保持。因此，本文選擇10 m寬水平采伐改造帶(長100 m)為樣地，對2008—2015年低質林更新階段林地內的營養(yǎng)元素積累與循環(huán)進行模擬分析。

2 營養(yǎng)元素物質流模型的構建

低質林改造更新?lián)嵊到y(tǒng)的養(yǎng)分循環(huán)過程，包括林木的吸收、存留、凋落歸還、淋溶歸還、大氣降雨、徑流輸入和輸出等途徑。低質林區(qū)更新?lián)嵊镔|流分析系統(tǒng)如圖1所示。本研究在人工純林營養(yǎng)元素循環(huán)模型的基礎上[14]，考慮低質林更新階段中栽植不同種類苗木的情況，提出了適用于低質林更新系統(tǒng)的營養(yǎng)元素循環(huán)模型。

圖1 低質林更新?lián)嵊到y(tǒng)營養(yǎng)元素流動分析圖

2.1 營養(yǎng)元素存量計算

土壤中營養(yǎng)元素量：假設土壤中某種營養(yǎng)元素量用X1t表示(t=1、2、3、…、n)，低質林更新林分中更新樹種用s表示(s=1、2、3、…、m)，則第t年末土壤中該營養(yǎng)元素量為：

(1)

式中：F01t、F31t分別為第t年從降雨和枯枝落葉中分解輸入土壤中的某種營養(yǎng)元素量；F21ts為第t年第s個樹種經(jīng)淋溶作用進入土壤中的某種營養(yǎng)元素量；F12ts為第t年第s個樹種生長從土壤中吸收的某種營養(yǎng)元素量；F10t為第t年因徑流作用從土壤中流出的某種營養(yǎng)元素量；F001t為第t年土壤從施肥過程中吸收的某種營養(yǎng)元素量；F51t為撫育剩余物分解進入土壤中的營養(yǎng)元素量。

林木中營養(yǎng)元素量：在林木生長過程中，需要不斷地從土壤中吸收營養(yǎng)元素，并且在雨水的淋溶作用下將一部分營養(yǎng)元素歸還到土壤中；另外，還有一部分以凋落物的形式進入到枯枝落葉層[2]。林木中營養(yǎng)元素當前的量，是與上一年存留量、吸收量、淋溶輸出量、凋落物輸出量之間的動態(tài)平衡。第 個樹種中某種營養(yǎng)元素量，可以用公式(2)計算。當?shù)?年發(fā)生透光撫育或生長撫育活動時，一部分營養(yǎng)元素隨小徑材流出系統(tǒng)外，另一部分以撫育剩余物的形式進入枯枝落葉層，此時公式(2)變形為公式(3)。所有樹種中該種營養(yǎng)元素量，可用公式(4)計算。

X2ts=X2(t-1)s+F12ts-F21ts-F23ts，?t≥1、s；

(2)

X2ts=X2(t-1)s+F12ts-F21ts-F23ts-F24ts，?t≥1、s；

(3)

(4)

式中：X2ts為第t年末第s個樹種某種營養(yǎng)元素量；F23ts為第t年第s個樹種營養(yǎng)元素進入枯枝落葉層的某種營養(yǎng)元素量；F24ts為轉移到撫育作業(yè)中的營養(yǎng)元素量。

枯枝落葉層營養(yǎng)元素量：枯枝落葉層中的某種營養(yǎng)元素量，是在上一年留存量的基礎上，隨著枯枝落葉輸入量的增加而增加，隨著分解量的增加而減少[15]。它們之間的關系如公式(5)。

(5)

撫育作業(yè)中營養(yǎng)元素量：當?shù)趖年發(fā)生透光撫育或生長撫育活動時，轉移到撫育作業(yè)中的營養(yǎng)元素用公式(6)表示，其中轉移到系統(tǒng)外(F40t)和轉移到撫育剩余物中(F45t)的某種營養(yǎng)元素量分別用公式(7)和公式(8)表示。式中的p為撫育強度，a40、a45分別為轉移到系統(tǒng)外和轉移到剩余物中的比例，總和為1(公式(9))。撫育剩余物留在林地中，最終通過分解將營養(yǎng)元素回流到土壤中(公式(10))。

(6)

F40t=a40·X4t；

(7)

F45t=a45·X4t；

(8)

a40+a45=1；

(9)

X5t=F45t=F51t。

(10)

2.2 營養(yǎng)元素流計算

由于降雨的作用而進入到土壤中的營養(yǎng)元素，可根據(jù)大氣降雨中的營養(yǎng)元素確定。因此，可以令

F01t=S1t，?t≥1。

(11)

式中：S1t為第t年降雨中某種營養(yǎng)元素量。

由于低質林更新階段，主要對幼齡林進行扶正、施肥、修枝、除草等，本研究不考慮撫育作業(yè)對徑流的影響，僅考慮因降雨產(chǎn)生的系統(tǒng)徑流，可用公式(12)表示，公式中的R為施肥量隨徑流流出的比例。

F10t=a10·F01t+R·F001t。

(12)

林木在生長過程中從土壤中吸收的營養(yǎng)元素，與其自身的元素量和所處的生長階段相關。由于進行施肥作業(yè)，可認為土壤中營養(yǎng)元素可供低質林更新林分生長吸收。根據(jù)吸收比例系數(shù)(f12s(t))和自身的元素量可計算出第t年的吸收量，如公式(13)所示，式中的f12s(t)為隨生長時間變化的函數(shù)。

F12ts=f12s(t)·X2(t-1)s，?t≥1、s。

(13)

淋溶和凋落物轉移的營養(yǎng)元素，與林木吸收的元素量呈一定的比例關系[12]，令

F21ts=a21·F12ts，?t≥1、s；

(14)

F23ts=(1-f23s(t))F12ts-F21ts，?t≥1、s。

(15)

式中：a21為淋溶比例系數(shù)；f23s(t)為流出比例系數(shù)，是隨林木生長時間變化的階段函數(shù)。

枯枝落葉當年的分解量，與上一年積累的某種營養(yǎng)元素量以及當年產(chǎn)生的某種營養(yǎng)元素量之和呈正比例關系，可用公式(16)計算。

(16)

3 結果與分析

3.1 模型參數(shù)的確定

結合研究區(qū)氣象資料，本研究大氣降雨中N、P、K量的取值，分別為40.0、3.5、15.0 kg/(hm2·a)。林區(qū)土壤中每年隨徑流流出的營養(yǎng)元素，約占降雨中營養(yǎng)元素量的70%[16]。不同樹種的枯枝落葉層中各元素之間的分解速率有所差異，但基本在3 a內可全部分解；因此，本研究假設各元素分解速率為100%。由于研究期內(2008—2015年)的苗木仍處于幼齡林階段，且未達到撫育間伐的年齡，因此，對于更新林分主要進行扶正苗木、松土除草等工作。另外，考慮到作業(yè)區(qū)的實際情況，改造項目中沒有采取施肥措施。

低質林更新林分的生長是一個動態(tài)過程，不同種類林分在不同的生長階段，其吸收和存留的營養(yǎng)元素也有一定的差異。通過查閱相關文獻可知，在幼齡林階段，云杉存量為95%以上，落葉松、紅松存量分別為80%、90%左右，而N、P、K三種營養(yǎng)元素之間的比例變化不大。因此，本研究在相關研究的基礎上，確定低質林改造更新?lián)嵊到y(tǒng)中的相關模型參數(shù)[13,17-23](見表1)。造林密度，與樹種、立地條件、造林技術、經(jīng)營水平密切相關。本研究中在10 m×100 m改造帶內分3段栽植紅皮云杉、興安落葉松、紅松幼苗，栽植比例為紅皮云杉∶興安落葉松∶紅松=6∶2∶2，造林時原則上與原有林分邊緣間隔1 m左右，株行距配置為1.5 m×1.5 m[24]。紅松、云杉、落葉松幼苗移栽到更新?lián)嵊到y(tǒng)中的林齡，分別為5、5、3 a；低質林采伐改造當年，土壤中的水解N、有效P、速效K，分別為215.0、9.4、98.0 kg/hm2。

表1 幼齡林木生長過程養(yǎng)分循環(huán)計算經(jīng)驗公式

注：t*為幼齡林階段的初始年限，t*=1 a。

3.2 模型精度檢驗

為了對營養(yǎng)元素循環(huán)模擬模型的精度進行檢驗，本研究將模擬得到的土壤中累積的營養(yǎng)元素量(水解N、有效P、速效K)與實際測量值進行對比分析(見圖2)。由圖2可見，實測值在造林后各測量年份之間存在一定的波動性，而預測值盡管波動性不大，但是能夠準確地捕捉到各指標的長期變化趨勢。各指標的預測值與實測值之間的相關系數(shù)，水解N為0.87、有效P為0.77、速效K為0.88。因此，該營養(yǎng)元素循環(huán)模型顯示了較高的精度，可以用于進一步的營養(yǎng)元素循環(huán)與積累分析。

3.3 營養(yǎng)元素循環(huán)與積累

以低質林改造更新系統(tǒng)內紅皮云杉營養(yǎng)元素循環(huán)為例，2015年，即幼苗栽植8 a后，云杉更新系統(tǒng)通過降雨過程從系統(tǒng)外吸收的N、P、K，分別為40.0、3.5、15.0 kg/hm2；同時，隨徑流流出的N、P、K，分別為28.00、2.45、10.50 kg/hm2。云杉幼林當年從土壤中吸收的N、P、K，分別為5.04、0.72、2.16 kg/hm2；通過淋溶作用歸還到土壤中的N、P、K，分別為0.10、0.02、0.11 kg/hm2；除此之外，通過凋落物回流至土壤中的N、P、K，分別為0.71、0.09、0.24 kg/hm2；土壤中營養(yǎng)元素N、P、K的量，達到286.85、14.35、123.65 kg/hm2。云杉幼林立木中營養(yǎng)元素N、P、K的累積量，達到27.65、3.95、11.85 kg/hm2。

圖2 土壤中營養(yǎng)元素積累實測值與預測值的對比

結合不同種類林分在不同生長階段營養(yǎng)元素存留、吸收的動態(tài)變化，分別計算不同林分養(yǎng)分循環(huán)過程，進一步按林分比例(紅皮云杉∶興安落葉松∶紅松=6∶2∶2)加和計算得到小興安嶺低質林水平帶狀改造帶營養(yǎng)元素循環(huán)計算結果(見表2)。低質林更新?lián)嵊到y(tǒng)造林8 a期間，通過降雨過程從系統(tǒng)外吸收的N、P、K，分別為320、28、120 kg/hm2；同時，隨徑流流出系統(tǒng)的N、P、K，分別為224.0、19.6、84.0 kg/hm2。更新林分從土壤中吸收的N、P、K，分別為26.45 、4.29、11.11 kg/hm2；通過淋溶作用歸還到土壤中的N、P、K，分別為0.36、0.13、0.56 kg/hm2。另外，幼齡林撫育措施及逐年累積的凋落物回流至土壤中的N、P、K，分別為4.17、0.63、1.31 kg/hm2。土壤中累積的營養(yǎng)元素N、P、K，分別達到289.08、14.27、124.75 kg/hm2。

表2 小興安嶺低質林改造更新8 a后養(yǎng)分元素循環(huán)計算結果 kg·hm-2

注：低質林改造帶內不同樹種的栽植比例為紅皮云杉∶興安落葉松∶紅松=6∶2∶2。

4 結論

通過養(yǎng)分元素流動模型得到的土壤養(yǎng)分指標預測值與實測值之間的相關系數(shù)，水解N為0.87、有效P為0.77、速效K為0.88，顯示了較高的模型精度，說明該模型可用于進一步的營養(yǎng)元素循環(huán)與積累分析。各土壤養(yǎng)分指標比改造當年均有較大幅度的增加，表明低質林改造對于林地內的養(yǎng)分循環(huán)具有一定的促進作用。

2008—2015年，更新系統(tǒng)通過降雨過程從系統(tǒng)外吸收的N、P、K，分別為320、28、120 kg/hm2；同時，隨徑流流出系統(tǒng)的N、P、K，分別為224.0、19.6、84.0 kg/hm2。林分從土壤中吸收的N、P、K，分別為26.45、4.29、11.11 kg/hm2；通過淋溶作用歸還到土壤中的N、P、K，分別為0.36、0.13、0.56 kg/hm2。另外，幼齡林撫育措施及逐年累積的凋落物回流至土壤中的N、P、K，分別為4.17、0.63、1.31 kg/hm2；土壤中累積的營養(yǎng)元素N、P、K，分別為289.08、14.27、124.75 kg/hm2。

盡管本研究提出的低質林更新階段營養(yǎng)元素流動模型顯示了很好的模型精度，但是，也存在一些不足，使得預測結果與實際值仍存在一定的偏差。在以后的研究中，將考慮結合光照、氣候等環(huán)境因素以及不同林分的養(yǎng)分含量及生物量預測模型對于預測結果進行進一步的修正。

[1] 聶道平.森林生態(tài)系統(tǒng)營養(yǎng)元素的生物循環(huán)[J].林業(yè)科學研究，1991，4(4)：435-440.

[2] 李俊清，宮偉光.東北主要造林樹種營養(yǎng)元素含量的分析[J].東北林業(yè)大學學報，1991，19(水胡黃椴?？?：228-235.

[3] 曹建華，李小波，趙春梅，等.森林生態(tài)系統(tǒng)養(yǎng)分循環(huán)研究進展[J].熱帶農業(yè)科學，2007，27(6)：68-79.

[4] 劉增文.森林生態(tài)系統(tǒng)的物質積累與循環(huán)[D].北京：中國林業(yè)出版社，2009.

[5] OVINGTON J D. Organic production， turnover and mineral cycling in woodlands[J]. Biological Reviews，1965，40(3)：295-336.

[6] BORMANN F H， LIKENS G E， FISHER D V， et al. Nutrient loss accelerated by clear-cutting of a forest ecosystem[J]. Science，1968，159：882-884.

[7] CLEVE K V， OLIVER L， SCHLENTNER R， et al. Productivity and nutrient cycling in taiga forest ecosystems[J]. Canadian Journal of Forest Research，1983，13(5)：747-766.

[8] ATTIWILL P M， ADAMS M A. Nutrient cycling in forests[J]. New Phytologist，1993，124(4)：561-582.

[9] 李茜，楊勝天，白曉輝，等.貴州典型森林植被群落營養(yǎng)元素生物循環(huán)空間信息模型構建[J].地理研究，2009，28(4)：1047-1058.

[10] 張泱，宋啟亮，董希斌.帶狀皆伐改造對小興安嶺低質林土壤養(yǎng)分的影響[J].東北林業(yè)大學學報，2013，41(6)：33-37.

[11] 呂海龍.小興安嶺低質林改造后生態(tài)系統(tǒng)恢復效果的評價[D].哈爾濱：東北林業(yè)大學，2012.

[12] 張?zhí)穑煊窠?，董希?小興安嶺用材林土壤肥力綜合評價及評價方法比較[J].東北林業(yè)大學學報，2016，44(12)：10-14，98.

[13] 張泱，宋啟亮，董希斌，等.小興安嶺林區(qū)低質林改造效果評價[J].東北林業(yè)大學學報，2009，37(12)：35-36，55.

[14] 陳俊松.基于生命周期評價理論的人工林作業(yè)資源—環(huán)境—經(jīng)濟影響研究[D].南京：南京林業(yè)大學，2011.

[15] 陳百靈，董希斌，崔莉，等.大興安嶺低質林生態(tài)改造后枯落物水文效應變化[J].東北林業(yè)大學學報，2015，43(6)：72-77，88.

[16] 張鼎華.皆伐煉山對林地養(yǎng)分元素資源遷移影響的研究[J].自然資源學報，1997，12(4)：343-348.

[17] 詹鴻振，任淑文，沈淑娟.紅松幼樹的營養(yǎng)特性研究[J].東北林業(yè)大學學報，1989，17(1)：9-14.

[18] 丁寶永，孫繼華.紅松人工林生態(tài)系統(tǒng)生物生產(chǎn)力及養(yǎng)分循環(huán)研究[J].東北林業(yè)大學學報，1989，17(專輯)：1-98.

[19] 曹國軍，周百娟，郭士偉，等.紅皮云杉苗期營養(yǎng)規(guī)律初探[J].吉林農業(yè)大學學報，1993，15(3)：43-48，103.

[20] 韋淑英，馬華文，黃青新，等.紅皮云杉苗木生物量群體多樣性的研究[J].北華大學學報(自然科學版)，2001，2(3)：246-250，256.

[21] 宿以明，劉興良，何飛，等.日本落葉松人工林營養(yǎng)元素含量、貯量與分配的研究[J].四川林業(yè)科技，2003，24(3)：35-39.

[22] 丁寶永，徐立英，張世英.興安落葉松人工林營養(yǎng)元素的分析[J].生態(tài)學報，1989，9(1)：71-76.

[23] 朱玉杰，董希斌，李祥.不同撫育強度對興安落葉松幼苗光合作用的影響[J].東北林業(yè)大學學報，2015，43(10)：51-54，67.

[24] 呂海龍，董希斌.基于主成分分析的小興安嶺低質林不同皆伐改造模式評價[J].林業(yè)科學，2011，47(12)：172-178.

SimulationonNutrientCyclingofALow-qualityForestStandRegenerationSystemafterLoggingTransformation//

WuJinzhuo,KongLinlin,DongXibin,ShenWei

(NortheastForestryUniversity,Harbin150040,P.R.China)

//JournalofNortheastForestryUniversity，2017,45(8)：25-29.

Accordingtothenutrientcyclingprocessoflow-qualityforeststandregenerationsystem,BulkInternalFlow(BIF)andMaterialFlowAnalysis(MFA)wereusedtoestablishamathematicalmodeltosimulatethenutrientcyclingofalow-qualityforeststandregenerationsystemafterloggingtransformation.Afterverifyingtheeffectivenessofthemodelbyfieldinvestigation,themodelwasusedtosimulateandanalyzethenutrientcyclingandaccumulationofalow-qualityforeststandintheregenerationstageaftertransformation(2008-2015)inthestudyarea.Thecorrelationcoefficientsbetweenthepredictedsoilnutrientindexesobtainedbythesimulationmodelandtheactualmeasuredvalueswereasfollows:hydrolysablenitrogencontentof0.87,availablephosphoruscontentof0.77,andavailablepotassiumcontentof0.88,indicatinghighermodelprecision.In2008-2015,theregenerationsystemabsorbed320kg/hm2N, 28kg/hm2P,and120kg/hm2Kfromrainfall,and224kg/hm2N, 19.6kg/hm2P, 84kg/hm2Kwereflowedoutofthesystem.Theregenerationstandsabsorbed26.45kg/hm2N, 4.29kg/hm2P, 11.11kg/hm2Kfromforestsoilsand0.36kg/hm2N, 0.13kg/hm2P, 0.56kg/hm2Kwerereturnedbacktosoilsbyeluviation.Inaddition, 4.17kg/hm2N, 0.63kg/hm2P, 1.31kg/hm2Kwerefinallyreturnedtosoilbydecompositionwithforestlitter.Theaccumulatednutrientinforestsoilswere289.08kg/hm2N, 14.27kg/hm2P,and124.75kg/hm2K,withasignificantincreasecomparedwiththatintheyearoflow-qualityforesttransformation,whichindicatedthatthelow-qualityforesttransformationhadacertainroleinpromotingthenutrientcyclingintheregeneratedforeststand.

Low-quality forest; Low-quality forest transformation; Low-quality forest regeneration; Forest nutrient

1)國家自然科學基金項目(31400539)。

吳金卓，女，1980年12月生，東北林業(yè)大學工程技術學院，副教授。E-mail：wjz@nefu.edu.cn。

董希斌，東北林業(yè)大學工程技術學院，教授。E-mail：xibindong@nefu.edu.cn。

2017年3月23日。

S756.5；S718.55

責任編輯：張 玉。