不同林型興安落葉松林土壤團聚體及其有機碳特征

王 冰，張鵬杰，張秋良

(內(nèi)蒙古農(nóng)業(yè)大學(xué)林學(xué)院，內(nèi)蒙古 呼和浩特 010019)

土壤團聚體是土壤結(jié)構(gòu)的基本單元[1]，其對有機碳的保護作用是穩(wěn)定土壤碳庫的重要機制[2]。由于不同粒徑土壤團聚體對有機碳的保護機制不同，有機碳在不同粒徑土壤團聚體的穩(wěn)定程度也有所差異。一般而言，粉黏粒通過化學(xué)或物理化學(xué)機制與有機碳緊密結(jié)合，比大團聚體和微團聚體更穩(wěn)定；微團聚體對有機碳的保護程度大于大團聚體[3]。土壤中的膠結(jié)物質(zhì)是土壤團聚體形成的核心和基礎(chǔ)[4]，不同生態(tài)系統(tǒng)下土壤膠結(jié)物質(zhì)的變化驅(qū)動著土壤團聚體組成發(fā)生分異，進而影響土壤有機碳在不同粒徑團聚體中的分配[5]。因此，土地利用方式[6-8]和林分類型[1,9-10]均會對土壤團聚體及其有機碳的分布產(chǎn)生影響。土壤膠結(jié)物質(zhì)可分為有機膠結(jié)物質(zhì)和無機膠結(jié)物質(zhì)，有機膠結(jié)物質(zhì)主要包括有機質(zhì)、微生物及其代謝產(chǎn)物、植物根系及其分泌物等，無機膠結(jié)物質(zhì)主要包括黏粒、多價金屬離子和氧化物等[11]，其中，表面帶有不同類型和不同數(shù)量電荷的鐵鋁氧化物可以吸附在黏粒表面形成氧化物膠膜，促使土粒更加緊密，形成穩(wěn)定的團聚體[12-13]。

興安落葉松是我國寒溫帶的地帶性植被，對呼倫貝爾乃至全國的生態(tài)環(huán)境發(fā)揮著極深遠的保護作用。有學(xué)者研究了凍土區(qū)、沼澤區(qū)的興安落葉松林土壤團聚體及其有機碳特征，分析了重度火燒和撫育間伐對興安落葉松林土壤團聚體及其有機碳的影響。結(jié)果表明，凍土區(qū)土壤以黏粒和粉粒為主，砂粒含量最少[14]；沼澤林0～40 cm土層均以粒徑>2.000 mm水穩(wěn)性團聚體含量最高[15]；重度火燒樣地0～20 cm土層中粒徑≥0.250～2.000 mm團聚體有機碳含量較對照均有不同程度下降[16]；撫育間伐對土壤團聚體穩(wěn)定性的影響主要集中在0～10 cm土層，30%強度間伐樣地土壤水穩(wěn)性最優(yōu)[17]。因立地條件差異及海拔差異，興安落葉松林呈現(xiàn)不同的森林類型，其中分布廣且具代表性的有杜鵑-興安落葉松林、杜香-興安落葉松林、草類-興安落葉松林等，其土壤有機碳含量大小依次為杜香-興安落葉松林>杜鵑-興安落葉松林>草類-興安落葉松林[18-19]。不同林型興安落葉松林土壤團聚體及其有機碳的分布是否存在差異，林下植被是否會對興安落葉松林土壤團聚體的穩(wěn)定性產(chǎn)生影響，是目前亟待探討的問題。本研究以內(nèi)蒙古大興安嶺地區(qū)不同林型的興安落葉松天然林為研究對象，分析其土壤團聚體有機碳的分布特征及影響因素，探明不同林型對土壤團聚體穩(wěn)定性及有機碳分布的影響，以期為我國興安落葉松林的可持續(xù)經(jīng)營、碳匯功能提升提供科學(xué)參考。

1 材料與方法

1.1 研究區(qū)概況及樣地調(diào)查

研究區(qū)位于內(nèi)蒙古大興安嶺森林生態(tài)系統(tǒng)國家野外科學(xué)觀測研究站的原始林試驗區(qū)(121°30′～121°31′E，50°49′～50°51′N)，地處寒溫帶濕潤氣候區(qū)，年均氣溫-5.4 ℃，年降水量 450～550 mm。該區(qū)土壤為棕色針葉林土，優(yōu)勢樹種為興安落葉松(Larixgmelinii)，并伴生有白樺(Betulaplatyphylla)、山楊(Populusdavidiana)等。林下植物主要有杜鵑(Rhododendronsimsii)、杜香(Ledumpalustre)、柴樺(Betulafruticosa)、越橘(Vacciniumvitis-idaea)、舞鶴草 (Maianthemumbifolium)、紅花鹿蹄草(Pyrolaasarifoliasubsp.incarnata)等[20]。

