土地利用變化對東北溫帶幼齡林土壤碳氮磷含量及其化學計量特征的影響

周正虎，王傳寬，張全智

東北林業大學生態研究中心, 哈爾濱 150040

土地利用變化對東北溫帶幼齡林土壤碳氮磷含量及其化學計量特征的影響

周正虎，王傳寬*，張全智

東北林業大學生態研究中心, 哈爾濱 150040

土地利用方式的改變打破森林生態系統原有的碳氮磷平衡，從而顯著地影響森林生態系統的生物地球化學循環過程。以地段相鄰、林齡相同(10年生)、原始植被一致但土地利用方式不同(無土壤翻動的天然次生林[NS]、間作大豆而土壤翻動中等的人工林[MS]、間作人參而土壤翻動嚴重的次生林[SS])的溫帶幼齡林為對象，探索土地利用變化對土壤碳、氮、磷含量及相互關系的影響。結果顯示：(1)土地利用方式顯著改變表層和深層土壤碳含量，各土壤層次碳含量均呈現NS>MS>SS；而氮含量僅在0—20 cm具有顯著性差異(P<0.05)；不同土地利用類型之間磷含量無顯著差異(P>0.05);表明碳氮磷對土地利用變化敏感程度不同。(2)SS土壤碳氮比(C/N)和碳磷比(C/P)低于NS和MS，而NS和MS之間C/N和C/P因土壤層次而異。不同土壤層次氮磷比(N/P)均隨土壤翻動強度的增加而顯著減小(NS>MS>SS，P<0.05)，且隨土層加深而降低;表明N/P相對于C/N和C/P可能對土地利用變化具有更優生態指示功能。(3)土地利用變化顯著改變土壤碳-氮、碳-磷、氮-磷的耦合關系。土壤碳-氮(C-N)之間存在極顯著(P<0.001)的線性關系，其中3種土地利用方式的土壤C-N關系的斜率差異不顯著(P=0.458，共同斜率為11.1)，但截距差異顯著(P<0.001)。結合本地區和全球文獻數據分析指出，森林土壤碳氮關系既有大尺度上的普適性，又有小尺度上對土地利用方式響應的局域分異性。

土地利用變化; 生物地球化學; 化學計量; 碳氮耦合

人為活動造成土地利用方式的改變勢必會打破陸地生態系統原有的碳氮磷平衡[1-3]。土壤翻動不但破壞了土壤團聚體，而且也加速了土壤中動物植物微生物殘體、排泄物、分泌物以及土壤腐殖質的破碎速率，使其更易被微生物分解、化學氧化、物理釋放進入大氣和水域生態系統[4-5]，從而使土壤成為大氣碳庫的重要來源[6]。與此同時，土地利用變化也會造成氮磷元素的丟失[7-9]，但丟失程度取決于干擾的強度[1, 10]。例如，Guo等[1]綜合分析全球數據表明：天然林轉變為人工林和農田過程中土壤碳分別減少了13%和42%，而農田轉變為天然次生林和人工林時土壤碳分別增加了18%和53%；Xu等[10]研究表明，在天然林到第一代人工林再到第二代人工林轉變過程中土壤碳氮磷含量均顯著減小。

生態化學計量受國內外學者的廣泛關注，特別是研究碳、氮、磷化學計量特征可以反映土壤內部碳氮磷循環，具有重要的生態指示作用，因而對認識生態系統碳匯潛力及其對氣候變化的響應具有重要意義[11-15]。Tian[16]、Cleveland等[17]分別從區域尺度和全球尺度上探討土壤C/N、N/P和C/P化學計量比時得出：雖然在大尺度上不同生態系統之間土壤碳、氮、磷含量具有顯著的變異性，但C/N的空間變異較??；而由于磷元素和碳氮元素來源的差異性，土壤C/P和N/P隨氣候、土壤類型、土壤深度和風化階段而呈現顯著的空間變異性。而且，天然生態系統土壤C/N在不同層次間基本保持恒定[16, 18]。Yang等[18]研究也發現，天然高寒草原土壤碳-氮之間具有顯著的共斜率線性關系，且這種共斜率關系也存在于不同天然生態系統中[17]及其不同演替階段[19]。

東北地區歷經半個多世紀的采伐、墾荒等人為干擾，林地面積顯著減小，原始林基本殆盡，森林年齡結構趨于單一化和低齡化。隨著退耕還林、天然林保護等林業生態工程的實施，出現大面積的天然次生林和人工林[20-21]。然而，這種土地利用方式轉變對土壤碳氮磷含量、生態化學計量特征及其相互關系的影響還不清楚。本研究以帽兒山生態站的相鄰地段、相同林齡(均為10年生)、原始植被一致但土地利用方式不同的3種溫帶幼齡林為對象，最大程度地消除土壤本身異質性，以便探索土壤碳、氮、磷含量及其生態化學計量關系對土地利用變化的響應。

