培肥措施對旱地農田土壤CO2排放和碳庫管理指數的影響

王曉嬌，蔡立群，齊鵬，王雅芝，陳曉龍，武均，2，4，張仁陟

（1.甘肅農業大學資源與環境學院，甘肅 蘭州730070；2.甘肅農業大學甘肅省干旱生境作物學重點實驗室，甘肅 蘭州730070；3.甘肅農業大學管理學院，甘肅 蘭州730070；4.甘肅省節水農業工程技術研究中心，甘肅 蘭州730070）

氮肥在農業發展中扮演著重要的角色，其使全球糧食產量增加了40%～50%，然而，氮肥過量施用對土壤造成的損害與日俱增，損失土壤有機碳、引起土壤板結化、加劇“溫室效應”、促進土壤面源污染［1］。在此背景下，我國提出化肥“零增長計劃”、“化肥減量、有機肥替代”等一系列政策，有機物料替代化肥的土壤培肥措施逐步開始在農田大面積應用，其目的是降低化肥施用量，減少環境問題。農田生態系統土壤CO2排放作為陸地生態系統土壤CO2排放的重要分支之一［2］，對陸地生態系統乃至全球碳排放的影響巨大［3］，預計到21 世紀末，全球平均地表溫度將升高1～4 ℃［4?5］。農田土壤CO2排放主要由土壤有機質的微生物分解驅動，有機物料培肥將攜帶大量碳輸入農田土壤提升土壤有機質［6］，從而影響土壤CO2排放。因此，在提高產量的基礎上，分析不同培肥模式下農田土壤CO2排放、作物碳排放效率和土壤碳庫管理指數成為評估環境友好培肥模式的重要依據。

已有研究表明，不同培肥措施下土壤CO2排放的差異很大，利用田間試驗［7?8］、模型［9］、meta 分析［10?11］等方式研究得出有機物料通過向土壤輸入外源碳提高土壤有機碳的積累和增強土壤微生物的活性刺激增加土壤CO2排放量，也有部分研究表明農田土壤CO2排放量在常規施肥、優化施肥和有機肥配合優化施肥等處理間差異不顯著［12］。有機物料投入與單施氮肥相比對產量的影響結果存在爭議，部分學者認為等氮添加下有機物料施肥與單施氮肥間無差異［13］，也有學者的研究表明有機物料施肥后產量會下降，從而也導致了碳排放效率的差異［14?15］。土壤碳庫管理指數是評價土壤質量和管理措施的重要指標［16］，現有在灰漠土、黑土、褐土、紅壤和棕壤等的研究認為含有機物料的培肥措施對土壤活性有機碳和碳庫管理指數、微生物特征都有提高作用，其中有機無機配施對土壤有機質的提高較大［13］，但研究大多集中在小麥（Triticum aestivum）、水稻（Oryza sativa）農田，多以短期試驗較多，對旱土高原玉米（Zea mays）農田的中長期定位試驗研究較少［17?19］；大多數值模擬和短期田間試驗的研究結果表明土壤CO2排放通量受溫度、水分、施肥措施和有機物料類型等因素影響［20］，在與土壤活性有機碳、脲酶、蔗糖酶和微生物數量的關系研究方面存在一定爭議［16，21?23］，陳哲［24］的研究認為土壤CO2排放通量與活性有機碳無相關關系，張文麗［25］的研究認為其與活性有機碳、微生物碳氮比有一定的關系。在黃土高原地區，外源有機物料的添加對土壤CO2排放量的影響因子研究大多停留在溫度和水分［26?29］，運用結構方程方法使相關的土壤碳庫、理化性質、微生物特征因子整體來解釋CO2排放機制缺乏文獻支持［20］。

