基于橡膠林生態系統養分循環指導施肥的效益分析

張以山+袁曉軍+曹建華

摘 要 以橡膠品種無性系PR107、RRIM600、熱研7-33-97為研究對象，借助已研究建立的橡膠樹精準施肥決策計算機系統和橡膠林生態系統養分循環試驗數據，比較分析了基于橡膠林生態系統養分循環施肥決策系統指導施肥的經濟效益和生態效益。結果表明，利用橡膠林生態系統養分循環施肥決策系統進行施肥，能促進天然橡膠增產，土壤有機質含量有所增加，膠園土壤酸化趨勢得到減緩，經濟效益和生態效益明顯。

關鍵詞 橡膠林 ；養分循環 ；施肥指導 ；經濟效益 ；生態效益

中圖分類號 S794.1 文獻標識碼 A Doi：10.12008/j.issn.1009-2196.2017.09.006

Benefit Analysis of Rubber Plantations Based on Nutrient Cycle

of Rubber Plantation Ecosystem

ZHANG Yishan1） YUAN Xiaojun2） CAO Jianhua2）

（1 Chinese Academy of Tropical Agricultural Sciences， Haikou 571737；

2 Rubber Research Institute， CATAS， Danzhou 571737）

Abstract On the basis of the available precision fertilizing decision system and experimental data， rubber clones PR107， RRIM600 and Reyan 7-33-97 were selected to analyze economic and ecological benefits from their rubber plantations dressed with fertilizer based on the fertilizer recommendation from the nutrient cycle of rubber plantation ecosystem. The results showed that with this decision system for fertilizer recommendation the rubber plantations increased in yield and soil organic content， and mitigated the trend of soil acidification， resulting in obvious economic and ecological benefits.

Keywords rubber plantation ； nutrient cycle ； fertilizer recommendation ； economic benefit ； ecological benefit

人工施肥、林下生物量和積累量以及其分解速率、降雨、養分有效性等均會影響橡膠園土壤養分循環。橡膠林下枯落物的分解速率較快，因此，林下生物量與積累量相對較低[1]，但各養分元素的分解釋放效率和歸還量亦不相同[2]，降雨對膠園養分的輸入、運轉、淋溶、利用效率等都有較大的影響，也是膠園土壤肥力的重要來源之一，其養分來自地表徑流、地下滲漏水、淋洗樹體分泌物以及枯枝、落葉、死根的分解釋放和溶解肥料，這對橡膠樹的生長產生了較大的影響[3]。土壤養分的有效性受pH[4]、施肥位置及深度[5]、肥料種類與形態[6]、土壤理化性狀[7]、養分元素間的相互作用[8]等因素的影響。另外，土壤微生物、砧木與接穗之間的親和關系、遺傳特性、年齡因素、林下間種[9-12]等都可影響養分的有效性和吸收量[13]。然而，膠園土壤肥力普遍偏低[14]，令人擔憂，這主要是長期割膠移走、土壤淋溶、水土流失[15-17]，而人為施肥不足造成的。橡膠研究所研究人員對海南、云南、廣東三大墾區膠園土壤進行了測評，結果表明：與20世紀七、八十年代相比，膠園土壤肥力急劇下降、酸化趨勢明顯。這一現象必須引起人們的警覺和足夠的重視。因此，開展膠園生態系統養分循環研究，指導生產科學合理施肥，對于維持膠園土壤地力、防止土壤酸化，維護農業生態平衡，保障橡膠產業可持續發展具有重要作用。

