傅 勇, 王 淘, 楊志平, 周 偉, 劉 琦, 任萬(wàn)軍, 陳 勇

四川農(nóng)業(yè)大學(xué)農(nóng)學(xué)院/作物生理生態(tài)及栽培四川省重點(diǎn)實(shí)驗(yàn)室, 溫江 611130

工業(yè)革命以來(lái),由于人類活動(dòng)的影響,大氣CO2濃度急劇升高,增加了約 32%,全球正經(jīng)歷著變暖的危機(jī)[1,2]。農(nóng)田生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)的重要組成部分,占全球陸地面積的10.5%,是全球碳循環(huán)中最活躍的部分,CO2占人為溫室氣體排放量的21%—25%,且受人類干擾頻繁[3]。土壤碳庫(kù)是大氣碳庫(kù)的2倍,土壤中的CO2濃度發(fā)生微量的變化,就會(huì)對(duì)大氣中CO2濃度造成巨大影響[4, 5]。土壤呼吸是農(nóng)田生態(tài)系統(tǒng)土壤碳排放的主要途徑,土壤通過(guò)呼吸作用向大氣釋放CO2的年通量為5.0×1016—7.5×1016g,是土壤碳素以CO2形式回歸大氣碳庫(kù)的主要路徑,占陸地生態(tài)系統(tǒng)與大氣之間碳交換總量的2/3,對(duì)氣候變化和碳循環(huán)起著不可忽視的重要作用[6, 7]。土壤呼吸作為碳循環(huán)的核心問(wèn)題已引起了國(guó)內(nèi)外學(xué)者的密切關(guān)注和廣泛研究,并取得了大量的研究成果。

近年來(lái),我國(guó)農(nóng)田生態(tài)系統(tǒng)土壤呼吸的研究主要集中在不同作物類型、土地利用方式或栽培措施下土壤呼吸的時(shí)空變異特征及調(diào)控因素的研究,研究對(duì)象涉及小麥、玉米、水稻等主要農(nóng)作物[8- 10]。很多研究表明,農(nóng)田土壤呼吸受環(huán)境因子共同調(diào)控,尤其對(duì)土壤溫度和土壤含水量的響應(yīng)十分敏感[11-12]。也有研究表明,與土壤溫度相比,土壤水分對(duì)土壤呼吸季節(jié)變化的影響更為重要[13]。但前人關(guān)于作物系統(tǒng)的研究對(duì)象多為旱地輪作系統(tǒng)或單一作物系統(tǒng),缺乏對(duì)稻田水旱輪作系統(tǒng)周年土壤呼吸的探索。水旱輪作是我國(guó)南方稻區(qū)十分重要的種植制度,但由于旱、水土壤環(huán)境轉(zhuǎn)換,土壤理化性質(zhì)、土壤有機(jī)質(zhì)礦化分解速率、微生物種類、微生物數(shù)量和活性、根際生物量等發(fā)生劇烈變化,土壤呼吸隨之受到極大的影響[14,15]。為了從周年種植模式角度探討不同水旱輪作系統(tǒng)周年土壤呼吸特征及其影響因素,本研究基于蒜-稻、麥-稻以及油-稻三種川西平原灌區(qū)重要的水旱輪作模式的定位試驗(yàn),探明不同輪作模式及其不同生長(zhǎng)季土壤呼吸的變化趨勢(shì)和差異特征;探明土壤呼吸與土壤溫度、土壤含水量等生態(tài)因子的關(guān)系,以期為水旱輪作系統(tǒng)減排降耗提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

試驗(yàn)設(shè)在四川省崇州市四川農(nóng)業(yè)大學(xué)現(xiàn)代農(nóng)業(yè)研發(fā)基地(30°33′N,103°38′E),位于川西平原都江堰灌區(qū),屬亞熱帶濕潤(rùn)季風(fēng)氣候區(qū)。該區(qū)年平均氣溫15.9℃,平均降雨量1012.4 mm,平均日照時(shí)數(shù)為1161.5 h。試驗(yàn)地基礎(chǔ)養(yǎng)分含量為土壤有機(jī)質(zhì)31.58 g/kg,全氮1.85 g/kg,速效磷11.81 mg/kg,速效鉀85.48 mg/kg。

1.2 試驗(yàn)設(shè)計(jì)與管理

