孫 濤，賈志清，馬全林，李銀科，張瑩花，王耀林,3

(1.中國(guó)林業(yè)科學(xué)研究院 荒漠化研究所，北京 100091；2.甘肅省治沙研究所/甘肅省荒漠化與風(fēng)沙災(zāi)害防治重點(diǎn)實(shí)驗(yàn)室-省部共建國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地，甘肅 蘭州 730070；3.甘肅省林業(yè)廳 GEF項(xiàng)目辦公室，甘肅 蘭州 730070)

景電灌區(qū)次生鹽堿地枸杞林土壤有機(jī)碳庫(kù)變化的情景模擬

孫 濤1,2，賈志清1，馬全林1，李銀科1，張瑩花1，王耀林1,3

(1.中國(guó)林業(yè)科學(xué)研究院 荒漠化研究所，北京 100091；2.甘肅省治沙研究所/甘肅省荒漠化與風(fēng)沙災(zāi)害防治重點(diǎn)實(shí)驗(yàn)室-省部共建國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地，甘肅 蘭州 730070；3.甘肅省林業(yè)廳 GEF項(xiàng)目辦公室，甘肅 蘭州 730070)

次生鹽堿地；枸杞；土壤有機(jī)碳；CENTURY模型；情景模擬

以景泰紅躍村典型次生鹽堿地枸杞林為研究對(duì)象，基于CENTURY模型模擬成熟枸杞林在實(shí)際田間管理?xiàng)l件和優(yōu)化灌溉模式條件兩種情景下生長(zhǎng)過(guò)程中表層0～20 cm土壤有機(jī)碳庫(kù)的變化，研究結(jié)果表明：次生鹽堿地開墾種植枸杞使得土壤總有機(jī)碳發(fā)生了顯著的變化。在當(dāng)前的田間管理水平情景模擬條件下，成熟枸杞林地土壤總有機(jī)碳逐漸降低，隨著枸杞生長(zhǎng)年限的增加，在未來(lái)20年后土壤有機(jī)碳含量將降低到土地開墾前的水平，未來(lái)將會(huì)變成一個(gè)弱的碳排放源；而在優(yōu)化灌溉模式情景模擬條件中，成熟枸杞林地土壤有機(jī)碳庫(kù)呈先緩慢下降后逐漸回升的狀態(tài)，在未來(lái)20年后土壤有機(jī)碳庫(kù)趨于穩(wěn)定，穩(wěn)定狀態(tài)下表層0～20 cm土壤有機(jī)碳含量相比次生鹽堿地增加約7.57%。因此，改變田間管理方式、優(yōu)化灌水模式對(duì)提高枸杞林土壤質(zhì)量、增強(qiáng)碳匯能力和減緩溫室效應(yīng)具有重要的現(xiàn)實(shí)意義。

土壤有機(jī)質(zhì)模型能夠模擬土壤有機(jī)質(zhì)的變化規(guī)律，已被廣泛用于模擬土壤有機(jī)質(zhì)的動(dòng)態(tài)變化。土壤有機(jī)質(zhì)模擬已經(jīng)在當(dāng)前大量的有關(guān)氣候變化和全球碳循環(huán)研究中起到了重要的作用[1-2]，成為更好地理解和管理陸地碳循環(huán)的一個(gè)重要組成部分。CENTURY模型是應(yīng)用比較廣泛的表征土壤有機(jī)質(zhì)動(dòng)態(tài)的模型，是美國(guó)科羅拉多州立大學(xué)的Parton等于20世紀(jì)80年代末建立的，起初用于模擬草地生態(tài)系統(tǒng)的碳、氮、磷、硫等元素的長(zhǎng)期演變過(guò)程，以后加以改進(jìn)，其應(yīng)用擴(kuò)展到稀樹草原、森林、農(nóng)業(yè)等生態(tài)系統(tǒng)中。它主要基于土壤的結(jié)構(gòu)功能，從碳、氮、磷的生物地球化學(xué)循環(huán)因子入手，結(jié)合溫度、降水等驅(qū)動(dòng)因子，對(duì)生態(tài)系統(tǒng)的生產(chǎn)力進(jìn)行模擬和預(yù)測(cè)[3-4]。國(guó)內(nèi)外眾多學(xué)者先后將CENTURY模型應(yīng)用于農(nóng)業(yè)或林業(yè)生態(tài)系統(tǒng)中，并對(duì)土壤有機(jī)碳進(jìn)行了模擬，其模擬值與實(shí)測(cè)值具有極高的相關(guān)性，并可用于不同的農(nóng)林業(yè)生態(tài)系統(tǒng)的模擬研究中[5-9]。