作物生長動態模擬通過適當的數學模型來定量描述各個子系統隨時間和相應生態環境變化而變化的狀態，并預測作物生長發育的未來過程，提供優化管理的決策建議。氮素對作物生長有著舉足輕重的影響，國外許多作物生長模型通過與氮素模型耦合，實現了氮素限制下的生產力動態模擬[18]，其中以CERES模型[19]、ORYZA模型[20]最具代表性。在人工林養分循環動態模擬方面，運用不同方法研究了養分循環的通量特征[21-23]，建立了養分循環的分室模型[24-25]、生理化學模型[26]、生物量模型NBM[27]、生態模型和施肥模型[28]、林木養分管理模型[29]、森林管理評估模型[30]和智能施肥決策系統[31]等；而FORCYTE（Forest nutrient cycling and Yield Trend Evaluator）[32]則是一個典型的森林生態系統養分循環模型研發與應用軟件。這對林業的科學經營管理與發展起到了巨大的科技支撐作用。橡膠養分循環的動態模擬與應用研究，是解決目前橡膠園營養診斷方法過時問題的重要而有效的新手段之一，對于完善橡膠平衡施肥和精準施肥理論具有重要的科學指導意義。

1 材料與方法

1.1 材料

以現行普遍推廣的乙烯利刺激割膠制度（d/3×ET）為背景，結合研究團隊前期研究建立的橡膠樹精準施肥決策計算機系統，開展以不同年齡段的海南主栽橡膠品系PR107、RRIM600、熱研7-33-97為研究對象，在海南儋州西聯、西慶、中國熱帶農業科學院試驗農場選擇樹齡6～29 a的10個林段，株行距為3 m×7 m，約476株/hm2。試驗共進行23個月。endprint

1.2 方法

1.2.1 試驗設計

考慮橡膠樹生物量的分配格局與樹齡有關，結合現有的位于試驗區正常生長的橡膠林段，選擇橡膠樹齡在6～29 a的林段，選擇橡膠無性系PR107樹齡為10、15、19、21、27 a，每種樹齡2個林段，共10個林段；橡膠無性系RRIM600樹齡為9、12、16、25、29 a，每種樹齡2個林段，共10個林段；橡膠無性系熱研7-33-97樹齡為6、10、14、20 a，每種樹齡2個林段，共8個林段。每個林段選擇一個樹位（300～400株），每個樹位為一個試驗小區，每個小區內設置3次重復，每一重復選擇橡膠樹70～80株。

1.2.2 橡膠樹施肥量效益計算分析

研究團隊在采樣數據建立的橡膠樹養分模型的基礎上，建立了橡膠樹精準施肥決策計算機系統，系統能夠綜合樹齡、季節、土壤等因素給出科學的精準施肥方案，本研究借助研究建立的系統進行無性系PR107、RRIM600、熱研7-33-97精準施肥診斷，計算分析橡膠樹施肥量效益。

1.2.3 產量分析

每個品種每種樹齡的2個林段分別設置原施肥標準和基于養分循環的施肥診斷決策系統新施肥標準，每個林段均連續測7～9個月，取月平均值。產量分析處理方法如表1所示。表1中的處理方式，“處理”代表“原施肥標準”，“對照”代表“新施肥標準”，下同。

1.2.4 投入產出分析

根據肥料尿素、過磷酸鈣、氯化鉀市場價格和干膠近10 a平均價格，分別計算出肥料投入成本、單株增產值、公頃增產值、單株凈增產值、公頃凈增產值。

1.2.5 土壤肥力分析

每個品種每個月均采集試驗小區的橡膠樹土壤，分別測定土壤中N、P、K、有機質含量。

1.2.6 膠園土壤pH分析

每個品種每個月均在試驗小區測定橡膠園土壤pH，取月平均值。

2 結果與分析

2.1 經濟效益分析

2.1.1 產量分析

2.1.1.1 PR107產量分析

PR107產量結果如表2所示。

