辛贊紅， 江 洪，2， 接程月， 魏曉華， BLANCO Juan， 周?chē)?guó)模

（1. 浙江農(nóng)林大學(xué) 國(guó)際空間生態(tài)與生態(tài)系統(tǒng)生態(tài)研究中心， 浙江 臨安311300； 2. 南京大學(xué) 國(guó)際地球系統(tǒng)科學(xué)研究所， 江蘇 南京210093； 3. 不列顛哥倫比亞大學(xué) 地球與環(huán)境科學(xué)系， 不列顛哥倫比亞 基隆那V1V1V7）

森林生態(tài)系統(tǒng)的養(yǎng)分循環(huán)是生態(tài)系統(tǒng)的重要功能過(guò)程之一， 直接影響著森林的生產(chǎn)力［1-2］。 作為生命必需元素之一的氮是影響森林生態(tài)系統(tǒng)生產(chǎn)力的重要限制因素［2］。 研究一定立地條件下森林生態(tài)系統(tǒng)中氮元素的分配、 積累和循環(huán)， 可以進(jìn)一步認(rèn)識(shí)森林生態(tài)系統(tǒng)的功能［3］。 這對(duì)于改進(jìn)森林經(jīng)營(yíng)策略， 促進(jìn)林木速生、 豐產(chǎn)、 優(yōu)質(zhì)也有重要意義。 杉木Cunninghamia lanceolata 是中國(guó)南方最重要的用材樹(shù)種， 生長(zhǎng)快， 材質(zhì)好， 單產(chǎn)高， 已有1 000 多年的栽培歷史， 在中國(guó)南方林業(yè)生產(chǎn)中占舉足輕重的地位［4］。 近年來(lái)由于杉木連栽造成生產(chǎn)力日趨下降， 已嚴(yán)重影響了杉木人工林的持續(xù)經(jīng)營(yíng)［5-8］。 長(zhǎng)期定位試驗(yàn)是研究森林連栽或輪伐對(duì)地力影響的重要方法， 但是， 其需時(shí)長(zhǎng)， 不確定性大， 操作起來(lái)有很大的難度。 現(xiàn)階段， 利用模型進(jìn)行相關(guān)的模擬， 已是一個(gè)重要的發(fā)展趨勢(shì)［9-10］。 森林生態(tài)系統(tǒng)生物地球化學(xué)過(guò)程模擬模型是了解森林生態(tài)系統(tǒng)生物組分與其無(wú)機(jī)環(huán)境之間相互作用的重要工具［9］。 目前已經(jīng)開(kāi)發(fā)出了許多種生態(tài)系統(tǒng)過(guò)程模擬模型用以森林生態(tài)系統(tǒng)的研究［10-13］， FORECAST 模型也是其中重要的一種［14］。 本研究通過(guò)運(yùn)用FORECAST 模型模擬杉木人工林在不同輪伐期的生長(zhǎng)過(guò)程， 揭示其氮循環(huán)的特點(diǎn)， 從而為深入了解杉木人工林在連栽條件下的養(yǎng)分循環(huán)機(jī)理和杉木人工林的合理經(jīng)營(yíng)管理提供基礎(chǔ)數(shù)據(jù)。

1 研究方法

1.1 FORECAST 及其氮循環(huán)的簡(jiǎn)單示意圖

FORECAST 模型是一個(gè)以林分群體特征及林地養(yǎng)分循環(huán)為基礎(chǔ)的森林生態(tài)系統(tǒng)管理模型［14］。 FORECAST 模型關(guān)注的是土壤有機(jī)物和營(yíng)養(yǎng)庫(kù)不斷變化的層面上樹(shù)木的生長(zhǎng)狀況和林材的收獲產(chǎn)量。 運(yùn)用該模型進(jìn)行模擬可以得到土壤有機(jī)物和營(yíng)養(yǎng)庫(kù)中各種營(yíng)養(yǎng)元素循環(huán)變化的情況。 FORECAST 模型使用的模擬方法是將生物產(chǎn)量表方法和過(guò)程模擬方法相結(jié)合的混合模擬方法。 模型的驅(qū)動(dòng)函數(shù)是葉氮同化率（FNE）。 FORECAST 模型氮循環(huán)的簡(jiǎn)單示意圖1 所示。