引水灌溉對(duì)荒漠、半荒漠區(qū)域的農(nóng)業(yè)發(fā)展起到了重要的推動(dòng)作用，但是在不適當(dāng)?shù)墓喔扰潘畻l件下造成了灌區(qū)土壤次生鹽堿化等生態(tài)退化問(wèn)題[10]。甘肅景泰川灌區(qū)多年來(lái)不合理的大水漫灌、串灌等使得景泰川這個(gè)半封閉型盆地的地下水位迅速上升，導(dǎo)致土地次生鹽堿化面積逐年增大，土地退化問(wèn)題尤為嚴(yán)重[11]，給農(nóng)業(yè)生產(chǎn)和人民生活造成了嚴(yán)重影響。

從20世紀(jì)90年代中期開始，當(dāng)?shù)馗鶕?jù)實(shí)際情況采取生物措施對(duì)次生鹽堿地進(jìn)行治理，引進(jìn)枸杞并種植成功，現(xiàn)在已經(jīng)在景泰川地區(qū)大面積地推廣種植開來(lái)，不光產(chǎn)生了經(jīng)濟(jì)效益，也充分利用和改良了鹽堿地。在枸杞種植過(guò)程中土壤性質(zhì)會(huì)發(fā)生變化，對(duì)土壤質(zhì)地產(chǎn)生一定的影響，根系生物量和微生物的增加，會(huì)提高土壤碳庫(kù)活度。隨著生長(zhǎng)年限的增加和農(nóng)田耕作措施的變化，土壤有機(jī)碳庫(kù)也發(fā)生相應(yīng)的變化。當(dāng)前，枸杞種植戶為了獲得較高的枸杞產(chǎn)量和經(jīng)濟(jì)效益，增加灌溉次數(shù)、采用大水漫灌等粗放的管理方式在灌區(qū)非常普遍。由于景泰川地區(qū)是一個(gè)淺層半封閉型的盆地，這種粗放的農(nóng)田管理方式直接導(dǎo)致局部地區(qū)地下水位進(jìn)一步抬升，地下水不能及時(shí)排出，使得枸杞林土壤結(jié)構(gòu)改變，土質(zhì)下降，根腐病較為嚴(yán)重，病死枯死植株也較多。在當(dāng)前這種灌溉方式和田間管理模式下，水鹽動(dòng)態(tài)被進(jìn)一步打破，土壤碳庫(kù)的變化特征將如何？如果改變當(dāng)前不合理的灌溉方式，采用合理的管理模式，在未來(lái)的管理過(guò)程中其特征變化和演變趨勢(shì)如何？是否能提高土壤的碳庫(kù)效益和增加固碳能力？運(yùn)用模型擬合研究在不同的情景模式下土壤碳庫(kù)的變化規(guī)律，并對(duì)其變化趨勢(shì)做出分析和預(yù)測(cè)是一個(gè)較好的研究方法。因此，利用CENTURY模型模擬在次生鹽漬化改良過(guò)程中土壤有機(jī)碳的變化特征，設(shè)置兩種不同的情景模式，研究成熟枸杞林對(duì)次生鹽堿地的改良作用和土壤碳的動(dòng)態(tài)變化，進(jìn)而預(yù)測(cè)在不同情景模式下未來(lái)土壤碳庫(kù)的演變趨勢(shì)，闡明土壤碳庫(kù)的變化特征和差異，可為次生鹽堿地土壤質(zhì)量和健康評(píng)價(jià)提供基礎(chǔ)數(shù)據(jù)，也可為當(dāng)前枸杞林地田間管理措施的科學(xué)評(píng)估提供科學(xué)依據(jù)。

1 研究地區(qū)與研究方法

1.1 自然概況

