劉慧霞， 孫宗玖， 石宇堃， 武文超， 鄭麗， 艾提劍

(新疆農(nóng)業(yè)大學草業(yè)學院，烏魯木齊830052)

荒漠生態(tài)系統(tǒng)是我國干旱區(qū)典型的生態(tài)系統(tǒng)之一，具有獨特的特征和功能[1]。梭梭(Haloxylonammodendron)和白梭梭(Haloxylonpersicum)是荒漠生態(tài)區(qū)重要的優(yōu)勢物種[2-3]，主要分布于我國西北地區(qū)，其中以準噶爾盆地古爾班通古特沙漠分布范圍最為廣泛，占我國梭梭荒漠總面積的68%，占新疆梭梭荒漠總面積的94%，對維護荒漠區(qū)生態(tài)環(huán)境穩(wěn)定具有十分重要的作用[4-6]。古爾班通古特沙漠生態(tài)系統(tǒng)極其脆弱，在梭梭荒漠的防護下，表層細粒物質(zhì)的積累使風沙土得到固定、土壤有機碳含量得到提高，風沙土的營養(yǎng)成分得到改善[7]，由此表明，梭梭荒漠能夠較好地起到防風固沙的作用。然而，研究表明，準噶爾盆地天然梭梭荒漠近年來出現(xiàn)大規(guī)模退化現(xiàn)象，嚴重影響了梭梭荒漠的固碳能力[6]。為更好地了解梭梭荒漠的固碳現(xiàn)狀，定量準確評估其土壤有機碳含量以及查明其影響因素至關重要。

土壤有機碳庫不僅影響土壤的物理、化學和生物特性，同時也會影響土壤肥力和植物的生產(chǎn)力[8]。目前，有關土壤有機碳的研究多集中在濕潤、半濕潤地區(qū)的森林、草地等生態(tài)系統(tǒng)[9-12]，而干旱區(qū)沙質(zhì)荒漠生態(tài)系統(tǒng)土壤有機碳的研究，往往因其植被類型單一、分布稀疏而被忽略。有學者認為沙質(zhì)荒漠土壤碳密度幾乎為零[13]；但也有一些學者認為天然梭梭荒漠群落覆蓋度雖較低，但其分布范圍較廣，而土壤有機碳儲量取決于植被生長狀況，因此不能被忽略[14-15]。吉小敏等[14]指出，新疆甘家湖自然保護區(qū)梭梭土壤有機碳平均含量為1.79～4.17 g·kg-1；牛攀新等[15]指出，準噶爾盆地梭梭0—100 cm各土層土壤有機碳含量為0.08～0.67 g·kg-1。前人研究多集中于樣點及樣線尺度，在區(qū)域尺度上研究不足，同時對梭梭荒漠生態(tài)系統(tǒng)土壤有機碳含量與植被、土壤、氣候、地理等因子間相互關系的研究相對缺乏。因此，研究荒漠生態(tài)系統(tǒng)中土壤有機碳含量與環(huán)境因子間的關系以及環(huán)境因子對梭梭荒漠土壤有機碳含量的影響具有重要科學意義。本研究以準噶爾盆地梭梭沙質(zhì)荒漠為對象，采用路線調(diào)查和典型樣地布設相結合方法，通過對不同區(qū)域典型梭梭群落土壤有機碳含量進行測定，運用統(tǒng)計學和數(shù)量生態(tài)學分析方法，研究土壤有機碳含量與環(huán)境因子之間的內(nèi)在聯(lián)系，以期為退化梭梭荒漠的治理及可持續(xù)利用提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于為準噶爾盆地，地理位置處于N 44.24°—46.70°，E 82.68°—91.23°，海拔224～964 m，年均降水小于200 mm，年均溫度6.09～31.79 ℃，≥10 ℃年積溫3 000～3 500 ℃，氣候干旱，屬于典型的溫帶荒漠氣候。試驗區(qū)植被群落蓋度小于30%，土壤質(zhì)地為沙質(zhì)，植被以梭梭(H.ammodendron)、白梭梭(H.persicum)為建群種，沙拐棗(Calligonummongolicum)、駝絨藜(Ceratoideslatens)、沙漠絹蒿(Seriphidiumsantolinum)、刺沙蓬(Salsolaruthenica)、琵琶柴(Reaumuriasongonica)、假木賊(Anabasissp.)等為伴生種，同時也混有一些短命、類短命植物。沙質(zhì)荒漠草地研究區(qū)樣點基本概況詳見表1。

續(xù)表Continued

1.2 研究方法