試驗(yàn)地采用水旱輪作方式進(jìn)行試驗(yàn),已連續(xù)進(jìn)行3年的水旱輪作模式研究。采用單因素隨機(jī)區(qū)組設(shè)計(jì),設(shè)置3個(gè)水旱輪作模式處理,分別為小麥-水稻(WR)、大蒜-水稻(GR)和油菜-水稻(RR),每年的處理保持一致, 3次重復(fù),共9個(gè)小區(qū),小區(qū)面積為26 m2(4 m×6.5 m)。除水稻季施肥采用氮肥后移技術(shù)外,其余各作物施肥及管理方案均按照農(nóng)民常規(guī)方法進(jìn)行。本研究的時(shí)間范圍為2017年5月—2018年5月,為1個(gè)完整的水旱輪作周期,其中2017年5月—9月為水稻季,2017年9月—2018年5月為旱作季,輪作周期的各作物播栽期及播栽密度見(jiàn)表1,耕作、施肥及管理措施見(jiàn)表2。

表1 各作物播栽期及播栽密度

表2 不同作物耕作、施肥及管理措施

1.3 測(cè)定指標(biāo)與方法

1.3.1土壤呼吸測(cè)定

土壤呼吸測(cè)定于2017年5月25日開(kāi)始,每隔30 d選擇晴天(遇雨適當(dāng)延后)進(jìn)行土壤呼吸測(cè)定[16,17]。其中,土壤呼吸季節(jié)動(dòng)態(tài)選擇取樣當(dāng)天9:00—11:00進(jìn)行取樣,土壤呼吸日動(dòng)態(tài)取樣從9:00至18:00每隔3 h取樣一次;18:00至次日6:00每隔4 h取樣一次。水稻移栽后,于各小區(qū)行間分別放置1個(gè)由PVC材料制成的土壤呼吸氣體收集裝置,該裝置參照王維鈺等[18]的設(shè)計(jì),并根據(jù)田間取樣環(huán)境進(jìn)行了一定的改進(jìn),形成腔室高25 cm,直徑為16 cm,插入土壤5 cm。腔室頂部為開(kāi)放結(jié)構(gòu),開(kāi)放部分直徑為11.5 cm,確保腔室內(nèi)外部環(huán)境條件一致。取樣時(shí)加蓋密封,蓋頂開(kāi)有兩個(gè)小孔,分別用于安裝溫度計(jì)和連接橡膠軟管(橡膠軟管連接至PVC腔室中心)。取樣時(shí)加蓋密封后,蓋頂用橡膠軟管連接注射器,采樣時(shí)間為加蓋后0、10、20和30 min,取60 mL氣體置于真空袋中待測(cè),取樣后立即用鐵夾夾住橡膠軟管保持密閉。樣品濃度用島津2010Plus氣相色譜儀檢測(cè)分析,通過(guò)對(duì)4個(gè)氣樣濃度進(jìn)行線性回歸,計(jì)算氣體排放速率。氣體排放速率計(jì)算公式如下[19]

式中,F為CO2排放速率,mg m-2h-1;ρ為標(biāo)準(zhǔn)大氣壓下的CO2密度,為 1.97 kg/m3;T為采樣過(guò)程中腔室內(nèi)的平均溫度,℃;H是腔室的凈高度,m;dC/dt是采樣氣體濃度的變化率。

土壤呼吸累積排放量計(jì)算公式如下：

式中,M為土壤CO2排放總量(kg/hm2);F為CO2排放速率(mg m-2h-1);i為第i次測(cè)定;(ti+1-ti)為相鄰兩次測(cè)定間隔天數(shù)。

1.3.2水熱生態(tài)因子測(cè)定

于每月土壤呼吸取樣時(shí),同步測(cè)定土壤溫度、土壤含水量,稻季測(cè)定水層高度。土壤溫度采用便攜式溫度測(cè)定儀測(cè)定,分0—10 cm和10—20 cm測(cè)定。各小區(qū)按五點(diǎn)取樣法采集0—20 cm土樣,采用烘干法測(cè)定土壤含水量;田間水層高度用直尺測(cè)量。

1.3.3土壤呼吸與水熱因子擬合

土壤呼吸速率和土壤溫度、土壤含水量和水層高度分別采用指數(shù)函數(shù)模型和一元二次模型擬合[20,21]：