由表2可知，利用橡膠樹精準施肥決策計算機系統診斷施肥后，不同年齡PR107平均干含提高了1.01%～2.34%，產量提高了0.72～1.06 kg/（株·年）、增幅21.72%～23.63%，其中樹齡10 a的平均增產比率為22.91%、單株增產0.09 kg，樹齡15 a的平均增產比率為22.94%、單株增產0.12 kg，樹齡19 a的平均增產比率為21.72%、單株增產0.13 kg，樹齡21 a的平均增產比率為23.63%、單株增產0.12 kg，樹齡27 a的平均增產比率為22.70%、單株增產0.14 kg。綜上，橡膠無性系PR107樹齡為10、15、19、21和27 a的新施肥標準的產量均高于原施肥標準的產量，增產效果明顯。

2.1.1.2 RRIM600產量分析

RRIM600產量結果如表3所示。

由表3可知，利用橡膠樹精準施肥決策計算機系統診斷施肥后，不同年齡RRIM600平均干含提高了1.93%～4.95%，產量提高了0.525～0.693 kg/（株·年）、增幅14.9%～28.6%，其中樹齡9 a的平均增產比率為26.65%、單株增產0.07 kg，樹齡12 a的平均增產比率為22.66%、單株增產0.07 kg，樹齡16 a的平均增產比率為23.70%、單株增產0.08 kg，樹齡25 a的平均增產比率為17.62%、單株增產0.08 kg，樹齡29 a的平均增產比率為13.82%、單株增產0.07 kg。綜上，橡膠無性系RRIM600樹齡為9、12、16、25和29 a的新施肥標準的產量均高于原施肥標準的產量，增產效果明顯。

2.1.1.3 熱研7-33-97產量分析

熱研7-33-97產量結果 如表4所示。

由表4可知，利用橡膠樹精準施肥決策計算機系統診斷施肥后，不同年齡熱研7-33-97橡膠樹干膠的單株產量在1.05～3.64 kg，對照（常規施肥）產量在0.86～2.85 kg，平均單株增產為0.45 kg，平均公頃增產為203.7 kg，平均增長率為21.01%，其中樹齡6 a的平均增產比率為22.09%、單株增產0.07 kg，樹齡10 a的平均增產比率為21.18%、單株增產0.50 kg，樹齡14 a的平均增產比率為27.72 %，單株增產0.79kg，樹齡20 a的平均增產比率為13.04%、單株增產0.33 kg。綜上，橡膠無性系熱研7-33-97樹齡為6、10、14、20 a的新施肥標準的產量均高于原施肥標準的產量，增產效果明顯。

2.1.2 投入產出分析

橡膠樹施肥量效益如表5所示。

表5的數據顯示， PR107單株肥力投入成本包括氮肥1.59元、磷肥0.34元、鉀肥0.53元，每株共計投入2.46元，單株增收17.20元，每公頃增收7 740元，單株凈增效益14.74元，每公頃凈增效益6 633元； RRIM600單株肥力投入成本包括氮肥2.04元、磷肥0.47元、鉀肥0.81元，每株共計投入3.32元，單株增收10.92元，每公頃增收4 911.9元，單株凈增效益7.60元，每公頃凈增效益3 420元；熱研7-33-97單株肥力投入成本包括氮肥0.87元、磷肥0.74元、鉀肥0.49元，每株共計投入2.10元，單株增收13.44元，每公頃增收6 074.95元，單株凈增效益11.34元，每公頃凈增效益5 169.75元。利用橡膠樹精準施肥決策計算機系統診斷，按時、按需科學施肥后，減少了肥料的浪費，直接經濟效益明顯增加。

2.2 生態效益分析endprint

2.2.1 土壤肥力分析

2.2.1.1 PR107膠園土壤肥力分析

PR107膠園土壤（0～40 cm）養分含量結果如表6所示。

表6的數據顯示，利用橡膠樹精準施肥決策計算機系統診斷施肥后，不同年齡橡膠無性系PR107 膠園土壤肥力有所提升，其中：硝態氮含量增加0.59～1.66 mg/kg、提高10.59%～32.8%；氨態氮含量增加0.65～1.85 mg/kg、提高10.96%～32.5%；速效磷含量增加0.54～0.76 mg/kg、提高15.74%～20.91%；速效鉀含量增加4.42～7.10 mg/kg、提高22.49%～23.08%；有機質含量增加0.02%～0.17%、提高1.21%～12.07%。表明通過新的營養診斷與施肥，補充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機質含量有所增加。