1.2 杉木FORECAST 模型的校準(zhǔn)和應(yīng)用

FORECAST 模型的應(yīng)用有4 個(gè)步驟： ①數(shù)據(jù)的收集、 輸入和驗(yàn)證； ②建立生態(tài)系統(tǒng)模擬運(yùn)行的初始條件； ③定義森林管理方式或者干擾方式； ④模擬各種管理方式并分析模型輸出結(jié)果。 前2 個(gè)步驟即為模型的參數(shù)校準(zhǔn)。 模型的校準(zhǔn)數(shù)據(jù)是由描述不同養(yǎng)分條件林地上林木（地上部分和地下部分）和小型植被按照林分生長(zhǎng)的時(shí)間序列而逐漸累積的生物量的數(shù)據(jù)組成。 林木生物量和林分自疏率通常由樹(shù)高、 胸徑和林分密度結(jié)合具有樹(shù)種特異性的組分生物量異速生長(zhǎng)方程得到。 在模型的養(yǎng)分模擬部分， 需要不同養(yǎng)分在各種生物量組分中的含量的數(shù)據(jù)。 模型同時(shí)還需要不同葉量所產(chǎn)生的遮光度的數(shù)據(jù)和樹(shù)葉對(duì)不同光照強(qiáng)度響應(yīng)的數(shù)據(jù)。 模型校準(zhǔn)所需的這些數(shù)據(jù)主要來(lái)自于文獻(xiàn)， 包括福建、 湖南等主要杉木中心產(chǎn)區(qū)的研究報(bào)告和江西、 廣西等氣候相似區(qū)域的杉木林的研究文獻(xiàn)。

第2 階段的校準(zhǔn)需要在“設(shè)置” 模式下運(yùn)行模型來(lái)建立立地初始條件。 在這個(gè)階段， 模型運(yùn)行時(shí)需關(guān)閉養(yǎng)分反饋， 從而使表示立地條件的植被、 凋落物和土壤有機(jī)質(zhì)得以積累， 這反映出歷史的積累模式。

杉木人工林的模型參數(shù)在借鑒Bi 等［15］的基礎(chǔ)上進(jìn)行了一些調(diào)整， 數(shù)據(jù)主要來(lái)源于田大倫等編著的《杉木林生態(tài)系統(tǒng)》以及盛煒彤等編著的《杉木人工林長(zhǎng)期生產(chǎn)力保持機(jī)制研究》。 模型中的降水營(yíng)養(yǎng)輸入、 坡面滲流、 礦質(zhì)土壤陽(yáng)離子和陰離子交換能力、 腐殖質(zhì)質(zhì)量、 凋落物分解速率等數(shù)據(jù)來(lái)源于廖利平等［16-17］、黃志群等［18］、 丁應(yīng)祥等［19］、 楊玉盛等［20］的研究結(jié)果。 樹(shù)木生物量、 死亡率和林分密度、 樹(shù)高以及林冠高度、 活組織的營(yíng)養(yǎng)濃度的數(shù)據(jù)以楊玉盛等［21］、 肖文發(fā)等［22］、 林開(kāi)敏等［23］、 劉茜等［24］、 周后盛等［25］發(fā)表的數(shù)據(jù)為基礎(chǔ)。

1.4 模擬方案

圖1 FORECAST 模型氮循環(huán)的簡(jiǎn)單示意圖Figure 1 Diagram of nitrogen cycling in FORECAST model

模型的校準(zhǔn)數(shù)據(jù)來(lái)自中國(guó)可觀測(cè)的杉木人工林范圍內(nèi)的所有立地， 立地按照質(zhì)量情況分為較差、 中等、 較好， 并被量化為17， 21 和27（從差到好， 基于優(yōu)勢(shì)木50 a 的生長(zhǎng)高度）來(lái)代表模型推演林木生長(zhǎng)的相關(guān)指數(shù)。 FORECAST 模型中數(shù)量化的立地質(zhì)量可以是但并不是必須基于立地指數(shù)即某一指數(shù)年齡優(yōu)勢(shì)木的平均高度［26-27］。 模擬的年限為200 a； 3 種不同的輪伐期為10， 25 和50 a； 模擬的收獲方式為莖干采伐法（每次輪伐期末除去90%的樹(shù)干生物量， 殘留部分均遺留在原地）。 除采伐收獲外沒(méi)有其他干擾。