Y=aebT

Y=X2+X+c

溫度敏感性系數(shù)Q10代表土壤呼吸的溫度敏感性,計(jì)算公式如下[17-18]：

Q10=e10b

式中,Y為土壤呼吸速率(mg m-2h-1),T為平均溫度(℃),a為0℃時(shí)土壤呼吸速率,b為溫度反應(yīng)系數(shù),X為土壤含水量(%),c為截距,Q10指溫度升高10℃時(shí)土壤呼吸速率變化的倍數(shù)。

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

采用SPSS13.0進(jìn)行數(shù)據(jù)處理,并用Microsoft Excel 2010作圖。

2 結(jié)果分析

2.1 不同水旱輪作模式土壤呼吸速率季節(jié)變化特征

圖1表明,三種模式土壤呼吸速率季節(jié)變化規(guī)律趨于一致,周年呈現(xiàn)雙峰型變化趨勢(shì),6—9月土壤呼吸速率顯著高于其余月份。三種模式土壤呼吸速率分別在6月和9月達(dá)到峰值,其中GR和RR模式在6月達(dá)到全年呼吸速率最大值,分別為63.1 mg m-2h-1和61.2 mg m-2h-1,WR模式則在9月達(dá)到全年呼吸速率最大值,為36.82 mg m-2h-1;而三種模式土壤呼吸速率最低值出現(xiàn)在1月,GR、WR、RR模式土壤呼吸速率最低值分別為3.89 mg m-2h-1、4.49 mg m-2h-1和1.42 mg m-2h-1。不同模式之間,GR模式全年平均土壤呼吸速率顯著高于其余兩種模式,而以WR模式最低。

2.2 不同作物生長(zhǎng)季節(jié)土壤呼吸累積排放量

圖2表明, GR、RR和WR三種模式周年土壤呼吸累積排放量分別為193.36 kg/hm2、157.28 kg/hm2和160.27 kg/hm2,其中累積排放最高的GR模式分別比RR和WR模式高17.11%和18.5%。季間空閑期,WR模式土壤呼吸累積排放量為41.84 kg/hm2,顯著高于其余兩種模式;從不同作物生長(zhǎng)季來(lái)看,GR、WR、RR三種模式均在稻季累積排放量最高,分別占周年累積排放量的49.28%、48.11%和58.33%; GR、WR、RR三種模式稻季土壤呼吸累積排放量較旱季分別高16.43 kg/hm2、35.88 kg/hm2、53.68 kg/hm2。在稻季,GR模式和RR模式的土壤呼吸累計(jì)排放量顯著高于WR模式,兩種模式分別比WR模式高出25.93%和23.56%;而旱季GR模式土壤呼吸累計(jì)排放量達(dá)78.5 kg/hm2,顯著高于WR模式和RR模式,分別比WR模式和RR模式分別高出98.23%和98.08%。

圖1 不同水旱輪作模式周年土壤呼吸速率特征 Fig.1 Characteristics of soil respiration rate in different paddy-upland rotation

圖2 不同水旱輪作模式土壤呼吸累積排放量 Fig.2 Cumulative emission of soil respiration in different paddy-upland crop rotation system同組中不同英文字母分別代表多重比較下在P<0.01水平下差異達(dá)顯著水平

2.3 不同水旱輪作模式土壤呼吸速率日動(dòng)態(tài)變化特征

如圖3所示,三種模式土壤呼吸速率晝夜變化主要呈現(xiàn)單峰型變化,總體表現(xiàn)為白天高,夜晚低,具有一定波動(dòng)性。各模式土壤呼吸速率在12:00—15:00達(dá)到最大值,其中WR模式在2月、10月和12月最大值出現(xiàn)的時(shí)間滯后于其余兩種模式;三種模式土壤呼吸速率在2:00—6:00出現(xiàn)最低值,不同季節(jié)各模式最低值出現(xiàn)的時(shí)間略有差異。GR、RR模式土壤呼吸速率日變化幅度在6月出現(xiàn)最大值,分別為31.11 mg m-2h-1和52.52 mg m-2h-1,WR模式土壤呼吸速率日變化幅度在8月出現(xiàn)最大值,為28.29 mg m-2h-1;GR和RR模式在2月出現(xiàn)日變化幅度最小值,分別為6.52 mg m-2h-1和3.56 mg m-2h-1,而WR模式日變化幅度最小值出現(xiàn)在12月,為7.63 mg m-2h-1。GR模式各月份土壤呼吸速率日動(dòng)態(tài)波動(dòng)較小,而WR模式和RR模式土壤呼吸速率日動(dòng)態(tài)變化相對(duì)劇烈;WR模式土壤呼吸速率在2月和4月的15:00—18:00之間分別下降54.34%和58.07%;RR模式土壤呼吸速率在4月和6月的15:00—18:00之間分別下降65.21%和30.22%。