1.2.1樣地選取及野外取樣 2018年8—9月，在阿勒泰地區(qū)、昌吉回族自治州、克拉瑪依市、博爾塔拉蒙古自治州及塔城地區(qū)布設梭梭荒漠典型樣點20個(表1和圖1)。每個樣點內(nèi)設置1個100 m×100 m的樣地，樣地內(nèi)等距離布設5個10 m×10 m的灌木樣方，分物種收獲當年新生枝條的生物量，裝袋標記并稱重；同時等距離布設3條平行樣線，每條樣線等距離布設3個1 m×1 m的草本樣方，齊地分種刈割樣方內(nèi)所有物種的生物量，同時收集樣方內(nèi)所有地上凋落物，裝袋標記并稱重。完成地上植被測定后，在樣地中心挖取面積為1 m×1 m的土壤剖面，采用切土塊法，在3個垂直剖面上采用固定深度法獲取0—5、5—10、10—20、20—30、30—50、50—70、70—100 cm土層的土壤樣品，用于土壤養(yǎng)分分析；同時采用環(huán)刀法進行土壤容重的測定，分裝帶回實驗室進行室內(nèi)分析。

圖1 研究區(qū)概況及樣點布設圖

1.2.2室內(nèi)分析 土壤樣品在剔除可見植物根系、動物糞便等雜物后，室內(nèi)陰干將土樣混勻，過0.25 mm篩用于土壤有機碳的測定；過1 mm篩用于土壤pH及電導率的測定。土壤有機碳含量采用重鉻酸鉀外加熱法進行測定，pH采用水土比5∶1酸度計法進行測定，電導率采用水土比5∶1電導法進行測定，容重采用105 ℃、24 h烘干法進行測定[16]。

1.3 環(huán)境因素數(shù)據(jù)獲取

通過國家氣象臺站點(www.cnern.org.cn)下載近60 a的氣象數(shù)據(jù)，以平均值為基準提取各氣象站點的氣象數(shù)據(jù)，利用ArcGIS 10.2中克里金插值法獲取各采樣點≥10 ℃的年積溫、年均溫度、年均降水量等氣象數(shù)據(jù)，以表示實際采樣點的氣候條件。經(jīng)緯度、海拔采用手持GPS測定儀進行實地測定。以2000—2018年GLEAM V3數(shù)據(jù)集(www.GLEAM.eu)衛(wèi)星觀測陸地蒸發(fā)量和根區(qū)土壤濕度數(shù)據(jù)為基準，根據(jù)研究區(qū)的矢量數(shù)據(jù)對影像進行裁剪，按照采樣點地理坐標反演得到土壤濕度和蒸散發(fā)量[17]。

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

利用Microsoft Excel 2013和 SPSS 20.0 對數(shù)據(jù)進行統(tǒng)計分析，采用單因素方差、獨立樣本t檢驗分析土壤有機碳含量間的差異，多重比較采用Duncan方法檢驗顯著性。采用Canoco 5.0進行冗余分析，利用偏冗余分析進一步得出各環(huán)境因子的解釋率。數(shù)據(jù)可視化在 Origin 9.0下完成。

2 結果與分析

2.1 梭梭荒漠土壤有機碳總體體征

沙質(zhì)荒漠草地不同土層土壤有機碳含量結果(表2)表明，0—100 cm土層有機碳含量平均為1.62 g·kg-1，最小值為0.00 g·kg-1，最大值為11.46 g·kg-1，變異系數(shù)為81.67%，呈中等變異水平[18]。0—5、5—10、10—20、20—30、30—50、50—70、70—100 cm土壤有機碳含量依次為2.32、1.76、1.62、1.58、1.43、1.35、1.26 g·kg-1，變異系數(shù)為64.33%～87.07%，均呈中等變異水平。其中，0—5 cm土壤有機碳含量顯著高于5—100 cm土層。對土壤有機碳含量與土層深度進行相關分析表明(圖2)，土壤有機碳含量與土層深度極顯著相關，即隨著土層深度的增加，土壤有機碳含量逐漸降低。

表2 不同土層的土壤有機碳含量

注：不同小寫字母代表不同土層間差異在P<0.05水平顯著。

2.2 兩種類型梭梭荒漠土壤有機碳比較

不同植被類型土壤有機碳含量及密度剖面分布特征結果(圖3)表明，白梭梭荒漠與梭梭荒漠間0—100 cm土壤有機碳含量呈極顯著差異(P<0.01)；隨土層深度增加，白梭梭荒漠與梭梭荒漠土壤的有機碳含量均呈降低趨勢，且0—5 cm土層有機碳含量均顯著高于10—100 cm土層。其中，梭梭荒漠的土壤有機碳含量(1.74 g·kg-1)是白梭梭荒漠土壤有機碳含量(1.25 g·kg-1)的1.40倍；具體分析兩種植被下不同土層的有機碳含量表明，僅10—20 cm和20—30 cm土層梭梭荒漠土壤的有機碳含量顯著高于白梭梭荒漠。與白梭梭荒漠相比，梭梭荒漠土壤有機碳含量從0—5 cm到5—10 cm土層的降幅更明顯。