2 模擬結(jié)果

2.1 FORECAST 模型模擬結(jié)果的檢驗(yàn)

模型驗(yàn)證有很多方法。 本研究采用模型模擬結(jié)果與實(shí)驗(yàn)測(cè)量數(shù)據(jù)對(duì)比的方法來(lái)進(jìn)行驗(yàn)證。 由于比較不同輪伐期（不同代數(shù)）養(yǎng)分氮循環(huán)的研究很少， 所以我們獲得的相關(guān)文獻(xiàn)數(shù)據(jù)非常有限。 根據(jù)劉愛(ài)琴等［1］的結(jié)果， 可以得到立地指數(shù)為16 的2 代11 年生和31 年生杉木人工林氮吸收量的數(shù)據(jù)。 由模擬結(jié)果可知， 3 種不同立地20 a 的優(yōu)勢(shì)木高分別為： 較差立地12.1 m， 中等立地15.3 m， 較好立地19.1 m。 對(duì)照文獻(xiàn)地位指數(shù)， 可以選用中等立地11 a 輪伐期末和31 a 輪伐期末的氮吸收量數(shù)據(jù)與文獻(xiàn)數(shù)據(jù)進(jìn)行比較。 模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)的對(duì)比如表1。 將模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)作回歸分析， 可得方程y = 1.895 6x-76.728， r2= 0.969 3**。 由此可見(jiàn)： 模型模擬第1 代和第2 代杉木人工林氮的吸收量有較高的可信度。

2.2 土壤有效氮的模擬結(jié)果

由圖2 可見(jiàn)： 3 種不同立地條件下，輪伐期長(zhǎng)短對(duì)土壤有效氮含量有顯著的影響。 在較差立地上， 10 a 輪伐期的前5個(gè)輪伐期土壤有效氮含量較大， 下降也較快； 從第6 個(gè)輪伐期開(kāi)始有小幅下降， 然后逐漸過(guò)渡到平穩(wěn)。 25 a 輪伐期的第1 個(gè)輪伐期， 其土壤有效氮含量小于10 a 輪伐期第1 個(gè)輪伐期的值， 但是25 a 輪伐期的下降幅度較小， 并且在第2 個(gè)輪伐期之后就逐漸平穩(wěn)。 50 a輪伐期的土壤有效氮含量高于25 a 輪伐期的土壤有效氮含量， 并且高于第2 個(gè)輪伐期之后的10 a 輪伐期的值； 50 a 輪伐期的土壤有效氮含量在各個(gè)輪伐期之間的變化不大。 中等立地和較好立地有相似的規(guī)律， 只是相比差等立地25 a 輪伐期的值要低于10 a 輪伐期的值。

表1 氮吸收量預(yù)測(cè)數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)對(duì)比Table 1 Comparison between predicted data and observed data

2.3 氮的吸收量和歸還量

圖3 不同輪伐期氮的年均吸收量Figure 3 Uptake N in different rotations of Chinese fir plantation

圖4 不同輪伐期凋落物形式氮的年歸還量Figure 4 Average quantity of N returned by litterfall at different rotations

由圖3 和圖4 可見(jiàn)： 3 種不同立地條件下， 氮的年吸收量和凋落物形式氮的年歸還量， 總體趨勢(shì)都是隨著輪伐期的增長(zhǎng)而增加， 隨著輪伐代數(shù)的增加而減少。 10 a 輪伐期的吸收量和歸還量在前4 個(gè)輪伐期有明顯的下降， 而且下降幅度相較25 a 輪伐期要大； 50 a 輪伐期的吸收量和歸還量， 在各輪伐期之間變化不大， 趨勢(shì)相對(duì)比較平穩(wěn)。 由于模擬采用的是莖干采伐的收獲方式， 枝葉和根均遺留在原地， 所以大量的氮以凋落物的形式歸還給土壤。 以中等立地為例， 不同輪伐期氮的循環(huán)系數(shù)［28］（即氮的年歸還量比氮的年吸收量）分別為： 0.544， 0.766， 0.847 和0.879， 可見(jiàn)隨著輪伐期增長(zhǎng)氮的循環(huán)系數(shù)增大。