圖3 不同水旱輪作模式土壤呼吸速率日動(dòng)態(tài)特征Fig.3 Daily dynamic characteristics of soil respiration rate in different paddy-upland crop rotation system

2.4 土壤呼吸速率與土壤含水量、水層高度及土壤溫度的關(guān)系

通過(guò)對(duì)土壤呼吸速率與旱季土壤含水量、稻季水層高度進(jìn)行回歸分析,結(jié)果表明,在旱季,土壤呼吸速率與土壤含水量沒(méi)有顯著相關(guān)性(P>0.05),但稻季受水層高度影響較大,其相關(guān)性達(dá)到顯著水平(P<0.05)。通過(guò)對(duì)土壤呼吸速率與土壤溫度的回歸分析發(fā)現(xiàn),不同模式土壤呼吸速率與0—10 cm、10—20 cm土壤溫度均存在極顯著的相關(guān)性(P<0.01),即土壤呼吸速率隨著土壤溫度的升高而增高。不同模式土壤呼吸速率對(duì)土壤溫度的敏感性有所差異,表現(xiàn)為GR>WR>RR;從不同土層來(lái)看,三種模式Q10均表現(xiàn)為0—10 cm>10—20 cm,表明隨土層的深入,土壤呼吸速率與土壤溫度的敏感性逐漸減小。不同模式不同土層深度Q10變化范圍為1.15—2.72,以GR模式0—10 cm最高,RR模式10—20 cm最低。

3 討論

3.1 不同水旱輪作模式對(duì)土壤呼吸速率影響

本研究結(jié)果表明,三種模式土壤周年呼吸累積排放量表現(xiàn)為GR>RR>WR。研究表明[22,23],秸稈還田能夠明顯提高農(nóng)田土壤碳固定,有效提高土壤有機(jī)碳及活性有機(jī)碳組分含量,而土壤有機(jī)質(zhì)是微生物進(jìn)行分解活動(dòng)排放CO2的物質(zhì)基礎(chǔ)[24]。前人研究表明[25],由于“啟動(dòng)效應(yīng)”的存在,土壤中碳的投入會(huì)導(dǎo)致土壤礦化的短期增加,從而導(dǎo)致碳流失。在本研究中,GR模式由于旱季進(jìn)行秸稈覆蓋還田,增加了外源有機(jī)質(zhì)的輸入,提高了土壤有機(jī)碳含量,為微生物活動(dòng)提供了大量底物,同時(shí)加速土壤礦化,增加了土壤呼吸速率,從而使GR模式土壤呼吸周年累積排放量顯著高于WR、RR模式。此外,王春新、張慶忠[26, 27]等研究表明,施氮量對(duì)土壤呼吸具有顯著影響,在本研究中,GR模式旱季施氮量遠(yuǎn)高于WR、RR模式,GR模式旱季較高的氮投入不僅促進(jìn)了大蒜地上和地下部分生物量的增加,同時(shí)殘留的氮素在稻季得到利用,增加了GR模式稻季生物量的累積,從而間接提高了GR模式土壤呼吸速率,這也可能是導(dǎo)致GR模式土壤呼吸累積排放量高于其余兩種模式的原因之一。

表3 土壤呼吸速率與含水量、水層的關(guān)系

*—P<0.05 ;**—P<0.01

圖4 土壤呼吸速率與土壤溫度的關(guān)系Fig.4 Relationship between soil respiration rate and soil temperature

