馬原，遲美靜，張玉玲, 范慶峰，虞娜，鄒洪濤

黑土旱地改稻田土壤水穩(wěn)性團聚體有機碳和全氮的變化特征

馬原，遲美靜，張玉玲, 范慶峰，虞娜，鄒洪濤

（沈陽農(nóng)業(yè)大學(xué)土地與環(huán)境學(xué)院/農(nóng)業(yè)農(nóng)村部東北耕地保育重點實驗室，沈陽 110866）

【目的】分析東北黑土旱地改稻田后土壤團聚體組成及其穩(wěn)定性、各粒級團聚體有機碳、全氮含量及其13C、15N自然豐度值的動態(tài)變化，探討旱地改稻田后土壤團聚體有機碳、全氮的賦存能力及穩(wěn)定性，揭示旱地改稻田后土壤團聚體及其有機碳、全氮的演變規(guī)律?！痉椒ā窟x擇東北典型黑土旱地土壤（種植大豆年限大于60年，作為對照）和改種不同年限的稻田土壤（3、5、10、17、20和25年，改稻田前種植作物均為大豆），利用土壤團聚體濕篩分離技術(shù)和穩(wěn)定同位素分析技術(shù)，研究旱地改稻田后土壤團聚體有機碳、全氮的動態(tài)變化特征?！窘Y(jié)果】在0—60 cm土層，與對照土壤相比，改種水稻各年限土壤中2—0.25 mm團聚體組成有所減少，0.25—0.053 mm和＜0.053 mm團聚體組成有所增加，＞2 mm團聚體組成的變化無明顯規(guī)律，但旱地改稻田不同年限均以2—0.053 mm團聚體為主；團聚體平均重量直徑（MWD）與＞2 mm團聚體組成之間呈顯著線性正相關(guān)關(guān)系（＜0.01），與0.25—0.053 mm、＜0.053 mm團聚體組成之間均呈顯著線性負相關(guān)關(guān)系（＜0.01或＜0.05）；水穩(wěn)性團聚體組成變化受水稻種植年限和土層深度的顯著影響，而MWD的變化則受土層深度的顯著影響。與對照土壤相比，在0—40 cm土層，2—0.25 mm、0.25—0.053 mm團聚體有機碳和全氮含量在改種水稻3年時均有所下降，在改種水稻3—25年間均隨水稻種植年限延長大體上呈增加趨勢。總體上，2—0.25 mm、0.25—0.053 mm團聚體是賦存有機碳和全氮的主要粒級；在0—60 cm土層，＞2 mm團聚體有機碳、全氮含量與其團聚體組成之間呈顯著正相關(guān)關(guān)系（＜0.01或＜0.05），在0—20 cm土層，2—0.25 mm團聚體有機碳、全氮含量與其團聚體組成之間也呈顯著正相關(guān)關(guān)系（＜0.01或＜0.05）；＜2 mm團聚體有機碳和全氮含量的變化受水稻種植年限影響顯著，而＞0.25 mm團聚體有機碳和全氮含量的變化則受土層深度影響顯著。與對照土壤相比，各粒級團聚體中δ13C在改種水稻3年時均明顯增加，在改種水稻5年時均明顯下降，在改種水稻5—25年間變化不明顯，各粒級團聚體中δ15N在改種水稻25年間均略有下降?？傮w上，在改稻田3—25年間，團聚體中δ13C、δ15N的變化受水稻種植年限和土層深度的顯著影響，其數(shù)值均隨粒級的減少而增加，相同年限各粒級團聚體δ13C隨著土層的加深而增大，δ15N無明顯變化規(guī)律?！窘Y(jié)論】東北典型黑土旱地改稻田25年間，土壤中非水穩(wěn)性大團聚體遭受破壞形成了粒徑較小的團聚體，2—0.053 mm水穩(wěn)性團聚體是有機碳、全氮固存的主要載體，較小粒級團聚體賦存的有機碳較為穩(wěn)定，其穩(wěn)定性隨水稻種植年限延長、土層加深而增強。

黑土；旱地；稻田；水穩(wěn)性團聚體；有機碳；全氮；13C和15N 自然豐度

0 引言