2.4 土壤中氮的年淋溶損失

由圖5 可知： 杉木人工林的土壤氮淋溶損失量總體是隨著立地條件的變好而增大。 3 種不同立地條件下， 10 a 輪伐期的年平均土壤氮淋溶損失量較大， 并且在前3 個(gè)輪伐期表現(xiàn)尤為明顯， 3 個(gè)輪伐期之后逐漸平穩(wěn)； 25 a 輪伐期的年均土壤氮淋溶損失量較小， 只在第1 個(gè)輪伐期和第2 個(gè)輪伐期之間變化比較大； 50 a 輪伐期的年均土壤氮淋溶損失量較小， 輪伐期之間的變化也不大。 總體趨勢(shì)上， 較差立地和中等立地， 土壤氮淋溶損失量都是隨著輪伐期的增長(zhǎng)而減??； 較好立地50 a 輪伐期的土壤氮淋溶損失量要大于25 a 輪伐期的土壤氮淋溶損失量， 但是小于10 a 輪伐期的淋溶量。

圖5 不同輪伐期土壤氮的淋溶損失量Figure 5 Chinese fir plantation soil N leaching loss at different rotations

2.5 葉片含氮量的動(dòng)態(tài)

葉片含氮量是FORECAST 模型一個(gè)重要的輸出參數(shù)。 通過(guò)這個(gè)參數(shù)， 我們可以看出伴隨著林分生長(zhǎng)， 氮元素在葉片中積累的量。 圖6 顯示了3 種立地條件下， 不同輪伐期葉片含氮量的變化趨勢(shì)。 由圖6 可知： 隨著輪伐期的增長(zhǎng)積累在葉片中的氮逐漸增加。 10 a 輪伐期的葉片含氮量隨著立地條件變好與25 a 輪伐期的差距逐漸減小， 這可能是因?yàn)殡S著立地條件的變好， 林木迅速生長(zhǎng)， 使得10 a 輪伐期的密度優(yōu)勢(shì)逐漸顯現(xiàn)了出來(lái)， 從而造成了葉片含氮量在總量上的增加。

圖6 杉木人工林不同輪伐期葉片氮含量的動(dòng)態(tài)Figure 6 Chinese fir plantation foliage N dynamics at different rotations

3 分析與討論

3.1 FORECAST 模型模擬的精度

FORECAST 模型通過(guò)大量森林生態(tài)系統(tǒng)的模擬， 已得到了較好的檢驗(yàn)［29-34］。 對(duì)于杉木森林生態(tài)系統(tǒng)，Bi 等［15］作了實(shí)測(cè)數(shù)據(jù)與模擬結(jié)果的測(cè)試， 證明模擬結(jié)果精度較高。 本模擬的結(jié)果也證明， FORECAST 模型模擬杉木林生態(tài)系統(tǒng)氮循環(huán)的關(guān)鍵變量， 具有較高的模擬精度， 可以廣泛應(yīng)用在有關(guān)的模擬實(shí)踐中。

3.2 不同輪伐期杉木人工林生態(tài)系統(tǒng)氮元素的變化趨勢(shì)

