祁連山區(qū)成土過程驅(qū)動的土壤地形序列土壤有機碳固存特征①

谷 俊，楊 飛，李德成，宋效東，張甘霖,4*

祁連山區(qū)成土過程驅(qū)動的土壤地形序列土壤有機碳固存特征①

谷 俊1,2，楊 飛1,3，李德成1，宋效東1，張甘霖1,2,4*

(1土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點實驗室(中國科學(xué)院南京土壤研究所)，南京 210008；2 中國科學(xué)院大學(xué)，北京 100049；3 青海師范大學(xué)/青海省自然地理與環(huán)境過程重點實驗室，西寧 810008；4 中國科學(xué)院南京地理與湖泊研究所/流域地理學(xué)重點實驗室，南京 210008)

土壤有機碳(SOC)含量和動態(tài)是有機碳輸入輸出平衡和土壤固持能力共同作用的結(jié)果，前者主要受生物氣候條件控制，而后者則主要受黏粒及無機礦物等土壤理化屬性的影響。本研究通過沿海拔梯度選擇表土地形序列作為“自然試驗場”，探究祁連山區(qū)SOC含量及組分變化的控制因素。研究發(fā)現(xiàn)：該地形序列土壤母質(zhì)主要來源于黃土沉積，土壤黏土礦物以綠泥石和水云母為主，指示該區(qū)域的整體弱風化特征。在巨大的水熱梯度影響下，地形序列內(nèi)土壤成土強度差異明顯，低海拔地區(qū)土壤含有碳酸鹽，隨著海拔上升，碳酸鹽物質(zhì)逐漸淋失，SOC和鐵鋁氧化物含量增加。進一步分析發(fā)現(xiàn)，鐵鋁氧化物是SOC含量及組分方差的主要解釋變量。偏相關(guān)分析顯示，當控制鐵鋁氧化物后，氣候?qū)OC含量及組分的影響不顯著。這表明氣候?qū)OC的影響可能主要通過影響土壤屬性，造成鐵鋁氧化物等屬性的差異，間接影響SOC的長期固存，且該機制主要作用于礦物結(jié)合態(tài)有機碳(MOC)組分。本研究對理解SOC固存及其對氣候變化的動態(tài)響應(yīng)有重要啟示。

土壤地形序列；土壤有機碳；鐵鋁氧化物；祁連山

土壤有機碳(SOC)庫作為最大的陸地碳庫，其微小變化會引起大氣中二氧化碳濃度的明顯變化[1]，并可能導(dǎo)致全球變暖與陸地碳循環(huán)之間的正反饋[2]。高寒草甸、高寒草原等草地是青藏高原主要的生態(tài)系統(tǒng)，占總面積的60% 以上，共儲存了約2.32×1010t有機碳[3]。巨大的碳儲量及對氣候變化的敏感性，使青藏高原在全球陸地生態(tài)系統(tǒng)碳循環(huán)中扮演重要角色[4]。關(guān)于氣候和生物因素對SOC含量的單一和交互影響已有大量研究[5-7]。但目前對SOC動態(tài)變化的預(yù)測仍存在不確定性，主要原因是對SOC含量的控制因素，尤其是土壤因素了解不全面[8]。

目前，土壤增溫試驗和不同階段碳庫調(diào)查結(jié)果被用于探究青藏高原SOC對氣候變化的響應(yīng)。Yang等[9]通過對比1980—2004年間青藏高原SOC含量變化發(fā)現(xiàn)，氣候變暖導(dǎo)致生物量輸入增加，引起SOC含量增加。通過增溫試驗發(fā)現(xiàn)，升溫導(dǎo)致土壤中不穩(wěn)定SOC含量降低[10-12]。盡管增溫試驗和區(qū)域時序?qū)Ρ仍囼灋槔斫釹OC在短時間尺度的氣候響應(yīng)提供了有力證據(jù)[13]，但目前關(guān)于氣候在長時間尺度對SOC的影響仍缺乏系統(tǒng)認識。通過自然氣候梯度研究可以彌補這一缺陷，海拔梯度是研究生態(tài)系統(tǒng)對氣候變化響應(yīng)的“自然試驗(natural experiment)”之一[14]。隨海拔升高，氣候、植被發(fā)生改變，這不僅直接造成SOC的來源差異，還導(dǎo)致成土過程及土壤理化屬性發(fā)生變化[15-16]，并進一步對SOC的分解與固存產(chǎn)生影響。如海拔上升導(dǎo)致降水增加，加速礦物鹽基離子流失，土壤pH下降，酸性至中性且濕潤的成土環(huán)境加速土壤礦物風化，促進土壤發(fā)育[15]。Torn等[17]發(fā)現(xiàn)，隨著土壤發(fā)育程度增加，土壤中非晶質(zhì)礦物會先增加后減少，而非晶質(zhì)礦物具有較強吸附力和較大比表面積是SOC重要的固定機制之一[18]?？傮w而言，目前關(guān)于土壤發(fā)生過程如何影響SOC固存的系統(tǒng)研究十分匱乏。