1 研究方法

1.1 研究地概況與試驗設計

研究地設在帽兒山森林生態站(45°20′N，127°30′E)。平均海拔400 m，平均坡度10°—15°， 地帶性土壤為暗棕色森林土。氣候屬大陸性季風氣候，夏季溫暖濕潤，冬季寒冷干燥，平均年降水量約629 mm，50%的降水集中在6—8月，平均年蒸發量約864 mm，年均氣溫3.1 ℃，1月份平均氣溫-18.5 ℃，7月份平均氣溫 22.0 ℃。無霜期約為120—140 d[22]。現有植被是原地帶性植被闊葉-紅松林經干擾后形成的天然次生林和人工林。

本研究選擇帽兒山生態站相鄰地段、相同林齡(土壤擾動后自然更新10a，林分現已郁閉)但土地利用方式不同的3種溫帶森林類型，分別為：無土壤翻動的天然次生林 [NS]、間作大豆而土壤翻動中等的人工林[MS]、間作人參而土壤翻動嚴重的次生林[SS]。NS為皆伐后天然更新次生林，土壤無任何擾動，其中優勢樹種為白樺(BetulaplatyphyllaSuk.)，伴生有春榆(Ulmusjaponica(Rehc.) Sarg.)、水曲柳(FraxinusmandshuricaRupr.)、大青楊(PopulousussuriensisKom.)、黃菠蘿(PhellodendronamurenseRupr.)和五角槭(AcermonoMaxim.)。MS為皆伐后人工栽植水曲柳(株行距1.5 m × 1.5 m)，同時翻土深度15 cm左右，間作大豆2a之后混入天然更新樹種白樺、春榆等。SS為皆伐后翻土深度40 cm左右，作壟種植人參10a后天然更新的次生林，其樹種組成與NS相似。于2013年9月在每林型中隨機設立3塊20 m × 30 m的固定樣地，并對樣地進行植被調查。

1.2 樣品采集與分析

2013年10月，在每塊樣地內隨機選取5個樣點，采用內徑為50 mm的土鉆分層取樣，取樣深度分別為：0—10 cm、10—20 cm、20—30 cm、30—40 cm、40—50 cm、50—60 cm。土壤樣品風干后全部過2 mm篩，除去所有根系和石礫，將土樣粉碎后作為分析土樣。稱取65 ℃烘干粉碎土樣約0.100 g，采用multi N/C 3000分析儀和HT 1500 Solids Module分析儀(Analytik Jena AG, Germany)燃燒法測定土壤全碳含量。另稱取烘干的粉碎土樣約0.200 g，加入5 mL濃硫酸，3 mL過氧化氫420 ℃消煮90 min(消煮完全至白色)，待冷卻后定容至100 mL，靜置12 h后吸取上清液利用連續流動分析儀(BRAN+LUEBBE-AA3, Germany)測定氮、磷含量。

1.3 數據分析

以樣地為統計分析試驗單位，采用單因素方差分析(one-way ANOVA)和最小顯著差異法(LSD)比較不同林型、土壤層次的碳、氮、磷含量及其生態化學計量特征；回歸分析法建立不同土地利用方式下土壤碳、氮、磷含量之間的最優耦合模型(以每樣方各土層平均土壤碳、氮、磷含量為數據分析單位)，并采用SMATR檢驗耦合模型的斜率同質性[23]。直線回歸系數差異性檢驗由SMATR Version 2.0軟件完成，其它統計分析均由SPSS19.0軟件完成。

2 結果

2.1 不同土地利用方式下土壤碳、氮、磷含量的垂直分布

圖1 不同土地利用類型下土壤碳、氮、磷含量的垂直分布Fig.1 Vertical distribution of soil carbon (C), nitrogen (N) and phosphorus (P) contents for different land use types誤差線為標準誤(n=3)；誤差線上的不同字母表示同一土壤層次不同土地利用間顯著性差異分組(α=0.05)；NS: 無土壤翻動的天然次生林Naturally regenerated stands; MS: 間作大豆而土壤翻動中等的人工林Moderately disturbed stands;SS: 間作人參而土壤翻動嚴重的次生林Severely disturbed stands