甘肅隴中黃土高原旱作區環境惡劣，發展旱作糧飼兼用玉米是緩解糧食供需矛盾，解決資源性缺水雙重壓力的必然選擇［30］，選擇增產減排的旱作玉米培肥模式尤為重要。本研究以甘肅隴中黃土高原旱作玉米農田長期定位試驗（2012 年布設）為研究對象，研究旱作玉米農田土壤碳排放特征、產量、土壤理化指標及微生物特征對不同培肥措施的響應及土壤碳排放的影響機制。旨在明確：1）不同培肥措施下土壤CO2排放特征、產量、作物碳排放效率和碳庫管理指數，尋找增產減排的培肥模式；2）不同培肥措施下水分、溫度、有機碳、微生物量碳氮、脲酶、蔗糖酶等因子的變化特征；3）環境因子對土壤CO2排放通量的影響機制。以期為黃土高原旱作農業區選擇提升作物產量和環境友好的培肥模式以及為農田生態系統的碳減排提供科學依據。

1 材料與方法

1.1 試驗區概況

試驗地點設于甘肅省定西市安定區李家堡鎮麻子川村的甘肅農業大學旱作農業綜合試驗站，屬黃土高原西部甘肅隴中半干旱丘陵溝壑區。該區是典型的旱作雨養農業區，海拔2000 m 左右，日照時數2476.6 h，年 均 太 陽 輻 射594.7 kJ·cm?2，年 均 氣 溫6.4 ℃，≥10 ℃年積溫為2239.1 ℃，年均無霜期140 d，多年平均降水量390.9 mm，2018 年降水量472.05 mm（圖1），年蒸發量達到1531.0 mm，干燥度2.53；土壤為典型的黃綿土，土質綿軟，土層深厚，質地均勻，儲水性能良好，耕層肥力低（表1）。

圖1 2018 年研究區月降水量分布Fig.1 Monthly precipitation distribution of study area in 2018

表1 試驗地土壤主要理化性質Table 1 Soil physical and chemical properties of experimental field

1.2 試驗設計

試驗采用單因素完全隨機設計，按照當地常規氮肥用量（N 200 kg·hm?2）為標準，采用等氮投入的方法設置5個處理：空白對照（CK）、氮肥（nitrogen fertilizer，NF，N 200 kg·hm?2）、有機肥（organic manure，OM，6000 kg·hm?2）、玉米秸稈（straw，ST，28500 kg·hm?2）、有機肥結合氮肥（organic manure combined with nitrogen fertilizer，OMNF，N 100 kg·hm?2+3000 kg·hm?2），3 次重復。氮肥為尿素（含氮量46%），有機肥、秸稈均采用等氮投入的方法，通過其含氮量（有機肥3.3%、秸稈0.7%）折算投入量，磷肥為過磷酸鈣，含P2O516%，氮肥處理用量為150 kg·hm?2，有機肥結合氮肥處理用量為75 kg·hm?2。施肥方式為播種前一次性投入，耕作方式為旋耕，玉米秸稈粉碎還田（切割為5 cm 左右）后旋耕，與土壤均勻混合，起壟，地膜覆蓋，小區面積為42.9 m2（3.3 m×13.0 m），小區間走道50 cm。試驗從2012 年開始，本研究2018 年開始，參試玉米品種為“先玉335”，各處理播種量均為5.25 萬株·hm?2，玉米于4 月25 日用點播器進行播種，10 月15 日收獲。

1.3 測定指標與方法

1.3.1 土壤CO2排放速率 在玉米的生育期內，選擇無雨天氣，用LI-8100 A 開路式土壤碳通量測量系統（LICOR，USA）進行田間土壤CO2排放通量測定。每個小區安放1 個土壤呼吸測量腔室（直徑20 cm）。2018 年1 月7 日開始第一次測定，以后每隔15 d 左右測定一次，每次測定時間為北京時間上午09：00?11：00，每個小區重復測定3 次，每次測定時間為90 s［31］。為保證測量腔室內的土壤與覆膜土壤的均一性，用塑料地膜將腔室內覆蓋，在每次測定前1 d，去掉地膜，去除基座腔室內土壤表層的一切活體及凋落物，在整個生育期基座埋設位置保持不變。

