不同林齡橡膠凋落物葉分解特性與有機碳動態(tài)研究

金龍 吳志祥 楊川 管利民 賴華英

摘 要 通過網袋法，研究一年內6 a、12 a、20 a、28 a共4組不同林齡橡膠凋落物葉分解及有機碳動態(tài)變化的過程。結果表明：（1）不同林齡橡膠凋落物葉均表現(xiàn)為慢-快-慢的分解趨勢，經1 a的分解，4組林齡組橡膠凋落葉殘留率分別為28.33%、30.53%、31.59%和33.89%。（2）不同林齡橡膠凋落物葉分解系數(shù)為1.376 2～1.608 8，平均值為1.498 2，大小順序為6 a>12 a>18 a>28 a，分解50%和95%所需時間分別為0.43、0.45、0.47、0.50 a和1.86、1.96、2.03、2.18 a。（3）林齡顯著影響橡膠凋落物葉初始N含量、P含量、木質素/N比、C/P比與木質素/P比，但對C含量、木質素含量及C/N比影響不明顯。（4）初始C/N比極顯著影響橡膠凋落物葉分解速率，而初始P含量對凋落物分解影響不顯著（5）4組林齡組橡膠凋落物葉有機碳濃度是一個先上升后下降的過程，研究結束時，4組林齡橡膠凋落物葉有機碳分別釋放了77.58%、79.41%、75.89%和75.02%，其動態(tài)變化過程顯著符合一元多項式模型，決定系數(shù)（R2）在0.88以上。

關鍵詞 橡膠林；不同林齡；凋落物；分解速率；有機碳

中圖分類號 S718 文獻標識碼 A

森林凋落物是森林生態(tài)系統(tǒng)的重要組成部分，同時也是地表及地下分解者的主要食物來源，其生產和分解過程是森林生態(tài)系統(tǒng)物質循環(huán)和能量流動的重要環(huán)節(jié)之一[1]。另外，森林凋落物還在改善土壤理化性質、提升林地生產力、減少水土流失及保溫增濕等方面起著重要作用[2-4]。從廣義上講，森林凋落物主要包括枯葉、枯枝、死亡根系和倒木等部分[5]，其中葉是主要成分，約占凋落物總量的60%～80%，甚至更高[1-6]。森林凋落物及其分解一直都是土壤學、生態(tài)學、環(huán)境科學及林學等學科的研究重點，至今已有100多年的歷史，自德國生物學家Ebermayer于1876年首次報道森林凋落物及其在養(yǎng)分循環(huán)中的作用以來，國外就陸續(xù)開展了生態(tài)系統(tǒng)內凋落物的總量、組成成分、凋落節(jié)律、凋落物分解與養(yǎng)分釋放及分解模型等方面的研究[7-9]。而中國有關凋落物研究的時間較短，直到20世紀80年代，中國學者才開始對凋落物的分解進行研究，主要涉及生物因素和非生物因素對凋落物分解的影響、養(yǎng)分動態(tài)變化研究[10-11]。20世紀90年代后，伴隨對全球氣候變化的廣泛關注，有關凋落物分解的研究被結合在了對全球碳收支平衡的調節(jié)控制上[12-13]。近10年來，隨地下生態(tài)學研究不斷深入及地下生態(tài)過程在響應全球變暖過程中的重要性，凋落物對地下生態(tài)過程調控機制也越來越受到關注[13]。但由于影響凋落物分解的因素復雜眾多，凋落物中各種化學組分及其相互作用對凋落物分解的影響目前仍不十分不清楚[6， 13]，尤其是在氣候變化的背景下，C、N養(yǎng)分及其交互作用對凋落物分解的調控機制還需更深入的研究[13]。因而，盡管已經有大量森林凋落物研究報道的資料[14-15]，但全球仍在持續(xù)深入研究各類型森林凋落物的特征及其分解對氣候的響應機制。