【研究意義】土壤團聚體是土壤結(jié)構(gòu)的基礎(chǔ)，是由土壤顆粒和各種有機質(zhì)（微生物、植物、動物殘渣及其分泌物）膠結(jié)而成，不同粒級團聚體的分布影響著土壤養(yǎng)分的儲存、土壤生物學(xué)特性和土壤結(jié)構(gòu)[1-2]。因此，土壤團聚體的穩(wěn)定性和分布與土壤有機碳和全氮的固存密切相關(guān)[3]，對于研究土壤結(jié)構(gòu)和有機碳、全氮的固存和轉(zhuǎn)化有重要意義[4]。東北黑土有機質(zhì)含量豐富、結(jié)構(gòu)良好。近年來，隨著農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整，黑土旱地改稻田面積逐年增加，這一土地利用方式的改變及土壤環(huán)境條件的變化，影響了土壤有機碳、全氮的數(shù)量[5]，進而影響土壤中各粒級水穩(wěn)性團聚體的數(shù)量、分配及其穩(wěn)定性[6]。不同粒級團聚體及其有機碳、全氮的賦存數(shù)量的變化可以用來表征或反映不同環(huán)境條件下土壤質(zhì)量的變化[7-8]，因此，研究東北典型黑土旱地改稻田后土壤水穩(wěn)性團聚體及其有機碳、全氮的動態(tài)變化對于揭示黑土有機碳、全氮的演變特征具有重要意義?！厩叭搜芯窟M展】水稻土是長期人為水耕熟化的一類特殊類型土壤，具有較高的固碳潛力。許多研究表明，長期植稻有利于土壤有機碳的固存[9-10]，且與旱地土壤相比有機碳含量較高[5,11-14]，水稻種植過程主要通過物理保護和腐殖化過程來富集和固存土壤有機碳[10]，通過土壤團聚體對有機碳的保護機制可以解釋土壤有機碳的固存和分解[15]。研究發(fā)現(xiàn)，太湖地區(qū)水稻土中主要以2—0.25 mm、0.25—0.02 mm團聚體為主[15]；浙江平原區(qū)植稻17年土壤中2—0.25 mm團聚體是有機碳的主要載體[16]；紅壤荒地開墾種植水稻20年后，水穩(wěn)性大團聚體是有機碳、全氮的主要載體[17]；慈溪和上虞地區(qū)植稻300年間，大團聚體有所增多，微團聚體減少，大團聚體是土壤有機碳的主要載體[8]。但也有研究發(fā)現(xiàn)，慈溪地區(qū)植稻30—2 000年間，0.25—0.053 mm團聚體逐漸減少，＜0.053 mm團聚體增加，土壤有機碳主要分布在0.25—0.053 mm團聚體中[18]。稻田轉(zhuǎn)變?yōu)槭卟说?0年后＞0.25 mm團聚體含量及其固碳能力會下降[19]；本課題組前期研究發(fā)現(xiàn)，東北黑土植稻大于50年的水稻土中主要以2—0.25 mm團聚體為主[20]。土壤13C、15N自然豐度值可以從內(nèi)在反映出土壤有機碳的來源和周轉(zhuǎn)速率[21-24]，植被的轉(zhuǎn)變或者新碳的介入會對土壤有機碳、氮的周轉(zhuǎn)產(chǎn)生影響，從而導(dǎo)致各級團聚體δ13C、δ15N發(fā)生改變[25-27]。由此可見，不同土壤、土地利用方式及種植年限對土壤團聚體的影響不盡相同，多數(shù)研究集中在我國亞熱帶地區(qū)。【本研究切入點】課題組在前期研究中已明確東北黑土旱地改稻田后土壤有機碳、全氮的變化特征[11]，但關(guān)于東北黑土旱田改稻田后土壤各粒級團聚體有機碳、全氮的變化特征如何？還尚不清楚?！緮M解決的關(guān)鍵問題】利用土壤團聚體濕篩分離技術(shù)和穩(wěn)定同位素分析技術(shù)，研究東北黑土區(qū)旱地改稻田后土壤各粒級團聚體組成、團聚體中有機碳、全氮含量及其δ13C、δ15N的變化，闡明黑土區(qū)旱地改稻田后土壤團聚體有機碳、全氮的賦存能力及其穩(wěn)定性，為揭示黑土區(qū)旱地改稻田后土壤團聚體有機碳、全氮的演變規(guī)律提供重要的理論依據(jù)。

1 材料與方法

1.1 供試土壤

供試土壤采自黑龍江省綏化市慶安縣勤勞鎮(zhèn)勤樸村，屬于典型黑土，為黃土性沉積物發(fā)育的土壤。2015 年10月進行實地田間生產(chǎn)調(diào)研，確定連續(xù)種植大豆60年以上的旱地土壤（作為對照土壤）和旱地改稻田年限分別為3、5、10、17、20 和 25年的稻田土壤。供試土壤區(qū)域的氣候和地形（漫崗丘陵區(qū)，坡度小于3°）大致相同，旱地改種水稻前種植作物均為大豆，改種水稻前土壤性狀與對照土壤性狀基本一致。旱地改種水稻后每年施用氮、磷、鉀化肥量大致相同，但年限間化學(xué)肥料的品種、用量存在差異，供試土壤不施用有機肥料。

2015年10月水稻收獲后進行供試田塊土壤樣本采集。將每個相同年限的田塊作為一個采樣區(qū)域，并劃分為面積大約相同的3個采樣單元作為3次采樣重復(fù)，每個采樣單元以“S”法布設(shè)5—7點，用土鉆采集布設(shè)點樣本后均勻混合作為一個重復(fù)樣本，混合后土壤樣本裝在塑料箱中（避免破壞原土狀態(tài)）?？紤]黑土自身腐質(zhì)層比較深厚，同時旱地改稻田年限相對較短，稻田土壤犁底層形成相對較弱，因此確定采樣深度為0—20、20—40和40—60 cm。土樣運回室內(nèi)后，在室溫下自然風(fēng)干，在風(fēng)干過程中用手輕輕地把土塊沿自然結(jié)構(gòu)面掰成直徑約5 cm的小塊，并除去肉眼可以識別的粗根和石塊，風(fēng)干后混合均勻，其中，一部分土樣磨碎過0.15 mm篩備用，一部分土樣過8 mm篩用于制備各粒級團聚體。供試土壤的地理信息及理化性質(zhì)見表1。

1.2 各級團聚體的分離

土壤團聚體的分離依據(jù)SIX等[28]提供的土壤團聚體濕篩法，采用團粒分析儀（DIK-2012，日本）進行篩分。稱取50 g自然風(fēng)干土樣，放在2 mm、0.25 mm和0.053 mm組成的套篩上，調(diào)整水面高度使水面剛好沒過土樣，浸泡5 min后，以速度為30 次/min振蕩30 min，篩分結(jié)束后分別將套篩上的團聚體進行收集并轉(zhuǎn)移至蒸發(fā)皿中，并于60 ℃烘干至恒重（約48 h）并稱重，獲得＞2 mm、2—0.25 mm、0.25—0.053 mm和＜0.053 mm團聚體。并將制備好各粒級團聚體樣品磨碎過0.15 mm篩備用。

1.3 有機碳、全氮含量及其δ13C、δ15N自然豐度值測定

原土和各粒級團聚體樣品中的有機碳、全氮含量及其δ13C、δ15N均采用元素分析儀—穩(wěn)定同位素比例質(zhì)譜儀（EA–IRMS，Elementanalysis-Stable100 Isotope Ratio Mass Spectrometer，德國）聯(lián)用測定。