1.3.2 土壤理化性質及微生物指標 玉米收獲后在整個小區采用S 型5 點取樣方法采集0～5 cm、5～10 cm和10～30 cm 層土樣，混勻并挑除動植物殘體，立即過2 mm 篩，一部分土樣低溫（4 ℃）冷藏，用于測定土壤微生物量、酶活性等指標，剩余土樣風干后用于測定土壤有機碳與活性有機碳；每個小區隨機選取20 株地上部分帶回實驗室后，105 ℃殺青30 min，然后80 ℃烘干至恒重后，稱量計算產量和秸稈干物質重，根據種植密度計算并按照籽粒含水率12%折算公頃產量；根系生物量占地上部分生物量的系數為0.15［32］；每次采集氣體時測定腔室旁0～5 cm，5～10 cm，10～30 cm 土壤水分和溫度。

采用重鉻酸鉀?濃硫酸外加熱法測定有機碳含量（soil organic carbon，SOC）［33］；采用碘化鈉重液分組法進行分離，分離后碳氮聯合分析儀測定易氧化有機碳含量（readily oxidized organic carbon，ROOC）［33］；采用氯仿熏蒸，0.5 mol·L?1的硫酸鉀溶液浸提（土液比1∶4），浸提后碳氮聯合分析儀測定土壤微生物量碳氮含量（microbial bio?mass carbon，microbial biomass nitrogen，MBC，MBN）（Jena multi N/C 2100 s，Germany），換算 系數為0.38、0.45［34］；采用次氯酸鈉?苯酚鈉比色法測定脲酶（urease，UR）活性，采用3，5?二硝基水楊酸比色法測定蔗糖酶（invertase）活性［35］。土壤水分通過烘干法（105±2）℃測定；土壤溫度采用地溫計測定；土壤容重和土壤孔隙度采用環刀法測定［36］。

1.4 分析方法

1.4.1 土壤CO2累計排放量 利用公式（1）計算土壤CO2碳排放量（carbon emission，CE）［31］：

式中：R為CO2排放速率；i+1與i表示兩次測量之間相隔時間，t表示播種后的天數；系數0.1584 表示將C 排放數值單位μmol CO2·m?2·s?1轉換為g CO2·m?2·h?1，即CO2摩爾質量與3600 s 相乘后將單位轉化為g；0.2727 表示將單位g CO2·m?2·h?1轉換為g C·m?2·h?1，即12/44；24 與10 表示將C 排放數值單位由g C·m?2·h?1轉換為kg C·hm?2。

1.4.2 作物碳排放效率［37］

式中：Y為作物籽粒產量（yield），CEE為單位產量釋放的千克碳（carbon emission efficiency）。

1.4.3 土壤碳庫管理指數［16］碳庫指數（carbon pool management，CPI）=農田土壤總有機碳/參考農田土壤總有機碳；碳庫活度（activity，A）=易氧化有機碳/穩態碳；碳庫活度指數（activity index，AI）=碳庫活度/參考土壤碳庫活度；碳庫管理指數（carbon pool management index，CPMI）=CPI×AI×100；土壤碳素效率（soil carbon efficiency，SCE）=易氧化有機碳/總有機碳×100%；參考土壤選擇不施肥處理。

1.5 數據分析

使用Microsoft Excel 2010 軟件進行數據處理，采用SPSS 24.0 軟件進行方差分析，采用Duncan 法檢驗其差異顯著性（P<0.05），利用SigmaPlot 14.0 和R 3.6.2 編程軟件繪制圖表，用R 3.6.2 編程軟件內的corrplot 包和lavaan 包［38］進行結構方程模型分析。

2 結果與分析

2.1 土壤CO2排放特征