橡膠林生態(tài)系統(tǒng)是中國熱區(qū)一種最典型的人工林生態(tài)系統(tǒng)，從第九屆中國橡膠基礎研究研討會上獲悉，截至2012年，中國植膠面積達到了111萬hm2。在全球廣泛開展森林凋落物研究背景下，中國關于橡膠林凋落物的研究也取得了一定的進展。楊曾獎等[16]研究了27 a橡膠林間作砂仁模式下凋落物的組成及分解狀況；任泳紅等[17]研究了橡膠多層林凋落物的節(jié)律及分解變化；曹建華[18]報道了PR107品種橡膠3.0 m×7.0 m 栽培模式下9個年齡段（4～28 a）凋落物量及養(yǎng)分的釋放狀況；趙春梅等[19]研究了不同林齡橡膠枯落物的分解及氮素釋放過程等。Philip等[20]研究了15 a橡膠枯落物的分解與養(yǎng)分的釋放狀況。但對橡膠這類人工林凋落物的研究范圍仍較窄，主要集中在橡膠林凋落物量及養(yǎng)分釋放動態(tài)的研究，而有關養(yǎng)分的研究也主要集中在N、P、K、Ca、Mg等營養(yǎng)元素的釋放規(guī)律上，而對橡膠凋落葉基質質量與其分解速率的關系，分解過程中有機碳動態(tài)及其釋放量的研究還鮮有報道?；诖?，本研究通過野外分解袋法研究了不同林齡（6 a、12 a、20 a 和28 a）橡膠凋落葉分解過程及其有機碳動態(tài)變化的過程。旨在為準確評估橡膠林碳源/匯功能提供參考依據(jù)，同時也為林地養(yǎng)分管理、提高林地生產力等方面提供了理論依據(jù)。

1 材料與方法

1.1 試驗地概況

不同林齡橡膠凋落物葉分解實驗的試驗樣地設在海南儋州中國熱帶農業(yè)科學院試驗農場3隊內，地點位于19°31′47"N，109°29′30"E，平均海拔144 m，屬于熱帶季風氣候，年平均氣溫為20.5～28.5 ℃，最冷月平均氣溫為16.5～17.6 ℃，全年日平均氣溫≥10 ℃的積溫為8 500～9 100℃；全年旱雨兩季分明，5～10月為雨季，11月～翌年4月為旱季，年平均降雨量為1 607～2 000 mm，其中7、8、9共3個月降雨量占全年降雨量的70%以上；年平均相對濕度為83%；實驗區(qū)地形為緩坡丘陵（相對高度差<10 m），土壤為花崗巖母質化所形成的磚紅壤，土層厚度約為100 cm，pH值4.52～5.86。

該實驗區(qū)為第二代膠園，其中，分解實驗樣地林下隔行間作有用于壓青肥的天堂鳥，其他植被有從當年生草本到多年生小草本，包括弓果黍（Cyrtococcum patens）、飛機草（Chromolaene odorata）、地膽頭（Elephantopus scaber）、 含羞草（Mimosa pudica）等。

1.2 凋落物分解實驗

凋落物分解采用網袋法，同時為保證土壤性質和氣候環(huán)境的均勻一致性，本研究以試驗場內6 a、12 a、20 a 和28 a共4個林段（栽培品系分別為RR525、熱研7-33-97、熱研7-33-97和PR107）橡膠樹的新鮮枯落葉作為實驗材料，于2013年12月橡膠樹開始落葉后，分別收集4個不同年齡段新鮮凋落物葉，待樣品風干后，準確稱?。?0.00±0.50）g直接裝入大小為45 mm×35 mm的網袋（孔徑大小為1.2 mm×1.2 mm）內，于2014年1月10日分別隨機放入6～12林段（林齡為12 a，3 m×7 m常規(guī)種植模式）的地表，每林齡15袋，3個重復，共計60袋，并用枯落物將其覆蓋。2014年3月10日開始取第一次樣，然后每隔2個月取一次樣，共計5次，每次每個林齡各取回3袋，輕輕除去分解袋上附著的泥沙，并去除袋內根系等雜物，烘干至恒重后稱重，計算凋落物的分解速率；同時實驗室采用重鉻酸鉀法測定樣品的有機碳含量[21]，采用靛酚藍比色法測定全N[22]，采用鉬銻抗比色法測定全P[22]，采用溴乙酰-乙酸比色法測定木質素[23]。

1.3 質量殘留率與有機碳釋放速率的測定

凋落物分解過程中凋落物殘留率和有機碳分釋放速率采分別采用以下公式進行計算：