在試驗區(qū)(根河林業(yè)局境內(nèi))，按照草類-興安落葉松林(簡稱草類松林)、杜香-興安落葉松林(簡稱杜香松林)、杜鵑-興安落葉松林(簡稱杜鵑松林)3種林型設(shè)置樣地28塊，每塊樣地規(guī)格為30 m×30 m。記錄各樣地的經(jīng)緯度、海拔、坡度、坡向、坡位等地形信息，并對各樣地進行每木調(diào)查和灌草調(diào)查，樣地基本情況見表1。

表1 內(nèi)蒙古不同興安落葉松林型基本情況

1.2 土壤樣品采集及測定

土壤樣品采集時間為2017年7—8月，在每個樣地內(nèi)，沿對角線挖取3個土壤剖面，去除表面凋落物，以距地面0～10 cm、≥10～20 cm、≥20～40 cm和≥40～60 cm分層取土樣，然后將樣地內(nèi)各點同一土層的土壤混合均勻，取約1 kg待測土樣裝入塑封袋內(nèi)；并取環(huán)刀土用于土壤含水量和容重的測定。將采集的土樣帶回實驗室，去除土樣表面的植物殘體及石礫，自然風(fēng)干，過2 mm篩，用于土壤各指標的測定。

參照Six等[22]和Elliott[23]的濕篩法在微團聚體分離裝置中對各粒徑土壤團聚體進行提取與測定：在0.250 mm濾膜上放置50個玻璃珠，利用水的流動分散土壤顆粒，分離出>0.250 mm的土壤顆粒組分；在0.053 mm篩上收集微團聚體，用濕篩法從水穩(wěn)性微團聚體中分離出易分散的黏粉級組分(>0.053 mm)；離心分離懸濁液，獲得易分散的黏粉級組分(<0.053 mm)，并于60 ℃烘干，稱質(zhì)量。不同粒徑土壤團聚體內(nèi)的有機碳含量均利用TOC分析儀(SHIMADZU TOC-V CPH，日本)測定。

1.3 數(shù)據(jù)處理及統(tǒng)計

各粒徑團聚體的質(zhì)量分數(shù)、團聚體有機碳含量及其貢獻率的計算方法如下：

各粒徑土壤團聚體質(zhì)量分數(shù)=該粒徑團聚體質(zhì)量(g)/土壤樣品總質(zhì)量(g)×100%；

(1)

各粒徑土壤團聚體含碳量(g/kg)=該粒徑團聚體中有機碳質(zhì)量(g)/該粒徑團聚體質(zhì)量(kg)；

(2)

各粒徑土壤團聚體有機碳含量(g/kg)=該粒徑團聚體含碳量(g/kg)×該粒徑團聚體質(zhì)量分數(shù)(%)；

(3)

各粒徑土壤團聚體有機碳貢獻率=該粒徑土壤團聚體有機碳含量(g/kg)/土壤總有機碳含量(g/kg)×100%。

(4)

平均質(zhì)量直徑(MWD，公式中以MMWD表示)、幾何平均直徑(GMD，公式中以GGMD表示)和分形維數(shù)(D)等用于表征團聚體穩(wěn)定性的指標計算方法如下[24]：

(5)

(6)

(7)

采用單因素方差分析(One-way ANOVA)和多重比較(LSD)進行不同數(shù)據(jù)組間差異顯著性檢驗，采用Pearson相關(guān)分析和Linear逐步回歸方法分析各土壤指標對土壤團聚體及其有機碳含量的影響。數(shù)據(jù)處理及繪圖在Excel 2016、SPSS 22.0 軟件中完成。

2 結(jié)果與分析

2.1 不同林型土壤團聚體特征

2.1.1 土壤團聚體組成特征