將SOC分成不同碳池有助于深刻理解SOC的動態(tài)變化。Lavallee等[19]研究表明，將SOC分成顆粒態(tài)有機碳(POC)和礦物結(jié)合態(tài)有機碳(MOC)能簡單有效地概念化土壤有機質(zhì)。兩種組分代表不同周轉(zhuǎn)時間和穩(wěn)定機制的碳池，通常POC是游離的碳，屬于不穩(wěn)定、周轉(zhuǎn)速率快的碳池，而MOC是與黏粉粒級礦物顆粒結(jié)合的碳，相對穩(wěn)定[20]。區(qū)分不同SOC組分對各影響因子的響應(yīng)有助于更好地評估SOC動態(tài)變化及其對環(huán)境變化的響應(yīng)。

本研究在祁連山地區(qū)沿海拔梯度采集地形序列土壤，分析區(qū)域成土特征及其對環(huán)境梯度變化的響應(yīng)規(guī)律，并進一步通過研究SOC及其組分(POC和MOC)與土壤要素及氣候要素的關(guān)系，揭示成土環(huán)境及成土過程對SOC的影響與機制，以更好地理解寒旱區(qū)SOC固存特征及其對氣候變化的潛在響應(yīng)機制。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于青藏高原東北緣的祁連山區(qū)，該區(qū)具有典型大陸性氣候和高原氣候特征，地形是造成該地區(qū)氣候及景觀差異的主要因素。受高原隆升影響，該地區(qū)存在巨大海拔梯度，海拔由5 500 m以上的高山區(qū)向北迅速降至2 000 m以下的山前地帶。海拔是導(dǎo)致區(qū)域氣候變化的主要驅(qū)動因素，隨海拔升高，年均溫下降，年均降水增加。同時該區(qū)域?qū)儆跂|亞季風影響區(qū)的末端，降水呈現(xiàn)由東向西減少的趨勢，年平均降水量為300 mm，年均溫小于4 ℃。自然景觀為從干旱地區(qū)的荒漠草原到亞高山地區(qū)的草原或針葉林以及高山環(huán)境中的草甸或灌木林地。

1.2 樣品采集與測定

在研究區(qū)布設(shè)29個采樣點(圖1A)，海拔跨度為1 800 ～ 3 800 m，年均降水量沿海拔從200 mm增加到600 mm，年均溫從6.5 ℃降至–3 ℃。景觀類型包括荒漠草原和亞高山草原、高山草甸和高山灌叢草甸(圖1B)。每個采樣點取3個位置0 ～ 15 cm的表層土樣混合成1個土壤樣品。

樣品前處理及實驗室土壤理化屬性測定參考《土壤調(diào)查實驗室分析方法》[21]。將土壤樣品風干，除去根系和礫石，研磨過10目篩后，再分別研磨過60目、100目篩備用。使用激光衍射粒度儀(LS230，Beckman Coulter)分析土壤粒度分布；重鉻酸鉀–硫酸消化法測定總有機碳含量(SOC)；水土質(zhì)量比為2.5∶1的混合液使用pH計(pHS-3C， Shanghai REX Instrument Factory，China)測定土壤pH；氣量法測定土壤碳酸鈣相當物含量；檸檬酸三鈉、碳酸氫鈉和連二硫酸鈉(DCB法)提取游離態(tài)鐵鋁氧化物(Ald、Fed)，草酸–草酸銨溶液提取結(jié)晶度較差的非晶質(zhì)鐵鋁氧化物(Feo、Alo)，以上步驟所得到的上清液使用ICP-AES(Optima 8000，PerkinElmer)測定溶液中鐵鋁含量，鐵氧化物的結(jié)晶度用Feo/Fed衡量； X射線衍射(XRD)儀(日本理學(xué) Ultima Ⅳ全自動 X射線衍射儀)測定礦物組成。