式中，R、R分別為樣品13C/12C 原子比值和15N/14N原子比值；R和R值為0.0112372（以美國南卡羅來納洲白堊紀皮狄組層中的擬箭石化石（Pee Dee Belemnite, PDB）為標準物質(zhì)）和0.0036765（以純凈大氣氮（Air Nitrogen, AN）為標準物質(zhì)）。

1.4 數(shù)據(jù)計算

（1）團聚體平均重量直徑（MWD）

（2）某粒級團聚體有機碳、全氮含量

aggregate（aggregate）=con-aggregate）con-aggregate）×aggregate（4）

式中，aggregate（aggregate）為某粒級團聚體有機碳（全氮）含量（g·kg-1soil），con-aggregate（con-aggregate）為某粒組團聚體有機碳（全氮）測試濃度（g·kg-1aggregate）。

1.5 數(shù)據(jù)處理

采用 Excel 2010 和SPSS 19.0 （SPSS Inc., Chicago, IL, USA）進行數(shù)據(jù)處理，采用 Origin9.0 軟件對數(shù)據(jù)進行繪圖。以年限和土層為因素進行雙因素方差分析（two-ways ANOVA），采用Tukey法進行多重比較，若數(shù)據(jù)誤差方差為非齊性，則使用非參數(shù)檢驗（non-parametric test），Kruskal-Wallis單因素方差分析（one-way ANOVA）。若兩個因素存在顯著的交互作用，則進行單因素方差分析（one-way ANOVA），采用Tukey法或Dunnett T3法進行多重比較，并以小寫字母表示差異顯著性（＜0.05）。數(shù)據(jù)均為3次重復(fù)的平均值。

表1 供試土壤的地理信息及理化性質(zhì)

2 結(jié)果

2.1 旱地改稻田土壤水穩(wěn)性團聚體組成及平均重量直徑的變化

由圖1可以看出，旱地改稻田3—25年間，在0—20 cm土層，＞2 mm、2—0.25 mm、0.25—0.053 mm和＜0.053 mm團聚體組成分別為8.96%—27.11%、27.68%—38.58%、22.29%—37.12%和14.65%—26.05%；在20—40 cm土層，＞2 mm、2—0.25 mm、0.25—0.053 mm和＜0.053 mm團聚體組成分別為2.87%— 24.37%、24.45%—34.40%、27.38%—38.12%和18.25% —32.32%；在40—60 cm土層，＞2 mm、2—0.25 mm、0.25—0.053 mm和＜0.053 mm團聚體組成分別為1.84%—7.54%、25.69%—34.41%、33.35%—38.91%和23.00%—33.42%。總體上，在0—60 cm土層，與種植大豆的旱地土壤相比，改種水稻各年限土壤中，2—0.25 mm團聚體組成有所減少，0.25—0.053 mm和＜0.053mm團聚體組成有所增加，＞2 mm團聚體組成的變化無明顯規(guī)律。另外，水稻種植年限和土層深度之間的交互作用對各粒級團聚體組成均無顯著影響，但水稻種植年限對2—0.25 mm、0.25— 0.053 mm和＜0.053 mm團聚體組成影響顯著，土層深度則對＞2 mm、0.25—0.053 mm和＜0.053 mm團聚體組成影響顯著。由此可見，旱地改稻田，由于每年在水稻種植期間淹水，使土壤中非水穩(wěn)性大團聚體遭受破壞，形成了粒徑較小的團聚體，但均以2—0.25 mm和0.25—0.053 mm團聚體為主；水稻種植年限和土層深度是影響水穩(wěn)性團聚體組成變化的重要因素。

數(shù)值為平均值±標準誤。* ，**，***，ns分別表示為0.05，0.01，0.001水平顯著和不顯著。下同

由圖2所示，土層深度是顯著影響土壤團聚體平均重量直徑（MWD）變化的重要因素；總體上，MWD隨土層深度的加深而減少，0—20 cm土層的MWD顯著（＜0.01或＜0.05）大于20—40 cm和40—60 cm土層的MWD（＜0.05）。相關(guān)分析顯示（表2），在3個土層中，MWD與＞2 mm團聚體組成之間呈顯著線性正相關(guān)關(guān)系（＜0.01），與0.25—0.053 mm、＜0.053 mm團聚體組成之間均呈顯著線性負相關(guān)關(guān)系（＜0.01或＜0.05）。這表明旱地改稻田后，＞2 mm團聚體數(shù)量的增加提高了土壤水穩(wěn)性團聚體的穩(wěn)定性，而較小粒級團聚體（＜0.25 mm）的增加則降低了土壤水穩(wěn)性團聚體的穩(wěn)定性。

***為0.001水平顯著，ns為不顯著

2.2 旱地改稻田土壤水穩(wěn)性團聚體有機碳（氮）含量的變化

由表3可以看出，總體上，2—0.25 mm、0.25— 0.053 mm團聚體是賦存有機碳和全氮的主要粒級。與對照土壤相比，在0—20 cm和20—40 cm土層，2—0.25 mm、0.25—0.053 mm團聚體有機碳和全氮含量在改種水稻3年時均有所下降；在改種水稻3—25年間均隨水稻種植年限延長大體上呈增加趨勢，且在改種水稻25年時，兩個粒級團聚體有機碳和全氮含量與對照土壤的相接近。在40—60 cm土層，2—0.25 mm、0.25—0.053 mm團聚體有機碳和全氮含量在改種水稻3年時均有所下降；兩個粒級團聚體有機碳含量在改種水稻3—17年間大體呈逐年增加趨勢，而在改種水稻17—25年間則呈逐年降低趨勢；2—0.25 mm團聚體全氮含量在改種水稻3—25年間大體呈逐年增加趨勢，0.25—0.053 mm團聚體全氮含量在改種水稻3—17年間有所增加，在改種水稻17—25年間有所下降。＞2 mm和＜0.053 mm團聚體有機碳和全氮含量在改種水稻后無顯明規(guī)律性。另外，水稻種植年限和土層深度之間的交互作用對各粒級團聚體有機碳和全氮含量均無顯著影響，但水稻種植年限和土層深度對不同粒級團聚體有機碳和全氮含量有重要影響。總體上，水稻種植年限對2—0.25 mm、0.25—0.053 mm和＜0.053 mm團聚體有機碳和全氮含量影響顯著，土層深度對＞2 mm和2—0.25 mm團聚體有機碳和全氮含量影響顯著。