紅躍村地處景泰縣中部，景電一期引黃提灌灌區(qū)中心地帶，位于景泰縣草窩灘鎮(zhèn)以東3 km，縣城以北8 km處，是由西南向東北傾斜的微型盆地，海拔1 565 m。氣候干燥、干旱少雨、風(fēng)沙大，年均降水量185 mm，降水年際變化不大，雨水集中在7—9月，年蒸發(fā)量3 040 mm，年均氣溫8.5 ℃，雨熱同期，地表水資源較為貧乏，地下水礦化度高、水質(zhì)差，20世紀(jì)70—80年代建成了景電一、二期高揚(yáng)程電力引黃提灌工程[12]。全村現(xiàn)有耕地面積330.7 hm2，約80%的耕地為次生鹽堿化土地。自2000年以來(lái)，該村開始推廣引種耐鹽堿的經(jīng)濟(jì)樹種枸杞，取得了良好的經(jīng)濟(jì)和生態(tài)效益[13]。本研究地區(qū)土壤母質(zhì)類型一致，地帶性土壤為荒漠灰鈣土，土壤母質(zhì)包括黃土、風(fēng)積物和洪積物。地下水位均為1 m左右。次生鹽堿地主要植被為鹽爪爪〔Kalidiumfoliatum(Pall.) Moq.〕，蓋度8%～15%。

1.2 研究方法與數(shù)據(jù)采集

應(yīng)用模型前要經(jīng)過(guò)校正與驗(yàn)證，調(diào)整模型參數(shù)，反復(fù)檢驗(yàn)?zāi)M結(jié)果，檢驗(yàn)?zāi)M值與觀測(cè)值之間的吻合程度，使之達(dá)到最好的模擬效果?；谡{(diào)試后的參數(shù)，再通過(guò)對(duì)比分析其他年份的模擬與觀測(cè)結(jié)果進(jìn)行模型的模擬研究與驗(yàn)證平衡[1,14]。

本研究于2012年7月中旬開始進(jìn)行野外調(diào)查、土樣采集，并對(duì)農(nóng)戶進(jìn)行訪談式調(diào)查，了解枸杞種植和生長(zhǎng)情況。采用空間代替時(shí)間的方法，分別選擇4、7、11年的枸杞林地采集土樣，選擇臨近的次生鹽堿化土地作為對(duì)照。土壤采樣深度40 cm，每10 cm分層采樣，每塊樣地5個(gè)重復(fù)，同層次土樣混合作為1個(gè)混合樣，帶回實(shí)驗(yàn)室風(fēng)干后測(cè)定土壤有機(jī)質(zhì)、有機(jī)碳、pH值、鹽分等理化指標(biāo)。同時(shí)，用環(huán)刀法取樣測(cè)土壤容重，鋁盒裝少量新鮮土壤用于土壤水分測(cè)定。土壤粒度采用激光粒度分析儀測(cè)定，總有機(jī)碳用重鉻酸鉀氧化-外加熱法測(cè)定，活性有機(jī)碳用高錳酸鉀氧化法測(cè)定[15]，全鹽量測(cè)定水土比為5 ∶1，pH值用酸度計(jì)法(水土比為2.5 ∶1)測(cè)定[16]。經(jīng)測(cè)定，研究區(qū)次生鹽堿裸地的土壤有機(jī)質(zhì)、全氮、全磷、全鉀含量分別為10.90、0.76、0.95、22.85 g/kg，速效磷、速效鉀含量分別為3.04、212 mg/kg。

1.3 CENTURY模型的參數(shù)化設(shè)置

根據(jù)土壤有機(jī)質(zhì)的分解速率，CENTURY模型將土壤總有機(jī)碳(TOC)分成了3個(gè)碳庫(kù)，分別是活性、緩性和惰性土壤有機(jī)碳庫(kù)[1,3]。本研究所用的是CENTURY4.0版本，初始輸入數(shù)據(jù)包括基本的氣象數(shù)據(jù)、土壤質(zhì)地等，見表1。

表1 模型初始化所需主要參數(shù)

模型所運(yùn)用的氣象數(shù)據(jù)由景泰縣氣象站提供的1955—2012年氣象數(shù)據(jù)統(tǒng)計(jì)分析而得，由于研究區(qū)距離縣城僅8 km，地勢(shì)平坦，因此可以將此數(shù)據(jù)直接作為研究區(qū)的數(shù)據(jù)進(jìn)行模擬運(yùn)用。模型要求至少10年以上的氣象數(shù)據(jù)，所需資料包括：月均降水量(PREC)、月最低溫度(TMIN)、月最高溫度(TMAX)。CENTURY模型在運(yùn)行調(diào)試過(guò)程中首先要設(shè)定主要的植被生長(zhǎng)參數(shù)，包括植被生長(zhǎng)的開始月份、 衰老月份、休眠月份。本研究采用模型中自帶的GRASSLAND/CROP模塊，其適用于北方典型草地、低矮草原和荒漠草地。選擇CROP.100文件中的CPR為模型的初始植被，4月開始生長(zhǎng)，10月為衰老季節(jié)，11月份開始休眠，時(shí)間步長(zhǎng)為一年，進(jìn)行初始態(tài)的建立。部分參數(shù)應(yīng)用模型提供的缺省值。

據(jù)調(diào)查，當(dāng)?shù)卮紊}堿化土地植被以鹽爪爪為主，蓋度8%～15%，2000年開墾種植枸杞。枸杞在不同生長(zhǎng)階段田間管理措施有所不同：枸杞生長(zhǎng)初期(相當(dāng)于2000—2004年)，主要保證枸杞的正常生長(zhǎng)和成林，沒有修剪枝條，沒有果實(shí)產(chǎn)生，因而施肥澆水次數(shù)較少，均為每年3～4次；枸杞林成熟穩(wěn)產(chǎn)階段(相當(dāng)于2005—2013年)，最大特點(diǎn)是有果實(shí)產(chǎn)出和植被的移除，每年春季3、4月修剪枝條，近1/3枝條被移除，從6月份開始持續(xù)4～5個(gè)月采摘果實(shí)，為了獲得穩(wěn)定的果實(shí)產(chǎn)量，施肥和澆水次數(shù)增至每年5～6次。