圖1 土壤樣點分布(A)與代表性樣點地形地貌及對應(yīng)的土壤剖面(B)

SOC分組采用濕篩法[22]，具體為：10 g風干土加入30 mL 5 g/L六偏磷酸鈉溶液中，振蕩15 h，將得到的混合液過300目篩，將未通過篩的顆粒(POC組分)和通過篩(MOC組分)的泥漿在45 ℃下烘干并稱重，得到不同組分的質(zhì)量。干燥后的樣品研磨并分析其中有機碳含量，不同組分測得的有機碳含量分別乘以其質(zhì)量，并除以總土樣質(zhì)量(10 g)，得到顆粒態(tài)有機碳含量(POC，g/kg)和礦物結(jié)合態(tài)有機碳含量(MOC，g/kg)。

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

本研究獲取的氣候數(shù)據(jù)包括年均氣溫(MAT)和年均降水量(MAP)，數(shù)據(jù)來自中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所中國生態(tài)環(huán)境背景層面建造項目完成的柵格數(shù)據(jù)(分辨率1 000 m)，為1980—1999年的逐月平均值計算生成。歸一化植被指數(shù)(NDVI)來自于青藏高原科學(xué)數(shù)據(jù)中心，為2013—2018年多年NDVI平均值。

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

采用R軟件進行數(shù)據(jù)處理、統(tǒng)計分析和圖表制作，使用偏相關(guān)分析和方差分解量化氣候、土壤因素對SOC含量及組分的影響。

2 結(jié)果與分析

2.1 土壤粒度分布與礦物組成

土壤粒度分布曲線顯示，研究區(qū)地形序列內(nèi)土壤樣品具有較高的粒度分布一致性，大部分樣品在30 ～ 40 μm處存在主要的粒徑峰值，個別樣品峰值出現(xiàn)在70 ～ 80 μm(圖2)。研究區(qū)土壤粉粒含量高，質(zhì)地大多為粉壤質(zhì)，個別為壤質(zhì)，土壤質(zhì)地較均一。在土壤地形序列內(nèi)按海拔梯度選擇10個土壤樣品進行土壤礦物鑒定，XRD圖譜(圖3) 顯示，序列內(nèi)土壤原生礦物以石英和長石為主，黏土礦物以水云母和綠泥石為主，表明研究區(qū)土壤化學(xué)風化程度較弱。序列內(nèi)僅碳酸鹽礦物表現(xiàn)出較高的差異性，在海拔2 600 m及以下的表土中，方解石與白云石普遍大量存在(圖3)。

圖2 全部樣品粒度分布 (n=29)

圖3 地形序列土壤的X射線衍射圖譜(n=10)

2.2 土壤理化屬性變異特征及其水熱驅(qū)動

不同于地形序列內(nèi)土壤質(zhì)地及礦物組成的相對一致性，研究區(qū)土壤pH、碳酸鈣含量、鐵鋁氧化物含量、SOC含量及組分在地形序列內(nèi)表現(xiàn)出較高的空間變異性，其土壤基本理化屬性見表1。其中，SOC含量為(40.0 ± 22.7) g/kg，POC和MOC含量分別為(12.8 ± 9.0) g/kg和(26.9 ± 14.4) g/kg，可見，在本研究區(qū)MOC是SOC主要存在形式。土壤pH呈中性至石灰性，最高為8.9。土壤中Fed含量變化范圍為5.5 ～ 14.8 g/kg，F(xiàn)eo含量變化范圍為0.7 ～ 4.1 g/kg，鐵活化度Feo/Fed變化范圍為9% ～ 32%，表明本地區(qū)鐵的活化度較高。主成分分析可以進一步表征土壤屬性的整體變異特征，結(jié)果顯示，研究區(qū)地形序列內(nèi)氣候尤其是MAP、碳酸鈣含量、海拔與鐵鋁氧化物含量和SOC含量及組分表現(xiàn)出較高的協(xié)變關(guān)系(圖4)，表明隨著海拔上升，降水增加，土壤中碳酸鈣被淋洗，鐵鋁氧化物含量增加，土壤整體發(fā)育強度增加(圖5)。

表1 土樣基本理化屬性

圖4 主成分分析(PCA)：16個變量特征值的荷載圖