枯落物殘留率：Di=（Wi/W0）×100%，式中，Di為第i月的殘留率（%），Wi為第i月所取樣品的剩余重，W0為落葉初始質量。

有機碳釋放速率=[（初始濃度×初始質量）-（t時刻濃度

×t時刻質量）]/（初始濃度×初始質量）×100%。

1.4 基于啞變量的Olson指數(shù)衰減模型

目前研究凋落物分解的模型有多種，其中以Olson指數(shù)衰減模型應用最為廣泛，形式為：

Mt=M0e-kt，式中，M0是初始質量，Mt為某一時間t時的質量，k為分解速率常數(shù)。

目前針對不同林齡生長模型的預測時都是單獨建立模型，這樣建立的模型會造成整體和分量之間的不相容，增大了模型誤差。而采用啞變量（dummy variable）可把各林齡用定性代碼來表示，從而整合成一個模型來構建，這樣既減少了工作量又提高了模型精度。啞變量的定義為： 對于等級性（定性）數(shù)據(jù)x，用變量δ（x，i）表示成：

這種方法叫做定性因子（0，1）化展開，因此稱變量δ（x，i）為啞變量。一個定性變量（m個等級）對應一個向量δ（x，#）=[δ（x，1），…，δ（x，m）]。啞變量只取0或1，于是一個定性變量就變成取0或1的數(shù)值向量，便可用數(shù)值方法處理[24-25]。

本研究所收集的數(shù)據(jù)來自4個林齡，第i個林齡編號為Si，把定性數(shù)據(jù)Si轉化為（0，1）數(shù)據(jù)：

則可將Olson指數(shù)衰減模型改寫成啞變量形式如下：

y=（a1S1+a2S2+a3S3+a4S4）e[-（k1S1+k2S2+k3S3+k4S4）t]

式中，a1、a2、a3、a4是待定參數(shù)，S1、S2、S3和S4分別是6 a、12 a、20 a和28 a共4個林齡的代碼，k1、k2、k3和k4分別為4個林齡的分解系數(shù)，t為分解時間。

1.5 數(shù)據(jù)分析

采用Microsoft office excel 2007、origin9.0等軟件對所有取樣點數(shù)據(jù)進行整理、分析、計算、作圖。

2 結果與分析

2.1 不同林齡橡膠凋落物葉初始質量參數(shù)

凋落物葉化學組分是影響其分解的重要內在因子，不同林齡橡膠林凋落物葉初始質量參數(shù)見表1。不同林齡（未考慮品系因素，另有說明，下同）橡膠凋落葉有機碳含量介于49.29%～51.86%，平均含量為51.58%。橡膠凋落物葉中N、P濃度隨林齡的增長而呈波動變化，其中N濃度表現(xiàn)為28 a、12 a顯著低于6 a和20 a；P濃度表現(xiàn)為6 a、28 a顯著低于12 a和20 a。木質素含量、C/N比、木質素/N比和木質素/P比隨林齡變化均呈升高-降低趨勢，C/P比呈降低-升高趨勢，平均值分別為23.69%、24.52、10.97和154.62；方差分析結果表明，木質素濃度、C/N比在不同林齡間差異不顯著；木質素/N比為6 a顯著低于12 a、20 a、28 a；木質素/P比為6 a、12 a 顯著低于20 a、28 a；C/P比為6 a、28 a顯著低于12 a和20 a。綜上所述，橡膠凋落物葉的品質隨林齡的增加有所降低。

2.2 凋落物葉分解殘留率的月變化

通過分解袋法，計算樣品的殘留率來研究橡膠凋落葉的分解速度。由表2可知，4組不同林齡橡膠凋落物葉的殘留量隨時間的推移不斷降低，月質量損失率存在明顯的季節(jié)變化，表現(xiàn)為先升高后降低的趨勢，分解高峰出現(xiàn)在3～7月，分解高峰平均月分解速率為12.16%。經過將近1 a的分解，4組林齡橡膠凋落物葉最終殘留率隨林齡的增加表現(xiàn)出升高趨勢，分別為28.33%、30.53%、31.59%和33.89%，但不同林齡間分解速率差異不顯著（p>0.05）。說明林齡對橡膠凋落物葉分解速率有一定的影響，隨林齡的增加分解難度增加。