注：不同小寫字母代表不同土層土壤有機碳含量差異在P<0.05水平顯著；*表示不同植被類型間差異在P<0.05水平顯著。

2.3 海拔對土壤有機碳的影響

土壤有機碳含量隨海拔高度變化結果(表3)表明，不同海拔梯度0—100 cm土層土壤有機碳含量的平均值間存在差異。500～700 m海拔間的土壤有機碳含量最低，為1.09 g·kg-1，顯著低于其他海拔(43.12%～79.82%)，且各土層均表現(xiàn)為海拔500～700 m土壤有機碳含量最低；而海拔<300 m或>700 m 以及300～500 m的土壤有機碳含量間差異不顯著。對同一海拔下不同土層的土壤有機碳含量進行比較發(fā)現(xiàn)，在海拔<300 m和>700 m條件下，各土層間土壤有機碳含量差異不顯著；海拔300～500 m下，0—5 cm土層土壤的有機碳含量較10—100 cm各土層顯著增加55.67%～114.18%；海拔500～700 m下，0—5 cm土層較其余各土層土壤有機碳含量顯著增加72.27%～200.00%。

表3 不同海拔梯度不同土層的土壤有機碳含量

2.4 環(huán)境因素對土壤有機碳含量的影響

對土壤有機碳含量與環(huán)境因子進行冗余分析，如圖4所示，實心箭頭代表響應變量，空心箭頭代表環(huán)境因子，實心箭頭與空心箭頭夾角的余弦值代表土壤有機碳含量及密度與環(huán)境因子間的相關性，箭頭長度反映了環(huán)境因子對響應變量的解釋量[19]。從箭頭夾角的余弦值可以看出(圖4A)，土壤有機碳含量與土壤pH、土壤含水率、根部土壤濕度、表層土壤濕度、緯度、植被覆蓋度、凋落物生物量、土壤電導率、土壤容重、海拔、年均降水量及≥10 ℃年積溫均呈正相關關系，與年均溫度、地上生物量、經(jīng)度、土壤蒸散發(fā)量呈負相關關系。根據(jù)偏冗余進一步分析結果(圖4B)表明，單個環(huán)境因子對土壤有機碳含量的解釋能力依次為土壤pH(28.1%)>植被覆蓋度(20.5%)>土壤含水率(14.8%)>年均降水量(10.3%)>土壤蒸散發(fā)量(7.7%)>表層土壤濕度(5.7%)>年均溫度(3.9%)>土壤電導率(2.3%)>海拔(1.8%)>≥10 ℃年積溫(1.4%)>緯度(1.1%)>經(jīng)度(1.0%)>凋落物生物量(0.9%)>地上生物量(0.3%)>根部土壤濕度(0.2%)。其中，土壤pH、植被覆蓋度和土壤含水率達顯著水平(P<0.05)。綜上所述，引起土壤有機碳含量空間變異的主要影響因素是土壤pH、植被覆蓋度和土壤含水率，其他環(huán)境因子對土壤有機碳含量的影響相對較弱。

注：A：冗余分析；B：偏冗余分析。SOC—土壤有機碳含量；SOCD—土壤有機碳密度； Long—經(jīng)度；Lat—緯度；Ele—海拔；≥10 ℃AT—≥10 ℃年積溫；MAP—年均降水；MAT—年均溫度；BD—土壤容重；SWC—土壤含水率；EC—電導率；Cov—植被覆蓋度；AGB—地上生物量；Litter—凋落物生物量；Evap—蒸散發(fā)；SSM—表層土壤濕度；RSM—根部土壤濕度。

3 討論

梭梭荒漠生態(tài)系統(tǒng)雖然土壤肥力較低，植被稀疏，但植被在防風固沙方面的作用不可忽視。本研究中0—30 cm土層土壤有機碳含量為1.82 g·kg-1，低于世界干旱區(qū)平均水平[20]。一方面與古爾班通古特沙漠氣候極端干旱，蒸發(fā)大、降水少、植被類型及凋落物分解產(chǎn)生的差異有強烈關系；另一方面與人類頻繁的干擾行為(如過牧、采挖、開墾等)有關。有研究指出，梭梭保護區(qū)土壤有機碳含量在高蓋度下為7.70 g·kg-1，而古爾班通古特沙漠南緣(退化梭梭荒漠)土壤有機碳含量僅為0.31 g·kg-1[14-15]。因此，必需對古爾班通古特沙漠南緣退化的天然梭梭荒漠實施保護措施，提高梭梭荒漠土壤有機碳含量，對梭梭沙質(zhì)荒漠土壤恢復及土壤固碳能力的提升具有重要意義。