2.3 環(huán)境因子對土壤有機碳含量及組分的影響

相關(guān)分析結(jié)果顯示，研究區(qū)SOC及其組分碳含量均與MAP呈極顯著正相關(guān)，與MAT呈極顯著負相關(guān)(表2)。SOC及MOC與氣候因子的相關(guān)系數(shù)接近，均高于POC與氣候因子的相關(guān)性。SOC及其組分碳含量與海拔、NDVI呈極顯著正相關(guān)。SOC及其組分碳含量與土壤pH呈顯著負相關(guān)，與土壤黏粒含量呈極顯著正相關(guān)，而與砂粒、粉粒含量的相關(guān)性不顯著。不同形態(tài)的鐵鋁氧化物均與SOC、POC和MOC呈顯著或極顯著相關(guān)性。MOC與不同形態(tài)鐵鋁氧化物相關(guān)性均高于SOC，而POC與不同鐵鋁氧化物相關(guān)性均低于SOC。對于各形態(tài)鐵鋁氧化物，Ald和Alo與SOC、MOC和POC相關(guān)系數(shù)最高。

偏相關(guān)分析顯示，當控制土壤pH、質(zhì)地和鐵鋁氧化物之后，MAT、MAP與SOC及其組分之間的相關(guān)系數(shù)大幅降低(表3)。在所有變量中Ald、Alo最為重要，具體表現(xiàn)為控制Ald、Alo后，SOC及其組分與氣候因子之間偏相關(guān)系數(shù)下降最為劇烈。POC、MOC與MAT呈極顯著負相關(guān)，但在控制Ald后，相關(guān)系數(shù)降低為–0.34和–0.27。POC、MOC與MAP呈極顯著正相關(guān)，但在控制Alo后，相關(guān)系數(shù)分別降低為0.35和0.26。

方差分解(圖6)結(jié)果顯示，氣候、質(zhì)地和鐵鋁氧化物對于SOC及其組分的重要性有明顯差異。對于SOC及其組分，質(zhì)地的重要性最低，與僅有氣候和質(zhì)地變量的模型相比，加入鐵鋁氧化物變量之后，模型對SOC、POC和MOC含量的解釋能力均增加，分別達到71.2%、47.5%、78.6%，其中對MOC的方差解釋力增加最大；鐵鋁氧化物的單獨解釋能力在SOC和POC、MOC中分別為5.8%、3.1%、15.5%，在MOC中解釋力最高，表明該部分碳庫受到土壤礦物影響更大。雖然鐵鋁氧化物的單獨解釋力大于氣候因子，但氣候因子與鐵鋁氧化物的共同作用不可忽視，鐵鋁氧化物與氣候因子共同解釋力在SOC、POC和MOC分別達到了21.3%、21.5%、17.7%，意味著鐵鋁氧化物與氣候因子對SOC具有強烈的共同作用。

3 討論

3.1 土壤母質(zhì)來源一致性與區(qū)域弱風化特征

土壤粒度特征對判斷土壤母質(zhì)來源有較好的指示作用，尤其對于風積物。本研究中，土壤粒徑峰值大多在30 ～ 40 μm，與典型黃土接近，少量樣品粒徑峰值在70 ～ 80 μm，為砂質(zhì)黃土特征[23]。粒度分布特征指示本研究涉及的表層土壤均為風積物發(fā)育，這與前人在該地區(qū)的研究結(jié)論相一致[24]。本地區(qū)黏土礦物以水云母和綠泥石為主，揭示區(qū)域土壤整體風化程度較弱的特點，這主要是由于高山區(qū)冰凍期長，年均溫低，而祁連山北麓山前地區(qū)氣候干旱，土壤化學(xué)風化受干冷氣候限制。另外，本地區(qū)土壤(表層)均為全新世以來的黃土沉積[25]，土壤發(fā)育時間較短。成土母質(zhì)、氣候與成土?xí)r間共同決定了序列土壤質(zhì)地均一、黏粒含量較低且黏土礦物類型相對一致的特點。

圖5 土壤屬性與MAP關(guān)系的散點圖

表2 土壤有機碳不同組分與環(huán)境變量的相關(guān)性(n=29)

注：*、**、***分別代表在<0.05、<0.01、<0.001水平顯著相關(guān)；下同。

表3 土壤有機碳不同組分與氣候因子之間的偏相關(guān)系數(shù)(n=29)