2.1.1 土壤CO2排放速率變化特征 土壤CO2排放速率在全年和生育期均基本呈先增后降的趨勢，非生育期呈逐漸增大（1 月7 日?4 月12 日）和逐漸減小（10 月15 日?12 月10 日）的趨勢（圖2a）。生育期內，6 月23 日土壤CO2排放速率最大，最大為（1083.67±109.91）mg·m?2·h?1（ST 處理），ST、OM、NF 和OMNF 分別比CK 高93.31%、76.97%、75.04%和40.33%；10 月15 日（收獲期）各處理土壤CO2排放速率出現谷值，由大到小依次為ST、NF、OM、OMNF、CK 處理，各處理分別比CK 處理高85.48%、79.07%、78.40%、74.81%；生育期土壤CO2排放平均速率最大值為（592.89±252.89）mg·m?2·h?1（ST 處理），ST、OMNF、OM、NF 處理分別比CK 處理高42.72%、30.82%、29.79%、27.28%，中位數集中在423.63～564.14 mg·m?2·h?1（圖2b），依次為OMNF（574.67 mg·m?2·h?1）、ST（564.14 mg·m?2·h?1）、OM（548.15 mg·m?2·h?1）、NF（501.11 mg·m?2·h?1）、CK（423.63 mg·m?2·h?1）；非生育期各處理土壤CO2排放速率與土壤溫度的增減趨勢一致（圖2 和圖3a），平均值由大到小依次為ST、OMNF、OM、NF、CK 處理，各處理分別比CK 高12.99%、8.53%、7.31%、1.34%；全年變化的平均值依次為ST［（428.89±301.91）mg·m?2·h?1］、OMNF［（395.53±271.42）mg·m?2·h?1］、OM［（392.22±264.37）mg·m?2·h?1］、NF［（382.77±273.34）mg·m?2·h?1］、CK［（310.55±208.77）mg·m?2·h?1］，中位數集 中在366.26～470.62 mg·m?2·h?1（圖2b），依 次 為OMNF（470.62 mg·m?2·h?1）、OM（467.66 mg·m?2·h?1）、NF（416.59 mg·m?2·h?1）、ST（396.32 mg·m?2·h?1）、CK（336.62 mg·m?2·h?1）。通過生育期、非生育期、全年的土壤CO2排放通量的最大值、最小值、平均數和中位數分析表明，有機物料處理下，全年土壤CO2排放平均速率高于氮肥處理，其中ST 處理下土壤CO2排放平均速率最大且變幅最大。

圖2 不同時期土壤CO2排放通量Fig.2 Soil CO2 emission flux in different periods

2.1.2 作物碳排放效率 各處理的籽粒產量結果表明（表2），NF、OMNF 處理籽粒產量較高，兩者間差異不顯著（P>0.05），與OM、ST、CK 間差異顯著（P<0.05），較OM、ST、CK 分別提高了75.38%、54.33%、85.74%、54.33%和151.41%、121.23%；生育期和全年CE 在ST、OM、OMNF 處理間差異不顯著（P>0.05），但ST 和OM 顯著高于NF、CK 處理（P<0.05），生育期ST、OM 處理分別比CK 提高了36.35%、32.45%，OMNF 與NF 處理間差異不顯著（P>0.05），非生育期處理間CE 差異均不顯著（P>0.05）；全年CEE 表現為NF、OMNF 處理間差異不顯著（P>0.05），但顯著低于其他處理（P<0.05），ST 處理最高。

表2 不同培肥處理對農田作物產量和碳排放效率的影響Table 2 Crop yield and carbon emission efficiency under different fertilization treatments

2.2 土壤理化性質、微生物特征和碳庫管理指數