2.3 含有啞變量Olson指數(shù)衰減模型擬合

本研究利用含有啞變量的Olson指數(shù)衰減模型進行擬合，擬合結果為：

Y=（93.469S1+95.348S2+95.014S3+92.042S4）e[-（1.608 8S1+1.532 3S2+1.475 4S3+1.376 2S4）t]

式中S1、S2、S3和S4分別是6 a、12 a、20 a和28 a共4個林齡的代碼；t為分解時間。6 a、12 a、20 a和28 a共4個林齡的擬合參數(shù)為93.469、95.348、95.014、92.042；分解系數(shù)（K）分別為1.608 8、1.532 3、1.475 4和1.376 2。4個林齡的擬合結果均達到極顯著水平（p<0.01），決定系數(shù)在95%以上，擬合效果較好。比較不同林齡橡膠凋落物葉的分解系數(shù)，不同林齡之間K值的大小順序為6 a>12 a>20 a>28 a（表3）。根據(jù)該模型計算分解50%、95%所需的時間，結果顯示4林齡分解50%所需時間分別為0.43、0.45、0.47、0.50 a，分解95%所需時間為1.86、1.96、2.03、2.18 a，該值要略大于任泳紅等[17]和曹建華等[18]所報道的1.33 a和12個月（表3）。

2.4 凋落物葉初始質量與分解系數(shù)的相關性

在相同分解條件下，凋落物分解速率主要受凋落物化學成分的影響，且不同分解階段影響的因子不同，前期主要受養(yǎng)分含量、水溶性碳化合物和結構碳化合物含量的強烈影響，而后期則更多地受木質素及纖維素/木質素比的支配[26]。通過分析不同林齡橡膠凋落物葉分解系數(shù)與各主要化學成分之間的關系（表4），分解系數(shù)與初始C/N比呈極顯著負相關，與C、N、C/P比含量呈正相關，與木質素、P、木質素/N比、木質素/P比呈負相關，其中與木質素、木質素/N的相關性在0.7以上?？梢姡谕粴夂颦h(huán)境下，C/N比是引起不同林齡橡膠凋落物葉分解速率差異的最為主要原因。

2.5 凋落物葉分解過程中有機碳濃度的變化

由圖1可知，橡膠凋落物葉分解過程中，有機碳濃度表現(xiàn)出一定的變化趨勢：4組不同林齡橡膠凋落葉有機碳濃度在分解初期（1～3月）有一個上升的過程，且上升幅度在不同林齡間存在差異，其中，上升幅度以12 a最高，以20 a最低，相對初始濃度分別增加了10.46%和4.81%，6 a和28 a則分別增加了7.27%和6.18%；分解開始后的3～7月，4組林齡橡膠凋落葉有機碳濃度均表現(xiàn)為明顯的快速下降過程；7月后，除12 a在9～11月有一個突降外，其余3林齡組有機碳濃度趨于平穩(wěn)，取樣結束時，4組林齡橡膠凋落物葉有機碳濃度分別為39.66%、35.01%、38.24%和36.33%，較初始濃度分別下降了10.46%、16.85%、11.87%和12.96%。

2.6 凋落物葉分解過程中有機碳釋放速率的動態(tài)變化

由圖2可知，4組不同林齡橡膠凋落物葉有機碳釋放速率呈現(xiàn)以下特征：在分解初期（1～3月），有機碳釋放速率慢，甚至表現(xiàn)為凈積累過程，如6 a和12 a橡膠凋落葉，凈積累量分別為0.84%和7.48%；分解開始后的3～7月，為有機碳快速釋放階段；7月后，釋放速率開始放緩，整個過程中，4林齡橡膠凋落物葉有機碳釋放速率較為一致；至11月最后一次取樣時，總有機碳量分別釋放了77.58%、79.41%、75.89%和75.02%。

本研究運用6種不同模型對各林齡橡膠凋落物葉分解過程中有機碳釋放速率隨時間變化過程進行擬合，結果以一元多項式擬合方程最為理想（表5）。由表5可知，4組不同林齡橡膠凋落物葉在分解過程中，有機碳釋放速率隨時間的變化均顯著（p<0.05）符合一元多項式模型，決定系數(shù)在0.88以上。因此，用一元多項式方程對橡膠凋落葉有機碳釋放速率進行模擬可取得較好的效果。