2 結(jié)果與分析

2.1 模型參數(shù)初始化設(shè)置及穩(wěn)態(tài)的建立

用CENTURY模型模擬研究土壤變化時(shí)首先要進(jìn)行模型參數(shù)的初始化，在參數(shù)平衡狀態(tài)下才能進(jìn)行真正的模擬運(yùn)算[1-2]。主要是模擬研究區(qū)在自然狀態(tài)下土壤活性碳庫(kù)、緩性碳庫(kù)、惰性碳庫(kù)和土壤總碳庫(kù)(包含結(jié)構(gòu)性碳庫(kù)和代謝性碳庫(kù))逐漸積累到穩(wěn)定平衡階段的過(guò)程。如圖1所示土壤碳庫(kù)快速積累變化所需時(shí)間約為2 000 a，之后緩慢積累到2 400 a時(shí)，土壤各碳庫(kù)逐步建立平衡，各碳庫(kù)變化穩(wěn)定，基本在土壤碳含量平均水平上下波動(dòng)(圖中水平虛線所示)，達(dá)到了試驗(yàn)區(qū)土壤碳水平，此時(shí)活性、緩性和惰性有機(jī)碳庫(kù)有機(jī)碳含量分別占總有機(jī)碳庫(kù)的3.28%、68.67%和28.05%。從圖1中土壤有機(jī)碳含量的年際變化可看出：土壤總碳庫(kù)和緩性碳庫(kù)有機(jī)碳含量年際變化較大，但在一定的水平上波動(dòng)；惰性有機(jī)碳庫(kù)變化穩(wěn)定、曲線平滑；活性有機(jī)碳庫(kù)有機(jī)碳含量年際變化敏感，但是其值較小，所占比例小。研究地為次生鹽堿地，其土壤有機(jī)碳含量為3 300 g/m2左右，而CENTURY模型在平衡狀態(tài)下模擬的土壤總有機(jī)碳含量在3 300～3 400 g/m2間波動(dòng)，兩者相差很小，說(shuō)明CENTURY模型正確模擬了研究區(qū)土壤有機(jī)碳的動(dòng)態(tài)積累過(guò)程。

圖1 CENTURY模型參數(shù)初始化過(guò)程

2.2 不同生長(zhǎng)階段枸杞林對(duì)次生鹽漬化土壤有機(jī)碳庫(kù)的影響