注：零階代表零階偏相關(guān)系數(shù)；All代表控制所有土壤變量后的偏相關(guān)系數(shù)。

3.2 環(huán)境梯度下的成土強度分異

雖然區(qū)域內(nèi)序列表土主要起源于黃土母質(zhì)，但不同于干旱環(huán)境中成土強度低的典型黃土。本區(qū)域跨巨大的環(huán)境梯度，高程差超過2 000 m，水熱條件差異顯著，土壤屬性在空間上有較強的分異。隨海拔上升，MAP升高，MAT下降，有利于SOC積累和碳酸鹽礦物的溶解及淋洗。山前地區(qū)SOC含量低，碳酸鈣含量高，而隨著海拔上升，土壤類型由山前鈣積土向高山草甸土轉(zhuǎn)變，高山草甸土已經(jīng)完成淋溶脫鈣過程，表層土壤碳酸鈣含量低，并開始土壤交換性鹽基的淋溶酸化過程，土壤pH由弱堿性降至中性或弱酸性。游離態(tài)鐵鋁氧化物與非晶質(zhì)態(tài)鐵鋁氧化物的豐度是成土強度的良好指標[26]。本研究中，游離態(tài)鐵鋁(Fed、Ald)和非晶質(zhì)鐵鋁(Feo、Alo)均與MAP表現(xiàn)出顯著相關(guān)性(圖5)，表明水分條件是序列內(nèi)驅(qū)動土壤母質(zhì)風化與成土強度分異的主要因素。

(氣候因子包括MAT、MAP；質(zhì)地因子包括砂粒、粉粒、黏粒；鐵鋁氧化物因子包括Fed、Ald、Feo、Alo、Feo/Fed)

3.3 鐵鋁氧化物主控有機碳變異

盡管相關(guān)分析顯示氣候因子與SOC及其組分有較高相關(guān)性(表2)，但根據(jù)偏相關(guān)分析的結(jié)果，發(fā)現(xiàn)控制土壤變量后，氣候因子與SOC及其組分間的相關(guān)性有不同程度下降，尤其是控制鐵鋁氧化物時，相關(guān)系數(shù)下降最為顯著(表3)。這表明氣候因子對SOC及其組分的影響更多是通過影響土壤屬性間接影響，尤其是受水熱條件控制的鐵鋁氧化物。這與Doetterl等[27]跨氣候帶研究發(fā)現(xiàn)地球化學(xué)特征與氣候因子交互作用控制SOC儲量的結(jié)論一致。方差分解也進一步表明，氣候因子和鐵鋁氧化物的共同作用在解釋SOC的方差變化中占據(jù)重要地位。

鐵鋁氧化物具有較高的比表面積和帶電表面，能夠吸附SOC的官能團，并形成有機–黏土礦物復(fù)合物[28]，減緩SOC的礦化分解[29]。此外，鐵鋁氧化物參與團聚體的形成[30]，通過物理保護減少SOC與微生物的接觸，間接增強SOC的穩(wěn)定性。這些機制共同作用抑制了SOC的分解，使得SOC更加穩(wěn)定。von Fromm等[31]也發(fā)現(xiàn)，土壤中的鐵鋁氧化物是預(yù)測大陸尺度SOC儲量的最重要因子，其中非晶質(zhì)鐵鋁氧化物尤為重要。同樣，其他學(xué)者在青藏高原的研究也強調(diào)了土壤鐵鋁氧化物對草地SOC儲量的重要影響[32-33]。

與鐵鋁氧化物相比，本研究中土壤質(zhì)地對SOC的影響較弱。這主要是由于序列表土均發(fā)育于黃土母質(zhì)，顆粒組成相對均一，且受氣候及成土?xí)r間的制約，整體處于較弱的風化階段，黏土礦物類型以水云母和綠泥石為主。母質(zhì)的同源性及黏土礦物的一致性決定了質(zhì)地對SOC及其組分含量的影響有限。Six等[34]的研究表明，當黏粒含量變化幅度足夠大時，黏粒在預(yù)測SOC方面具有重要作用。而本研究的黏粒含量整體水平較低，且變異程度低(8.1% ～ 16.5%)，因此方差分解結(jié)果顯示質(zhì)地對解釋SOC的方差變異不起作用。該結(jié)果與Percival等[35]在新西蘭草原的研究結(jié)果一致，黏土含量與長期SOC積累的關(guān)系不大，而其他地球化學(xué)特征(如鐵鋁等)具有更高的解釋力。