2.2.1 土壤理化性質和微生物特征 土壤溫度（0～30 cm）呈拋物線趨勢（圖3a），各處理土壤溫度變化趨勢基本一致，全年土壤溫度和生育期土壤溫度分別在?0.099～25.800 ℃和7.80～25.80 ℃內波動，均表現為CK>NF>OM>OMNF>ST，各處理均在6 月10 日出現峰值，變化范圍在27.41～29.51 ℃；非生育期ST 處理土壤溫度高于其他處理。各處理土壤含水量（0～30 cm）均呈波動趨勢（圖3b），土壤含水量在雨季相對較高（圖1）。全年各處理平均值在17.98%～19.87%間波動，土壤平均含水量變化表現為ST>OMNF>CK>NF>OM，生育期土壤含水量在17.21%～18.34%間波動，平均含水量變化均表現為ST>OMNF>NF>CK>OM；非生育期土壤含水量在16.71%～24.67%間波動，平均含水量變化表現為ST>OM>CK>NF>OMNF，非生育期ST 處理土壤含水量高于其他處理，最高達34.12%。各處理土壤有機碳含量均隨土層加深而降低（圖3c），在7.13～9.45 g·kg?1變化，0～5 cm、5～10 cm 土層處理間有機碳差異顯著（P<0.05），均表現為ST≥OM≥OMNF、OM?NF≥NF>CK，10～30 cm 土層ST、OM 處理有機碳含量顯著高于其他處理（P<0.05）；0～30 cm 土層各施肥處理有機碳差異不顯著（P>0.05），與CK 處理差異顯著（P<0.05），0～30 cm 土層OM、ST、NF、OMNF 處理有機碳分別比CK 增加了18.84%、17.31%、11.14%、10.07%。各處理土壤易氧化有機碳含量均隨土層加深而降低（圖3d），在1.17～2.64 g·kg?1變化，0～30 cm 間各土層處理間易氧化有機碳差異顯著（P<0.05），0～5 cm、5～10 cm 土層均表現為ST 最大，OM、OMNF 次之；0～30 cm 土層表現為OM、OMNF、ST 最大，OM、ST、OMNF、NF分別比CK 增加了95.08%、90.16%、89.34%、24.59%。

各處理土壤微生物量碳含量均隨土層加深而降低（圖3e），在38.50～268.85 mg·kg?1變化，0～30 cm 各土層處理間土壤微生物量碳差異顯著（P<0.05），均表現為ST 最大，OM 次之，在0～30 cm 土層ST、OM、OMNF、NF分別比CK 增加了355.08%、251.18%、216.85%、192.16%。各處理土壤微生物量氮含量均隨土層加深而降低（圖3f），在10.43～37.08 mg·kg?1變化，0～30 cm 各土層處理間土壤微生物量氮差異顯著（P<0.05），均表現為ST 最 大，OMNF、OM 次 之，在0～30cm 土層ST、OMNF、OM、NF 分別 比CK 增加 了106.34%、34.61%、31.23%、13.53%。各處理土壤脲酶活性均隨土層加深而降低（圖3g），在0.96～1.48 mg·g?1變化，0～30 cm 各土層處理間土壤脲酶活性差異顯著（P<0.05），在0～30 cm 土層NF、ST、OMNF、OM 分別比CK 增加了41.51%、37.87%、27.56%、25.05%。各處理土壤蔗糖酶活性均隨土層加深而降低（圖3h），在19.38～34.35 mg·g?1變化，0～30 cm 各土層處理間土壤蔗糖酶活性差異顯著（P<0.05），均表現為ST 最大，OM 和OMNF 次之，在0～30 cm土層ST、OM、OMNF、NF 分別比CK 增加了48.57%、40.41%、21.73%、17.90%。

圖3 不同培肥處理下土壤理化性質和微生物指標Fig.3 Soil physical chemistry properties and microbial indexes under different fertilization treatments

2.2.2 土壤碳庫管理指數 施肥處理下各土層碳庫活度、碳庫活度指數、碳庫指數、碳庫管理指數和土壤碳素效率與CK 差異顯著（P<0.05），ST、OM、OMNF 處理間差異不顯著（P>0.05），在0～5 cm、5～10 cm 土層，ST、OM 和OMNF 處理均能提高碳庫管理指數和土壤碳素效率；在10～30 cm 土層，OMNF 處理的碳庫管理指數和土壤碳素效率最大；在0～30 cm 土層，OMNF、OM 處理的碳庫管理指數和土壤碳素效率均最大（表3）。綜上，OM?NF、ST、OM、NF 均能提高土壤0～30 cm 碳庫管理指數和土壤碳素效率，OMNF 提高幅度最大。