研究區(qū)自2000年開始開墾次生鹽堿地種植枸杞，土地類型發(fā)生了變化，自然植被轉(zhuǎn)變成了人工栽植的經(jīng)濟(jì)作物，土壤有機(jī)碳隨之發(fā)生了顯著的變化(圖2)。從次生鹽堿裸地到開墾種植枸杞，再到枸杞成熟穩(wěn)產(chǎn)，在這過(guò)程中因?yàn)橹脖活愋筒煌?、生長(zhǎng)年限逐漸增加，故所對(duì)應(yīng)的措施也不同。在CENTURY模型模擬過(guò)程中，不同的生長(zhǎng)階段對(duì)應(yīng)不同的管理措施。由于每個(gè)階段的土地耕作方式和田間管理模式不同，因此土壤活性、緩性、惰性和土壤總碳庫(kù)也隨著生長(zhǎng)年限的增加發(fā)生了顯著的變化(圖2)。

圖2 枸杞生長(zhǎng)不同階段土壤各個(gè)碳庫(kù)隨時(shí)間的變化

初期的1—4年為枸杞幼苗生長(zhǎng)期，也是各個(gè)碳庫(kù)快速增加時(shí)期，如圖2中B時(shí)間段。土壤活性、緩性、惰性和總有機(jī)碳庫(kù)在這一時(shí)間段內(nèi)均呈快速增加的態(tài)勢(shì)，土壤活性碳含量增幅最大，相比次生鹽堿裸地增加了137%，其次是土壤總碳、 土壤緩性碳、土壤惰性碳含量，其增幅分別為22.59%、7.69%和0.33%。4年后為枸杞林成熟穩(wěn)產(chǎn)期，各個(gè)碳庫(kù)變化逐漸穩(wěn)定，如圖2中C時(shí)間段。緩性碳和惰性碳在這個(gè)階段繼續(xù)積累，但是累積速率相比枸杞生長(zhǎng)初期低，因此變化相對(duì)緩慢，有滯后效應(yīng)，在枸杞生長(zhǎng)到第9—10年(相當(dāng)于2009—2010年)達(dá)到最高值；活性碳變化較為敏感，下降幅度大，2008年以后下降速度趨緩，逐漸穩(wěn)定，其值低于開墾前活性碳的含量；土壤總有機(jī)碳含量隨著活性有機(jī)碳含量的快速降低而逐漸減小，從開始種植時(shí)隨著枸杞生長(zhǎng)年限的增加而增加變?yōu)橹饾u降低，呈倒“V”形變化。造成總有機(jī)碳逐漸降低的原因主要是春季修剪枝條，近1/3的枝條被剪除，同時(shí)枸杞果實(shí)也被大量采摘，這些生物量均從枸杞林移除，導(dǎo)致枯落物碳儲(chǔ)量逐漸下降，從而也導(dǎo)致土壤總碳庫(kù)含量的變化；另一個(gè)原因是植株衰老，更新能力減弱。

2.3 枸杞林實(shí)測(cè)值與模擬值相關(guān)性

以生長(zhǎng)年限為0(次生鹽堿裸地)、4、7、11年的枸杞林土壤有機(jī)碳實(shí)測(cè)數(shù)據(jù)與模型模擬數(shù)據(jù)進(jìn)行數(shù)據(jù)擬合(結(jié)果如圖3)，可知土壤實(shí)測(cè)值與模擬值具有較好的相關(guān)性，變化趨勢(shì)一致，R2達(dá)到了0.947，統(tǒng)計(jì)分析P=0.27>0.05，表明模型值與實(shí)測(cè)值之間差異不顯著，說(shuō)明運(yùn)用CENTURY模型來(lái)模擬次生鹽漬化土地枸杞林有機(jī)碳的動(dòng)態(tài)變化基本符合實(shí)際情況[17]，可用來(lái)模擬研究試驗(yàn)區(qū)的土壤動(dòng)態(tài)變化。

圖3 枸杞林土壤有機(jī)碳實(shí)測(cè)值與模擬值的關(guān)系

2.4 未來(lái)20年枸杞林土壤有機(jī)碳庫(kù)情景模擬