表2 各粒級團聚體組成與土壤團聚體平均重量直徑的相關(guān)系數(shù)

*，**分別表示為0.05，0.01水平顯著。下同

With＜ 0.05(*),＜ 0.01(**)The same as below

表3 旱地改稻田后不同年限土壤團聚體有機碳和全氮含量

相關(guān)分析顯示（表4），在3個土層，＞2 mm團聚體有機碳、全氮含量與其團聚體組成之間均呈顯著正相關(guān)關(guān)系（＜0.01或＜0.05）。0—20 cm土層，2—0.25 mm團聚體有機碳、全氮含量與其團聚體組成之間也呈顯著正相關(guān)關(guān)系（＜0.01或＜0.05）。這說明旱地改種水稻田后，有機碳、全氮對于土壤水穩(wěn)性大團聚體（＞0.25 mm），尤其是表層水穩(wěn)性大團聚體（＞2 mm）形成具有重要作用，同時大團聚體的形成與穩(wěn)定也有利于有機碳、全氮的賦存。

表4 各粒級團聚體有機碳、全氮含量與各粒級土壤團聚體的相關(guān)系數(shù)

2.3 旱地改稻田土壤水穩(wěn)性團聚體中δ13C、δ15N的變化

由表5可見，與對照土壤相比，各粒級團聚體中δ13C在改種水稻3年時均明顯增加，在改種水稻3—5年間明顯下降，在改種水稻5—25年間變化不明顯；在改種水稻25年間，各粒級團聚體中δ15N值均略有下降?？傮w上，在改種水稻3—25年間，團聚體中δ13C、δ15N隨粒級的減小呈增加趨勢，相同年限各粒級團聚體δ13C隨著土層的加深而增大，δ15N的變化無明顯規(guī)律。水稻種植年限顯著影響4個粒級團聚體中δ13C和δ15N，土層深度顯著影響＜2 mm粒級團聚體中δ13C和δ15N，二者的交互作用僅顯著影響0.25—0.053 mm團聚體中δ13C和δ15N。這說明隨著團聚體的粒級變小，土壤有機碳（氮）易分解性越小，相對越穩(wěn)定；其穩(wěn)定性隨水稻種植年限延長、土層加深度而增強。

3 討論

3.1 旱地改稻田對土壤團聚體組成及其穩(wěn)定性的影響

土地利用方式和耕作管理（作物類型、施肥、排水和灌溉）是影響土壤結(jié)構(gòu)和特性的最重要因素，并對土壤團聚體的形成和破裂產(chǎn)生影響[29]。Huang等[14]對東北旱地和稻田土壤的研究發(fā)現(xiàn)，土壤團聚體主要以2—0.25 mm、0.25—0.053 mm粒級為主，旱地改種稻田27年時稻田土壤中＞0.25 mm團聚體顯著高于旱地土壤。Zou等[9]研究表明，隨著種稻年限的延長大團聚體隨之增加，微團聚體隨之減少。李昌新等[11]研究表明，長期玉米和水稻種植下，稻田土壤中＞2 mm、＜0.053 mm團聚體均顯著高于旱地，而2—0.25 mm團聚體則相反。Zheng等[8]和Pan等[15]研究表明，耕作會引起2—0.25 mm團聚體發(fā)生改變。本研究中，旱地改稻田3年時，2—0.25 mm團聚體下降，＜0.25 mm團聚體組成上升；旱地改稻田3—25年間，2—0.25 mm團聚體有所上升，＜0.25 mm團聚體有所下降；旱地改稻田25年時2—0.25 mm團聚體低于旱地土壤，＜0.25 mm團聚體高于旱地土壤（圖1），這說明黑土在旱地改稻田年限較短時，干濕交替過程[30]和淹水條件[31]會導(dǎo)致大團聚體破裂形成微團聚體，其中2—0.25 mm團聚體對于耕作環(huán)境的變化比較敏感，但隨著水稻種植年限的延長，由于長期淹水會使微團聚體進一步團聚，進而促進大團聚體的形成[15]。

表5 旱地改稻田后不同年限土壤團聚體中的δ13C和δ15N

不同小寫字母表示相同粒級團聚體，年限與土層之間存在顯著交互作用達0.05顯著水平

Indicate the significant (＜0.05) differences between years and soil depth by different lowercase letters

土壤團聚體平均重量直徑（MWD）可以反映出不同粒級團聚體的分布狀況和穩(wěn)定性[1]，MWD值越大，團聚體穩(wěn)定性越好[32]。Wang等[19]研究表明，大團聚體在干濕交替過程中易于破裂為微團聚體，并導(dǎo)致團聚體水穩(wěn)性顯著下降。本研究中，MWD值與＞2 mm團聚體組成呈顯著正相關(guān)，與＜0.25 mm團聚體組成呈顯著負相關(guān)（表2）；水稻種植年限對MWD影響不顯著，而土層深度則對MWD影響顯著，0—20 cm土層團聚體穩(wěn)定性顯著高于20—60 cm土層（圖2）。這說明旱地改稻田后，由于干濕交替和淹水過程，影響了0—60 cm土層各粒級團聚體的變化，水稻種植年限為17年時，由于0—40 cm土層中＞2 mm和2—0.25 mm大團聚體的增加，0.25—0.053 mm和＜0.053 mm團聚體的減少（圖1），使0—40 cm土層團聚體水穩(wěn)性有所上升（圖2）。