不同林型興安落葉松林土壤團聚體的組成與剖面分布特征見表2和圖1。總體而言，不同粒徑土壤團聚體含量在各林型間的規(guī)律不同。粒徑≥0.250～2.000 mm團聚體在林型間的含量大小依次為杜鵑松林>草類松林>杜香松林；粒徑≥0.053～0.250 mm團聚體在林型間的含量大小依次為草類松林>杜香松林>杜鵑松林；粒徑<0.053 mm團聚體在林型間的含量大小排序與粒徑≥0.250～2.000 mm團聚體恰好相反，即杜香松林>草類松林>杜鵑松林。其中，草類松林粒徑≥0.053～0.250 mm團聚體含量顯著高于杜香松林和杜鵑松林(P<0.05，F(xiàn)=5.829，df=2)。各林型興安落葉松林土壤團聚體基本以大團聚體(粒徑≥0.250～2.000 mm)為主，占整個團聚體含量的40.86%～48.02%；草類松林和杜鵑松林各粒徑團聚體含量從大到小均為≥0.250～2.000 mm團聚體、<0.053 mm團聚體、≥0.053～0.250 mm團聚體，而杜香松林各粒徑團聚體含量從大到小均為<0.053 mm團聚體、≥0.250～2.000 mm團聚體、≥0.053～0.250 mm團聚體。

表2 不同林型興安落葉松林各粒徑土壤團聚體指標值

各林型粒徑≥0.250～2.000 mm團聚體均具有明顯的表聚特征，土壤表層(0～10 cm)含量顯著大于其他各層；而表層粒徑<0.053 mm和粒徑≥0.053～0.250 mm團聚體含量明顯低于其他各土層。粒徑≥0.250～2.000 mm團聚體含量在各林型間的差異主要表現(xiàn)在≥40～60 cm土層，草類松林和杜鵑松林顯著高于杜香松林(P<0.05，F(xiàn)=10.225，df=2)。粒徑≥0.053～0.250 mm土壤團聚體含量在各林型間的差異主要表現(xiàn)在0～10 cm和≥20～40 cm土層，均為草類松林顯著大于杜鵑松林(P<0.05，F(xiàn)0～10 cm=10.659，F(xiàn)20～40 cm=6.293，df=1)。粒徑<0.053 mm土壤團聚體含量在各林型間的差異特征與粒徑≥0.250～2.000 mm土壤團聚體較一致，也主要表現(xiàn)在≥40～60 cm土層，但大小關(guān)系恰好相反，即草類松林和杜鵑松林顯著低于杜香松林(P<0.05，F(xiàn)=11.997，df=2)。

圖1 不同林型興安落葉松林土壤團聚體剖面特征Fig.1 Profile characteristics of soil aggregates in different Larix gmelinii forest types

2.1.2 土壤團聚體穩(wěn)定性特征

不同林型興安落葉松林土壤團聚體穩(wěn)定性指標與剖面分布特征見表2和圖2。總體而言，平均質(zhì)量直徑(MWD)和幾何平均直徑(GMD)在各林型間無顯著差異，杜鵑松林中較大，草類松林的分形維數(shù)(D)顯著小于杜香松林(P<0.05，F(xiàn)=5.103，df=1)。由圖2可知，不同林型各穩(wěn)定性指標在各土層間均無顯著性差異，不同土層各穩(wěn)定性指標在各林型間的顯著性差異主要發(fā)生在≥40～60 cm土層。≥40～60 cm土層，GMD值在3種林型間存在顯著差異(P<0.05，F(xiàn)=8.020，df=2)，杜香松林的MWD值顯著小于草類松林和杜鵑松林的(P<0.05，F(xiàn)=10.483，df=2)。隨土層深度增加，草類松林和杜鵑松林的變化趨勢較一致，即MWD和GMD呈先減小后增大趨勢，D呈先增大后減小趨勢；杜香松林的MWD和GMD呈逐漸減小趨勢，D呈逐漸增大趨勢。

圖2 不同林型興安落葉松林土壤團聚體穩(wěn)定性 指標的剖面特征Fig.2 Profile characteristics of aggregate stability indexes in different Larix gmelinii forest types

圖3 不同林型興安落葉松林土壤團聚體有機碳剖面特征Fig.3 Profile characteristics of soil aggregate organic carbon in different Larix gmelinii forest types

2.2 不同林型土壤團聚體有機碳特征

2.2.1 土壤團聚體有機碳含量