2.3 土壤CO2排放速率與土壤環境因子之間的關系

土壤CO2排放速率與土壤環境因子（0～30 cm）之間的結構方程模型表明（圖4）：環境因子對土壤CO2排放通量的總解釋度為53%，影響總效應為5.65，其中STP（2.36）、MBC（1.59）、SWC（1.18）影響較大，Invertase （0.16），ROOC（0.14）、MBN（0.13）的影響較小，脲酶無顯著影響。STP、MBC 和SMC 均間接地影響著土壤CO2排放通量，其中，STP 通過作用MBC（0.93）、蔗糖酶（0.56）、ROOC（0.01）影響土壤CO2排放通量，MBC 通過作用MBN（1.19）、ROOC（0.81）影響土壤CO2排放通量，SMC 通過作用MBC（0.39）影響土壤CO2排放通量。MBC、STP 均 對ROOC 有顯著正影響（P<0.05），MBC（0.81）影響較大，MBC 對MBN 有顯著正影響（P<0.05）。

表3 不同施肥措施對土壤碳庫管理指數的影響Table 3 Soil carbon pool management index under different fertilization treatments

圖4 土壤理化性質和微生物指標對土壤CO2排放通量的影響Fig.4 Effects of soil physical chemistry properties and microbial indexes on soil CO2 emission flux

3 討論

3.1 不同培肥措施對土壤CO2排放速率和排放效率的影響

土壤CO2排放是一個復雜的生物化學過程，在玉米、小麥農田中的諸多研究［16，23，39?40］均得出，土壤CO2平均排放速率在全年內呈拋物線規律且在生育期內先增強后減弱，土壤CO2排放速率最大值一般出現在大喇叭口期［10，41］，本試驗結果與此結論相符，主要是由于大喇叭口期為營養生長的最旺盛時期，該時期溫度較高，降水適中，加之地膜覆蓋，土壤濕度變化不大，土壤微生物活躍，土壤微生物呼吸和根系呼吸強烈。

針對土壤CO2排放速率變化劇烈的情況，本研究引入了平均值和中位數來共同分析，研究發現無論是生育期還是非生育期，有機物料處理均能顯著提高土壤CO2平均排放速率，其中ST 處理最大，OM、OMNF 處理次之，因為有機物料的投入增加了土壤的有機碳庫（圖3c），為土壤微生物提供了營養，增加了微生物活性，提高了微生物呼吸、有機碳庫和有機物料的分解；其次，秸稈保溫保水性好（圖3a，b），產生的促進作用更大；再次，黃土高原干旱少雨，較少的降水為微生物活動提供了水分，促進微生物產生較多有機質分解胞內、胞外酶［42］，進而加劇秸稈和有機肥料的礦化，產生的營養物質促進植物根系生長，同時為土壤微生物和微小動物提供營養，這些協同作用使根系和微生物呼吸加強，土壤溫度、水分和土壤微生物量碳對土壤CO2排放速率的效應也證實了此結論（圖4）。已有研究表明少量降水對土壤呼吸的影響可能存在激發效應［43］，隨著降水的增加，土壤中水分含量增大，會填充土壤的孔隙，從而抑制土壤呼吸。本研究不排除降水激發效應的存在，但本研究未在每次降水前后增加采樣次數，采樣頻率不足，不能很好地解釋降水的激發效應，還需要進一步開展降水模擬試驗區去驗證，特別是結合降水強度、頻度、地溫等因素開展相關研究。最后，有機物料可以顯著降低土壤容重、提高孔隙度、有機質和養分含量，增加根系分泌物，促進根系呼吸［44］。而土壤CO2平均排放速率和中位數差異可能主要是由于水分、溫度、養分和微生物的交互作用影響（圖4），劉曉雨［45］的研究也證實了這一點。本研究發現有機物料培肥處理的土壤CO2排放量顯著高于氮肥處理，有機物料處理間無顯著差異，但ST 處理有明顯增高的趨勢，無差異的原因可能是由于地膜覆蓋增加了土壤溫度，減少了土壤蒸發，降低了土壤水分和溫度的差異性（圖3a，b），減緩了不同處理的效應［46］。本研究還得出有機肥結合無機肥處理較其他有機物料處理顯著提高了作物產量，降低了土壤碳排放效率，這與馬德帝［47］的研究結論一致，產量增加的主要原因是有機無機配施初期，有機肥因分解慢，無機肥可滿足作物早期的養分需求，有利于作物的生長，到作物生長后期，有機肥分解的營養物質又滿足了氮肥的虧缺，其作用等同于緩釋肥；另一方面，有機無機配施投入后，同時增加了碳素和氮素投入，可調整土壤C/N 處于合理的范圍內，使微生物活性增大（圖4），促進有機肥的礦化，為作物提供了充足的養分，而秸稈處理和單純的有機肥處理無機營養物質釋放緩慢，還可能會導致作物和土壤微生物競爭氮素，從而降低作物產量，進而增大作物的碳排放效率，也說明了秸稈和無機肥配施比單施秸稈或有機肥的效果好［48］。