輪伐期的長(zhǎng)短會(huì)影響杉木人工林氮元素的循環(huán)和利用。 彭長(zhǎng)輝等［10］、 江洪等［11］運(yùn)用Century 4.0 模型對(duì)加拿大北方森林和中國(guó)北方森林所做的研究顯示， 在采伐強(qiáng)度一定的情況下， 土壤有效氮含量在長(zhǎng)輪伐期下最大， 并且隨著輪伐期的縮短而減小。 極短的輪伐期會(huì)造成北方森林減產(chǎn)達(dá)65%。 通過(guò)FORECAST 模型的模擬結(jié)果我們可以看出： 在采伐方式一定的情況下， 短輪伐期（10 a）土壤有效氮含量下降，氮吸收量降低， 長(zhǎng)輪伐期（50 a）土壤有效氮和氮的吸收都保持較大的量并且輪伐期之間的變化較平穩(wěn)。 隨著輪伐期的增長(zhǎng)， 凋落物形式氮的年歸還量也逐漸增大。 以中等立地為例， 10， 25， 50 a 輪伐期土壤有效氮200 a 的平均值分別為96.24， 86.87 和96.07 kg·hm-2； 25 a 輪伐期相對(duì)10 a 輪伐期減少了9.7%；50 a 輪伐期與10 a 輪伐期基本持平， 相對(duì)25 a 輪伐期增加10.6%。 吸收氮的平均值分別為49.22， 58.44和69.07 kg·hm-2； 25 a 輪伐期相對(duì)10 a 年輪伐期增加了18.7%； 50 a 輪伐期相對(duì)10 a 輪伐期增加了40.3%， 相對(duì)25 a 輪伐期增加18.2%。 凋落物形式歸還氮的平均值分別為26.75， 44.76 和60.74 kg·hm-2；25 a 輪伐期相對(duì)10 a 年輪伐期增加了67.3%； 50 a 輪伐期相對(duì)10 a 輪伐期增加了127.1%， 相對(duì)25 a 輪伐期增加35.7%。 從模擬結(jié)果看： 10 a 短輪伐期下土壤氮的淋溶損失總量， 較差立地為1 167.57 kg·hm-2，中等立地為1 075.27 kg·hm-2， 較 好立 地為1 085.70 kg·hm-2； 25 a 輪伐期分別為895.05， 893.60 和924.28 kg·hm-2； 50 a 輪伐期分別為707.33， 765.98 和959.05 kg·hm-2。 10 a 短輪伐期土壤氮的淋溶損失總量均高于同等立地條件下25 a 輪伐期和50 a 輪伐期的總量。 從以上模擬結(jié)果可以看出： 輪伐期的長(zhǎng)短對(duì)杉木人工林生態(tài)系統(tǒng)氮元素的循環(huán)和利用具有重要影響， 過(guò)短的輪伐期會(huì)導(dǎo)致土壤有效氮含量下降， 氮吸收量降低， 淋溶損失增加， 從而破壞土壤肥力， 造成林木生產(chǎn)力的下降。

3.3 連栽代數(shù)對(duì)杉木人工林氮元素循環(huán)和利用的影響

從模擬結(jié)果可以看出： 相同輪伐期不同代數(shù)杉木林的氮循環(huán)和利用存在著顯著差異。 隨著代數(shù)的增加， 土壤有效氮含量、 年吸收量和年歸還量均呈遞減趨勢(shì)， 而且這種趨勢(shì)隨著輪伐周期的變短而愈加明顯。 以中等立地前4 代（即前4 個(gè)輪伐期）為例， 10 a 輪伐期的第2 代、 第3 代和第4 代土壤有效氮較第1 代分別減少了17.2%， 23.6%和27.2%； 年吸收量分別減少了10.1%， 16.9%和20.4%； 年歸還量分別減少了8.1%， 14.1%和17.9%。 25 a 輪伐期的第2 代、 第3 代和第4 代土壤有效氮相較第1 代分別減少了10.6%， 13.1%和14.5%； 年吸收量分別減少了9.2%， 12.0%和13.9%； 年歸還量分別減少了8.9%，11.8%和13.4%。 可見(jiàn)， 隨著杉木連栽代數(shù)的增加， 土壤有效氮減少， 土壤供肥和保肥的能力下降， 吸收氮和凋落物歸還氮的量也在減少， 從而可能導(dǎo)致林分產(chǎn)量的降低。 本研究的模擬結(jié)果中， 第1 代、 第2 代和第3 代土壤有效氮的變化幅度較大， 這與盛煒彤［6］、 楊玉盛等［22］研究者的結(jié)論相似。

3.4 展望