不同林型興安落葉松林土壤團聚體有機碳含量無顯著差異，變化范圍為5.50～50.10 g/kg(表2)。杜香松林各粒徑團聚體有機碳含量均表現(xiàn)為最高；粒徑≥0.250～2.000 mm土壤團聚體有機碳在林型間的含量大小依次為杜香松林>杜鵑松林>草類松林；粒徑≥0.053～0.250 mm和<0.053 mm土壤團聚體有機碳含量在林型間的大小均表現(xiàn)為杜香松林>草類松林>杜鵑松林。經(jīng)分析可知(圖3)，不同林型各粒徑團聚體有機碳含量均呈現(xiàn)明顯的垂直分布特征，即隨土層深度增加而遞減；遞減速率各粒徑大小依次為≥0.250～2.000 mm團聚體、≥0.053～0.250 mm團聚體、0.053 mm團聚體；且具有明顯的表聚特征。不同土層各粒徑土壤團聚體有機碳含量在各林型間的規(guī)律存在一定差異。在≥20 cm土層中，粒徑≥0.250～2.000 mm土壤團聚體有機碳在各林型間的含量大小排序依次為杜香松林>杜鵑松林>草類松林，20 cm以下土層則表現(xiàn)為杜鵑松林>草類松林>杜香松林，表明杜香松林的粒徑≥0.250～2.000 mm土壤團聚體有機碳隨剖面下降速率較快。20 cm以上土層，杜鵑松林粒徑≥0.053～0.250 mm土壤團聚體有機碳含量最低；20 cm以下土層，則杜香松林為最低，表明杜香松林的粒徑≥0.053～0.250 mm土壤團聚體有機碳隨剖面下降速率較快。≥10 cm以上土層，杜香松林的粒徑< 0.053 mm土壤團聚體有機碳含量最低；10 cm以下土層，則杜鵑松林為最低，表明杜香松林粒徑<0.053 mm土壤團聚體有機碳隨剖面下降速率較慢。

2.2.2 土壤團聚體有機碳貢獻率

分析發(fā)現(xiàn)，不同林型興安落葉松林各粒徑土壤團聚體碳貢獻率為14.65%～58.01%(表2)。草類松林粒徑≥0.053～0.250 mm土壤團聚體有機碳占比顯著高于杜香松林(P<0.05，F(xiàn)=3.977，df=1)，粒徑<0.053 mm土壤團聚體有機碳占比顯著高于杜鵑松林(P<0.05，F(xiàn)=5.165，df=1)。

各林型不同粒徑土壤團聚體有機碳貢獻率呈現(xiàn)不同的垂直分布規(guī)律(圖4)：總體上，隨土層深度增加，粒徑≥0.250～2.000 mm土壤團聚體有機碳貢獻率呈下降趨勢，而粒徑≥0.053～0.250 mm和<0.053 mm土壤團聚體則呈升高趨勢，其中粒徑≥0.053～0.250 mm的貢獻率較穩(wěn)定；各林型間的顯著差異主要表現(xiàn)在20 cm以下土層，草類松林和杜香松林粒徑≥0.250～2.000 mm土壤團聚體有機碳貢獻率顯著小于杜鵑松林(P<0.05，F(xiàn)=10.535，df=2)，而粒徑<0.053 mm土壤團聚體有機碳貢獻率顯著大于杜鵑松林(P<0.05，F(xiàn)=12.967，df=2)。隨著土層深度增加，不同林型各粒徑土壤團聚體有機碳貢獻率發(fā)生明顯變化，≥20 cm土層中，草類松林和杜香松林均以粒徑≥0.250～2.000 mm土壤團聚體有機碳占主導(dǎo)，20 cm以下土層，則均以粒徑<0.053 mm土壤團聚體有機碳占主導(dǎo)；杜鵑松林則一直以粒徑≥0.250～2.000 mm土壤團聚體有機碳占主導(dǎo)。粒徑≥0.250～2.000 mm土壤團聚體有機碳貢獻率在各土層均以杜鵑松林為最高；≥40 cm土層中，粒徑≥0.250～2.000 mm土壤團聚體有機碳貢獻率在各林型間的大小排序依次為杜鵑松林>杜香松林>草類松林。≥20 cm土層中，粒徑≥0.053～0.250 mm土壤團聚體貢獻率以杜鵑松林最小，20 cm以下土層，則以杜香松林最小。粒徑<0.053 mm土壤團聚體有機碳貢獻率在各土層均以杜鵑松林為最低，而杜香松林較高。

圖4 不同林型興安落葉松林土壤團聚體有機碳貢 獻率的剖面特征Fig.4 Profile characteristics of contribution rates of soil aggregate organic carbon in different Larix gmelinii forest types

2.3 不同林型土壤團聚體及其有機碳影響因素分析