2.2.1.2 RRIM600膠園土壤肥力分析

RRIM600膠園土壤（0～40 cm）養分含量結果如表7所示。

表7的數據顯示，利用橡膠樹精準施肥決策計算機系統診斷施肥后，不同年齡橡膠無性系RRIM 600膠園土壤肥力有所提升，其中：硝態氮含量增加0.79～1.05 mg/kg、提高18.23%～47.33%；氨態氮含量增加1.02～2.30 mg/kg、提高19.97%～48.44%；速效磷含量增加1.04～4.22 mg/kg、提高36.61%～57.28%；速效鉀含量增加10.51～32.12 mg/kg、提高21.61%～56.73%；有機質含量增加0.18%～0.37%、提高23.70%～41.95%。該試驗地砂質土較重，土壤肥力相對貧瘠，故施肥后肥效迅速。通過新的營養診斷與施肥，補充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機質含量有所增加。

2.2.1.3 熱研7-33-97膠園土壤肥力分析

熱研7-33-97膠園土壤（0～40 cm）養分含量結果如表8所示。

表8的數據顯示，利用橡膠樹精準施肥決策計算機系統診斷施肥后，不同年齡橡膠無性系熱研7-33-97膠園土壤肥力有所提升，其中：硝態氮含量增加0.77～1.01 mg/kg、提高17.08%～25.21%；氨態氮含量增加1.33～1.46 mg/kg、提高23.07%～30.60%；速效磷含量增加1.52～2.33 mg/kg、提高35.27%～41.35%；速效鉀含量增加7.42～17.57 mg/kg、提高22.19%～42.83%；有機質含量增加0.14%～0.20%、提高13.26%～20.52%。試驗結果表明，通過新的營養診斷與施肥，補充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機質含量有所增加。

2.2.2 膠園土壤pH分析

膠園土壤pH結果如表9所示。

由表9可知，PR107膠園土壤pH增加0.03～0.09、提高0.58%～2.21%，其中樹齡10、15、19 a的pH均增加了0.03，樹齡21 a的pH增加了0.07，樹齡27 a的pH增加了0.09。RRIM600膠園土壤pH增加0.16～0.34、提高3.34%～6.99%，其中樹齡10 a的pH增加了0.16，樹齡15 a的pH增加了0.26，樹齡19 a的pH增加了0.18，樹齡21 a的pH增加了0.34，樹齡27 a的pH增加了0.21。熱研7-33-97膠園土壤 pH增加0.09～0.19、提高2.07%～4.5%，其中樹齡6 a的pH增加了0.09，樹齡10 a的pH增加了0.14，樹齡14 a的pH增加了0.10，樹齡20 a的pH增加了0.19。按照新施肥標準施肥，土壤pH均增加，膠園土壤酸化趨勢得到減緩。

3 結論與討論

3.1 結論

（1）借助已研究建立的橡膠樹精準施肥決策計算機系統和橡膠林生態系統養分循環試驗數據，PR107、RRIM600、熱研7-33-97三種品系的總體效果為：新施肥標準的產量均高于原施肥標準的產量，增產效果明顯；按時、按需科學施肥后，減少了肥料的浪費，直接經濟效益明顯增加；均補充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機質含量有所增加；土壤pH均增加，膠園土壤酸化趨勢得到減緩。

（2）采用新的施肥標準，PR107、RRIM600、熱研7-33-97三種品系的增產比率分別為：21.72%～23.63%、14.9%～28.6%、13.04%～27.72%，其中PR107增產效果比較穩定，增產差異不大；3種品系的投入產出比分別為16.69%、43.68%、18.52%，說明采用新標準施肥后，PR107和熱研7-33-97的效益較明顯；3種品系的土壤肥力提高分別10.59%～32.8%、18.23%～47.33%、17.08%～25.21%，其中RRIM600的土壤肥力提高幅度最大。

