基于橡膠林生態(tài)系統(tǒng)養(yǎng)分循環(huán)指導(dǎo)施肥的效益分析

張以山+袁曉軍+曹建華

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

關(guān)鍵詞 橡膠林 ；養(yǎng)分循環(huán) ；施肥指導(dǎo) ；經(jīng)濟(jì)效益 ；生態(tài)效益

中圖分類號 S794.1 文獻(xiàn)標(biāo)識碼 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

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

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

1 材料與方法

1.1 材料

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

1.2 方法

1.2.1 試驗(yàn)設(shè)計(jì)

考慮橡膠樹生物量的分配格局與樹齡有關(guān)，結(jié)合現(xiàn)有的位于試驗(yàn)區(qū)正常生長的橡膠林段，選擇橡膠樹齡在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株），每個樹位為一個試驗(yàn)小區(qū)，每個小區(qū)內(nèi)設(shè)置3次重復(fù)，每一重復(fù)選擇橡膠樹70～80株。

1.2.2 橡膠樹施肥量效益計(jì)算分析

研究團(tuán)隊(duì)在采樣數(shù)據(jù)建立的橡膠樹養(yǎng)分模型的基礎(chǔ)上，建立了橡膠樹精準(zhǔn)施肥決策計(jì)算機(jī)系統(tǒng)，系統(tǒng)能夠綜合樹齡、季節(jié)、土壤等因素給出科學(xué)的精準(zhǔn)施肥方案，本研究借助研究建立的系統(tǒng)進(jìn)行無性系PR107、RRIM600、熱研7-33-97精準(zhǔn)施肥診斷，計(jì)算分析橡膠樹施肥量效益。

1.2.3 產(chǎn)量分析

每個品種每種樹齡的2個林段分別設(shè)置原施肥標(biāo)準(zhǔn)和基于養(yǎng)分循環(huán)的施肥診斷決策系統(tǒng)新施肥標(biāo)準(zhǔn)，每個林段均連續(xù)測7～9個月，取月平均值。產(chǎn)量分析處理方法如表1所示。表1中的處理方式，“處理”代表“原施肥標(biāo)準(zhǔn)”，“對照”代表“新施肥標(biāo)準(zhǔn)”，下同。

1.2.4 投入產(chǎn)出分析

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

1.2.5 土壤肥力分析

每個品種每個月均采集試驗(yàn)小區(qū)的橡膠樹土壤，分別測定土壤中N、P、K、有機(jī)質(zhì)含量。

1.2.6 膠園土壤pH分析

每個品種每個月均在試驗(yàn)小區(qū)測定橡膠園土壤pH，取月平均值。

2 結(jié)果與分析

2.1 經(jīng)濟(jì)效益分析

2.1.1 產(chǎn)量分析

2.1.1.1 PR107產(chǎn)量分析

PR107產(chǎn)量結(jié)果如表2所示。

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

2.1.1.2 RRIM600產(chǎn)量分析

RRIM600產(chǎn)量結(jié)果如表3所示。

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

2.1.1.3 熱研7-33-97產(chǎn)量分析

熱研7-33-97產(chǎn)量結(jié)果 如表4所示。

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

2.1.2 投入產(chǎn)出分析

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

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

2.2 生態(tài)效益分析endprint

2.2.1 土壤肥力分析

2.2.1.1 PR107膠園土壤肥力分析

PR107膠園土壤（0～40 cm）養(yǎng)分含量結(jié)果如表6所示。

表6的數(shù)據(jù)顯示，利用橡膠樹精準(zhǔn)施肥決策計(jì)算機(jī)系統(tǒng)診斷施肥后，不同年齡橡膠無性系PR107 膠園土壤肥力有所提升，其中：硝態(tài)氮含量增加0.59～1.66 mg/kg、提高10.59%～32.8%；氨態(tài)氮含量增加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%；有機(jī)質(zhì)含量增加0.02%～0.17%、提高1.21%～12.07%。表明通過新的營養(yǎng)診斷與施肥，補(bǔ)充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機(jī)質(zhì)含量有所增加。

2.2.1.2 RRIM600膠園土壤肥力分析

RRIM600膠園土壤（0～40 cm）養(yǎng)分含量結(jié)果如表7所示。

表7的數(shù)據(jù)顯示，利用橡膠樹精準(zhǔn)施肥決策計(jì)算機(jī)系統(tǒng)診斷施肥后，不同年齡橡膠無性系RRIM 600膠園土壤肥力有所提升，其中：硝態(tài)氮含量增加0.79～1.05 mg/kg、提高18.23%～47.33%；氨態(tài)氮含量增加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%；有機(jī)質(zhì)含量增加0.18%～0.37%、提高23.70%～41.95%。該試驗(yàn)地砂質(zhì)土較重，土壤肥力相對貧瘠，故施肥后肥效迅速。通過新的營養(yǎng)診斷與施肥，補(bǔ)充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機(jī)質(zhì)含量有所增加。

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

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