從土壤呼吸速率季節(jié)變化特征來(lái)看,水稻生長(zhǎng)季土壤呼吸速率要顯著高于其余季節(jié)。在水稻生長(zhǎng)季中,6月底處于水稻分蘗末期曬田階段,植株根系生長(zhǎng)旺盛,根際微生物活躍,土壤呼吸速率增強(qiáng);此外由于沒(méi)有了水層的阻隔,土壤中閉蓄態(tài)CO2得以釋放[28, 29],也使得土壤呼吸速率增大。9月水稻收獲后,田間含水量降低且不受植株遮陰影響,土壤呼吸速率又出現(xiàn)一個(gè)峰值。GR模式大蒜播種后,受秸稈覆蓋還田和氮肥投入影響,GR模式土壤呼吸速率高于其余兩種模式。旱季作物生長(zhǎng)過(guò)程中,隨著氣溫的逐漸降低,三種模式土壤呼吸速率逐漸下降,至1月降至最低,之后土壤呼吸速率緩慢回升。但三種模式土壤呼吸速率在12月份有所上升,這可能是由于作物生育進(jìn)程和地下生物量的增加從而增加了土壤微生物的活性,進(jìn)而提高了土壤呼吸速率。

三種水旱輪作模式土壤呼吸日動(dòng)態(tài)主要表現(xiàn)為單峰形式,土壤呼吸速率總體白天高于夜間,峰值出現(xiàn)在12:00—15:00之間,最低值則出現(xiàn)在凌晨2:00—6:00之間,這與高會(huì)議、韓廣軒等[30,31]研究結(jié)果一致。但4月WR模式22:00土壤呼吸速率高于18:00土壤呼吸速率,2月GR模式呼吸土壤速率最低值出現(xiàn)在22:00,與氣溫和土壤溫度變化規(guī)律不一致,這與陳亮等[32]在黃河三角洲濱海濕地的研究結(jié)果類似,這可能與溫度降低后地表受霜降形成的水霧有關(guān)。GR模式土壤呼吸速率日均值在8月、10月和12月高于其余兩種模式,這主要是由于秸稈覆蓋還田增加了土壤的有機(jī)質(zhì)含量,為微生物活動(dòng)提供了大量底物,加速了土壤礦化,從而提高了土壤呼吸速率。此外,各模式不同季節(jié)土壤呼吸速率日動(dòng)態(tài)的平均值接近9:00的土壤呼吸速率,與季節(jié)動(dòng)態(tài)的取樣時(shí)間相吻合,這與前人研究結(jié)果相一致[33]。

3.2 土壤呼吸速率與溫度、土壤含水量及水層高度的關(guān)系

土壤溫度是影響土壤呼吸速率的主要因素之一,研究表明[8,34],土壤溫度在一定范圍內(nèi)增加可以提高土壤微生物活性,從而影響土壤呼吸,但當(dāng)土壤溫度超過(guò)一定范圍,抑制了土壤微生物活性,使土壤呼吸速率降低。本研究結(jié)果表明,土壤呼吸速率與土壤溫度呈指數(shù)正相關(guān)關(guān)系,這與崔海[35]、張俊麗[36]等研究結(jié)果一致。研究者們通常用Q10來(lái)表示土壤呼吸速率與土壤溫度之間的關(guān)系,前人研究結(jié)果表明[37],Q10值變化范圍一般在1.3—3.3之間。在本研究中,三種模式Q10值變化范圍在1.15—2.72之間,與Raich等[37]的研究結(jié)果類似。不同模式間,GR模式Q10值最高,表明GR模式土壤呼吸速率與土壤溫度敏感性最高,這可能是與GR模式采用秸稈覆蓋還田有關(guān)。