3.2 旱地改稻田對土壤團聚體有機碳、全氮含量及其穩(wěn)定性的影響

土壤團聚體是土壤有機碳、全氮固定[33]和物質(zhì)能量轉(zhuǎn)化[17]的重要場所，不僅可以保護土壤有機質(zhì)也可以延遲其礦化[34]，由于各粒級團聚體對土壤有機碳、全氮的保護機制不同，因此對土壤有機碳、全氮的固持也存在差異[33]。大團聚體主要是通過根系、真菌菌絲等膠結(jié)而成。在大多數(shù)情況下，有機物先與其結(jié)合，碳濃度最高且周轉(zhuǎn)率最快[2]，而微團聚體受芳香族化合物束縛，碳周轉(zhuǎn)速率較低[35]。李昌新等[10]對種植20多年的南方稻田和旱地土壤的研究表明，稻田和旱地土壤有機碳含量均主要分配在2—0.25 mm和0.25—0.053 mm團聚體中，與旱地土壤相比，稻田土壤降低了2—0.25 mm團聚體有機碳的比例，但提高了其他粒級團聚體有機碳的比例。HUANG等[14]研究表明，旱地轉(zhuǎn)變?yōu)榈咎?7年后顯著提高了2—0.25 mm、0.25—0.053 mm團聚體有機碳含量，但其他粒級團聚體有機碳含量則無顯著變化。徐文靜等[20]對東北黑土區(qū)種稻大于50年的水稻土研究表明，2—0.25 mm團聚體是賦存有機碳的主要載體，具有明顯的固碳能力。許多研究表明，稻田土壤中有機碳含量的變化與大團聚體緊密相關(guān)[10,15,36]，各粒級團聚體有機碳含量隨種稻年限的延長而增加[9]，但也有研究表明，稻田土壤有機碳可向微團聚體和粉黏粒中轉(zhuǎn)移以利于有機碳的長期賦存[11,14,18]。本研究中，旱地改種稻田后各粒級團聚體中有機碳、全氮含量的變化大致相同，土壤有機碳、全氮含量與土壤團聚體組成之間緊密相關(guān)（表4），在改種水稻年限較短時（3年間），各粒級團聚體中有機碳、全氮含量下降（表3），這主要是由于土地利用方式的改變使土壤中各粒級團聚體遭受不同程度破壞（圖1），進而加速了有機碳（氮）分解和周轉(zhuǎn)速率[37]；在改稻田3—25年間，2—0.25 mm、0.25—0.05 mm團聚體有機碳、全氮含量則有所增加（圖3），這主要是由于隨著改種年限延長，淹水條件可延緩?fù)寥烙袡C質(zhì)的分解速率，降低土壤有機碳、氮的礦化速率[38]，同時由于大團聚體破裂為微團聚體，使其微團聚體易與有機質(zhì)粘結(jié)[39]。由此可見，東北黑土旱地改稻田25間，2—0.25 mm、0.25—0.053 mm團聚體是土壤有機碳、全氮賦存的主要粒級，對有機碳（氮）的物理保護發(fā)揮重要作用。

采用δ13C可以研究土壤有機碳的動態(tài)變化，土壤的δ13C隨時間的變化主要受植被輸入的有機碳控制，植物殘體進入土壤后在微生物作用下可發(fā)生礦化分解過程和腐殖化過程，致使土壤的δ13C會與植物本身的δ13C有所不同[24]；其次由生物衰變過程控制[22]，土壤有機碳分解過程中微生物對12C的優(yōu)先分解會引起土壤δ13C增高，δ13C越大，其有機碳的降解程度越低[25]。大豆和水稻均屬于C3作物，其δ13C介于-29.1‰—-24.2‰[21]，水稻的平均δ13C為-26.34‰[4]。本研究中，長期連續(xù)種植大豆的旱地，在改種水稻年限較短時（3年），由于土地利用方式及植被類型的改變，新鮮有機質(zhì)會不同程度進入到不同粒級團聚體中，但又由于水稻單產(chǎn)顯著高于大豆單產(chǎn)，對土壤養(yǎng)分的需求大幅度增加，同時大團聚體遭受破壞，加速土壤中易分解有機碳的分解使碳同位素分餾效應(yīng)增加，因此，致使旱地改稻田3年間各粒級團聚體中δ13C均有所增加（表5），這與課題組前期關(guān)于旱地改稻田土壤δ13C的變化結(jié)果相一致[5]；在旱地改稻田5—25年，由于大量水稻根系和地上殘體歸還土壤使有機碳得到補充和更新，又由于稻田土壤環(huán)境逐漸趨于穩(wěn)定可使有機碳分解緩慢，進而使各級團聚體中δ13C均有所降低，但隨種稻年限延長變化較為平緩（表5）。LIU等[27]研究發(fā)現(xiàn)，新鮮有機碳的賦存一般從大團聚體（＞0.25 mm）開始，經(jīng)過分解過程和微生物消耗后，降解的有機碳則被螯合在微團聚體（＜0.25 mm）中，植物殘留物的腐殖化作用會伴隨著13C的輕微富集。本研究中，團聚體的粒徑越小，其δ13C越大（表5），這也說明團聚體粒徑越大，賦存新鮮有機碳越多，這一結(jié)果與LIU等[27]研究結(jié)果相類似。另外，相同年限各粒級團聚體中δ13C隨土層加深而增大（表5），這進一步說明，外源新鮮有機碳不能及時地進入補充到深層土壤，同時深層土壤有機碳相對較“老”[23,27]，礦化分解速度較慢。