3 討論與結論

3.1 不同林齡橡膠凋落物葉初始養(yǎng)分分析

N、P是植物體的重要組成部分，參與生命物質的合成。隨林齡的增加，橡膠凋落葉初始N含量、P含量均呈現(xiàn)出顯著性波動變化，其中N表現(xiàn)為6 a>20 a>12 a>28 a，P表現(xiàn)為12 a>20 a>28 a>6 a，總體養(yǎng)分及化學組分別隨林齡的增加有所降低。N、P等養(yǎng)分在不同林齡橡膠間的波動變化，不僅與土壤養(yǎng)分供應有關，還與養(yǎng)分轉移及在各器官之間的分配有關[27]。橡膠樹在幼齡期（6 a）由于生物量小，對N、P等養(yǎng)分總需求量低，且主要用于生物量的增長，因而單位體積土壤能提供相對充足的養(yǎng)分；進入開割期后（12 a、 20 a），雖然有一部分養(yǎng)分伴隨膠乳流失，但這種刺激可能會激發(fā)對養(yǎng)分的吸收，同時由于凋落物量也基本進入穩(wěn)定期，再加上人為施肥等措施能夠及時補充橡膠正常生長所需。對于老齡（28 a）即將更新樹，地上生物量大、人工扶管措施減弱，前期土壤肥力流失嚴重，因此導致對N、P等養(yǎng)分吸收競爭激烈，難度增大。植物體的元素化學組分主要來自于根系對土壤養(yǎng)分的吸收，因此，土壤質量的不同必然引起橡膠生長狀況的差異，而這種差異可能首先就反映在植物中心部位-葉上。另外，元素轉移在植物體內普遍存在，橡膠葉在掉落前，養(yǎng)分在不同林齡橡膠植物體內的轉移差異，也會影響到凋落葉的養(yǎng)分含量。而對于P為何6 a顯著性最低，是因為其余3林齡均為割膠林段，而割膠刺激可能激發(fā)了橡膠樹對P的吸收。研究結果還表明，林齡間的N、P養(yǎng)分顯著性差異還導致了木質素/N、木質素/P的顯著性差異。

3.2 不同林齡橡膠凋落物葉分解速率分析

4組林齡橡膠凋落物葉隨時間變化均呈慢-快-慢的分解過程，周轉時間隨林齡增大逐漸增長。在1 a的分解期內，橡膠凋落物葉分解速率是雨季分解快，旱季分解慢，前期分解快，后期分解慢，分解高峰出現(xiàn)在3～7月，這與前人研究的結論基本一致[19]。分解前期（1～2月）為一年中最冷月份，溫度低，降雨量少，微生物活性低，因淋溶作用和微生物代謝作用引起的凋落物質量損失少；夏季（雨季）高溫多雨，生物代謝和淋溶作用增強，因而凋落物分解速度快速增大；分解后期（7～11月），隨易分解物質的消耗殆盡，難分解物質大量積累，微生物可食用物質減少，凋落物葉上的生物開始出逃或則死亡，分解速度也開始減慢。在同一分解環(huán)境下，橡膠凋落葉品質是影響其分解速度的一個最為重要的因素。橡膠凋落葉品質隨林齡增加總體有所降低，其中，作為預測凋落物分解速度的C/N比[26]，隨林齡的增大逐漸增大。因此，不同林齡橡膠凋落葉品質差異是導致其分解快慢最為重要的因素。

經含有啞變量Olson指數(shù)衰減模型預測，4組林齡橡膠凋落葉徹底分解所需時間在1.86 a以上，高于任泳紅等[17]和趙春梅等[19]同類報道的1.33 a和12個月。這可能與林下微氣候環(huán)境、土壤生物及土壤肥力不同有關。本實驗林下由于間作用于膠園壓青的天堂鳥，其寬大密集的葉片對地表水熱通量起到一定的攔截，相比較凋落物質量，溫度和降水在凋落物分解中起主導作用[8]，因此，水、熱等環(huán)境因子的微小變化，必將對凋落物的分解產生重大影響。另外，地表間作天堂鳥可能會引起地表生物數(shù)量、群落結構及生物活性的改變，這也是影響凋落物分解速率差異的一個重要原因。