土地利用方式顯著影響土壤碳含量(圖1)。各層次土壤碳含量均呈現NS最高、MS居中、SS最低。NS各土壤層次碳含量平均比MS高19.4%(變化范圍為7.4%—32.4%)，比SS高77.4%(變化范圍為69.5%—83.1%)。其中，除0—10 cm外，其它各層次碳含量在NS和MS之間差異不顯著(P>0.05)，但兩者均顯著高于SS(P<0.05)。NS(F5,12=19.52,P<0.001)、MS(F5,12=38.71,P<0.001)和SS(F5,12=31.55,P<0.001)土壤碳含量均隨土層加深而顯著降低(圖1A)，但同一林型20 cm以下不同土層之間碳含量均無顯著性差異(P>0.05)。

與碳含量相同，各土壤層次氮含量均呈現NS>MS>SS，且具有相同垂直分布格局(圖1)。但是，除了0—10 cm(F2,6=7.18,P=0.026)、10—20 cm(F2,6=5.77,P=0.040)外，其它各土層氮含量在3種土地利用類型之間差異均不顯著(P>0.05)。NS各土壤層次氮含量平均比MS高38.0%(變化范圍為29.5% — 47.3%)，比SS高56.8%(變化范圍為51.8% — 60.8%)。

然而，與碳、氮含量不同，不同土地利用類型之間各土層的磷含量均無顯著差異(P>0.05)。NS和MS土壤磷含量隨著土層的加深而遞減，而SS的土壤磷含量呈現先遞減、20 cm以下保持基本穩定的趨勢(圖1)。

2.2 土地利用變化對土壤碳、氮、磷生態化學計量特征的影響

土地利用變化顯著影響土壤C/N, C/P和N/P(圖2)。其中，除0—10 cm外，其它各層C/N在不同土地利用方式之間具有顯著差異(P<0.05)，并呈現MS>NS>SS。NS不同土層之間C/N基本保持恒定(F5,12=0.923,P=0.499)，平均值為15.5；而MS的C/N隨著土壤深度的增加而顯著增加(F5,12=7.64,P=0.002)，變化范圍為13.4—28.6；SS的C/N隨著土壤深度的增加而減小，但不同土層之間無顯著差異(F5,12=1.744,P=0.199)，平均C/N為7.7。NS和MS各土層C/P均顯著高于SS(P<0.05)，而NS和MS之間僅在0—10 cm層具有顯著差異(P=0. 006)。土壤翻動顯著降低了土壤N/P，不同土地利用方式下N/P呈現：NS>MS>SS，且隨土層加深而顯著減小(P<0.05)。

圖2 不同土地利用類型下土壤C/N、C/P和N/P的垂直格局Fig.2 Vertical patterns of C/N,C/P and N/P for different land use types誤差線為標準誤(n=3)；誤差線上的不同字母表示同一土壤層次不同土地利用間顯著性差異分組(α=0.05)

土壤碳氮磷含量的耦合模型因元素和土地利用方式而異(圖3)。碳-氮(C-N)之間具有極顯著(P<0.001)的線性關系。進一步的SMATR斜率同質性檢驗表明，不同土地利用方式之間C-N化學計量關系斜率無顯著性差異(t=2.048，P=0.458)，共同斜率為11.1，NS和MS之間截距無顯著差異(P>0.05)，但兩者均顯著大于SS(P<0.001)(表1)。土壤碳磷(C-P)、氮磷(N-P)相互關系因土地利用方式而異，NS和MS呈現極顯著(P<0.001)的指數關系，而SS卻為直線關系(圖3)。

3 討論

3.1 土地利用方式對土壤碳、氮、磷含量及其垂直分布的影響

本研究中相鄰地段、相同林齡的不同土地利用方式下3種溫帶幼齡林，因土壤翻動歷史、程度和類型不同，其土壤碳、氮含量產生了顯著的差異。無論是表層還是深層，土壤翻動均促使土壤碳、氮含量的丟失(圖1)。其中，MS平均碳丟失率為19.4%，與Guo等[1]報道的全球平均值(13%)相當。這種中等程度土壤擾動導致的碳丟失主要發生在土壤表層(0—10 cm)，而對深層土壤的碳含量影響不顯著(圖1)，這是因為間作大豆2a過程中只翻動了0—15 cm，深層土壤結構沒有受到破壞，但對表層土壤的翻動過程中無疑增加了土壤的通透性，一定程度上促進了深層土壤碳的釋放，從而呈現出輕微的碳丟失，這與以往的研究相符[24-25]。然而，翻動強度大(作壟深達40 cm左右)、時間長(10a)的SS平均碳丟失率高達77.4%，遠高于Guo等[1]報道的全球森林向農田轉變過程中42%的碳丟失率，也高于Ellert等[8]報道的長期耕作導致34%的碳丟失率。這種重度土壤翻動顯著促使了表層和深層土壤的碳丟失，這是因為種植人參過程中作壟可達40 cm處，一方面土壤結構遭到嚴重破壞，另一方面也加大了深層土壤和外界環境的接觸機會與面積，從而促進了微生物對土壤有機質的分解導致的碳釋放以及隨地表徑流的丟失。此外，Don等[26]整合分析也表明土地利用變化過程中土地利用變化的類型、土壤水熱狀況以及土壤屬性(如土壤類型，粘土含量等)只解釋了55%的土壤碳變化，而將近一半的變異很可能是由土壤翻動強度的差異性所致。