3.2 不同培肥措施對土壤碳庫和微生物特征的影響

本研究得出，有機物料處理顯著提高了各土層土壤微生物碳氮和蔗糖酶活性，且均以秸稈處理提高最大，有機肥和有機無機配施次之（圖3）。國內外關于施肥體系下土壤微生物特征的研究結果一致表明施用有機物料后，較不施肥或單施化肥顯著提高土壤微生物碳氮和蔗糖酶活性以及土壤碳庫［15，49?54］，形成原因主要是有機物料可以為土壤微生物增加額外的營養源和能源，土壤酶產生于土壤微生物的生命活動，與微生物量成正比。土壤碳庫管理指數是土壤有機碳庫變化的靈敏指標，土壤碳素效率高低能反映土壤有機碳被微生物分解利用的速度［16］。本研究得出施用有機物料處理均能提高各研究土層碳庫管理指數和土壤碳素效率，有機無機配施提升幅度最大，也均能提高0～30 cm 有機碳和易氧化有機碳含量，有機肥處理提高幅度最大，與已有研究基本一致［50，54?55］。主要是由于有機物料的投入增加了土壤碳，為微生物活動提供了營養物質，促進了微生物的活性，加速了有機殘余物質的分解、某些無機化合物的轉化等生物化學過程，特別是有機無機配施，使土壤C/N 處于合理的范圍之內，更有助于酶促反應的進行，從而提高易氧化有機碳含量及碳庫管理指數。本研究得出與不施肥處理相比，氮肥處理的土壤有機碳、易氧化有機碳含量也有所升高，這與徐陽春等［56］的研究結論不一致，主要是由于黃土高原旱地土壤鉀肥充足，氮肥和磷肥施入等同于氮磷鉀肥同施，所以提高了土壤有機碳和易氧化有機碳，與Dai 等［57］的研究結果一致，但提高幅度不大。本試驗中，玉米采用旋耕，耕層在0～20 cm 左右，有機物料施入土壤埋藏深度較淺，這也說明了土壤化學和微生物特性隨土壤深度逐漸減低的現象。秸稈處理的水分含量在部分非生育期較高，進一步分析發現1 月7 日、2 月6 日、3 月4 日水分含量較高是由于秸稈處理10～30 cm 水分含量較高，分析其原因可能是秸稈處理土壤孔隙度較其他處理大，玉米收獲后，降水（2017 年10?12 月降水量17.1 mm、2018 年2 月6 日前降水3 mm、2 月6 日?3 月4 日降水4.7 mm）因無植株截留而大部分滲入土壤空隙，同時地膜對水分也有一定的保護作用，加之該時間段氣溫和土壤溫度在全年內均很低，從而減少了水分的散失，已有研究也表明秸稈和地膜交互作用呈現雙重的蓄水保墑效應［58］。