3.3 橡膠凋落物葉質量與分解速率的關系

橡膠凋落物葉初始C/N比極顯著影響分解速率，而初始P含量及C/P比對分解速率影響不顯著，均表現(xiàn)為弱相關性。已有研究結果發(fā)現(xiàn)，凋落物分解速率隨初始N、P濃度增加而增加，隨初始木質素濃度、C/N比和木質素/N比的增大而減小[26]。還有研究結果表明，受氮沉降的影響，P在熱帶地區(qū)凋落物分解中的重要性要大于N[28]。本研究結果發(fā)現(xiàn)，N、C/N比對橡膠凋落葉分解的重要性要強于P及C/P比，隨C/N比增加，分解速度逐漸降低，且C/N比與分解系數(shù)呈極顯著相關，可用來作為預測橡膠凋落物葉分解的指標。至于P及C/P比為何在橡膠凋落葉分解中沒有表現(xiàn)出相當?shù)闹匾?，甚至P表現(xiàn)為弱相關性。分析原因，可能是因為橡膠凋落物葉初始C/P比值（318.84～380.22）遠低于大多數(shù)凋落物分解的臨界C/P比值（1 000～1 300），而凋落物分解的臨界C/P比值就是分解微生物的臨界C/P比值[26]，因此，橡膠凋落物葉中較為豐富的P可能不是微生物分解的限制因子。另外，P的行為還與凋落物類型有關[29]。唐仕姍等[30]總結了中國不同森林類型凋落物研究成果后發(fā)現(xiàn)，中國主要森林類型的凋落物初始P濃度與分解系數(shù)的關系不顯著，其中，初始P濃度介于0.1～2.7 g/kg，且主要分布在小于1 g/kg的范圍內。本研究結果還發(fā)現(xiàn)，碳組分/養(yǎng)分與分解系數(shù)的相關性優(yōu)于單獨養(yǎng)分，相關性絕對值大小前4位為C/N>木質素/N>木質素>木質素/P，說明凋落物碳組分與N、P等養(yǎng)分的相互作用在調控橡膠凋落物葉分解過程中起著關鍵性作用。

3.4 橡膠凋落葉有機碳動態(tài)分析

4組林齡橡膠凋落物葉有機碳相對含量總體呈先上升后下降，最后趨于平穩(wěn)的態(tài)勢，其變化過程與凋落物的分解速率基本一致。有機碳相對含量是指有機碳絕對量與其干物質量的比值，有機碳濃度的上升下降過程，不僅體現(xiàn)了有機碳絕對量的釋放過程，同時也反映出了干物質的損失過程。凋落物分解前期主要為可溶性物質的淋溶損失過程，后期則為有機物的代謝損失過程[26]。前期分解過程中，由于作為有機碳主要成分的纖維素、木質素還未分解，有機碳損失量小于干物質損失量，因此有機碳濃度表現(xiàn)為上升過程；隨著纖維素、半纖維素等有機物大量分解，有機碳濃度也開始快速降低；最后殘留下的木質素等難分解的物質，限制了微生物的利用，導致干物質損失速率開始降低，有機碳濃度變化也基本趨于平穩(wěn)。同時，通過對比前人的研究結果，發(fā)現(xiàn)橡膠凋落葉有機碳動態(tài)變化過程與李正才等[31]報道的北亞熱帶6種人工林有機碳動態(tài)變化比較一致，與劉強等[32]報道的變化趨勢有所不同。另外，由于熱帶的水熱條件優(yōu)于亞熱帶，橡膠凋落葉分解1 a后的有機碳濃度變化量高于亞熱帶人工林樹種。同時，由于凋落物類型、分解環(huán)境及參與分解的生物類群、數(shù)量及活性差異，其釋放量低于尖峰嶺熱帶雨林等樣地。

3.5 橡膠凋落物葉有機碳釋放模型分析