3.4 不同有機碳組分影響因素差異揭示固碳機制

通過進一步分析SOC及兩種碳組分(POC和MOC)的影響因素發(fā)現(xiàn)，MOC與SOC行為表現(xiàn)較為接近，即各影響要素與SOC、MOC的相關(guān)系數(shù)和偏相關(guān)系數(shù)較接近，而POC與上述二者存在較大差別，這可能指示各要素對SOC的影響主要是通過影響MOC的途徑實現(xiàn)。由于POC是未受到保護或僅僅被保護在大團聚體中，其含量主要取決于生物量的輸入與微生物分解之間的平衡；而MOC經(jīng)微生物加工后具有更多的吸附點位且富有易與金屬氧化物結(jié)合的官能團，這使得MOC與鐵鋁氧化物更容易結(jié)合，形成有機–無機聚合物或被吸附在礦物表面，或被保護在團聚體內(nèi)部[30]，這解釋了土壤礦物對POC的影響要弱于MOC。POC一般穩(wěn)定性較差，易受氣候、生物等要素影響[36]，不易在土壤中長期存在。與POC相比，MOC與鐵鋁氧化物(尤其是Ald、Alo)的相關(guān)性更高，且地形序列土壤中MOC占比超過一半，因此鐵鋁氧化物與MOC的結(jié)合可能是SOC固存的主要機制。氣候?qū)OC水平的直接影響有限，更多是通過影響鐵鋁氧化物等土壤屬性影響MOC水平，進而在長時間尺度影響SOC存量。

4 結(jié)論

通過對祁連山地區(qū)發(fā)育于全新世黃土母質(zhì)的土壤地形序列進行研究，發(fā)現(xiàn)受成土環(huán)境與成土?xí)r間影響，序列內(nèi)表土質(zhì)地及黏土礦物類型保持相對一致，而土壤碳酸鹽、pH、鐵鋁氧化物等土壤理化屬性隨海拔/降水表現(xiàn)出高度相關(guān)性，指示巨大的高程差造成的水熱梯度驅(qū)動成土過程差異。偏相關(guān)分析及方差分解分析揭示，氣候?qū)OC的作用主要通過影響鐵鋁氧化物等土壤屬性介導(dǎo)，且主要通過影響其中的MOC組分影響SOC存量。因此，本研究結(jié)果揭示土壤發(fā)生過程/強度驅(qū)動了祁連山地區(qū)SOC的長期固存。

盡管本研究利用均一母質(zhì)(黃土)與海拔–氣候梯度構(gòu)建的地形序列屬于控制變量的自然試驗，但考慮黃土母質(zhì)在全球分布的廣泛性及組成均一性，本研究對理解不同地區(qū)黃土母質(zhì)發(fā)育土壤SOC長期固存及動態(tài)周轉(zhuǎn)有重要啟示意義，即不同發(fā)育階段的土壤SOC固存機制存在差異，其對氣候因子變化的響應(yīng)很大程度上依賴于SOC的存在形態(tài)。

[1] Schlesinger W H. Evidence from chronosequence studies for a low carbon-storage potential of soils[J]. Nature, 1990, 348(6298): 232–234.

[2] Post W M, Emanuel W R, Zinke P J, et al. Soil carbon pools and world life zones[J]. Nature, 1982, 298(5870): 156–159.

[3] 王根緒, 程國棟, 沈永平. 青藏高原草地土壤有機碳庫及其全球意義[J]. 冰川凍土, 2002, 24(6): 693–700.

[4] Mu C, Zhang T J, Peng X, et al. The organic carbon pool of permafrost regions on the Qinghai–Xizang (Tibetan) Plateau[J]. The Cryosphere Discussions, 2014, 8: 5015–5033.

[5] Johnston C A, Groffman P, Breshears D D, et al. Carbon cycling in soil[J]. Frontiers in Ecology and the Environment, 2004, 2(10): 522–528.

[6] 馬薔, 金會軍. 氣候變暖對多年凍土區(qū)土壤有機碳庫的影響[J]. 冰川凍土, 2020, 42(1): 91–103.

[7] 王蓓, 孫庚, 羅鵬, 等. 模擬升溫和放牧對高寒草甸土壤有機碳氮組分和微生物生物量的影響[J]. 生態(tài)學(xué)報, 2011, 31(6): 1506–1514.