3.3 土壤CO2排放與環境因子之間的關系

已有研究普遍認為農田土壤CO2排放的最重要的影響因素是土壤溫度、降水量、土壤水分［10，25，28］，與土壤溫度正相關，與土壤水分呈Birch 效應［42］，即水分含量在一定的范圍內是激發效應，隨著含水量的增大，土壤CO2排放量被抑制。本研究通過結構方程模型得出土壤溫度和水分是影響土壤CO2排放通量的主要因素，土壤溫度對土壤CO2排放通量呈顯著正效應（P<0.05），這是因為溫度是決定陸地生態系統碳循環過程的關鍵因素，溫度影響著土壤呼吸過程的所有環節，溫度升高可能會改變土壤微生物群落并有助于土壤團聚體形成的微生物次生化合物的分解［59］，在一定范圍內增加溫度可以提高土壤微生物活性，進而影響有機碳的礦化，從而使氣溫對土壤CO2排放具有極顯著正效應［60?61］；本研究得出土壤含水量對土壤CO2排放通量呈顯著正效應（P<0.05），與已有研究［10，25，28］略有區別，但這也符合黃土高原旱區的實際情況，黃土高原屬于雨養農業，旱地農田土壤含水量低，降水促進了土壤微生物、土壤微小動物和根系的活性，相應的呼吸作用增加了土壤CO2排放量［42］。本研究還得出微生物量碳和蔗糖酶與土壤CO2排放通量顯著正相關（P<0.05），與已有研究結果一致［62?63］。土壤呼吸速率的重要組分土壤微生物呼吸與土壤微生物量碳存在高度相關關系［10］，有機碳的礦化也需要微生物高度參與［44，61］，蔗糖酶能促進土壤中蔗糖水解為葡萄糖和果糖，隨后會進一步分解產生CO2，因此蔗糖酶與土壤呼吸強度也有關［59］，同時土壤溫度也促進了微生物量碳和蔗糖酶活性的提高（圖4）。本研究得出，微生物量碳對微生物量氮有顯著影響，可能是輸入外源碳后出現爭氮現象［64］。本研究得出土壤理化、微生物因子對土壤CO2排放通量的總體解釋度為53%，土壤CO2排放可能還存在其他的解釋因子，還需進一步開展糖苷酶、木糖苷酶等土壤胞外酶指標的研究［62?65］。

4 結論

不同培肥措施對黃土高原旱地土壤CO2排放速率、有機碳、易氧化有機碳、土壤碳庫管理指數、微生物量碳氮、蔗糖酶、脲酶、玉米產量和作物碳排放效率均產生顯著影響。秸稈、有機無機培肥、有機肥處理均能提高土壤CO2排放平均速率，其中秸稈處理提高幅度最大，秸稈、有機肥處理土壤CO2排放量顯著高于氮肥和不施肥處理。有機物料處理均能顯著提高0～5 cm 土壤易氧化有機碳、微生物碳氮、蔗糖酶活性和碳庫管理指數，其中，有機無機肥結合處理對土壤碳庫管理指數的提高幅度最大。有機無機配施處理和單施氮肥處理能顯著提高作物產量，降低碳排放效率。影響土壤CO2排放通量的重要因子包括土壤溫度、微生物量碳、土壤水分，均對土壤CO2排放通量呈顯著正效應，土壤溫度促進了微生物量碳和蔗糖酶活性的提高，微生物量碳促進了微生物量氮和易氧化有機碳的增加。綜合來看，有機無機配施處理可以提升土壤碳庫管理指數，保持微生物活性，增加作物產量，降低土壤碳排放效率，是比較適宜的農田培肥措施。在農村環境治理和改善農村人居環境的大背景下，若秸稈還田建議配施無機肥。