（3）采用新的施肥標準，樹齡分別為10、15、19、21、27 a的PR107的效果為：隨著樹齡的增長，增產比率呈下降趨；樹齡分別為9、12、16、25、29 a的RRIM600的效果為：隨著樹齡的增長，增產比率平穩；樹齡分別為6、10、14、20 a的熱研7-33-97的效果為：隨著樹齡的增長，樹齡14 a以下增產比率呈穩定上升趨勢，樹齡20 a的增產效果明顯下降。樹齡分別為10、15、19、21、27 a的PR107膠園土壤（0～40 cm）養分含量均有增加，樹齡越小，養分增加越明顯；樹齡分別為10、15、19、21、27 a的RRIM600膠園土壤（0～40 cm）養分含量均有增加，養分增加較穩定；樹齡分別為6、10、14、20 a的熱研7-33-97膠園土壤（0～40 cm）養分含量均有增加，養分增加較穩定。3種品系的膠園土壤pH均隨樹齡的增加而增加。endprint

3.2 討論

關于橡膠林生態系統養分循環效益分析方面的研究還鮮有報道，但在橡膠生物量數學模型估測方面，卻建立了較為成熟的數學模型。Shorrocks[33]首次建立了橡膠樹地上部分生物量與樹圍的冪函數回歸方程；胡耀華[34]采用標準木法伐樹測定了不同品種膠樹各部分的生物量，利用統計學中的最小二乘法求出膠樹1.5 m處莖圍與橡膠各部分生物量的關系方程式；周再知[35]采用FI（Furnivals index）指數，建立起了橡膠樹單木帶皮、去皮材積擬合方程，利用模型編制了橡膠樹立木材積表，并利用英國J. S. Huxley的“生物體各部分器官與測樹因子之間普遍存在著相對生長規律”理論，建立起了橡膠生物量測定的數學表達式及數學模型；趙春梅[36]利用數學模型對橡膠林生態系統養分循環通量及特征進行了研究，計算了要達到收支平衡各分室養分元素的盈虧數量及需要補充的數量。這些研究成果，為橡膠林養分循環與動態模擬研究打下了良好的基礎。

（1）新的施肥標準與傳統的標準相比，傳統生產上的肥料施用標準為，① 未開割幼齡樹：年施尿素190.4 kg/（hm2·a）[平均0.4 kg/（株·a）]；② 開割樹：尿素142.8 kg/（hm2·a）[平均0.3 kg/（株·a）]，專用復合肥（N 13%，P2O5 7%，K2O 10%）952 kg/（hm2·a）[平均2 kg/（株·a）]，有機肥2 kg/（株·a）。經折算后，尿素施用量約為408.4 kg/（hm2·a），P2O5約為66.4 kg/（hm2·a），K2O 95.21 kg/（hm2·a），這些施肥量與橡膠樹正常生長所需的養分數量還有相當的差距，造成了養分元素的不足。生產上這種固定配方、固定施肥量的管理措施，與不同年齡橡膠樹在一年中不同季節對養分的實際需求尚有差距，加上養分隨膠乳及更新后生物體外移出膠園生態系統，最終導致膠園土壤肥力狀況下降、地力衰退。利用橡膠樹精準施肥決策計算機系統診斷施肥，及時、足量補充了養分需求，使膠園土壤肥力有所提升，pH值有不同程度的增加、土壤酸化趨勢有所減緩，這對于恢復膠園土壤肥力、保障橡膠樹高產、穩產及可持續發展具有重要的生態意義。

（2）本研究分析了基于橡膠林生態系統養分循環指導施肥的效益，說明采用新的施肥標準是可以幫助提高土壤肥力，增產效果明顯，具有穩定的效益。