圖3 不同土地利用類型下土壤碳、氮、磷相互關系Fig.3 Relationships between soil carbon (C), nitrogen (N) and phosphorus (P) contents for different land use types誤差線為標準誤(n=18)

表1 不同土地利用方式下土壤碳-氮關系斜率、截距同質性檢驗

將本研究結果與本地區的頂極群落和成熟林相比：NS的碳含量相對于天然次生成熟林(蒙古櫟林、硬闊葉林、楊樺林[27])的表層土壤(0—10 cm)碳丟失率為6.6%；而3種天然次生成熟林表層土壤平均碳含量比本地區闊葉紅松原始林[3]低47.9%。與碳含量相同，各土壤層次氮含量也均呈現NS>MS>SS(圖1)。NS各土壤層次氮含量平均比MS高38.0%，比SS高56.8%，可見土壤氮含量對土地利用方式的響應與碳含量相似。土壤翻動造成碳氮含量的同時丟失已得到眾多學者的證實[7-9]，這是因為土壤有機質在分解過程中每釋放25份碳素就會消耗1份氮素[28]。綜上所述，在溫帶森林從原始林-天然次生林-人工林-農田轉變過程中，土壤碳氮含量不斷降低。因此，采用適當的退耕還林和森林更新措施可以逆轉土壤的碳氮丟失，從而充分發揮森林土壤的碳氮固持潛能[26, 29]。

土壤碳和氮主要來源于植物地上地下凋落物的分解，其含量受植被、氣候、人為活動等影響較大[28]。因此，土壤對土壤碳氮的垂直分布的影響會因土地利用而發生變化，但其主導能力隨著土壤深度的增加而降低；而深層土壤碳氮含量主要決定于土壤本身質地，如粘土含量、土壤類型等[18, 30-32]。然而，長期的重度土壤翻動也會導致深層碳氮含量減小(圖1)。這表明量化在退耕還林過程中長期耕作農田的深層土壤碳氮匯功能也是必要的[26, 33]。與土壤碳氮不同，土壤磷主要來源于巖石風化，其含量大小更大程度上取決于立地的地質特征，而受人為影響較小[3, 34]。本研究中0—30 cm土層磷含量呈現NS>MS>SS，這可能是土壤翻動促使土壤有機磷的分解，并隨地表徑流而流失的緣故[32]。

3.2 土壤碳、氮、磷生態化學計量特征對土地利用變化的響應

Sistla等[35]指出，從局部到整體研究生態化學計量的可塑性和適應性是理解生態系統碳和養分循環對全球變化響應的有效途徑。與Tian[16]、Cleveland等[17]大尺度研究結果不同，本研究中NS和MS土壤C/N隨土層加深而呈現出增加的趨勢。Hamer等[36]也報道了C/N的這種垂直分布格局。進一步分析劉爽等[27]報道的當地5種溫帶成熟林結果也顯示，10—20 cm土壤C/N(20.1)高于0—10 cm(16.7)，而C/N與土壤有機質分解速率成反比[14]，這說明本地區深層土壤具有較高的碳固存潛力[37]，且因本地區深層土壤溫度較低、土壤微生物分解速率較慢、粘粒含量較高而相對穩定。NS和MS之間的差異性可能是因為間作大豆2a中對土壤的翻動促進了表層有機質向深層的輸入而補充了翻動導致的有機質輸出，而在這些動態變化過程中碳對外界干擾較氮更敏感[38]；另一方面本研究中NS和MS具有不同的植被組成，而植被特征一定程度上會影響土壤C/N的垂直分布。例如，最近王建林等[39]發現，不同高寒草地類型下土壤C/N具有5種不同的剖面分布型。此外，也可能是由于深層土壤碳氮的小尺度空間異質性所導致[40]。種植人參過程中的重度土壤翻動導致了C/N的顯著減小，經過10a的恢復，SS僅表層土壤C/N和NS與MS無顯著差異，這說明盡管本地區森林深層土壤具有穩定而又較高的碳含量[37]，但是在土地利用變化過程中一旦發生碳丟失，便需要相當長的時間才能恢復。