土壤中δ15N的變異可反映出土壤氮素的轉(zhuǎn)化和遷移特征，其中礦化和硝化過程是導(dǎo)致氮同位素分餾從而影響δ15N的主要原因[40]。研究表明，土壤中δ15N在很大程度受土壤水分的影響[40-41]，非淹水條件下土壤中硝化作用較快，而在淹水條件下有利于反硝化作用，礦化作用相對較弱，微生物對氮的固定作用較快。本研究中，旱地改稻田后，由于稻田土壤多處于淹水或含水量較高的狀態(tài)，土壤硝化作用受到抑制，大量15N貧化的水稻植株根系或殘體進入土壤[26,35]，致使稻田土壤各粒級團聚體的δ15N低于旱地土壤（表5）。

4 結(jié)論

東北黑土旱地改種稻田25年間，稻田土壤中2—0.25 mm團聚體有所減少，0.25—0.053 mm和＜0.053 mm團聚體有所增加，＞2 mm團聚體則無明顯變化規(guī)律，但均以2—0.053 mm團聚體為主；各粒級團聚體有機碳、全氮含量隨改種年限延長呈大致相同的變化趨勢，總體上，2—0.053 mm團聚體是賦存有機碳和全氮的主要粒級；團聚體中δ13C、δ15N均隨粒級的減小呈增加趨勢，相同年限各粒級團聚體δ13C隨著土層的加深而增大，δ15N無明顯變化規(guī)律。因此，東北典型黑土旱地改稻田后，土壤中非水穩(wěn)性大團聚體遭受破壞形成了粒徑較小的團聚體，2—0.053 mm水穩(wěn)性團聚體是有機碳、全氮固存的主要載體，較小粒級團聚體賦存的有機碳較為穩(wěn)定，其穩(wěn)定性隨水稻種植年限延長、土層加深而增大。

[1] WANG Y, ZHANG J H, ZHANG Z H. Influences of intensive tillage on water-stable aggregate distribution on a steep hillslope., 2015, 151: 82-92.

[2] HUANG R, LAN M L, LIU J, GAO M. Soil aggregate and organic carbon distribution at dry land soil and paddy soil: the role of different straws returning., 2017, 24(36): 1-11.

[3] SIX J, BOSSUYT H, DEGRYZE S, DENEF K. A history of research on the link between (Micro) aggregates, soil biota, and soil organic matter dynamics., 2004, 79(1): 7-31.

[4] LIU Y, HU C, HU W, WANG L, LI Z G, PAN J F, CHEN F. Stable isotope fractionation provides information on carbon dynamics in soil aggregates subjected to different long-term fertilization practices., 2018, 177: 54-60.

[5] 賈樹海, 張佳楠, 張玉玲. 東北黑土區(qū)旱田改稻田后土壤有機碳、全氮的變化特征, 中國農(nóng)業(yè)科學(xué), 2017, 50(7): 1252-1262.

JIA S H, ZHANG J N, ZHANG Y L. Changes of the characteristics of soil organic carbon and total nitrogen after conversation from upland to paddy field in black soil region of Northeast China., 2017, 50(7): 1252-1262. (in Chinese)

[6] ZHAO J S, CHEN S, HU R G, LI Y Y. Aggregate stability and size distribution of red soils under different land uses integrally regulated by soil organic matter, and iron and aluminum oxides., 2017, 167: 73-79.

[7] 徐香茹, 汪景寬. 土壤團聚體與有機碳穩(wěn)定機制的研究進展. 土壤通報, 2017, 48(6): 1523-1529.

XU X R, WANG J K. A review on different stabilized mechanisms of soil aggregates and organic carbon.2017, 48(6): 1523- 1529. (in Chinese)

[8] ZHENG H B, LIU W R, ZHENG J Y, LUO Y, LI R P, WANG H, QI H. Effect of long-term tillage on soil aggregates and aggregate-associated carbon in black soil of Northeast China,, 2018, 13(6): e199523.

[9] ZOU P, FU J R, CAO Z H, YE J, YU Q G. Aggregate dynamics and associated soil organic matter in topsoils of two 2,000-year paddy soil chronosequences., 2015, 15(3): 510-522.

[10] WANG P, LIU Y L, LI L Q, CHENG K, ZHENG J F, ZHANG X H, ZHENG J W, JOSEPH S, PAN G X. Long-term rice cultivation stabilizes soil organic carbon and promotes soil microbial activity in a salt marsh derived soil chronosequence., 2015, 5: 15704.

[11] 李昌新, 黃山, 彭現(xiàn)憲, 黃欠如, 張衛(wèi)建. 南方紅壤稻田與旱地土壤有機碳及其組分的特征差異. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2009, 28(3): 606-611.

LI C X, HUANG S, PENG X X, HUANG Q R, ZHANG W J. Differences in soil organic carbon fractions between paddy field and upland field in red soil region of South China.2009, 28(3): 606-611. (in Chinese)

[12] YAN X, ZHOU H, ZHU Q H, WANG X F, ZHANG Y Z, YU X C, PENG X. Carbon sequestration efficiency in paddy soil and upland soil under long-term fertilization in Southern China., 2013, 130: 42-51.

[13] SUN Y N, HUANG S, YU X C, ZHANG W J. Differences in fertilization impacts on organic carbon content and stability in a paddy and an upland soil in Subtropical China., 2015, 397(1): 189-200.

[14] HUANG S, PAN X H, GUO J, QIAN C R, ZHANG W J. Differences in soil organic carbon stocks and fraction distributions between rice paddies and upland cropping systems in China., 2014, 14(1): 89-98.

[15] PAN G X, WU L S, LI L Q, ZHANG X H, GONG W, WOOD Y. Organic carbon stratification and size distribution of three typical paddy soils from Taihu lake region, China., 2008, 20(4): 456-463.

[16] 毛霞麗, 陸扣萍, 何麗芝, 宋照亮, 徐祖祥, 楊文葉, 徐進, 王海龍. 長期施肥對浙江稻田土壤團聚體及其有機碳分布的影響. 土壤學(xué)報, 2015(4): 828-838.