本研究表明，梭梭沙質(zhì)荒漠土壤有機碳集中于表層土壤，垂直分布表現(xiàn)為隨土層深度的增加逐漸降低，且深層土壤有機碳含量趨于穩(wěn)定。這可能是由于表層土壤通氣狀況較好，植物殘體中各種有機物都能被快速氧化，微生物在通氣良好的土壤中也能加速對腐殖質(zhì)的分解消耗[21]；而深層土壤緊實、孔隙少，凋落物等外源物質(zhì)的輸入較少，土壤溫度低，微生物數(shù)量及活性也較低，進而造成深層土壤有機碳含量顯著低于表層[22]。梭梭群落土壤的有機碳含量顯著高于白梭梭群落土壤，這可能是由于梭梭普遍生長于沙丘間低地，而白梭梭普遍生長于半流動沙丘和固定沙丘[23]。與沙丘相比，丘間低地受到風蝕的現(xiàn)象較弱[24]，風力作用對地表凋落物的再分配影響也較小[25]，進而導致沙質(zhì)草地土壤有機碳含量出現(xiàn)較大的異質(zhì)性。

此外，本研究還發(fā)現(xiàn)，不同土層土壤有機碳含量在不同海拔范圍間差異較小，但最低值均出現(xiàn)在500～700 m海拔范圍內(nèi)。經(jīng)調(diào)查發(fā)現(xiàn)，該海拔范圍的區(qū)域降水較少，梭梭群落總蓋度僅為11.27%，有的地區(qū)梭梭蓋度僅為2%，因此，大面積的地表裸露使風蝕作用加劇、土壤粗?；?，導致土壤中砂粒含量驟增，造成對土壤有機碳的累積量逐漸減少；草本稀疏使得降水不能有效儲存，水土流失嚴重；且500～700 m海拔范圍內(nèi)地上凋落物的固存量也顯著低于平均水平，導致有機質(zhì)的輸入量較少，因此，土壤有機碳含量較低[26]。除>700 m海拔外，其他海拔梯度土壤的有機碳含量均隨土層深度的增加呈現(xiàn)降低趨勢，且變化浮動較小，這可能是由于土壤中新的有機質(zhì)輸入較少、微生物降解基質(zhì)不足、土壤基質(zhì)較為均勻所致。當海拔>700 m時，土壤有機碳含量在10—20 cm土層出現(xiàn)最大值，這可能是凋落物受風力作用的掩埋、穩(wěn)定堆積后的結果[27]，隨著時間的延長，凋落物在土壤微生物的作用下緩慢分解，將有機碳固定在土壤中。雖然不同海拔梯度土壤有機碳含量隨土層的分布狀況不盡相同，但均表現(xiàn)為0—30 cm土層土壤的有機碳含量較高，與前人研究結果相一致[28]。

土壤有機碳的輸入與輸出受多種環(huán)境因子調(diào)控，且有機碳在土壤中的積累存在較大的空間異質(zhì)性[29]。冗余分析表明，土壤pH、土壤含水率和植被覆蓋度是驅(qū)動土壤有機碳含量產(chǎn)生變異的主要因素。本研究中植被覆蓋度與土壤有機碳含量呈顯著正相關。植被覆蓋度不僅可以反映生物量的高低，還可以減緩土壤水分的蒸發(fā)速度，有利于土壤有機碳的積累。除植被覆蓋度外，土壤pH與土壤有機碳含量也呈正相關，與張磊等[30]研究結果一致，可能是由于適度的堿性土壤有利于土壤中硝化細菌的硝化作用，進而有利于荒漠植物對土壤中氮素的利用，促進植物生長[31]。土壤含水率與土壤有機碳含量顯著正相關，荒漠區(qū)植被對水分的需求遠高于其他生境中的植被[30]，光合作用的產(chǎn)物大部分被用于地下根系的生長，根系不斷向下以獲取更多的水資源用于維持植株的生長，進而促進土壤有機碳的積累。但本研究發(fā)現(xiàn)，地上生物量與土壤有機碳含量表現(xiàn)為負相關關系，這可能是由于地上生物量雖然較多，但地表的固定凋落物較少，地表凋落物受強烈風力作用再分配的影響[25]，導致地上生物量對土壤有機碳含量的變化響應不敏感。本研究區(qū)域梭梭群落總植被的覆蓋度均小于30%，土壤有機碳含量平均值為1.62 g·kg-1，均低于甘家湖梭梭自然保護區(qū)[14]。高植被覆蓋度可以減輕風力作用，使凋落物能夠更好地被植物截獲，促進風沙土區(qū)域植被的正向演替[32]，由此表明，天然梭梭荒漠草地具有巨大的碳匯潛力。