[8] Davidson E A, Janssens I A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change[J]. Nature, 2006, 440(7081): 165–173.

[9] Yang Y H, Fang J Y, Smith P, et al. Changes in topsoil carbon stock in the Tibetan grasslands between the 1980s and 2004[J]. Global Change Biology, 2009, 15(11): 2723–2729.

[10] Jia J, Feng X J, He J S, et al. Comparing microbial carbon sequestration and priming in the subsoil versus topsoil of a Qinghai-Tibetan alpine grassland[J]. Soil Biology and Biochemistry, 2017, 104: 141–151.

[11] 徐滿厚, 薛嫻. 青藏高原高寒草甸夏季植被特征及對模擬增溫的短期響應(yīng)[J]. 生態(tài)學(xué)報, 2013, 33(7): 2071–2083.

[12] 李巖, 干珠扎布, 胡國錚, 等. 增溫對青藏高原高寒草甸呼吸作用的影響[J]. 生態(tài)學(xué)報, 2020, 40(1): 266–273.

[13] 李娜, 王根緒, 高永恒, 等. 青藏高原生態(tài)系統(tǒng)土壤有機碳研究進展[J]. 土壤, 2009, 41(4): 512–519.

[14] K?rner C. The use of ‘a(chǎn)ltitude’ in ecological research[J]. Trends in Ecology & Evolution, 2007, 22(11): 569–574.

[15] 黃昌勇. 土壤學(xué)[M]. 北京: 中國農(nóng)業(yè)出版社, 2000.

[16] Yang F, Huang L M, Rossiter D G, et al. Evolution of loess-derived soil along a climatic toposequence in the Qilian Mountains, NE Tibetan Plateau[J]. European Journal of Soil Science, 2017, 68(3): 270–280.

[17] Torn M S, Trumbore S E, Chadwick O A, et al. Mineral control of soil organic carbon storage and turnover[J]. Nature, 1997, 389(6647): 170–173.

[18] Keiluweit M, Bougoure J J, Nico P S, et al. Mineral protection of soil carbon counteracted by root exudates[J]. Nature Climate Change, 2015, 5(6): 588–595.

[19] Lavallee J M, Soong J L, Cotrufo M F. Conceptualizing soil organic matter into particulate and mineral- associated forms to address global change in the 21st century[J]. Global Change Biology, 2020, 26(1): 261–273.

[20] Lugato E, Lavallee J M, Haddix M L, et al. Different climate sensitivity of particulate and mineral-associated soil organic matter[J]. Nature Geoscience, 2021, 14(5): 295–300.

[21] 張甘霖, 龔子同. 土壤調(diào)查實驗室分析方法[M]. 北京: 科學(xué)出版社, 2012.

[22] Cambardella C A, Elliott E T. Particulate soil organic-matter changes across a grassland cultivation sequence[J]. Soil Science Society of America Journal, 1992, 56(3): 777–783.

[23] Nottebaum V, Lehmkuhl F, Stauch G, et al. Regional grain size variations in aeolian sediments along the transition between Tibetan Highlands and north-western Chinese Deserts - the influence of geomorphological settings on aeolian transport pathways[J]. Earth Surface Processes and Landforms, 2014, 39(14): 1960–1978.

[24] Yang F, Zhang G L, Yang F, et al. Pedogenetic interpretations of particle-size distribution curves for an alpine environment[J]. Geoderma, 2016, 282: 9–15.

[25] Stauch G. Geomorphological and palaeoclimate dynamics recorded by the formation of aeolian archives on the Tibetan Plateau[J]. Earth-Science Reviews, 2015, 150: 393–408.

[26] 龔子同, 張甘霖, 陳志誠, 等. 土壤發(fā)生與系統(tǒng)分類[M]. 北京: 科學(xué)出版社, 2007.

[27] Doetterl S, Stevens A, Six J, et al. Soil carbon storage controlled by interactions between geochemistry and climate[J]. Nature Geoscience, 2015, 8(10): 780–783.

[28] 周萍, 宋國菡, 潘根興, 等. 三種南方典型水稻土長期試驗下有機碳積累機制研究Ⅱ.團聚體內(nèi)有機碳的化學(xué)結(jié)合機制[J]. 土壤學(xué)報, 2009, 46(2): 263–273.