以2010年CENTURY模型輸出結(jié)果為土壤有機(jī)碳的初始值，氣象參數(shù)運(yùn)用模型自動(dòng)統(tǒng)計(jì)分析1955—2012年的歷史氣象資料得到的氣象數(shù)據(jù)，輸出時(shí)間步長(zhǎng)為月份，分析預(yù)測(cè)成熟枸杞林在不同的灌溉次數(shù)和用水量條件下土壤有機(jī)碳未來(lái)20年的模擬變化趨勢(shì)。設(shè)計(jì)兩種情景進(jìn)行模擬：一種情景表示在當(dāng)前實(shí)際田間管理方式下的土壤碳變化(曲線A50)，即為了讓成熟枸杞林獲得穩(wěn)定產(chǎn)量，枸杞林地每年澆水5～6次，澆水方式為大水漫灌；另一種情景為優(yōu)化灌溉模式狀態(tài)下土壤碳的變化(曲線A25)，即減少灌溉次數(shù)為3次，用水量為田間持水量的25%，其他措施均不變。相應(yīng)的在CENTURY模型IRRI.100文件中輸入對(duì)應(yīng)參數(shù)進(jìn)行兩種情景下的模擬，結(jié)果見圖4。由圖4可知，隨著植株生長(zhǎng)年限的逐漸增加，在未來(lái)15—20年間，土壤有機(jī)總碳庫(kù)均趨于穩(wěn)定，但其含量不同：在當(dāng)前灌水條件下(A50)枸杞林的土壤有機(jī)碳庫(kù)呈逐漸下降趨勢(shì)，未來(lái)20年后表層0～20 cm土壤有機(jī)碳含量平均約3 200 g/m2，與枸杞種植前的次生鹽堿裸地含量相當(dāng)；而在A25情景模擬條件下，未來(lái)20年后土壤有機(jī)碳含量約為3 550 g/m2，相對(duì)枸杞種植前的次生鹽堿地土壤有機(jī)碳含量提高了約7.57%?；钚?、緩性和惰性碳庫(kù)均隨之發(fā)生變化，活性碳庫(kù)在兩種情景模擬下均逐漸降低，但是在A25情景下降低速度要小于A50，而緩性和惰性碳庫(kù)在當(dāng)前的灌溉管理方式下均逐漸降低，在減少灌溉的優(yōu)化模式下均逐漸增加。各個(gè)碳庫(kù)年內(nèi)變化波動(dòng)小，主要是受季節(jié)性的氣溫、降水等變化影響，但是對(duì)整體變化趨勢(shì)影響不大。

圖4 成熟枸杞林未來(lái)土壤各碳庫(kù)變化趨勢(shì)

3 討論與結(jié)論

研究區(qū)枸杞林種植前是以鹽爪爪為主的次生鹽堿地，經(jīng)過(guò)大規(guī)模的開墾、施肥、排鹽等措施建成了枸杞林地，種植枸杞在產(chǎn)生經(jīng)濟(jì)效益的同時(shí)又能改良土壤，起到固碳的作用。枸杞林對(duì)土壤有機(jī)碳的影響有兩個(gè)階段，即幼苗生長(zhǎng)階段和成熟穩(wěn)定階段。開墾導(dǎo)致土壤溫度和濕度的改變，將極大地促進(jìn)土壤的呼吸作用，加速表層土壤有機(jī)質(zhì)的分解[18-19]，因此在枸杞生長(zhǎng)初期(前1—4年)土壤活性、緩性、惰性碳含量均呈快速增加的態(tài)勢(shì)。隨著枸杞幼苗地上和地下生物量逐漸增加，土壤碳的累積速率也在栽植4—5年間達(dá)到最大，并超過(guò)自然狀態(tài)下的總碳庫(kù)，2004年枸杞林土壤總碳庫(kù)的碳含量比開墾前增加了22.59%，表明這段時(shí)間枸杞林地土壤碳總體呈碳匯趨勢(shì)。但是在成熟枸杞林內(nèi)田間管理措施發(fā)生了變化，因而土壤碳庫(kù)也發(fā)生了變化。地上生物量的剪除和采摘果實(shí)使返回地表的枯落物減少、腐殖質(zhì)含量降低，活性碳庫(kù)隨之減小[20]。灌水次數(shù)的增加和水量的無(wú)序利用，使得土壤結(jié)構(gòu)和性狀被逐漸改變，土壤緩性碳和惰性碳在緩慢積累之后迅速降低，土壤總碳庫(kù)隨之降低，土壤的碳匯功能也被大大降低或被逆轉(zhuǎn)[21-22]。