土壤含水量是除土壤溫度外,影響土壤呼吸速率的又一關(guān)鍵因素。土壤含水量通過(guò)調(diào)節(jié)土壤通氣狀況,影響微生物活性以及調(diào)節(jié)土壤氧化還原狀況來(lái)調(diào)節(jié)土壤呼吸;土壤水分變化對(duì)溫度的影響(通過(guò)水分有效性驅(qū)動(dòng)的感熱通量和潛熱通量的變化)也間接影響土壤呼吸。在本研究旱季作物生長(zhǎng)過(guò)程中,土壤呼吸速率與土壤含水量沒(méi)有顯著相關(guān)性, GR、WR、RR模式分別在含水量達(dá)到22.57%、17.98%、19.44%時(shí),各模式土壤呼吸速率達(dá)到最大值。土壤含水量對(duì)土壤呼吸的影響十分復(fù)雜,往往同時(shí)受與土壤溫度的協(xié)調(diào)情況影響,目前關(guān)于土壤含水量與土壤呼吸速率關(guān)系的研究尚沒(méi)有定論。多數(shù)研究者認(rèn)為,土壤含水量在一定范圍時(shí),隨著土壤含水量的增加,土壤呼吸速率相應(yīng)上升,當(dāng)土壤含水量超出這個(gè)范圍時(shí),則會(huì)抑制土壤呼吸速率。王雪等[38]研究表明,土壤水分與土壤呼吸速率呈倒“U”型相關(guān)關(guān)系,即土壤水分含量較低時(shí),土壤水分的增加可以促進(jìn)土壤呼吸作用,但土壤含水量達(dá)到一定程度后,土壤水分的增加則抑制土壤呼吸作用,這可能是由于土壤在干旱階段,微生物活動(dòng)由于受到水分的限制;而在土壤處于濕潤(rùn)階段時(shí),有機(jī)物的分解過(guò)程從好氧轉(zhuǎn)變?yōu)閰捬?而厭氧分解速率僅為好氧分解速率的30%—40%[39]。丁新宇等[40]研究指出,在不同土壤條件下,土壤含水量呈現(xiàn)完全相反的相關(guān)性,在干旱階段,土壤呼吸速率與土壤含水量正相關(guān),但在濕潤(rùn)階段土壤呼吸速率與土壤含水量則呈負(fù)相關(guān)。

在水稻生長(zhǎng)季,淹水狀態(tài)下土壤呼吸速率受水層深度制約[41]。本研究結(jié)果表明,土壤呼吸速率與稻田水層高度呈顯著負(fù)相關(guān),這與武文明等[28]在西雙版納地區(qū)的研究結(jié)果一致。 由于水層的存在,土壤呼吸受到阻隔,同時(shí)水層能夠溶解部分CO2;此外,水層深度的變化還會(huì)影響土壤溫度、有機(jī)質(zhì)分解、微生物活性和土壤動(dòng)物的數(shù)量等其他因子,進(jìn)而影響土壤呼吸速率。因此水層深度對(duì)土壤呼吸的影響還具有不確定性和復(fù)雜性[42]。

研究表明[43],土壤溫度和含水量的協(xié)同作用能夠更好地解釋土壤呼吸變化情況,解釋能力高于50%,一般比單一因子的解釋能力更強(qiáng),但由于農(nóng)田生態(tài)系統(tǒng)復(fù)雜的環(huán)境因子,目前水熱因子同土壤呼吸之間的關(guān)系仍然難以定量描述。

4 結(jié)論

通過(guò)對(duì)川西平原灌區(qū)三種不同水旱輪作模式的土壤呼吸特征及其水熱影響因素研究,結(jié)果表明,土壤呼吸周年累積排放量以GR模式最高,從不同作物生長(zhǎng)季來(lái)看,稻季>旱季>空閑季。土壤呼吸速率季節(jié)動(dòng)態(tài)呈雙峰曲線,在6月和9月達(dá)到峰值。土壤呼吸速率日動(dòng)態(tài)呈單峰曲線,土壤呼吸速率在12:00—15:00間出現(xiàn)最大值,在2:00—6:00之間出現(xiàn)最低值。三種模式土壤呼吸速率與0—10 cm和10—20 cm土壤溫度呈指數(shù)正相關(guān)且對(duì)溫度變化敏感,Q10值變化范圍在1.15—2.72之間。不同模式間相比,GR模式Q10值最高,土壤呼吸速率與土壤溫度敏感性最高。三種模式土壤呼吸與土壤水分沒(méi)有顯著相關(guān)性,但稻季與水層高度表現(xiàn)為顯著負(fù)相關(guān)。然而,水熱生態(tài)因子對(duì)水旱輪作土壤呼吸的影響是一個(gè)十分復(fù)雜的過(guò)程,既會(huì)直接影響土壤呼吸的排放,也會(huì)通過(guò)改變土壤微生物區(qū)系組成、外源碳的進(jìn)入量等間接影響土壤呼吸,其作用機(jī)制還需要進(jìn)一步論證研究。總體來(lái)看,川西平原灌區(qū)水旱輪作栽培管理過(guò)程中,應(yīng)注意控制蒜-稻輪作模式的氮肥投入,稻季合理管理田間水層,從而有效管理土壤呼吸碳排放。