（3）新的施肥標準，能幫助減少因施肥不當造成的肥料浪費，減少對環境的污染，這對于保護生態環境具有積極意義。

參考文獻

[1] Krishnakumar A K， Chandra Gupta， Sinha R R， et al. Ecological impact of rubber （Hevea brasiliensis） plantations in North East India： 2. Soil properties and biomass recycling[J]. Indian Journal of Rubber Research. 1991， 4（2）： 134-141.

[2] Murbach M R， Boaretto A E， Muraoka T， et al. Nutrient cycling in a RRIM 600 clone rubber plantation[J]. Scientia Agricola， 2003， 60（2）： 353-357.

[3] 周再知，鄭海水，楊曾獎，等. 橡膠-砂仁復合系統生物產量、營養元素空間格局的研究[J]. 生態學報，1997，17（3）：225-233.

[4] Choudhury M， Sarma A C， Pal T K， et al. Available nutrient status of the rubber （Hevea brasiliensis） growing soils of Meghalaya， Mizoram， Manipur and southern parts of Assam[J]. Indian Journal of Hill Farming， 2001， 14（2）： 76-80.

[5] Hardjono A. Response of GT1 Hevea seedlings in the nursery on red-yellow podzolic soils to N rate and P placement[J]. Menara Perkebunan Indonesia， 1981， 49（6）： 153-158.

[6] George E S， Nair R B， Chitra M， et al. Influence of soil inorganic phosphorus fractions on leaf phosphorus status in Hevea brasiliensis[J]. Indian Journal of Natural Rubber Research， 2001， 14（1）： 75-78.

[7] Prasannakumari P， Joseph M， George S， et al. Movement of applied potassium in a sandy clay loam soil under rubber （Hevea brasiliensis） plantation[J]. Natural Rubber Research，2004， 17（2）： 168-171.

[8] Viegas I de J M， Sampaio M do C T， Costacurta C R， et al. Effect of calcium levels on dry matter production and macronutrients concentration in the tissue of young plants of rubber （Hevea spp）[J]. Revista de Ciencias Agrarias， 2001， 17 （35）： 25-40.endprint

[9] 陳 青，姚行成，周 珺，等. 不同覆蓋材料對幼齡膠園土壤養分及溫度的調節作用[J]. 廣東農業科學，2016（09）：72-76.

[10] 張志揚，譚 昕，陳 建，等. 不同郁閉度橡膠林間作豆科植物對土壤肥力的影響[J]. 湖北農業科學，2016（15）： 3 863-3 866.

[11] 李 娟，林位夫，周立軍. 長期間作咖啡對膠園土壤養分與土壤酶的影響[J]. 西南農業學報，2016（06）：1 371-1 374.

[12] 賀軍軍，姚艷麗，羅 萍，等. 幼齡膠園間種菠蘿對土壤肥力和種植效益的影響[J]. 西南農業學報，2016（02）：337-341.

[13] 曹建華，陶忠良，蔣菊生，等. 不同年齡橡膠樹PR107各器官養分利用效率研究[J]. 熱帶作物學報，2010，31（12）：2 091-2 097.

[14] 王巧環，劉志崴，王甲因. 海南農墾膠園土壤養分現狀及平衡施肥[J]. 熱帶農業科學，1999，19（5）：8-14.

[15] 趙春梅，蔣菊生，曹建華，等. 橡膠人工林養分循環通量及特征[J]. 生態學報，2009，29（7）：3 782-3 789.

[16] 陳桂良，王家銀，黎小清，等. 江城橡膠公司橡膠園土壤養分時空變異特征分析[J]. 西南農業學報，2017（01）：134-140.

[17] 曹建華，蔣菊生，陶忠良，等. 不同年齡橡膠樹膠乳礦質養分的流失[J]. 熱帶作物學報，2010，31（1）：1-5.