同樣，土地利用變化也顯著影響土壤C/P和N/P。3種不同土地利用方式下的C/P的剖面分布與C/N類似。相對于C/N、C/P，N/P垂直分布，與Tian等[16]研究結果相符，且不同土壤層次之間均呈現NS>MS>SS。這說明在消除了土壤本身異質性的小尺度上，土壤N/P可能比C/N和C/P與土地利用變化具有更高的同步性(隨土壤翻動強度的增加而顯著減小，圖2)。林麗等[13]也發現，土壤N/P對草地退化演替的敏感性較高。因此，N/P比C/N和C/P對外界環境變化可能具有更優的生態指示功能。表層土壤C/N、C/P和N/P均呈現SS < MS < NS，說明在森林恢復過程中會受到氮和磷限制。然而，盡管Koerselman[41]、Tessier[42]等先后通過綜述植被N/P來分析其對養分限制的指示作用，但由于受土壤本身碳氮磷絕對含量的影響，目前有關土壤碳氮磷化學計量對土壤養分限制的指示作用還無定論[14]。

本研究結果表明，3種林型土壤碳-氮-磷之間存在極顯著的耦聯關系(圖3)。Cleveland等[17]在綜述陸地生態系統土壤和微生物的“Redfield ratio[43]”時、Yang等[19]在綜述全球森林演替過程中土壤碳氮動態變化以及Tian等[16]在綜合我國土壤碳氮磷化學計量特征時均發現土壤碳氮之間存在著顯著的耦合關系。進一步的回歸系數差異性分析顯示，土壤碳-氮之間均存在顯著的共斜率線性關系，但土地利用變化顯著影響其截距(P<0.001，表2)。這表明土地利用變化并沒有顯著改變土壤碳-氮關系，而是通過改變其直線關系的截距來影響土壤碳氮的丟失以及C/N，這與最近研究相符[34, 44-45]。

圖4 普適性土壤碳-氮關系 Fig.4 Generalized relationships between soil carbon (C) and nitrogen (N)NS、MS、SS分別表示本研究中不同土地利用方式幼齡林；MF表示本地區的成熟林(mature forest)土壤碳-氮關系，土壤碳氮數據來自劉爽等[27]；Global表示全球森林生態系統土壤碳-氮關系，數據來源于Cleveland等[17]；方框中數據表示異常值

為了進一步驗證森林土壤碳氮化學計量關系的普適性，將本研究中的幼齡林和同一地區的5種典型森林類型的成熟林(即紅松人工林、落葉松人工林、硬闊葉林、蒙古櫟和楊樺林)[27]以及Cleveland等[17]的全球不同森林生態系統土壤碳氮數據作綜合分析。結果發現，不同尺度的土壤碳-氮之間線性關系的斜率均無顯著性差異(t=2.336,P=0.749)，其共同斜率為11.2(95%置信區間為[10.3, 11.8]，P<0.001，圖4)。全球尺度上碳-氮之間線性關系的截距與0無顯著性差異(95%置信區間為[-4.5, 15.9])。這些結果反映了土壤C/N在大尺度上恒定性，這種森林土壤碳氮關系的普適性在陸地生態系統生物地球化學循環研究中具有重要意義[18]。然而，大尺度上的異質性又掩蓋了C-N關系截距的差異性，這從另一方面說明小尺度上研究碳-氮之間線性關系截距的差異性也具有重要的生態指示作用。

與碳-氮關系不同，在最大程度消除土壤本身異質性的情況下，土壤碳-磷、氮-磷之間的相互關系因土地利用方式而異。NS和MS呈現極顯著的指數關系，而SS卻為直線關系(圖3)。Chen等[3]發現土壤磷含量對土地利用變化的響應與碳和氮并不相同，土壤全磷含量變化的不確定性同樣在落演替過程中發現[13]。這說明磷元素和碳氮元素不僅在來源上不同，且對外界干擾的響應機制可能也不同[16]，從而導致碳-磷、氮-磷關系的局域性。

4 結論