2.3.1 土壤團聚體特征值與土壤各因子間的相關(guān)分析

表3 興安落葉松林土壤團聚體特征值與土壤理化因子間的相關(guān)性

2.3.2 土壤團聚體與土壤各因子間的逐步回歸分析

為進一步明確不同林型興安落葉松林各土壤因子對團聚體的影響，利用逐步回歸分析方法篩選出對不同林型有顯著影響的主導(dǎo)因子，并建立土壤團聚體與土壤因子間的回歸方程。以3種林型各粒徑的土壤，團聚體含量為因變量Y，分別記為Y11(草類-興安落葉松林中粒徑≥0.250～2.000 mm土壤團聚體)、Y12(草類-興安落葉松林中粒徑≥0.053～0.250 mm土壤團聚體)、Y13(草類-興安落葉松林中粒徑<0.053 mm土壤團聚體)、Y21(杜香-興安落葉松林中粒徑≥0.250～2.000 mm土壤團聚體)、Y22(杜香-興安落葉松林中粒徑≥0.053～0.250 mm土壤團聚體)、Y23(杜香-興安落葉松林中粒徑<0.053 mm土壤團聚體)、Y31(杜鵑-興安落葉松林中粒徑≥0.250～2.000 mm土壤團聚體)、Y32(杜鵑-興安落葉松林中粒徑≥0.053～0.250 mm土壤團聚體)、Y33(杜鵑-興安落葉松林中粒徑<0.053 mm土壤團聚體)，以各土壤因子為自變量，進行逐步線性回歸，結(jié)果見表4。

表4 不同林型興安落葉松林各粒徑土壤團聚體回歸方程

3 討 論

由于有機質(zhì)輸入的差異和人為擾動的影響，不同土地利用類型、森林類型的土壤團聚體含量和穩(wěn)定性存在差異[7-8]。在本研究中，各林型興安落葉松林土壤團聚體的主要存在形式為粒徑>0.250 mm大團聚體。粒徑>0.250 mm團聚體的含量可在一定程度上表征土壤質(zhì)量的優(yōu)劣，其含量越高，土壤團聚體越穩(wěn)定，土壤結(jié)構(gòu)越好，質(zhì)量越佳[25]。林地植物根系通過纏繞和聯(lián)結(jié)土壤顆粒并釋放分泌物，促進了土壤大團聚體的形成與穩(wěn)定[5]。由于3種林型興安落葉松林地上、地下凋落物的輸入數(shù)量和質(zhì)量及其分解速率不同，改變了土壤有機碳的積累和含量，使形成大團聚體的膠結(jié)物質(zhì)發(fā)生差異，最終表現(xiàn)出不同的團聚體組成特征[26-27]。杜鵑-興安落葉松林的土壤大團聚體(粒徑≥0.250～2.000 mm)含量在3種林型中最高，而草類-興安落葉松林和杜香-興安落葉松林的粒徑<0.250 mm土壤團聚體占比較杜鵑松林高，表明草類-興安落葉松林和杜香-興安落葉松林土壤團聚體在向小粒徑方向演變。

團聚體組成及有機碳含量對穩(wěn)定土壤結(jié)構(gòu)和維持肥力有重要作用[28]。已有研究發(fā)現(xiàn)，大團聚體有機碳含量比微團聚體高[29]。3種林型興安落葉松林土壤團聚體有機碳含量均表現(xiàn)為粒徑≥0.250～2.000 mm大團聚體的有機碳含量最高，這與王富華等[24]、Mikha等[30]的研究結(jié)果一致。林地在促進黏粉粒和微團聚體形成大團聚體的同時，使更多的土壤有機碳向大團聚體富集，使大團聚體成為碳賦存的主體。各粒徑均表現(xiàn)為杜香-興安落葉松林的土壤團聚體有機碳含量最高。相同氣候類型和土壤母質(zhì)條件下，植物自身的生物學(xué)和生態(tài)學(xué)特性對土壤團聚體有機碳的形成和分布有著極大貢獻[23,31]。3種林型均以興安落葉松為絕對優(yōu)勢種，林分結(jié)構(gòu)簡單，灌木與草本的分布存在明顯垂直地帶性[18]。研究區(qū)杜香-興安落葉松林和杜鵑-興安落葉松林林分密度較大，杜香主要分布于土壤水分含量高、土壤溫度低的陰坡、半陰坡，群落發(fā)展穩(wěn)定，植被蓋度大，土壤環(huán)境得到改善，較適宜微生物活動，同時地表凋落物多，增加了土壤有機質(zhì)的輸入量[18,32]，所以杜香-興安落葉松林的土壤有機碳含量最高，為土壤團聚體提供了更多的碳源；而興安杜鵑雖與杜香同屬杜鵑花科灌木，但葉片近革質(zhì)，不易于分解。草類-興安落葉松林地處山坡，土層薄，胸徑小、根系不發(fā)達的樹木極易倒折，經(jīng)過長時間的自然更新后導(dǎo)致林分密度小，所以有機碳的輸入量較少[18]。研究區(qū)地處大興安嶺北部原始林區(qū)，未受人為干擾，植被生長旺盛，凋落物層較厚，加之興安落葉松的根系多分布于土壤表層，分解后形成的腐殖質(zhì)主要在土壤表層積累[26]，土壤有機碳含量呈現(xiàn)明顯的表聚特征，而土壤有機碳是大團聚體的主要膠結(jié)物質(zhì)，所以土壤大團聚體有機碳含量的表聚特征更加明顯，且隨深度增加下降明顯。