MAO X L, LU K P, HE L Z, SONG Z L, XU Z X, YANG W Y, XU J, WANG H L. Effect of long-term fertilizer application on distribution of aggregate-associated organic carbon in paddy doil., 2015(4): 828-838. (in Chinese)

[17] 陳曉芬, 李忠佩, 劉明, 江春玉. 不同施肥處理對紅壤水稻土團聚體有機碳、氮分布和微生物生物量的影響. 中國農(nóng)業(yè)科學(xué), 2013, 46(5): 950-960.

CHEN X F, LI Z P, LIU M, JIANG C Y. Effects of different fertilizations on organic carbon and nitrogen contents in water-stable aggregates and microbial biomass content in paddy soil of Subtropical China., 2013, 46(5): 950-960. (in Chinese)

[18] 王欣欣, 符建榮, 鄒平, 陳維, 葉靜, 俞巧鋼, 姜麗娜, 王強. 長期植稻年限序列水稻土團聚體有機碳分布特征. 應(yīng)用生態(tài)學(xué)報, 2013, 24(3): 719-724.

WANG X X, FU J R, ZOU P, CHEN W, YE J, YU Q G, JIANG L N, WANG Q. Distribution characteristics of aggregates organic carbon in a paddy soil chronosequence., 2013, 24(3): 719-724. (in Chinese)

[19] WANG H, GUAN D S, ZHANG R D, CHEN Y J, HU Y T, XIAO L. Soil aggregates and organic carbon affected by the land use change from rice paddy to vegetable field., 2014, 70: 206-211.

[20] 徐文靜, 叢耀輝, 張玉玲, 段鵬鵬, 范慶鋒, 張玉龍. 黑土區(qū)水稻土水穩(wěn)性團聚體有機碳及其顆粒有機碳的分布特征. 水土保持學(xué)報, 2016, 30(4): 210-215.

XU W J, CONG Y H, ZHANG Y L, DUAN P P, FAN Q F, ZHANG Y L. Distribution of organic carbon and partivulate organic carbon in water-stable aggregates of paddy soil in black soil area., 2016, 30(4): 210-215. (in Chinese)

[21] 竇森, 張晉京. 用δ13C值研究土壤有機質(zhì)周轉(zhuǎn)的方法及其評價. 吉林農(nóng)業(yè)大學(xué)學(xué)報, 2001, 23(2): 64-67.

DOU S, ZHANG J J. Introduction of a method for studying turnover of soil organic matter., 2001, 23(2): 64-67. (in Chinese)

[22] BAI E, BOUTTON T W, LIU F, WU X B, HALLMARK C T, ARCHER S R. Spatial variation of soil δ13C and its relation to carbon input and soil texture in a subtropical lowland woodland., 2012, 44(1): 102-112.

[23] PERI P L, LADD B, PEPPER D A, BONSER S P, LAFFAN S W, AMELUNG W. Carbon (δ13C) and nitrogen (δ15N) stable isotope composition in plant and soil in Southern Patagonia's Native Forests., 2012, 18(1): 311-321.

[24] KULSAWAT W, PORNTEPKASEMSAN B, NOCHIT P. Paddy soil profile dstribution of δ13C subjected to rice straw amendment and burning., 2019, 886: 3-7.

[25] GUNINA A, KUZYAKOV Y. Pathways of Litter C by Formation of aggregates and SOM density fractions: implications from13C natural abundance., 2014, 71: 95-104.

[26] LIM S S, KWAK J H, LEE K S, CHANG S X, YOON K S, KIM H Y, CHOI W J. Soil and plant nitrogen pools in paddy and upland ecosystems have contrasting δ15N., 2015, 51(2): 231-239.

[27] LIU Y, LIU W Z, WU L H, LIU C, WANG L, CHEN F, LI Z G. Soil aggregate-associated organic carbon dynamics subjected to different types of land use: evidence from13C natural abundance., 2018, 122: 295-302.

[28] SIX J, ELLIOTT E T, PAUSTIAN K, DORAN J W. Aggregation and soil organic matter accumulation in cultivated and native grassland soils., 1998, 62(5): 1367-1377.

[29] ZHAO J S, CHEN S, HU R G, LI Y L. Aggregate stability and size distribution of red soils under different land uses integrally regulated by soil organic matter, and iron and aluminum oxides., 2017, 167: 73-79.

[30] VAN VEEN A J, KUIKMAN P J. Soil structural aspects of decomposition of organic matter by micro-organisms., 1990, 11(3): 213-223.

[31] JIANG X J, XIE D T. Combining ridge with no-tillage in lowland rice-based cropping system: long-term effect on soil and rice yield., 2009, 19(4): 515-522.

[32] 郝翔翔, 楊春葆, 苑亞茹. 連續(xù)秸稈還田對黑土團聚體中有機碳含量及土壤肥力的影響. 中國農(nóng)學(xué)通報, 2013, 29(35): 263-269.

HAO X X, YANG C B, YUAN Y R. Effects of continuous straw returning on organic carbon content in aggregates and fertility of black soil., 2013, 29(35): 263-269. (in Chinese)

[33] 徐虎, 張敬業(yè), 蔡岸冬, 王小利, 張文菊. 外源有機物料碳氮在紅壤團聚體中的殘留特征. 中國農(nóng)業(yè)科學(xué), 2015, 48(23): 4660-4668.

Xu H, ZHANG J Y, CAI A D，WANG X L, ZHANG W J. Retention characteristic of carbon and nitrogen from amendments in different size aggregates of red soil., 2015, 48(23): 4660-4668. (in Chinese)

[34] RABBI S M F, WILSON B R, LOCKWOOD P V, DANIEL H, YOUNG I M. Aggregate hierarchy and carbon mineralization in two oxisols of new south wales, Australia., 2015, 146: 193-203.