4組林齡橡膠凋落物葉有機碳釋放速率顯著符合一元多項式模型，與凋落物分解指數(shù)衰減模型不同。橡膠凋落物葉有機碳釋放速率模型與干物質損失模型不同，說明有機碳釋放過程與凋落物分解過程不完全相同，這是因為對每種凋落物而言， 凋落物有機碳的分解與凋落物干物質降解的過程不同[31]；但對于橡膠凋落葉這一特定類型的凋落物，其有機碳分解快慢的順序與其干物質分解快慢的順序總體上較為一致， 說明橡膠凋落物有機碳的釋放和其干物質的分解機制存在一定的關系。和同類報道相比，與李正才等[31]、陳莉莎等[33]報道的有機碳指數(shù)釋放模型也不同，這可能是因為區(qū)域氣候環(huán)境及凋落物種類差異，導致凋落物分解速率及有機碳釋放過程中的調控機制不同有關。

本研究的不足之處主要是在分析林齡對橡膠凋落物葉分解及養(yǎng)分影響時，未考慮橡膠品系間的差異可能引起分解速率及化學組分含量的變化。但根據(jù)前人的研究結果，相同林齡不同品系的橡膠枯落物分解速率一致，養(yǎng)分含量也沒有顯著性差異變化[34-35]；相反，同一品系不同林齡橡膠枯枝、枯葉及鮮枝、鮮葉養(yǎng)分含量存在顯著差異[18]，因此，可推定不同林齡橡膠凋落物葉品質及分解速率差異主要由林齡效應引起。

參考文獻

[1] 于恩娜， 王金貴， 初寶順. 凋落物及其在森林生態(tài)中的作用[J]. 現(xiàn)代農業(yè)科技， 2009， 23（2）： 286-288.

[2] 薛 立， 何躍君， 屈 明，等. 華南典型人工林凋落物的持水特性[J]. 植物生態(tài)學報， 2005， 29（3）： 415-421.

[3] 肖 洋， 陳麗華，余新曉. 北京密云麻櫟人工混交林凋落物養(yǎng)分歸還特征[J]. 東北林業(yè)大學學報， 2010， 38（7）： 13-15.

[4] 楊麗萍，蘇海鵬. 橡膠林凋落物研究進展[J]. 熱帶農業(yè)工程，2011， 35（4）： 40-43.

[5] 趙鵬武， 宋彩玲， 蘇日娜，等. 森林生態(tài)系統(tǒng)凋落物研究綜述[J]. 內蒙古農業(yè)大學學報（自然科學版）， 2009， 30（2）： 292-299.

[6] 宋曰欽， 翟明普，賈黎明. 三倍體毛白楊地上凋落物對年齡的響應[J]. 東北林業(yè)大學學報， 2010， 38（3）： 17-19， 22.

[7] Bray J R， Gorham E. Litter production in forests of the word[J]. Advances in Ecological Research， 1964， 2（1）： 101-157.

[8] Aerts R. Climate， leaf litter chemistry and leaf litter decomposition in terrestrial ecosystems：a triangular relationship[J]. Oikos， 1997， 79（3）： 439-449.

[9] 林成芳， 高 人， 陳光水，等. 凋落物分解模型研究進展[J]. 福建林業(yè)科技， 2007， 34（3）： 227-233.

[10] 王 斌， 楊校生. 不同氣候區(qū)4種典型地帶性植被凋落物比較研究[J]. 世界林業(yè)研究， 2009， 29（S）： 44-48.

[11] 賈開心， 鄭 征，張一平. 西雙版納橡膠林生物量隨海拔梯度的變化[J]. 生態(tài)學雜志， 2006， 25（9）： 1 028-1 032.

[12] 彭少麟，劉 強. 森林凋落物動態(tài)及其對全球變暖的響應[J]. 生態(tài)學報， 2002， 22（9）： 1 534-1 544.

[13] 宋 飄， 張乃莉， 馬克平，等. 全球氣候變暖對凋落物分解的影響[J]. 生態(tài)學報， 2014， 34（6）： 1 327-1 339.

[14] Talbot J M， Lehtonen A. Needle and branch biomass turnover rates of Norway spruce（Picea abies）[J]. Canadian Journal of Forest Research， 2004， 34（12）： 2 517-2 527.