興安落葉松林TOC含量與MWD、GMD、D的相關(guān)性顯著，說明土壤有機碳含量影響興安落葉松林土壤團聚體的穩(wěn)定性，這與趙友朋等[4]、王心怡等[33]對不同林分類型土壤團聚體及其穩(wěn)定性的研究結(jié)果較一致。TOC含量與粒徑≥0.250～2.000 mm土壤團聚體含量的相關(guān)性顯著說明，土壤有機碳含量對興安落葉松林土壤團聚體的形成有顯著影響，其促進了大團聚體(粒徑≥0.250～2.000 mm)的形成，是大團聚體的主要膠結(jié)物質(zhì)。這一研究結(jié)果與劉艷等[34]對北京地區(qū)栓皮櫟和油松人工林的研究結(jié)果較一致。微團聚體通過有機質(zhì)黏合形成大團聚體[2]，土壤有機碳含量增加為土壤中較小顆粒膠結(jié)成較大團聚體創(chuàng)造了有利條件。粒徑≥0.250～2.000 mm粒徑團聚體含量、MWD、GMD均與pH、SWC、AK、AP、TP存在顯著相關(guān)性，表明土壤的含水量、養(yǎng)分含量和酸堿水平對興安落葉松林土壤大團聚體的形成和土壤結(jié)構(gòu)的穩(wěn)定具有重要影響，其影響主要通過改變土壤有機碳的積累來實現(xiàn)。土壤中氮、磷、鉀的增加可以促進植物生長，從而增加凋落物歸還量、根系分泌物和提高微生物活性，有利于土壤有機質(zhì)的重要組成成分多糖的累積[35-36]。pH主要通過影響土壤微生物的種類、數(shù)量和活性而對土壤有機質(zhì)的分解產(chǎn)生影響[37]。土壤各金屬氧化物含量與微團聚體(粒徑<0.250 mm)含量間的顯著相關(guān)性表明，F(xiàn)e、Al等金屬氧化物對團聚體的膠結(jié)作用更多體現(xiàn)在微團聚體中。金屬氧化物表面常帶正電荷，可以吸附在黏粒表面而形成氧化物膠膜[12]，在土壤團聚過程中充當“橋鍵”作用，促進土粒更加緊密，形成黏粒-多價金屬-有機質(zhì)復(fù)合體，改變土壤結(jié)構(gòu)[38]，這與金屬氧化物在土壤中的形態(tài)、表面電荷密度等有關(guān)[39]。

綜上可知，興安落葉松林土壤團聚體以大團聚體(粒徑≥0.250～2.000 mm)為主，各林型大小表現(xiàn)為杜鵑-興安落葉松林>草類-興安落葉松林>杜香-興安落葉松林。各粒徑土壤團聚體有機碳含量均為杜香-興安落葉松林最高。TOC是3種林型粒徑≥0.250～2.000 mm、<0.053 mm土壤團聚體的共同主導(dǎo)因子，而粒徑≥0.053～0.250 mm土壤團聚體的主導(dǎo)因子因林型不同而不同。土壤理化指標對土壤團聚體的形成和穩(wěn)定具有重要影響，土壤含水量、養(yǎng)分含量較高的酸性土壤更有利于興安落葉松林土壤大團聚體的形成和土壤結(jié)構(gòu)的穩(wěn)定，具體影響機制還有待深入研究。