[35] 高崇升, 王建國. 黑土農(nóng)田土壤有機碳演變研究進展. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2011, 19(6): 1468-1474.

Gao C S, Wang J G. A review of researches on evolution of soil organic carbon in mollisols farmland.2011, 19(6): 1468-1474. (in Chinese)

[36] CHEN Z D, TI J S, CHEN F. Soil aggregates response to tillage and residue management in a double paddy rice soil of the Southern China., 2017, 109(9): 1-12.

[37] WANG B S, GAO L L, YU W S, WEI X Q, LI J, LI S P, SONG X J, LIANG G P, CAI D X, WU X P. Distribution of soil aggregates and organic carbon in deep soil under long-term conservation tillage with residual retention in dryland., 2019, 11(2): 241-254.

[38] K?GEL-KNABNER I, AMELUNG W, CAO Z H, FIEDLER S, FRENZEL P, JAHN R, KALBITZ K, K?IBI A, SCHLOTER M. Biogeochemistry of paddy soils., 2010, 157(1): 1-14.

[39] QIAN J, LIU J J, WANG P F, WANG C, HU J, LI K, LU B, TIAN X, GUAN W Y. Effects of riparian land use changes on soil aggregates and organic carbon., 2018, 112: 82-88.

[40] 慈恩, 楊林章, 倪九派, 高明, 謝德體. 不同區(qū)域水稻土的氮素分配及δ15N特征. 水土保持學(xué)報, 2009, 23(2): 103-108.

CI E, YANG L Z, NI J P, GAO M, XIE D T. Distribution and δ15N character of nitrogen in paddy soils located in different regions., 2009, 23(2): 103-108. (in Chinese)

[41] CHOI W J, RO H M. Differences in isotopic fractionation of nitrogen in water-saturated and unsaturated soils., 2003, 35(3): 483-486.

Change Characteristics of Organic Carbon and Total Nitrogen in Water-stable Aggregate After Conversion from Upland to Paddy Field in Black Soil

MA Yuan, CHI MeiJing, ZHANG YuLing, FAN QingFeng, YU Na, ZOU HongTao

(College of Land and Environment, Shenyang Agricultural University/Key Laboratory of Northeast Arable Land Conservation, Ministry of Agriculture and Rural Affairs, Shenyang 110866)

【Objective】 The objectives of this study wereto analyze the composition and stability of soil aggregate, the changes of organic carbon (OC), total nitrogen (TN) content,natural abundance of13C and 15N in different-sized aggregates, to explore the sequestration and stability of soil aggregate organic C, TN, and to reveal the evolution of soil aggregate organic C, TN changes after the conversion from upland to paddy field in black soil region of Northeast China.【Method】 Soil samples were collected from upland (soybean planted for over 60 years) in typical black soil and paddy soil with different years (3, 5, 10, 17, 20 and 25 years, soybean was planted in all the fields before conversion to paddy field). The dynamic characteristics of OC and TN in soil aggregates were studied by using wet-sieving method and stable isotope analysis technology.【Result】In the 0-60 cm soil layers, compared with the control treatment, the composition of 2-0.25 mm aggregates in the soil of different years after rice planting was decreased, which of 0.25-0.053 mm and ＜0.053 mm aggregates was increased. There was no obvious change in the composition of ＞2 mm aggregates, but the different years of dry land change to paddy fields were dominated by 2-0.053 mm aggregates; the mean weight diameter (MWD) of aggregates was significantly positive correlated with the proportion of ＞2 mm aggregates (＜0.01), and significantly negative correlated with the proportion of 0.25-0.053 mm and ＜0.053 mm aggregates (＜0.01 or＜0.05). The change of aggregate composition was significantly affected by different rice planting years and soil depth, whereas the MWD was significantly affected by soil depth. Compared with the control soil, in the 0-40 cm soil layer, the OC and TN contents in the size of 2-0.25 mm and 0.25-0.053 mm aggregates were declined in the 3 years, however there showed increased trend with the extension of rice cultivation in 3-25 years. Generally, OC and TN were mainly accumulated in the 2-0.25 mm and 0.25-0.053 mm aggregates. There existed significant positive correlation between OC and TN contents and aggregate composition in ＞ 2 mm aggregates (＜0.01 or＜0.05) in the 0-60 cm soil layers, as well as 2-0.25 mm aggregates in the 0-20cm soil layer (＜0.01 or＜0.05). The OC and TN contents variation in ＜2 mm aggregates were significantly affected by rice cultivation time, while soil depth significantly affected ＞0.25 mm aggregates OC and TN contents. Compared with the control soil, the δ13C in each size of aggregates significantly increased in 3 rice planting years and decreased in 5 rice planting years, respectively, while there was no significant change in the 5-25 rice planting years, and the δ15N in all size of aggregates decreased slightly during the 25 years of rice replanting. In general, the δ13C and δ15N of soil in aggregates were significantly affected by rice cultivation time and soil depth, which increased with the decreasing of aggregate size. The δ13C increased with soil depth in the same year, while δ15N had no significant change.【Conclusion】After the conversion from dry land to paddy field for 25 years, non-water-stable macro-aggregates in the soil were damaged and formed into small sized aggregates. The 2-0.053 mm water-stable aggregates were the main carrier of OC and TN sequestration, while OC in small size aggregates more stable, and its stability was increased by the rice cultivation time and soil depth increased.

black soil; upland; paddy field; water-stable aggregate; organic carbon; total nitrogen;13C and15N natural abundance

2019-07-12；

2019-11-06

國家自然科學(xué)基金面上項目（41571280）、國家自然科學(xué)基金青年基金項目（41101276）

馬原，E-mail：mayuan275@163.com。通信作者張玉玲，E-mail：zhangyuling@syau.edu.cn

（責任編輯 李云霞）