[29] 江家彬, 祝貞科, 林森, 等. 針鐵礦吸附態(tài)和包裹態(tài)有機碳在稻田土壤中的礦化及其激發(fā)效應(yīng)[J]. 土壤學(xué)報, 2021, 58(6): 1530–1539.

[30] 王璐瑩, 秦雷, 呂憲國, 等. 鐵促進土壤有機碳累積作用研究進展[J]. 土壤學(xué)報, 2018, 55(5): 1041–1050.

[31] von Fromm S F, Hoyt A M, Lange M, et al. Continental-scale controls on soil organic carbon across sub-Saharan Africa[J]. Soil, 2021, 7(1): 305–332.

[32] Mu C C, Zhang T J, Zhao Q, et al. Soil organic carbon stabilization by iron in permafrost regions of the Qinghai-Tibet Plateau[J]. Geophysical Research Letters, 2016, 43(19): 10286–10294.

[33] Fang K, Qin S Q, Chen L Y, et al. Al/Fe mineral controls on soil organic carbon stock across Tibetan alpine grasslands[J]. Journal of Geophysical Research: Biogeosciences, 2019, 124(2): 247–259.

[34] Six J, Conant R T, Paul E A, et al. Stabilization mechanisms of soil organic matter: Implications for C-saturation of soils[J]. Plant and Soil, 2002, 241(2): 155–176.

[35] Percival H J, Parfitt R L, Scott N A. Factors controlling soil carbon levels in New Zealand grasslands is clay content important?[J]. Soil Science Society of America Journal, 2000, 64(5): 1623–1630.

[36] von Lützow M, K?gel-Knabner I, Ekschmitt K, et al. SOM fractionation methods: Relevance to functional pools and to stabilization mechanisms[J]. Soil Biology and Biochemistry, 2007, 39(9): 2183–2207.

Soil Organic Carbon Sequestration Driven by Pedogenesis Along Soil Toposequence in Qilian Mountains

GU Jun1,2, YANG Fei1,3, LI Decheng1, SONG Xiaodong1, ZHANG Ganlin1,2,4*

(1 State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China; 2 University of Chinese Academy of Sciences, Beijing 100049, China; 3 Qinghai Province Key Laboratory of Physical Geography and Environmental Process, Qinghai Normal University, Xining 810008, China; 4 Key Laboratory of Watershed Geographic Sciences, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, China)

Storage and dynamics of soil organic carbon (SOC) are the results of input-output balance and soil retention capacity for SOC, which are mainly controlled by bio-climatic conditions, and by soil physiochemical properties such as clay content and soil minerals, respectively. A soil toposequence along an altitudinal gradient was constructed as a natural experiment, aiming to explore the controlling factors for SOC and its fractions in the Qilian Mountains. It is found that soil parent material of this topographic sequence is mainly derived from loess deposits. Chlorite and hydromica are the dominant clay minerals, pointing to the relatively weak chemical weathering intensity of the study region. Meanwhile, the large hydrothermal gradient drives the differences in pedogenic intensities within the topographic sequence. Carbonates exist in soils from the low altitude area; with the rise of altitude, the carbonate is leached away while the contents of iron and aluminum oxide and SOC increase. Further analysis revealed that Al/Fe- oxide was the main explanatory variable for the variance of SOC and its fractions. Partial correlation analysis showed that once the effect of Al/Fe- oxides was controlled, the effects of climate on SOC and its fractions were insignificant. This indicated that the effect of climate on SOC took place mainly via its influence on soil properties, resulting in differences in properties such as Fe and Al oxides. This mechanism mainly occurs in mineral-associated organic carbon (MOC) fraction and make an influence on long-term sequestration of SOC. This study has great implications for understanding the preservation of SOC as well as its responses to climate changes.

Soil toposequence; Soil organic carbon; Al/Fe - oxides; Qilian Mountains

S159.2

A

10.13758/j.cnki.tr.2023.04.021

谷俊, 楊飛, 李德成, 等. 祁連山區(qū)成土過程驅(qū)動的土壤地形序列土壤有機碳固存特征. 土壤, 2023, 55(4): 871–879.

第二次青藏高原綜合科學(xué)考察研究項目(2019QZKK0306-02)和國家自然科學(xué)基金項目(42007001) 資助。

(glzhang@issas.ac.cn)

谷俊 (1995—)，男，安徽馬鞍山人，博士研究生，主要從事土壤有機碳方面研究。E-mail: gujun@issas.ac.cn