在這次模擬工作中， 只是研究不同立地條件下輪伐期長(zhǎng)短對(duì)杉木人工林生態(tài)系統(tǒng)氮?jiǎng)討B(tài)的影響， 而沒(méi)有考慮煉山、 整地和疏伐等育林干擾。 我們將進(jìn)一步完善模型， 加入煉山、 整地、 疏伐等育林措施，使模型模擬更貼近現(xiàn)實(shí)的生產(chǎn)管理， 為杉木人工林的合理經(jīng)營(yíng)提供更為準(zhǔn)確的數(shù)據(jù)。

［1］ 劉愛(ài)琴， 范少輝， 林開(kāi)敏， 等. 不同栽植代數(shù)杉木林養(yǎng)分循環(huán)的比較研究［J］. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào)， 2005， 11 （2）：273 - 278.LIU Aiqin， FAN Shaohui， LIN Kaimin， et al. Comparison on nutrient cycling in different generation plantations of Chinese-fir ［J］. Plant Nutr Fert Sci， 2005， 11 （2）： 273 - 278.

［2］ 諶小勇， 潘維儔. 杉木人工林生態(tài)系統(tǒng)中氮素的動(dòng)態(tài)特征［J］. 生態(tài)學(xué)報(bào)， 1989， 9 （3）： 201 - 206.SHEN Xiaoyong， PAN Weichou.Nitrogen dynamics of Chinese-fir plantation ecosystem ［J］.Acta Ecol Sin， 1989， 9 （3）：201 - 206.

［3］ 沈國(guó)舫， 董世仁， 聶道平. 油松人工林養(yǎng)分循環(huán)的研究營(yíng)養(yǎng)元素的含量及分布［J］. 北京林學(xué)院學(xué)報(bào)， 1985， 7 （4）：1- 14.SHEN Guofang， DONG Shiren， NIE Daoping. Study on nutrient cycling in Pinus tabulaeformis plantation ecosystem［J］. J Beijing For Coll， 1985， 7 （4）： 1 - 14.

［4］ 范少輝， 盛煒彤， 馬祥慶， 等. 多代連栽對(duì)不同發(fā)育階段杉木人工林生產(chǎn)力的影響［J］. 林業(yè)科學(xué)研究， 2003， 16（5）： 560 - 567.FAN Shaohui， SHENG Weitong， MA Xiangqing， et al. Effects of successive-generation planting on Chinese-fir plantation productivity at different development stages ［J］. For Res， 2003， 16 （5）： 560 - 567.

［5］ 徐化成. 關(guān)于人工林的地力下降問(wèn)題［J］. 世界林業(yè)研究， 1992， 5 （1）： 66 - 73.XU Huacheng. On the plantation soil fertility declination ［J］. World For Res， 1992， 5 （1）： 66 - 73.

［6］ 盛煒彤. 人工林地力衰退研究［M］. 北京： 中國(guó)科學(xué)技術(shù)出版社， 1992.

［7］ 俞新妥. 杉木林地持續(xù)利用問(wèn)題的研究和看法［J］. 世界林業(yè)研究， 1993， 6 （2）： 80 - 86.YU Xintuo. Study and suggestions on the sustainable utilization of Chinese fir land ［J］. World For Res， 1993， 6 （2）：80 - 86.

［8］ 馬祥慶. 杉木人工林連栽生產(chǎn)力下降研究進(jìn)展［J］. 福建林學(xué)院學(xué)報(bào)， 2001， 21 （4）： 380 - 384.MA Xiangqing.Advance in researches on productivity decline of repeating Chinese fir forests ［J］.J Fujian Coll For， 2001，21 （4）： 380 - 384.

［9］ KROS H， WARFVINGE P. Evaluation of model behaviour with respect to the biogeochemistry at the Solling spruce site ［J］. Ecol Model， 1995， 83： 255 - 262.

［10］ PENG Changhui， JIANG Hong， APPS M J， et al. Effects of harvesting regimes on carbon and nitrogen dynamics of boreal forests in central Canada： a process model simulation ［J］. Ecol Model， 2002， 155： 177 - 189.