對(duì)成熟枸杞林土壤有機(jī)碳含量的變化進(jìn)行預(yù)測(cè)分析表明，在當(dāng)前的農(nóng)田管理措施下枸杞林地土壤總有機(jī)碳含量逐漸地降低(圖4)，其含量約為3 200 g/m2，與次生鹽堿裸地含量相當(dāng)，這表明土壤有機(jī)碳庫(kù)在未來(lái)20年內(nèi)極有可能逐漸轉(zhuǎn)變?yōu)橐粋€(gè)弱的碳排放源。研究區(qū)地形為一個(gè)半封閉的小盆地，大量灌溉用水滲入地下不能及時(shí)排出，使水位上升，同時(shí)，在實(shí)際調(diào)查中發(fā)現(xiàn)附近的排堿渠年久失修，功能喪失，未能發(fā)揮排鹽排堿的功能，導(dǎo)致該區(qū)域次生鹽堿化現(xiàn)象進(jìn)一步加重，地下水位在0.8～1 m間，反鹽現(xiàn)象明顯，灌溉水不能快速下滲，土壤結(jié)構(gòu)和性狀發(fā)生變化，土壤菌類數(shù)量及微生物總量顯著降低[23]。在這種土壤狀態(tài)下植被自然生長(zhǎng)不好，土壤有機(jī)碳含量減少，土壤質(zhì)量呈惡化的趨勢(shì)[24]，這與情景模擬A50狀態(tài)下土壤有機(jī)碳含量逐漸降低的變化趨勢(shì)相似。這表明當(dāng)前對(duì)枸杞林采取的田間管理措施是不合理的，這種耕作方式導(dǎo)致了土壤肥力下降，土壤質(zhì)量退化，進(jìn)而決定了土壤有機(jī)碳含量逐漸降低；而轉(zhuǎn)變灌溉方式，減少灌溉次數(shù)，提高用水利用率可以增加土壤的有機(jī)碳含量，提高固碳效益，增強(qiáng)碳匯能力。因此，改變田間管理方式、采用科學(xué)灌水設(shè)備、及時(shí)更新苗木等對(duì)提高枸杞林土壤質(zhì)量、增加土壤碳儲(chǔ)量和減緩溫室效應(yīng)具有重要的現(xiàn)實(shí)意義。

[1] 高魯鵬,梁文舉,姜勇,等.利用CENTURY模型研究東北黑土有機(jī)碳的動(dòng)態(tài)變化Ⅰ.自然狀態(tài)下土壤有機(jī)碳的積累[J].應(yīng)用生態(tài)學(xué)報(bào),2004,15(5):772-775.

[2] Paustian K.Modelling soil organic matter dynamics-global challenges[G]//Rees R M,Ball B C,Campbell C D,et al.Sustainable Management of Soil Organic Matter.Wallingford,Oxon,UK:CAB International,2001:43-53.

[3] Parton W J,McKeown B,Kirchner V,et al.CENTURY User’s Manual[M].Fort Collins:Natural Resource Ecology Laboratory,Colorado State University,1992：10-50.

[4] Parton W J,Schimel D S,Ojima D S.Observation and modeling of biomass and soil organic matter dynamics for the grassland biome worldwide[J] .Global Biogeochemical Cycles,1993,7(4):785-809.

[5] 呂新苗,鄭度.氣候變化對(duì)長(zhǎng)江源地區(qū)高寒草甸生態(tài)系統(tǒng)的影響[J].長(zhǎng)江流域資源與環(huán)境,2006,15(5):603-607.

[6] 李東.基于CENTURY模型的高寒草甸土壤有機(jī)碳動(dòng)態(tài)模擬研究[D].南京:南京農(nóng)業(yè)大學(xué),2011：20-30.

[7] Mikhailova E A,Bryant R B,DeGloria S D,et al.Modeling soil organic matter dynamics after conversion of native grassland to long-term continuous fallow using the CENTURY model[J].Ecological Modeling,2000,132(3):247-257.

[8] Miko U F K,Keryn I P.Modelling C and N dynamics in forest soils with a modified version of the CENTURY model[J].Soil Biology and Biochemistry,2002,34(3):341-354.

[9] Burker I C,Youker C M,Parton W J.Texture,climate,and cultivation effects on soil organic matter content in US grassland soils[J].Soil Society of America Journal,1989,53(3):800-805.

[10] 邢大韋.西北地區(qū)灌溉的生態(tài)環(huán)境及對(duì)策研究[J].生態(tài)農(nóng)業(yè)研究,1996,4(3):22-26.