(1)土壤翻動會導致整個土壤剖面碳含量的顯著減小，而氮含量僅在0—20 cm顯著變化，而土壤磷含量保持恒定，說明土壤碳氮磷含量對土地利用變化響應機制可能不同。NS各土壤層次碳含量平均比MS和SS分別高19.4%和77.4%，氮含量平均比MS和SS分別高38%和56.8%，表明土壤翻動強度越大、持續時間越長，造成土壤碳氮丟失越嚴重。

(2)相比與NS和MS，各土層SS的C/N、C/P均最低，而NS和MS之間大小因土壤層次而異。3種土地利用方式下C/N和C/P具有不同的垂直分布。與C/N和C/P不同，各土層N/P與土地利用變化同步(NS>MS>SS)，且均為從高到低剖面分布。表明土壤碳氮磷學計量比具有可塑性與適應性，氮磷化學計量比作為土地利用變化的敏感指標優于碳氮和碳磷計量比。

(3)不同土地利用類型的C-P和N-P具有不同的耦合關系。其中C-N關系的斜率無顯著性差異，而其截距的差異顯著。這反映了森林土壤碳氮關系在大尺度上的普適性和小尺度上對土地利用方式的局域差異性。

致謝：帽兒山森林生態站提供野外基礎支持。

[1] Guo L B, Gifford R M. Soil carbon stocks and land use change:a meta analysis. Global Change Biology, 2002, 8(4):345-360.

[2] Post W M, Kwon K C. Soil carbon sequestration and land-use change:processes and potential. Global Change Biology, 2000, 6(3):317-327.

[3] Chen X W, Li B L. Change in soil carbon and nutrient storage after human disturbance of a primary Korean pine forest in Northeast China. Forest Ecology and Management, 2003, 186(1/3):197-206.

[4] Baker J M, Ochsner T E, Venterea R T, Griffis T J. Tillage and soil carbon sequestration—What do we really know? Agriculture, Ecosystems & Environment, 2007, 118(1/4):1-5.

[5] Cambardella C A, Elliott E T. Carbon and nitrogen distribution in aggregates from cultivated and native grassland soils. Soil Science Society of America Journal, 1993, 57(4):1071-1076.

[6] Lal R. World cropland soils as a source or sink for atmospheric carbon. Advances in Agronomy, 2001, 71:145-191.

[7] Ding F, Hu Y L, Li L J, Li A, Shi S W, Lian P Y, Zeng D H. Changes in soil organic carbon and total nitrogen stocks after conversion of meadow to cropland in Northeast China. Plant and Soil, 2013, 373(1/2):659-672.

[8] Ellert B H, Gregorich E G. Storage of carbon, nitrogen and phosphorus in cultivated and adjacent forested soils of Ontario. Soil Science, 1996, 161(9):587-603.

[9] Gál A, Vyn T J, Michéli E, Kladivko E J, McFee W W. Soil carbon and nitrogen accumulation with long-term no-till versus moldboard plowing overestimated with tilled-zone sampling depths. Soil and Tillage Research, 2007, 96(1/2):42-51.

[10] Xu Z H, Ward S, Chen C R, Blumfield T, Prasolova N, Liu J X. Soil carbon and nutrient pools, microbial properties and gross nitrogen transformations in adjacent natural forest and hoop pine plantations of subtropical Australia. Journal of Soils and Sediments, 2008, 8(2):99-105.

[11] 程濱, 趙永軍, 張文廣, 安樹青. 生態化學計量學研究進展. 生態學報, 2010, 30(6):1628-1637.

[12] 李從娟, 雷加強, 徐新文, 唐清亮, 高培, 王永東. 塔克拉瑪干沙漠腹地人工植被及土壤CNP的化學計量特征. 生態學報, 2013, 33(18):5760-5767.

[13] 林麗, 李以康, 張法偉, 杜巖功, 郭小偉, 李婧, 劉淑麗, 曹廣民. 高寒矮嵩草群落退化演替系列氮、磷生態化學計量學特征. 生態學報, 2013, 33(17):5245-5251.

[14] 王紹強, 于貴瑞. 生態系統碳氮磷元素的生態化學計量學特征. 生態學報, 2008, 28(8):3937-3947.

[15] 王維奇, 徐玲琳, 曾從盛, 仝川, 張林海. 河口濕地植物活體-枯落物-土壤的碳氮磷生態化學計量特征. 生態學報, 2012, 31(23):7119-7124.

[16] Tian H Q, Chen G S, Zhang C, Melillo J M, Hall C A S. Pattern and variation of C:N:P ratios in China′s soils:a synthesis of observational data. Biogeochemistry, 2010, 98(1/3):139-151.