表8的數(shù)據(jù)顯示，利用橡膠樹精準(zhǔn)施肥決策計(jì)算機(jī)系統(tǒng)診斷施肥后，不同年齡橡膠無性系熱研7-33-97膠園土壤肥力有所提升，其中：硝態(tài)氮含量增加0.77～1.01 mg/kg、提高17.08%～25.21%；氨態(tài)氮含量增加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%；有機(jī)質(zhì)含量增加0.14%～0.20%、提高13.26%～20.52%。試驗(yàn)結(jié)果表明，通過新的營養(yǎng)診斷與施肥，補(bǔ)充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機(jī)質(zhì)含量有所增加。

2.2.2 膠園土壤pH分析

膠園土壤pH結(jié)果如表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。按照新施肥標(biāo)準(zhǔn)施肥，土壤pH均增加，膠園土壤酸化趨勢得到減緩。

3 結(jié)論與討論

3.1 結(jié)論

（1）借助已研究建立的橡膠樹精準(zhǔn)施肥決策計(jì)算機(jī)系統(tǒng)和橡膠林生態(tài)系統(tǒng)養(yǎng)分循環(huán)試驗(yàn)數(shù)據(jù)，PR107、RRIM600、熱研7-33-97三種品系的總體效果為：新施肥標(biāo)準(zhǔn)的產(chǎn)量均高于原施肥標(biāo)準(zhǔn)的產(chǎn)量，增產(chǎn)效果明顯；按時、按需科學(xué)施肥后，減少了肥料的浪費(fèi)，直接經(jīng)濟(jì)效益明顯增加；均補(bǔ)充了原有施肥量的不足和施肥的不及時，減少了膠樹對膠園土壤肥力的掠奪，使土壤有機(jī)質(zhì)含量有所增加；土壤pH均增加，膠園土壤酸化趨勢得到減緩。

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

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

3.2 討論

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

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

（2）本研究分析了基于橡膠林生態(tài)系統(tǒng)養(yǎng)分循環(huán)指導(dǎo)施肥的效益，說明采用新的施肥標(biāo)準(zhǔn)是可以幫助提高土壤肥力，增產(chǎn)效果明顯，具有穩(wěn)定的效益。

（3）新的施肥標(biāo)準(zhǔn)，能幫助減少因施肥不當(dāng)造成的肥料浪費(fèi)，減少對環(huán)境的污染，這對于保護(hù)生態(tài)環(huán)境具有積極意義。

參考文獻(xiàn)

[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] 周再知，鄭海水，楊曾獎，等. 橡膠-砂仁復(fù)合系統(tǒng)生物產(chǎn)量、營養(yǎng)元素空間格局的研究[J]. 生態(tài)學(xué)報，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] 陳 青，姚行成，周 珺，等. 不同覆蓋材料對幼齡膠園土壤養(yǎng)分及溫度的調(diào)節(jié)作用[J]. 廣東農(nóng)業(yè)科學(xué)，2016（09）：72-76.

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

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

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

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

[14] 王巧環(huán)，劉志崴，王甲因. 海南農(nóng)墾膠園土壤養(yǎng)分現(xiàn)狀及平衡施肥[J]. 熱帶農(nóng)業(yè)科學(xué)，1999，19（5）：8-14.

[15] 趙春梅，蔣菊生，曹建華，等. 橡膠人工林養(yǎng)分循環(huán)通量及特征[J]. 生態(tài)學(xué)報，2009，29（7）：3 782-3 789.

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

[17] 曹建華，蔣菊生，陶忠良，等. 不同年齡橡膠樹膠乳礦質(zhì)養(yǎng)分的流失[J]. 熱帶作物學(xué)報，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] 劉廣路，范少輝，漆良華. 閩西北不同類型毛竹林養(yǎng)分分布及生物循環(huán)特征[J]. 生態(tài)學(xué)雜志，2010，29（11）：2 155-2 161.

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

[25] 陳長青，卞新民，何園球. 中國紅壤坡地不同林地養(yǎng)分動態(tài)變化與模擬研究[J]. 土壤學(xué)報，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] 徐應(yīng)德. 一個有用的森林生態(tài)系統(tǒng)、農(nóng)林復(fù)合系統(tǒng)的模擬模型-FORCYTE[A]. 林農(nóng)復(fù)合生態(tài)系統(tǒng)學(xué)術(shù)討論會論文集[C]. 東北林業(yè)人學(xué)出版社，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] 胡耀華，王 釗，舒宜通，等. 橡膠樹生物量分配及膠園生產(chǎn)率的研究[J]. 熱帶作物學(xué)報，1982，3（1）： 15-26.

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

[36] 趙春梅，蔣菊生，曹建華. 橡膠人工林生態(tài)系統(tǒng)氮素循環(huán)模型[J]. 林業(yè)資源管理，2009（3）：66-70.endprint