[11] 高志海.引黃灌區(qū)農(nóng)田鹽堿化與沙漠化治理的綜合途徑——以甘肅省景泰川為例[J].農(nóng)業(yè)環(huán)境與發(fā)展,2008,25(5):9-12.

[12] 徐向宏,王婷婷,雷文文.景電灌區(qū)綠洲綜合生態(tài)系統(tǒng)管理耦合技術(shù)評(píng)價(jià)分析[J].草業(yè)科學(xué),2010,27(6):38-44.

[13] 孫濤,王耀琳,馬全林,等.景電灌區(qū)次生鹽堿地枸杞林經(jīng)濟(jì)效益分析[J].草業(yè)科學(xué),2013,30(9):1454-1458.

[14] 高崇升,楊國(guó)亭,王建國(guó),等.CENTURY模型在農(nóng)田生態(tài)系統(tǒng)中的應(yīng)用及其參數(shù)確定[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究,2006,22(1):50-52.

[15] 蔡太義,黃耀威,黃會(huì)娟,等.不同年限免耕秸稈覆蓋對(duì)土壤活性有機(jī)碳和碳庫(kù)管理指數(shù)的影響[J].生態(tài)學(xué)雜志,2011,30(9):1962-1968.

[16] 王勇輝,馬蓓,海米提·依米提.艾比湖主要補(bǔ)給河流下游河岸帶土壤鹽分特征[J].干旱區(qū)研究,2013,30(2):196-202.

[17] 高崇升,楊國(guó)亭,王建國(guó),等.移耕農(nóng)業(yè)制度下黑土農(nóng)田土壤有機(jī)碳的演變[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究,2008,24(1):88-92.

[18] 戴全厚,劉國(guó)彬,薛萐,等.侵蝕環(huán)境退耕撂荒地土壤活性有機(jī)碳與碳庫(kù)管理指數(shù)演變[J].西北林學(xué)院學(xué)報(bào),2008,23 (6):24-28.

[19] John B,Yamashita T,Ludwig B,et al.Storage of organic carbon in aggregate and density fractions of silty soils under different types of land use[J].Geoderma,2005,128(1-2):63-79.

[20] 許文強(qiáng),陳曦,羅格平,等.基于CENTURY模型研究干旱區(qū)人工綠洲開發(fā)與管理模式變化對(duì)土壤碳動(dòng)態(tài)的影響[J].生態(tài)學(xué)報(bào),2010,30(14):3707-3716.

[21] 李月梅,曹廣民,徐仁海.土地利用方式對(duì)高寒草甸土壤有機(jī)碳及其組分的影響[J].安徽農(nóng)業(yè)科學(xué),2008,36(14):5951-5953.

[22] 李月梅,曹廣民,王躍思.開墾對(duì)海北高寒草甸土壤有機(jī)碳的影響[J].生態(tài)學(xué)雜志,2006,25(8):911-914.

[23] 李鳳霞,王學(xué)琴,郭永忠,等.寧夏不同類型鹽漬化土壤微生物區(qū)系及多樣性[J].水土保持學(xué)報(bào),2011,25(5):107-111.

[24] 陳辰,王靖,潘學(xué)標(biāo),等.CENTURY模型在內(nèi)蒙古草地生態(tài)系統(tǒng)的適用性評(píng)價(jià)[J].草地學(xué)報(bào),2012,20(6):1012-1019.

(責(zé)任編輯 徐素霞)

“十二五”農(nóng)村領(lǐng)域國(guó)家科技計(jì)劃項(xiàng)目(2012BAD16B0102)；國(guó)家自然科學(xué)基金項(xiàng)目(31300595)；甘肅省青年科學(xué)基金項(xiàng)目(1208RJYA093，1208RJYA067)；全球環(huán)境基金項(xiàng)目(GEF/UNEP)；碳效益研究項(xiàng)目(GEF/53-4280)

S157

A

1000-0941(2015)09-0059-05

孫濤(1978—)，男，甘肅永昌縣人，助理研究員，博士研究生，主要從事荒漠生態(tài)、荒漠植被恢復(fù)及荒漠化防治研究工作；通信作者賈志清(1968—)，女，北京市人，研究員，博士生導(dǎo)師，博士，主要從事荒漠化防治與水土保持研究。

2015-02-15