[17] Cleveland C C, Liptzin D. C:N:P stoichiometry in soil:is there a “Redfield ratio” for the microbial biomass? Biogeochemistry, 2007, 85(3):235-252.

[18] Yang Y H, Fang J Y, Guo D L, Ji C J, Ma W H. Vertical patterns of soil carbon, nitrogen and carbon:nitrogen stoichiometry in Tibetan grasslands. Biogeosciences Discussions, 2010, 7(1):1-24.

[19] Yang Y H, Luo Y Q. Carbon:nitrogen stoichiometry in forest ecosystems during stand development. Global Ecology and Biogeography, 2011, 20(2):354-361.

[20] Chen X W, Li B L, Lin Z S. The acceleration of succession for the restoration of the mixed-broadleaved Korean pine forests in Northeast China. Forest Ecology and Management, 2003, 177(1/3):503-514.

[21] Deng X Z, Jiang Q O, Zhan J Y, He S J, Lin Y Z. Simulation on the dynamics of forest area changes in Northeast China. Journal of Geographical Sciences, 2010, 20(4):495-509.

[22] Wang C K, Han Y, Chen J Q, Wang X C, Zhang Q Z, Bond-Lamberty B. Seasonality of soil CO2efflux in a temperate forest:Biophysical effects of snowpack and spring freeze-thaw cycles. Agricultural and Forest Meteorology, 2013, 177:83-92.

[23] Warton D I, Wright I J, Falster D S, Westoby M. Bivariate line-fitting methods for allometry. Biological Reviews, 2006, 81(2):259-291.

[24] Mishra U, Ussiri D A N, Lal R. Tillage effects on soil organic carbon storage and dynamics in Corn Belt of Ohio USA. Soil and Tillage Research, 2010, 107(2):88-96.

[25] Tiessen H, Stewart J W B, Bettany J R. Cultivation effects on the amounts and concentration of carbon, nitrogen, and phosphorus in grassland soils. Agronomy Journal, 1982, 74(5):831-835.

[26] Don A, Schumacher J, Freibauer A. Impact of tropical land-use change on soil organic carbon stocks-a meta-analysis. Global Change Biology, 2011, 17(4):1658-1670.

[27] 劉爽, 王傳寬. 五種溫帶森林土壤微生物生物量碳氮的時空格局. 生態學報, 2010, 30(12):3135-3143.

[28] Chapin III F S, Matson P A, Vitousek P M. Principles of Terrestrial Ecosystem Ecology. 2nd ed. New York:Springer, 2011.

[29] Meyfroidt P, Rudel T K, Lambin E F. Forest transitions, trade, and the global displacement of land use. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(49):20917-20922.

[30] Jobbágy E G, Jackson R B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecological Applications, 2000, 10(2):423-436.

[31] Jobbágy E G, Jackson R B. The distribution of soil nutrients with depth:global patterns and the imprint of plants. Biogeochemistry, 2001, 53(1):51-77.

[32] Meersmans J, Van Wesemael B, De Ridder F, Fallas Dotti M, De Baets S, Van Molle M. Changes in organic carbon distribution with depth in agricultural soils in northern Belgium, 1960-2006. Global Change Biology, 2009, 15(11):2739-2750.

[33] Rumpel C, K?gel-Knabner I. Deep soil organic matter—a key but poorly understood component of terrestrial C cycle. Plant and Soil, 2011, 338(1/2):143-158.

[34] Wang W J, Qiu L, Zu Y G, Su D X, An J, Wang H Y, Zheng G Y, Sun W, Chen X Q. Changes in soil organic carbon, nitrogen, pH and bulk density with the development of larch (Larixgmelinii) plantations in China. Global Change Biology, 2011, 17(8):2657-2676.

[35] Sistla S A, Schimel J P. Stoichiometric flexibility as a regulator of carbon and nutrient cycling in terrestrial ecosystems under change. New Phytologist, 2012, 196(1):68-78.

[36] Hamer U, Potthast K, Burneo J I, Makeschin F. Nutrient stocks and phosphorus fractions in mountain soils of Southern Ecuador after conversion of forest to pasture. Biogeochemistry, 2013, 112(1/3):495-510.

[37] 楊金艷, 王傳寬. 東北東部森林生態系統土壤碳貯量和碳通量. 生態學報, 2005, 25(11):2875-2882.

[38] 羅亞勇, 張宇, 張靜輝, 卡召加, 尚倫宇, 王少影. 不同退化階段高寒草甸土壤化學計量特征. 生態學雜志, 2012, 31(2):254-260.