［11］ JIANG Hong， APPS M J， PENG Changhui， et al. Modeling the influence of harvesting on Chinese boreal forest carbon dynamics ［J］. For Ecol Manage， 2002， 169： 65 - 68.

［12］ FERRIER R C， WRIGHT R F， COSBY B J， et al. Application of the MAGIC model to the Norway spruce stand at Solling， Germany ［J］. Ecol Model， 1995， 83： 77 - 84.

［13］ VERBURG P S J， JOHNSON D W. A spreadsheet-based biogeochemical model to simulate nutrient cycling processes in forest ecosystems ［J］. Ecol Model， 2001， 141， 185 - 200.

［14］ 接程月， 辛贊紅， 信曉穎， 等. FORECAST 模型的原理、 方法和應(yīng)用［J］. 浙江林學(xué)院學(xué)報(bào)， 2009， 26 （6）： 909－ 915.JIE Chengyue， XIN Zanhong， XIN Xiaoying， et al.Principle， method and application of FORECAST model ［J］.J Zhejiang For Coll， 2009， 26 （6）： 909－915.

［15］ BI J， BlANCO J A， SEELY B， et al. Yield decline in Chinese-fir plantations： a simulation investigation with implications for model complexity ［J］. Can J For Res， 2007， 37： 1615 - 1630.

［16］ 廖利平， 楊躍軍， 汪思龍， 等. 杉木Cunninghamia lanceolata， 火力楠Michelia macclurei 純林及其混交林細(xì)根分布、 分解與養(yǎng)分歸還［J］. 生態(tài)學(xué)報(bào)， 1999， 19 （3）： 246 - 250.LIAO Liping， YANG Yuejun， WANG Silong， et al. Distribution， decomposition and nutrient return of the fine root in pure Cunninghamia lanceolata， Michelia macclurei and the mixed plantations ［J］. Acta Ecol Sin， 1999， 19 （3）：246 - 250.

［17］ 廖利平， 馬越強(qiáng)， 汪思龍， 等. 杉木與主要闊葉造林樹(shù)種葉凋落物的混合分解［J］. 植物生態(tài)學(xué)報(bào)， 2000， 24 （1）：27 - 33.LIAO Liping， MA Yueqiang， WANG Silong， et al. Decomposition of leaf litter of Chinese fir in mixture with major associated broad-leaved plantation species ［J］. Acta Phytoecol Sin， 2000， 24 （1）： 27 - 33.

［18］ 黃志群， 廖利平， 汪思龍， 等. 杉木根樁和周?chē)寥婪雍康淖兓捌浠行?yīng)［J］. 應(yīng)用生態(tài)學(xué)報(bào)， 2000， 11（2）： 190 - 192.HUANG Zhiqun， LIAO Liping， WANG Silong， et al. Dynamics of phenolics content of Chinese fir stump-roots and the rhizosphere soil and its allelopathy ［J］. Chin J Appl Ecol， 2000， 11 （2）： 190 - 192.

［19］ DING Yingxiang， TIAN Ye. A testing simulation with FOORTOON on long-term productivity of Chinese fir plantations ［J］. For Study China， 1999， 1 （2）： 34 - 38.

［20］ 楊玉盛， 陳光水， 黃寶龍. 杉木多世代連栽的土壤水分和養(yǎng)分變化［J］. 南京林業(yè)大學(xué)學(xué)報(bào)， 2000， 24 （2）： 25 -28.YANG Yusheng， CHEN Guangshui， HUANG Baolong. Variation in the soil water and nutrients between different rotation stands of Chinese fir ［J］. J Nanjing For Univ， 2000， 24 （2）： 25 - 28.

［21］ 楊玉盛， 邱仁輝， 何宗明， 等. 不同栽杉代數(shù)29 年生杉木林凈生產(chǎn)力及營(yíng)養(yǎng)元素生物循環(huán)的研究 ［J］. 林業(yè)科學(xué)， 1998， 34 （6）： 3 - 11.YANG Yusheng， QIU Renhui， HE Zongming， et al. Study on the stand net productivity and biological cycle of nutrients in the 29-year-old plantations of Chinese fir on different rotations ［J］. Sci Silv Sin， 1998， 34 （6）： 3 - 11.