[18] Jongschaap R E E， ROTASK. In： Plentinger，M.C， Penning de Vries F W T. Graduate School for Production Ecology and Resource Conservation[C]. Wageningen University： The Netherlands. 1997： 81-86.

[19] Ritehie J. The CERES-MaizeModel[A]. In： Jones C A and Kiniry J R.CERES-Maize：A simulation model of maize growth and development[C]. Texas AM university Press. College station， 1986： 83-88.

[20] Bouman B A M， Kropff M J， Tuong T P， et al. Modeling lowdand rice[M]. Los Banos（Philippines）： Intemational Rice Research Institute， and Wageningen：Wageningen University and Research Center， 2001： 97-127.

[21] Tumer J， Lambert M J. Nut rient cycling in age sequences of two Eucalyptus plantation species[J]. Forest Ecology and Management， 2008， 255 （5/6）： 1 701-1 712.

[22] Juan A B， Bosco I. Nut rient return via lit terfall in two cont rasting Pinus sylvest ris forest s in the Pyrenees under different thinning intensities[J]. Forest Ecology and Management， 2008， 256： 1 840-1 852.

[23] 劉廣路，范少輝，漆良華. 閩西北不同類型毛竹林養分分布及生物循環特征[J]. 生態學雜志，2010，29（11）：2 155-2 161.

[24] 聶道平，沈國舫，董世仁. 油松人工林養分循環的研究Ⅲ. 養分元素生物循環和林分養分平衡[J].北京林業大學學報，1986，8（2）：8-13.

[25] 陳長青，卞新民，何園球. 中國紅壤坡地不同林地養分動態變化與模擬研究[J]. 土壤學報，2006，43（2）：240-246.

[26] Verburg P S J， Johnson D W A. spreadsheet-based biogeochemical model to simulate nutrient cycling processes in forest ecosystems[J]. Ecological Modelling， 2001， 141（1/3）： 185-200.

[27] Zhu Z X， Arp P A， Meng-FanRui， et al. A forest nutrient cycling and biomass model （For NBM） based on year-round， monthly weather conditions， Part I： Assumption， structure and processing[J]. Ecological-Modelling，2003， 169（2/3）： 347-360.

[28] Ingestad T. dFertilization moel. Ecosystem model of biomass production， nutrient cycling， and fertilization[J].Agro， 1986（5）： 1-2， 5-11.endprint

[29] Panpan Gao； Zhihong Qie； Xiangbin Kong. Simulation of Random dynamic changes of soil organic matter based on fuzzy Markov chain model[J]. Control and Decision Conference， 2010： 1 484-1 488.

[30] Blanco J A，Zavala M A， Bosco-Imbert J， et al. Sustainability of forest management practices： evaluation through a simulation model of nutrient cycling[J]. Forest Ecology and Management， 2005， 213（1/2/3）： 209-228.

[31] Yang Y s， Wang F L， Zhao J. Intelligent fertilization decision support system based on knowledge model and WebGIS： Decision for fertilization Computer Science and Information Technology[C]. 2nd IEEE International Conference on Digital Object Identifier： 2009， 232-235.

[32] 徐應德. 一個有用的森林生態系統、農林復合系統的模擬模型-FORCYTE[A]. 林農復合生態系統學術討論會論文集[C]. 東北林業人學出版社，1988：148-153.

[33] Hovind E L. Precipitation distribution around a windy mountain peak[J]. Journal of Geophysical Research， 1965， 70（14）： 3 271-3 278.

[34] 胡耀華，王 釗，舒宜通，等. 橡膠樹生物量分配及膠園生產率的研究[J]. 熱帶作物學報，1982，3（1）： 15-26.

[35] 周再知，鄭海水，尹光天，等. 橡膠樹生物量估測的數學模型[J]. 林業科學研究，1995，8（6）：624-629.

[36] 趙春梅，蔣菊生，曹建華. 橡膠人工林生態系統氮素循環模型[J]. 林業資源管理，2009（3）：66-70.endprint