[39] 王建林, 鐘志明, 王忠紅, 陳寶雄, 余成群, 胡興祥, 沈振西, 大次卓嘎, 張憲洲. 青藏高原高寒草原生態系統土壤氮磷比的分布特征. 應用生態學報, 2013, 24(12):3399-3406.

[40] 武小鋼, 郭晉平, 田旭平, 楊秀云. 蘆芽山亞高山草甸、云杉林土壤有機碳、全氮含量的小尺度空間異質性. 生態學報, 2013, 33(24):7756-7764.

[41] Koerselman W, Meuleman A F M. The vegetation N:P ratio:a new tool to detect the nature of nutrient limitation. Journal of Applied Ecology, 1996, 33(6):1441-1450.

[42] Tessier J T, Raynal D J. Use of nitrogen to phosphorus ratios in plant tissue as an indicator of nutrient limitation and nitrogen saturation. Journal of Applied Ecology, 2003, 40(3):523-534.

[43] Redfield A C. The biological control of chemical factors in the environment. American Scientist, 1958, 46(3):205-221.

[44] Li H, Li J, He Y L, Li S J, Liang Z S, Peng C H, Polle A, Luo Z B. Changes in carbon, nutrients and stoichiometric relations under different soil depths, plant tissues and ages in black locust plantations. Acta Physiologiae Plantarum, 2013, 35(10):2951-2964.

[45] Tischer A, Potthast K, Hamer U. Land-use and soil depth affect resource and microbial stoichiometry in a tropical mountain rainforest region of southern Ecuador. Oecologia, 2014, 175(1):375-393.

The effect of land use change on soil carbon, nitrogen, and phosphorus contents and their stoichiometry in temperate sapling stands in northeastern China

ZHOU Zhenghu, WANG Chuankuan*, ZHANG Quanzhi

CenterforEcologicalResearch,NortheastForestryUniversity,Harbin150040,China

Anthropogenic land use change alters the equilibrium among soil carbon (C), nitrogen (N), and phosphorus (P) in the forest ecosystem, likely affecting the biogeochemical cycles. In this study, conducted in northeastern China, we measured C, N, and P contents along soil profiles in adjacent sapling stands identical in age (10 years) and primary vegetation, but differing in land use type. Land use types included naturally regenerated stands without soil disturbance (NS), moderately disturbed stands (two-year interplanting of soybean after artificial reforestation, MS), and severely disturbed stands (natural reforestation after plantingPanaxginsengfor 10 years, SS). Our objective was to examine the effects of land use change on soil C, N, and P contents and their stoichiometry. Our results showed that:(1) Land use change significantly altered both top-and sub-soil C contents, which descended in the order NS>MS>SS for each soil horizon. However, land use markedly influenced only top-soil N contents and did not significantly influence P content along the entire soil profile (P> 0.05). These results suggest differing sensitivities of elements to land use change. (2) The C/N and C/P ratios among the three land use types were lowest for SS and differed between NS and MS depending on the soil horizon. Vertical patterns of C/N and C/P ratios differed among land use types, but N/P ratio was consistently constrained by land use change. N/P ratios descended in the order NS>MS>SS within each soil horizon and decreased with soil depth (P< 0.05), suggesting that N/P may be a better ecological indicator of response to land use change than C/N and C/P. (3) Land use change significantly affected coupled relationships between C-N, C-P, and N-P. Soil C and N contents were linearly correlated (P< 0.001), regardless of the land use type, with a common slope of 11.1 but significantly different intercepts (P< 0.001). This result, together with the synthesis of published local and global datasets, reflects a large-scale general C-N relationship and a small-scale heterogeneous response of C-N stoichiometry to local land use change.

land use change; biogeochemistry; stoichiometry; carbon-nitrogen coupling

林業公益性行業專項(201104009-05); 國家“十二五”科技支撐項目(2011BAD37B01); 教育部長江學者和創新團隊發展計劃(IRT1054)

2014-03-29; < class="emphasis_bold">網絡出版日期:

日期:2014-12-18

10.5846/stxb201403290589

*通訊作者Corresponding author.E-mail: wangck-cf@nefu.edu.cn

周正虎，王傳寬，張全智.土地利用變化對東北溫帶幼齡林土壤碳氮磷含量及其化學計量特征的影響.生態學報,2015,35(20):6694-6702.

Zhou Z H, Wang C K, Zhang Q Z.The effect of land use change on soil carbon, nitrogen, and phosphorus contents and their stoichiometry in temperate sapling stands in northeastern China.Acta Ecologica Sinica,2015,35(20):6694-6702.