［22］ 肖文發(fā)， 聶道平， 張家誠(chéng). 我國(guó)杉木林生物量與能量利用率的研究［J］. 林業(yè)科學(xué)研究， 1999， 12 （3）： 253 -258.XIAO Wenfa， NIE Daoping， ZHANG Jiacheng. Study on biomass and energy use efficiency of the stands of Cunninghamia lanceolata in China ［J］. For Res， 1999， 12 （3）： 253 - 258.

［23］ 林開(kāi)敏， 俞新妥， 何智英， 等. 不同密度杉木林分生物量結(jié)構(gòu)與土壤肥力差異研究［J］. 林業(yè)科學(xué)， 1996， 32（5）：385 - 392.LIN Kaimin， YU Xintuo， HE Zhiying， et al. Study on the difference of biomass structure and soil fertility in Chinese fir stands of different densities ［J］. Sci Silv Sin， 1996， 32 （5）： 385 - 392.

［24］ 劉茜. 杉木種子營(yíng)養(yǎng)化學(xué)成分的研究［J］. 林業(yè)科學(xué)， 1998， 34 （2）： 92 - 95.LIU Qian. A study on the chemical composition of Chinese fir seeds ［J］. Sci Silv Sin， 1998， 34 （2）： 92 - 95.

［25］ 周后盛. 杉木多代連栽地營(yíng)造混交林生產(chǎn)力的初步研究［J］. 福建林學(xué)院學(xué)報(bào)， 1999， 19 （2）： 165 - 169.ZHOU Housheng.Study on productivity of mixed forest planted on site of continuously planted Chinese fir ［J］.J Fujian Coll For， 1999， 19 （2）： 165 - 169.

［26］ KIMMINS J P, MAILLY D, SEELY B. Modeling forest ecosystem net primary production： the hybrid simulation approach used in FORECAST ［J］. Ecol Model， 1999， 122： 195 - 224.

［27］ KIMMINS J P. Scientific Foundations for the Simulation of Ecosystem Function and Management in FORCYTE-l1［R］. Vancouver： University of British Columbia. Faculty of Forestry. Department of Forest Sciences， 1993： 3 - 7.

［28］ 方晰， 田大倫， 蔡寶玉， 等. 第2 代杉木人工林氮、 磷、 鉀的吸收、 積累和遷移［J］. 中南林學(xué)院學(xué)報(bào)， 2002， 22（2）： 1 - 6.FANG Xi， TIAN Dalun， CAI Baoyu， et al. Absorption， accumulation and transportation of N， P， K elements in the second-generation Chinese-fir plantation ［J］. J Central South For Univ， 2002， 22 （2）： 1 - 6.

［29］ MORRIS D M， KIMMINS J P， DUCKERT D R. The use of soil organic matter as a criterion of the sustainability of forest management alternatives： a modeling approach using FORECAST ［J］. For Ecol Manage， 1997， 94： 61 - 78.

［30］ WEI Xiaohua， LIU W， WATERHOUSE J， et al. Simulations on impacts of different management strategies on longterm site productivity in lodgepole pine forests of the central interior of British Columbia ［J］.For Ecol Manage， 2000，133： 217 - 229.

［31］ WEI Xiaohua， KIMMINS J P， ZHOU Guoyi. Disturbance and the sustainability of long-term site productivity in lodgepole pine forests in the central interior of British Columbia： an ecosystem modeling approach ［J］. Ecol Model，2003， 164： 239 - 256.

［32］ SEELY B， WELHAM C， KIMMINS J P. Carbon sequestration in boreal forest ecosystem： results from the ecosystem simulation model， FORECAST ［J］. For Ecol Manage， 2002， 169： 123 - 135.

［33］ WELHAM C， SEELY B， KIMMINS J P. The utility of the two-pass harvesting system： an analysis using the ecosystem simulation model FORECAST ［J］. Can J For Res， 2002， 32： 1071 - 1079.

［34］ REICH P B， GRIGAL D F， ABER J D， et al. Nitrogen mineralization and productivity in 50 hardwood and conifer stands on diverse soils ［J］. Ecology， 1997， 78： 335 - 347.