[15] 徐李亞， 楊萬勤， 李 晗，等. 高山森林林窗對凋落物分解過程中水溶性氮和磷的影響[J]. 水土保持學報， 2014， 28（3）：214-221.

[16] 楊曾獎， 鄭海水， 周再知，等. 橡膠間種砂仁模式下凋落物的特征[J]. 廣東林業(yè)科技， 1997， 13（4）： 19-24.

[17] 任泳紅， 曹 敏， 唐建維，等. 西雙版納季節(jié)雨林與橡膠多層林凋落物動態(tài)的比較研究[J]. 植物生態(tài)學報， 1999， 23（5）：418-425.

[18] 曹建華， 陶忠良， 趙春梅，等. 不同樹齡橡膠樹枯落物養(yǎng)分歸還比較[J]. 熱帶作物學報， 2011， 32（1）： 1-6.

[19] 趙春梅， 曹建華， 李曉波， 等， 橡膠林枯落物分解及其氮素釋放規(guī)律研究[J]. 熱帶作物學報， 2012， 33（9）： 1 535-1 539.

[20] Philip A， Philip V， George E S， et a1. Leaf litter decomposition and nutrient release in a fifteen year old rubber plantation[J]. Indian Journal of Natural Rubber Research， 2003， 16（1/2）： 81-84.

[21] 中國科學院南京土壤研究所. 土壤理化分析[M]. 上海： 上?？茖W技術出版社， 1978.

[22] 鮑土旦. 土壤農化分析（第三版）[M]. 北京：中國農業(yè)出版社，2000： 30-34.

[23] 李 靖， 程 舟， 楊曉伶，等. 紫外分光光度法測定微量人參木質素的含量[J]. 中藥材， 2006， 29（3）： 239-241.

[24] 王 欣， 高明達， 楊 飛，等. 不同林齡華北落葉松人工林葉凋落物分解及養(yǎng)分動態(tài)比較[J]. 東北林業(yè)學報， 2012， 40（10）： 55-60， 66.

[25] 李忠國， 孫曉梅， 陳東升，等. 基于啞變量的日本落葉松生長模型研究[J]. 西北農林科技大學學報：自然科學版， 2011， 39（8）： 69-74.

[26] Xu X， Hirata E. Decomposition patterns of leaf litter of seven common canopy species in a subtropical forest： N and P dynamics [J]. Plant and Soil， 2005， 273（2）： 279-289.

[27] 鄒積鑫， 曹建華， 蔣菊生，等. 橡膠樹無性系PR107不同樹齡葉片養(yǎng)分含量比較研究[J]. 熱帶農業(yè)科學， 2009， 29（10）： 9-12.

[28] Hastwell G. Nutrient cycling and limitation： hawaii as a model system[J]. Austral Ecology， 2005， 30（3）： 609-610.

[29] Moro M J， Domingo F. Litter decomposition in four woody species in a mediterranean climate： weight loss， N and P dynamics[J]. Annals of Botany， 2000， 86（6）： 1 065-1 071.

[30] 唐仕姍， 楊萬勤， 殷 睿，等. 中國森林生態(tài)系統(tǒng)凋落葉分解速率的分布特征及其控制因子[J]. 植物生態(tài)學報， 2014， 38（6）： 529-539.

[31] 李正才， 徐德應， 楊校生， 等. 北亞熱帶6種森林類型凋落物分解過程中有機碳動態(tài)變化[J]. 林業(yè)科學研究， 2008， 21（5）：675-680.

[32] 劉 強， 彭少麟， 畢 華，等. 熱帶亞熱帶森林凋落物交互分解的養(yǎng)分動態(tài)[J]. 北京林業(yè)大學學報， 2005， 27（1）： 24-32.

[33] 陳莉莎， 張金池， 莊家堯， 等. 蘇南丘陵區(qū)3種林分凋落葉客置分解特征[J]. 浙江農林大學學報， 2013， 30（5）： 681-688.

[34] 謝學方. PR107和熱研7-33-97橡膠人工林膠園土壤養(yǎng)分平衡研究[D]. 海南： 海南大學， 2011.

[35] 王麗華，林釗沐. 不同產量水平橡膠品種幼樹的氮磷鉀營養(yǎng)差異研究[J]. 廣東農業(yè)科學， 2009， 36（3）： 31-35.