• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    免耕條件下秸稈還田對冬小麥-夏玉米輪作系統(tǒng)土壤呼吸及土壤水熱狀況的影響

    2016-07-18 07:44:20王維鈺KashifAKHTAR任廣鑫馮永忠
    中國農(nóng)業(yè)科學 2016年11期
    關鍵詞:免耕秸稈還田夏玉米

    王維鈺,喬 博,Kashif AKHTAR,袁 率,任廣鑫,馮永忠

    (1西北農(nóng)林科技大學農(nóng)學院,陜西楊凌 712100;2陜西省循環(huán)農(nóng)業(yè)工程技術研究中心,陜西楊凌 712100)

    ?

    免耕條件下秸稈還田對冬小麥-夏玉米輪作系統(tǒng)土壤呼吸及土壤水熱狀況的影響

    王維鈺1,2,喬博1,2,Kashif AKHTAR2,袁率1,2,任廣鑫1,2,馮永忠1,2

    (1西北農(nóng)林科技大學農(nóng)學院,陜西楊凌 712100;2陜西省循環(huán)農(nóng)業(yè)工程技術研究中心,陜西楊凌 712100)

    摘要:【目的】研究免耕條件下秸稈還田對旱地冬小麥-夏玉米輪作系統(tǒng)土壤呼吸及土壤水熱狀況的影響?!痉椒ā?011年10月至2014年9月,在陜西楊凌設置秸稈全量還田+施肥(S1F1)、秸稈全量還田+不施肥(S1F0)、秸稈半量還田+施肥(S1/2F1)、秸稈半量還田+不施肥(S1/2F0)、秸稈不還田+施肥(S0F1)、秸稈不還田+不施肥(S0F0)6 種不同耕作處理的3年定位試驗,測定并分析不同耕作處理下土壤呼吸、土壤水熱狀況、作物產(chǎn)量、土壤耕作層有機碳含量的差異。【結(jié)果】在冬小麥生育期內(nèi),各處理土壤呼吸速率均呈先下降后升高再下降的趨勢;在夏玉米生育期內(nèi),各處理土壤呼吸速率均表現(xiàn)為先升高后下降的趨勢。同一生育期內(nèi)各處理土壤呼吸平均速率及呼吸總量依次為S1F1>S1/2F1>S1/2F0>S0F1>S1F0>CK,同種作物不同生育期之間,各處理土壤生育期呼吸總量有逐年降低的趨勢。整個研究周期內(nèi),土壤溫度的變化趨勢與每月平均氣溫的變化趨勢相似,不同處理在同一生育期內(nèi)的土壤溫度變化趨勢相近,且各處理生育期土壤平均溫度均隨土壤深度的增加而降低;不同秸稈還田處理冬季土壤溫度均高于對照,但生育期土壤平均溫度均低于對照。土壤含水量隨土壤深度的增加而降低,但受降雨影響,不同輪作周期之間的土壤含水量波動較大,各處理同一生育期的土壤平均含水量均表現(xiàn)為S1F0>S1F1 >S1/2F0>S1/2F1>CK>S0F1,且不同秸稈還田處理的土壤含水量與對照間的差異均顯著(P<0.05);土壤溫度能夠解釋土壤呼吸速率變化的32.5%—60.4%,土壤含水量能夠解釋土壤呼吸速率變化的38.4%—82.5%,不同土層深度間,5 cm土層的溫度與土壤呼吸的擬合度性最高,而10—20 cm土層的含水量與土壤呼吸的擬合度最高。相同年份內(nèi),不同處理冬小麥和夏玉米產(chǎn)量均表現(xiàn)為S1F1>S1F0>S1/2F1>S0F1>S1/2F0>CK,這個研究周期內(nèi),冬小麥產(chǎn)量逐年增加,夏玉米在前兩季表現(xiàn)為增產(chǎn),但受極端炎熱天氣的影響,第三季的產(chǎn)量明顯降低。單季作物收獲后,各處理同一土層深度的有機碳含量均表現(xiàn)為S1F1>S1/2F1>S1F0>S1/2F0>S0F1>CK。且不同秸稈還田處理的土壤有機碳含量逐年升高?!窘Y(jié)論】長期免耕秸稈還田能夠有效降低農(nóng)田土壤碳排放、提高農(nóng)田土壤水分利用率及冬季土壤溫度、提高作物產(chǎn)量及土壤有機碳含量。不同秸稈還田處理間以S1F0處理的效果最優(yōu)。

    關鍵詞:免耕;冬小麥-夏玉米;輪作;秸稈還田;土壤呼吸;水熱狀況

    聯(lián)系方式:王維鈺,E-mail:lz2567@163.com。通信作者馮永忠,E-mail:fengyz@nwsuaf.edu.cn

    0 引言

    【研究意義】土壤呼吸是陸地碳循環(huán)的重要組成部分[1],全球每年因土壤呼吸作用向大氣中釋放的碳通量達75—120 Pg(1Pg=1×1015g)[2],即使土壤呼吸發(fā)生細微的變化都會顯著改變大氣中CO2濃度[3],進而影響氣候變化。農(nóng)田生態(tài)系統(tǒng)作為陸地生態(tài)系統(tǒng)的重要組成部分,每年通過土壤呼吸向大氣排放的CO2達640 g·m-2[4],且農(nóng)田生態(tài)系統(tǒng)本身易受到人類活動的干擾,不同的耕作方式、灌溉方式、施肥方法都會引起農(nóng)田土壤呼吸的變化[5-6],2009年中國因農(nóng)田利用所引起的碳排放量與 1999年相比增幅達93.9%[7]。因此,研究農(nóng)田生態(tài)系統(tǒng)土壤呼吸對農(nóng)業(yè)管理和氣候變化的響應具有十分重要的意義[8]?!厩叭搜芯窟M展】農(nóng)田管理措施的合理性直接影響農(nóng)田生態(tài)系統(tǒng)與大氣之間的碳循環(huán),合理的農(nóng)田管理措施能夠提高農(nóng)田土壤碳儲量,降低土壤碳排放[9-10]。已有的研究多集中于施肥方式,耕作方式等因素對土壤呼吸及土壤水熱條件的影響[11-14],或研究土壤呼吸的季節(jié)變化及日變化[15-16]。一般認為土壤呼吸由于受到土壤溫度及濕度的調(diào)控而存在一定的日變化和季節(jié)變化規(guī)律[17-18],不同的農(nóng)田管理措施對土壤呼吸產(chǎn)生明顯的影響,金皖豫等[19]、李銀坤等[20]研究表明農(nóng)田土壤呼吸速率隨施用氮肥量的增加而增加[19-20]。齊玉春等[21]認為一定范圍的灌溉能夠提高土壤碳排放,并影響土壤呼吸的溫度敏感性。YONEMURA等[22]、SCHWEN等[23]的研究結(jié)果表明免耕與傳統(tǒng)耕作相比,土壤的呼吸速率明顯降低。LIU等[24]、LAL等[25]的研究也指出免耕和秸稈還田能夠有效降低農(nóng)田土壤的碳排放,提高農(nóng)田生態(tài)的系統(tǒng)固碳減排效應?!颈狙芯壳腥朦c】前人對于免耕及秸稈還田對土壤呼吸影響的研究周期多為某種作物的 1—2個生長季或生育期內(nèi)的某一段時期[26-28],對長期免耕及秸稈還田對土壤呼吸影響的研究較少。輪作能夠最大限度地減緩土壤侵蝕,提高水分利用效率和土壤固碳潛力[29]?!緮M解決的關鍵問題】本研究以陜西關中地區(qū)常見的冬小麥-夏玉米輪作模式為研究對象。通過分析比較不同年份間免耕和秸稈還田條件下土壤呼吸的變化趨勢和土壤呼吸對土壤水熱因子的敏感性,結(jié)合作物產(chǎn)量,確定合理的農(nóng)田管理措施,為陜西關中地區(qū)制訂科學有效的農(nóng)田土壤減排制度提供理論基礎和技術支撐。

    1 材料與方法

    1.1 研究區(qū)概況

    試驗在西北農(nóng)林科技大學校內(nèi)試驗地進行。該地位于陜西省關中平原西部(E108°07′,N34°12′),海拔520 m,為大陸性暖溫帶季風區(qū)氣候,平均氣溫12.9℃,年降水量約為660 mm,主要集中在每年的7—9月,屬典型的一年兩熟農(nóng)作區(qū)。試驗土壤為土,試驗地土壤基本理化性狀見表 1。試驗期間當?shù)貧鉁嘏c降雨情況見圖1。

    表1 試驗地土壤養(yǎng)分特性Table 1 Soil chemical properties of the tested field

    圖1 試驗區(qū)氣溫與降雨量變化Fig. 1 Variation of air temperature and precipitation in the study area

    1.2 試驗設計

    長期定位試驗處理開始于2011年10月,在2011 —2014年內(nèi)的連續(xù) 3季冬小麥-夏玉米輪作周期內(nèi)進行。冬小麥于每年10月上旬播種,第二年6月中旬收獲,夏玉米在上茬冬小麥收獲后立即播種,每年 10月初收獲,冬小麥供試品種為西農(nóng) 889;夏玉米供試品種為漯單9號。

    試驗在免耕條件下共設置6個處理:秸稈全量還田+施肥(S1F1)、秸稈全量還田+不施肥(S1F0)、秸稈半量還田+施肥(S1/2F1)、秸稈半量還田+不施肥(S1/2F0)、秸稈不還田+施肥(S0F1)、秸稈不還田+不施肥(S0F0),其中S0F0處理作為對照。每個處理的面積均為68.8 m2,相互間隔0.5 m(圖2),每個處理共設3個重復,秸稈還田處理在前茬作物收獲后,利用秸稈粉碎機直接將作物秸稈粉碎后還田,秸稈全量還田量約為9 000 kg·hm-2,秸稈半量還田量約為4 500 kg·hm-2;秸稈不還田處理在前茬作物收獲后,人工拔出所有根茬。冬小麥施肥量為375 kg·hm-2尿素+375 kg·hm-2磷酸二銨;夏玉米施肥量為375 kg·hm-2尿素,均作為基肥施用。每一季冬小麥冬灌一次,每一季夏玉米根據(jù)當年天氣情況至多灌溉一次,其余時間依靠外界降水維持作物正常生長所需。

    圖2 試驗安排Fig. 2 Experiment arrangement

    1.3 測定項目

    1.3.1 氣體測定 土壤呼吸速率采用 GXH-3010E1型便攜式紅外線分析器(華云分析儀器研究所有限公司)測定,該儀器外形尺寸為248 mm×185 mm× 85 mm,測量范圍為0—1.000%,線性度≤±2% F·S,重復性≤1% F·S,響應時間≤10 s,采樣方式為泵吸式連續(xù)測量。在冬小麥和夏玉米播種后,在每個小區(qū)中部分別放置3個PVC 腔室,腔室高10 cm,直徑為16 cm,埋置深度為5 cm,腔室頂部為開放結(jié)構,開放部分直徑為11.5 cm,確保腔室內(nèi)部與外界環(huán)境條件的一致性,3個腔室呈“三角形”排列,相互間隔1 m。測量時用橡膠軟管連接便攜式紅外線分析器和腔室,測量過程中在腔室的頂部放置1個小風扇,使腔室內(nèi)的氣體混合均勻,持續(xù)3 min,然后測量,每一個處理重復測量3次,根據(jù)LUO等[30]描述,上午中段時間(大約10點)的土壤CO2通量很接近日平均值,9:00—11:00的測量結(jié)果也優(yōu)于整個白天的測量結(jié)果,以此作為依據(jù),結(jié)合當?shù)貙嶋H氣候情況,整個周期內(nèi)土壤CO2通量的測量均在9:00—10:00完成。

    土壤呼吸速率計算公式:R=k(X2-X1)H/Δt式中:R為土壤呼吸速率(μmol·m-2·s-1);k為換算系數(shù),k=1.80(25℃,1個標準大氣壓);X1、X2 分別為測定時 CO2初始質(zhì)量分數(shù)和測定結(jié)束后即時質(zhì)量分數(shù)(%);H為容器埋入地下的深度(m);Δt 為測定時間變化(s)。

    土壤呼吸速率自2011年10月開始測量,其中冬小麥分別于苗期、分蘗期、越冬期(冬灌前測定)、返青期各測定一次,返青期后每15天測定一次;夏玉米自6月中旬播種后,每隔15 d測定一次,為避免降水影響,如遇降水,適當后延3—5 d測定。

    1.3.2 土壤溫度測定 采用曲管地溫計測定土壤溫度。冬小麥及夏玉米播種后,將曲管地溫計埋至距PVC腔室10 cm處,每個處理3套。每次土壤呼吸速率測定的同時讀取5、10、15 cm層土壤溫度值。

    1.3.3 土壤含水量測定 采用烘干稱量法測定土壤含水量,每次土壤呼吸速率測定的當天,用土鉆取 0 —30 cm土層樣品(每10 cm為一層),3次重復。

    1.3.4 土壤有機碳測定 采用重鉻酸鉀容量法測定土壤耕作層有機碳含量,每季作物收獲后,分別采集0—10、10—20、20—30 cm的土層的土樣用于測定。

    1.3.5 作物產(chǎn)量測定 每季作物收獲期間,各處理以“V”字形隨機選樣,選樣面積為1 m2,重復3次,晾曬后,脫粒稱重。

    1.3.6 土壤呼吸累計釋放量測定 作物整個生育期內(nèi)的土壤呼吸累計釋放量的計算方法參考涂純等[14]、ZHAI等[31]研究。

    土壤呼吸累計釋放量計算公式:

    式中,Ra為土壤呼吸累計釋放量(t CO2·hm-2),Ri為第 i次測量的土壤呼吸速率(μmol·m-2·s-1),Ri+1為間隔 n天測量的土壤呼吸速率(μmol·m-2·s-1),n為相鄰兩次測量間隔天數(shù)(d),3600×24為s-1轉(zhuǎn)化為d-1,44為CO2的摩爾質(zhì)量,10-8為μg CO2·m-2轉(zhuǎn)化為t CO2·hm-2。

    1.4 數(shù)據(jù)分析方法及軟件

    土壤呼吸速率R對溫度T的敏感性采用指數(shù)函數(shù)進行非線性回歸分析,函數(shù)表達式為:R=aebT。

    土壤呼吸速率R對土壤含水量W的敏感性采用指數(shù)函數(shù)進行非線性回歸分析,函數(shù)表達式為:R=eaW2+bW+c。

    所有數(shù)據(jù)及圖表均采用Excel 2010、SPSS 20.0、Origin 9.1軟件進行分析處理,Duncan新復極差法進行方差分析。

    2 結(jié)果

    2.1 不同處理土壤呼吸的動態(tài)變化

    連續(xù)3年的測量結(jié)果表明,在整個輪作期間,不同作物土壤呼吸率的變化趨勢具有明顯的差異性,即在冬小麥生育期內(nèi),各處理土壤呼吸速率均呈先下降后升高再下降的趨勢;而在夏玉米生育期內(nèi),各處理土壤呼吸速率均表現(xiàn)為先升高后下降的趨勢(圖3)。自冬小麥苗期10月下旬(播種后10 d)至冬小麥分蘗期11月下旬(播種后50 d),隨著外界溫度的下降,土壤微生物活性降低,土壤呼吸速率不斷下降,在冬小麥越冬期(播種后95 d),各處理的土壤呼吸速率達到最低值。隨著外界溫度的持續(xù)升高,土壤呼吸速率逐漸升高,至冬小麥孕穗灌漿期間(播種后200—215 d),土壤呼吸速率達到峰值,隨后各處理的土壤呼吸速率又出現(xiàn)下降趨勢。夏玉米播種后,由于外界溫度較高,土壤呼吸速率自夏玉米苗期6月下旬逐漸升高,至8月上旬開花期(播種后60 d)達到峰值,隨后各處理的土壤呼吸速率開始下降。

    圖3 不同處理下冬小麥-夏玉米輪作系統(tǒng)土壤呼吸速率變化Fig. 3 Changes in soil respiration rate from the field of winter wheat-summer maize rotation in different treatments

    在整個研究周期內(nèi),冬小麥的土壤呼吸速率略低于夏玉米(P<0.01),2011—2014年間冬小麥全生育期內(nèi)的土壤呼吸速率范圍分別為 0.241—6.265、0.460—6.113、0.390—5.265 μmol·m-2·s-1;夏玉米全生育期內(nèi)的土壤呼吸速率范圍分別為 1.419—8.538、0.885—7.004、1.125—7.644 μmol·m-2·s-1。將不同作物同一生育期內(nèi)的土壤呼吸速率進行平均,各處理與對照間的差異均極顯著(P<0.01),各處理在3季冬小麥的土壤呼吸平均速率均表現(xiàn)為 S1F1>S1/2F1>S1/2F0>S0F1>S1F0>CK,說明秸稈還田和施肥均能夠提高土壤呼吸速率。不同處理的土壤呼吸平均速率第 1季冬小麥生育期內(nèi),較對照分別增加 21.7%、16.2%、11.9%、11.3%、7.3%;在第2季冬小麥生育期內(nèi),較對照分別增加28.4%、17.7%、11.7%、10.2%、5.9%;在第 3季冬小麥生育期內(nèi),較對照分別增加20%、14.4%、13.3%、8.8%、7.0%。3季夏玉米生育期內(nèi)各處理土壤呼吸平均速率間差異也均呈現(xiàn)出與冬小麥相似的規(guī)律。第1季夏玉米(2012年6月至2012 年10月),不同處理生育期土壤呼吸平均速率較對照分別增加 16.9%、13.9%、9.2%、7.5%、7.0%;第 2季夏玉米不同處理較對照分別增加 21.7%、16.7%、14.6%、4.5%、4.4%;第3季夏玉米不同處理較對照分別增加 9.0%、8.3%、7.4%、7.3%、3.1%??梢钥闯?,秸稈還田配合施肥處理的土壤呼吸速率增幅明顯高于單一施肥及單一秸稈還田的處理。

    冬小麥和夏玉米生育期內(nèi)各處理的土壤呼吸總量的變化規(guī)律與土壤呼吸平均速率相同,均表現(xiàn)為S1F1 >S1/2F1>S1/2F0>S0F1>S1F0>CK(表2),在施肥條件下,不同秸稈還田處理的土壤呼吸總量表現(xiàn)為秸稈全量還田>秸稈半量還田>秸稈不還田,但在不施肥條件下則表現(xiàn)為秸稈半量還田>秸稈全量還田>秸稈不還田。對比3個連續(xù)的生育期,發(fā)現(xiàn)各處理生育期的土壤呼吸總量有逐年下降的趨勢,冬小麥各處理的土壤呼吸總量均表現(xiàn)為第1季>第3季>第2季,后2季冬小麥的土壤呼吸總量均與第1季冬小麥的土壤呼吸總量差異極顯著(P<0.01),但后2季冬小麥之間的差異不顯著。各處理第2季冬小麥的土壤呼吸總量相比第1季冬小麥的降幅為14.2%—18.5%;第3季冬小麥的土壤呼吸總量相比第1季冬小麥的降幅為13.3%—16.3%。夏玉米各處理在不同年份間的土壤呼吸總量變化趨勢與冬小麥相同,且3季夏玉米相互間的差異均極顯著(P<0.01),各處理第2季夏玉米的土壤呼吸總量較第1季夏玉米降低22.8%—29.6%;第3季夏玉米的土壤呼吸總量較第 1季夏玉米降低11.0%—17.4%。

    表2 不同作物生育期內(nèi)各處理土壤呼吸總量Table 2 Different treatments cumulative soil respiration in different crops growth period

    2.2 不同作物土壤有機碳含量變化

    試驗結(jié)果顯示土壤有機碳含量隨土壤深度的增加而降低。CK處理的土壤有機碳含量與試驗開始之前(0 —10 cm 8.78 g·kg-1;10—20 cm 6.76 g·kg-1;20—30 cm 4.94 g·kg-1)相比有降低的趨勢,且這種趨勢隨著試驗周期的推移而不斷增大。這表明在當?shù)氐臍夂驐l件下,僅依靠免耕很難保持土壤有機碳的含量。單一施用氮肥盡管也能夠提高土壤有機碳的含量,但增幅十分有限。根據(jù)表3的結(jié)果顯示,每一季冬小麥收獲后,各處理同一土層深度的有機碳含量均依次為 S1F1>S1/2F1> S1F0>S1/2F0>S0F1>CK。表現(xiàn)為秸稈全量還田>秸稈半量還田>不還田,其中秸稈還田處理的土壤有機碳含量極顯著高于不還田處理(P<0.01)。不同年份間,相同處理的土壤有機碳含量逐年增加,但彼此間的差異不顯著。夏玉米收獲后的土壤有機碳含量低于冬小麥收獲后的含量(表4),但有機碳含量的變化趨勢與冬小麥相近,表現(xiàn)為秸稈還田處理高于不還田處理,且除對照外的各處理的有機碳含量逐年增加。綜合整個試驗周期內(nèi)土壤有機碳含量的變化趨勢,可以看出秸稈還田相比不還田更有利于土壤有機碳的積累。

    表3 冬小麥收獲后不同處理土壤有機碳含量變化Table 3 The change of soil organic carbon in difference treatment after winter wheat harvest

    表4 夏玉米收獲后不同處理土壤有機碳含量變化Table 4 The change of soil organic carbon in difference treatment after summer maize harvest

    2.3 土壤溫度與土壤呼吸速率的關系

    整個研究周期內(nèi)的土壤溫度具有明顯的差異性,其變化趨勢與每月平均氣溫的變化趨勢相似,不同處理在同一生育期內(nèi)的土壤溫度變化趨勢相近。冬小麥不同生育期內(nèi)的各處理土壤溫度均呈現(xiàn)先下降后升高的趨勢,而夏玉米不同生育期內(nèi)的各處理土壤溫度變化趨勢則表現(xiàn)為先升高后下降(圖 4)。各處理各生育期土壤平均溫度均隨土壤深度的增加而降低,其中,10和15 cm土層各生育期土壤平均溫度與5 cm的土層平均溫度均差異顯著(P<0.05),但二者之間的差異不顯著。4個秸稈還田處理在3季冬小麥越冬期5 cm深度的土壤溫度均高于2個不還田處理,說明一定量的秸稈還田有利于提高冬季表層土壤的溫度。夏玉米生育期內(nèi),土壤溫度持續(xù)升高,抽雄期前后(夏玉米播種45 d)土壤溫度達到峰值,之后隨外界氣溫的降低,逐漸下降。

    從連續(xù)3年的測定結(jié)果來看,除冬小麥越冬期外,秸稈還田處理對冬小麥及夏玉米其他生育時期的不同深度的土壤溫度均表現(xiàn)為降溫效應,這使得同一生育期內(nèi),不同深度下的各處理生育期土壤平均溫度均表現(xiàn)為CK>S0F1>S1/2F1>S1F1>S1F0>S1/2F0,CK處理對不同秸稈還田處理的差異均表現(xiàn)為極顯著(P <0.01),但不同秸稈還田處理間的差異不顯著,整個研究周期內(nèi),S1/2F1、S1F1、S1F0、S1/2F0處理的生育期土壤平均溫度與CK相比,最大溫差分別達1℃、1.1℃、1.4℃、2℃。不同冬小麥生育期之間,各處理0—15 cm的土層平均溫度均依次為第2季>第3季>第1季,第2季和第3季的土壤平均溫度均與第1季的差異極顯著(P<0.01);夏玉米不同生育期之間,各處理 0—15 cm的土層平均溫度表現(xiàn)為逐年上升,這一現(xiàn)象與研究周期內(nèi)每年6—9月平均氣溫的變化趨勢相近,第3季和第2季的土壤平均溫度均與第1季的差異極顯著(P<0.01),但第3季和第2季彼此間的土壤平均溫度差異不顯著。

    采用指數(shù)函數(shù)(R=aebT)擬合不同輪作周期內(nèi)土壤呼吸速率(R)和土壤溫度(T),回歸分析結(jié)果表明:不同處理土壤溫度可以依次解釋 2011.10—2012.10、2012.10—2013.10、2013.10—2014.9周期土壤呼吸速率變化的42.2%—59.8%、32.5%—56.7%、53.7%—60.4%(表5)。土壤溫度與土壤呼吸速率的擬合度隨土壤深度的增加而降低,5 cm土層溫度與土壤呼吸的相擬合度最高,可以依次解釋3個輪作周期土壤呼吸速率變化的48.5%—59.8%、37.8%—56.7%、54.0%—60.4%,而15 cm的土層溫度僅能依次解釋42.2%—49.0%、32.5% —46.7%、53.7%—57.3%的土壤呼吸速率。

    2.4 土壤濕度與土壤呼吸速率的關系

    不同處理在同一生育期內(nèi)的土壤含水量的變化趨勢相同,冬小麥不同生育期內(nèi)的各處理土壤含水量均呈現(xiàn)先下降后升高的趨勢,夏玉米生育期內(nèi)的各處理土壤含水量受到外界氣溫和降雨量的影響較大,不同年份的夏玉米土壤含水量變化趨勢略有不同,第1季和第2季夏玉米的土壤含水量變化趨勢近似于先升高后下降,而第3季夏玉米土壤含水量變化趨勢為先降低后升高(圖 5)。整個試驗周期內(nèi),耕作層土壤的含水量隨土壤深度的增加而降低,其中0—10 cm土層的土壤含水量與10—20 cm和20—30 cm土層的土壤含水量均極顯著差異(P<0.01),但 10—20 cm和20—30 cm土層的土壤含水量間的差異不顯著。在冬小麥越冬期間,土壤含水量降至最低,越冬期后,隨著外界土壤降雨的增多,各處理的土壤含水量也逐漸升高。夏玉米生育期正值當?shù)赜昙竞驮缕骄鶜鉁刈罡邥r期,受降雨和高溫的影響,各處理土壤含水量出現(xiàn)明顯的差異,其中第1季和第2季夏玉米生育期內(nèi)各處理土壤含水量在玉米播種后基本呈上升趨勢,第3季夏玉米生育期內(nèi)由于受到高溫和少雨的影響,各處理土壤含水量在播種后均出現(xiàn)明顯下降,直到抽雄期后,各處理土壤含水量才開始逐漸升高。

    圖4 不同輪作周期內(nèi)不同處理土壤溫度的變化Fig. 4 Change in soil temperature from different treatments in different rotation period

    夏玉米的生育期平均含水量顯著高于冬小麥(P <0.01)。同一生育期內(nèi),不同深度下的各處理生育期的土壤平均含水量均表現(xiàn)為 S1F0>S1F1>S1/2F0 >S1/2F1>CK>S0F1,土壤含水量均依次為秸稈全量還田>秸稈半量還田>秸稈不還田,且不同秸稈還田處理的土壤含水量與對照間的均差異顯著(P<0.05)。說明了秸稈還田具有一定保墑效應,且這種效應隨秸稈還田量的增加而增強。3季冬小麥之間,各處理在不同土層深度下的平均含水量均表現(xiàn)為第2季>第1季>第3季,相鄰兩季冬小麥之間不同深度土壤的平均含水量的差異均不顯著;3季夏玉米之間,各處理在各土層深度下的生育期平均含水量均表現(xiàn)為逐年下降,相同處理20—30 cm土層生育期的平均土壤含水量在3季夏玉米之間彼此差異極顯著(P<0.01)。

    采用指數(shù)函數(shù)(R=eaW2+bW+c)擬合不同輪作周期內(nèi)土壤呼吸速率(R)和土壤含水量(W),回歸分析結(jié)果表明,不同處理土壤含水量可以依次解釋2011.10—2012.10、2012.10—2013.10、2013.10—2014.9周期土壤呼吸速率變化的 46.7%—82.5%、38.4%—62.9%、44.2%—62.5%(表6),不同土壤深度下,10 —20 cm土層的含水量與土壤呼吸的擬合度最高,可以依次解釋3個輪作周期土壤呼吸速率變化的65.9% —82.5%、38.6%—62.9%、50.4%—62.5%。在第3個輪作周期內(nèi),所有處理的r2均低于對照,產(chǎn)生這種現(xiàn)象的可能是由于2014年夏季反常的高溫和降水,導致生育期土壤含水量差異較大。

    2.5 不同處理冬小麥及夏玉米產(chǎn)量的比較

    如表 7所示,5種處理的作物產(chǎn)量均高于對照(CK)。同一年份內(nèi),不同處理冬小麥和夏玉米產(chǎn)量均表現(xiàn)為 S1F1>S1F0>S1/2F1>S0F1>S1/2F0>S0F0(CK),且CK處理與其余5種處理差異極顯著(P<0.01),在整個試驗周期內(nèi),各處理的冬小麥產(chǎn)量逐年增加,且不同年份間的差異極顯著(P<0.01),其中 S1F1處理的冬小麥產(chǎn)量增幅最大,2014年較2012年相比,增產(chǎn)28.47%,S1F0處理次之。各處理2013年的夏玉米產(chǎn)量與2012年相比有所增加,增產(chǎn)10.63%—17.35%,且差異極顯著。但受極端炎熱天氣的影響,各處理 2014年夏玉米產(chǎn)量極顯著低于前兩年。

    圖5 不同輪作周期內(nèi)不同處理土壤含水量的變化Fig. 5 Change in soil moisture from different treatments in different rotation period

    表5 不同處理土壤呼吸速率(R)與土壤溫度(T)擬合方程及擬合度(r2)Table 5 Fitted equation of soil respiration rate (R) with soil temperature (T)

    表6 不同處理土壤呼吸速率(R)與土壤含水量(W)擬合方程及擬合度(r2)Table 6 Fitted equation of soil respiration rate (R) with soil moisture (W)

    表7 2012—2014 年冬小麥-夏玉米輪作體系內(nèi)不同處理的作物產(chǎn)量Table 7 Yield of crop in winter wheat-summer maize rotation system under different treatment in 2012-2014

    3 討論

    耕作方式能夠直接作用于土壤,通過改變土壤結(jié)構,進而影響土壤理化性質(zhì)。旋耕和翻耕對土壤擾動較大,而免耕可以極大程度的減小對土壤的擾動,有研究表明,免耕有效減少了土壤有機碳的損失,其土壤呼吸速率明顯低于翻耕和旋耕,從而降低了農(nóng)田土壤碳排放[32]。秸稈還田或施用氮肥作為常見的耕作管理措施,不僅能夠為土壤提供額外的養(yǎng)分,同時也會對土壤微生物的活性產(chǎn)生影響,促進土壤有機質(zhì)的礦化,進而引起土壤呼吸的變化。本試驗的結(jié)果表明在單季作物的生育期內(nèi),秸稈還田能夠提高土壤呼吸的平均速率和生育期呼吸總量,這與李瑋等[33]、于愛忠等[34]的研究結(jié)果一致,但比較連續(xù)多季作物生育期的觀測結(jié)果,發(fā)現(xiàn)免耕秸稈還田土壤的呼吸平均速率和呼吸總量呈現(xiàn)下降趨勢,這與涂純等[14]在陜西長武縣的研究結(jié)果相似,其中第2季冬小麥的土壤呼吸總量較第1季冬小麥降低14.2%—18.5%;第3季冬小麥的土壤呼吸總量相比第1季降低了13.3%—16.3%。第2季夏玉米的土壤呼吸總量較第1季夏玉米降低22.8% —29.6%;第3季夏玉米的土壤呼吸總量較第1季夏玉米降低11.0%—17.4%??赡艿脑蚴请S著免耕時間的推移,土壤團聚體由小變大,小團聚體逐漸減少,使土壤逐漸緊實,還田后的秸稈與土壤的接觸減少,分解速率減慢,土壤微生物活動所需的底物減少,進而導致土壤呼吸速率和呼吸總量的降低。在研究周期內(nèi),秸稈還田結(jié)合施肥的處理的生育期土壤平均速率和土壤呼吸總量均高于單一秸稈還田和單一施肥處理,這是由于向秸稈還田的基礎上額外施用氮肥后,土壤C/N降低,更有利于促進微生物的增殖及分解更多的有機質(zhì),進而增加了土壤有機質(zhì)中碳的分解與釋放[35]。本研究的結(jié)果還顯示增施氮肥會改變不同秸稈還田量對土壤呼吸總量的影響,在施肥條件下,不同秸稈還田處理的土壤呼吸總量表現(xiàn)為秸稈全量還田>秸稈半量還田>秸稈不還田,但在不施肥條件下則表現(xiàn)為秸稈半量還田>秸稈全量還田>秸稈不還田,這一方面可能是由于增施氮肥能在一定程度上影響秸稈的腐解狀況,加速秸稈的分解[36],加快了土壤中有機質(zhì)含量的積累,提高了有機質(zhì)分解釋放的CO2對土壤呼吸的貢獻。另一個方面可能是由于秸稈在分解過程中消耗一定的氮,改變了土壤原有的C/N。

    有研究表明,秸稈還田具有平衡和改善耕層土壤溫度狀況,當土壤溫度較低時具有保溫作用,當土壤溫度較高時具有降溫作用[37]。本研究的結(jié)果與前人研究結(jié)果相似,當土壤溫度低于 0℃,秸稈還田具有提高土壤溫度的效應,當土壤溫度高于 0℃。尤其是在夏玉米生育期內(nèi),秸稈還田具有明顯的降溫作用。就整個生育期的土壤平均溫度而言,相比對照,不同水平的秸稈還田表現(xiàn)為降溫效應,這與黃高寶等[38]、杜新艷等[39]研究結(jié)果相似。前人研究結(jié)果一致認為秸稈還田具有明顯的保墑效應,能夠抑制表層土壤水分的蒸發(fā),提高農(nóng)田土壤的水分利用率[40-41]。本研究結(jié)果表明:不同秸稈還田處理與裸地相比,生育期土壤含水量均有不同程度的提高,且生育期土壤含水量隨秸稈還田量的增加而增加,特別是在持續(xù)高溫的和葉面積指數(shù)較小的夏玉米生長前期,秸稈還田處理的土壤含水量明顯高于秸稈不還田處理,有效減小了株間土壤水分的無效蒸發(fā),這與前人的研究結(jié)果相一致[42]。

    土壤呼吸受到多種因子的調(diào)控,研究表明,土壤溫度與土壤濕度是影響土壤呼吸的主要因子[43-44]。土壤呼吸的季節(jié)變化主要受土壤溫度與土壤濕度的調(diào)控,表現(xiàn)為夏季最高冬季最低,在本研究期間,各生育期內(nèi)土壤呼吸的季節(jié)變化趨勢與土壤溫度的季節(jié)變化趨勢相近,但與土壤含水量的季節(jié)變化趨勢有所差異,這是因為土壤含水量的季節(jié)變化主要受外界降水影響,不同年份降雨量的不同導致土壤含水量季節(jié)變化的差異較大。土壤溫度可以依次解釋連續(xù)3個輪作周期土壤呼吸速率變化的 42.2%—59.8%、32.5%—56.7%、53.7%—60.4%,土壤濕度可以依次解釋連續(xù)3個輪作周期土壤呼吸速率變化的 46.7%—82.5%、38.4%—62.9%、44.2%—62.5%,其中土壤呼吸分別與地下5 cm深度的土壤溫度和地下10—20 cm深度的土壤含水量的相關性最好,這與前人的研究結(jié)果相近[45-46]。

    維持農(nóng)田土壤中較高的碳儲量目前已逐漸成為農(nóng)業(yè)生產(chǎn)和環(huán)境保護過程中重要的一環(huán)[47],但農(nóng)田土壤的有機碳累積受到氣候、土壤性質(zhì)及耕作管理措施等多種因素的影響[48],合理的耕作措施有利于提高土壤有機碳的輸入,減緩有機碳的分解,從而提高土壤有機碳的含量。本研究周期內(nèi)的結(jié)果顯示,在免耕條件下,秸稈還田處理的土壤有機碳含量高于秸稈不還田的處理,這與HAN等[49]及LENKA等[50]的研究結(jié)果一致。在整個試驗周期內(nèi),除對照外的其他處理不同深度的土壤有機碳含量均逐年增加,其中秸稈還田結(jié)合施肥處理的土壤有機碳含量增幅最明顯,單一秸稈還田處理的增幅次之,且增幅隨秸稈還田量的增多而增大,而單一施肥處理土壤有機碳含量的增幅最低。因此,說明合理的秸稈還田與施肥配比有利于土壤有機碳的積累。但由于同一處理不同年份間土壤有機碳含量的差異不顯著,因此,對于免耕秸稈還田下土壤有機碳含量的變化趨勢還需要進一步的研究。

    前人的研究結(jié)果表明,秸稈還田有助于提高作物產(chǎn)量[51-52],本文進一步驗證了前人的有關結(jié)論。整個輪作周期內(nèi),隨著免耕秸稈還田年限的增加,冬小麥產(chǎn)量逐年增加;第2季夏玉米的產(chǎn)量顯著高于第1季,但在第3季夏玉米生育期內(nèi)的連續(xù)高溫少雨天氣導致當年玉米減產(chǎn),這說明在無反常天氣的情況下,免耕秸稈還田對于作物產(chǎn)量的增加具有長期的效應。

    4 結(jié)論

    在研究周期內(nèi),不同作物生育期內(nèi)土壤呼吸速率的變化趨勢具有明顯的差異性。在單個生育期內(nèi),秸稈還田處理的土壤呼吸平均速率和呼吸總量均高于對照,而各處理在同種作物不同生育期間的土壤生育期呼吸總量有逐年下降的趨勢。土壤呼吸速率受到土壤溫度和土壤含水量的調(diào)控,整個研究周期內(nèi),土壤溫度能夠解釋土壤呼吸速率變化的32.5%—60.4%;土壤含水量能夠解釋土壤呼吸速率變化的38.4%—82.5%;土壤溫度及含水量受外界氣溫和降雨的影響,表現(xiàn)出較大的年際間差異,不同秸稈還田處理在冬小麥越冬期內(nèi)的土壤溫度高于對照,但生育期土壤平均溫度低于對照,不同秸稈還田處理的生育期土壤平均含水量高于對照,且秸稈還田量越高對土壤蒸發(fā)的抑制效果越好。長期免耕秸稈還田有利于提高土壤耕作層的有機碳含量,且有機碳含量的增幅隨秸稈還田量的增加而提高。總體而言,長期免耕秸稈還田有利于降低農(nóng)田土壤的呼吸速率和生育期碳排放,提高土壤水分的利用效率,降低作物株間土壤水分的無效蒸發(fā),提高冬季的土壤溫度,確保冬小麥安全越冬。另外,免耕條件下秸稈還田配合化肥的施用還有助于提高作物的產(chǎn)量。綜合評價各秸稈還田處理的土壤呼吸速率、土壤水熱狀況、作物產(chǎn)量等方面后可知,秸稈全量還田+不施肥(S1F0)處理在降低農(nóng)田土壤碳排放的同時提高了耕作層土壤有機碳含量及作物的產(chǎn)量,在提高農(nóng)田土壤水分利用率及冬季土壤溫度等方面也表現(xiàn)出較高的水平。

    References

    [1] BOND-LAMBERTY B, THOMSON A. A global database of soil respiration data. Biogeosciences, 2010, 7(6): 1915-1926.

    [2] HIBBARD K A, LAW B E, REICHSTEIN M. An analysis of soil respiration across northern hemisphere temperate ecosystems. Biogeochemistry, 2005, 73(1): 29-70.

    [3] SCHLESINGER W H, ANDREWS J A. Soil respiration and the global carbon cycle. Biogeochemistry, 2000, 48(1): 7-20.

    [4] 劉德輝, 陶于祥. 土壤、農(nóng)業(yè)與全球氣候變化. 火山地質(zhì)與礦產(chǎn),2000, 21(4): 290-295. LIU D H, TAO Y X. Soil, agriculture and global climate change. Volcanology & Mineral Resource, 2000, 21(4): 290-295. (in Chinese)

    [5] RAICH J W, TUFEKCIOGUL A. Vegetation and soil respiration:Correlations and controls. Biogeochemistry, 2000, 48(1): 71-90.

    [6] LAL R. Soil carbon sequestration impacts on global climate change and food security. Science, 2004, 304(5677): 1623-1627.

    [7] 段華平, 張悅, 趙建波, 卞新民. 中國農(nóng)田生態(tài)系統(tǒng)的碳足跡分析.水土保持學報, 2011, 25(5): 203-208. DUAN H P, ZHANG Y, ZHAO J B, BIAN X M. Carbon footprint analysis of farmland ecosystem in China. Journal of Soil and Water Conservation, 2011, 25(5): 203-208. (in Chinese)

    [8] WAN S, NORBY R J, LEDFORD J, WELTZIN F J. Responses of soil respiration to elevated CO2, air warming, and changing soil water availability in a model old-field grassland. Global Change Biology,2007, 13(11): 2411-2424.

    [9] WISE M, CALVIN K, THOMSON A, CLARKE L, BONDLAMBERTY B, SANDS R, SMITH S J, JANETOS A, EDMONDS J. Implications of limiting CO2concentrations for land use and energy. Science, 2009, 324(5931): 1183-1186.

    [10] 張前兵, 楊玲, 王進, 羅宏海, 張亞黎, 張旺鋒. 干旱區(qū)不同灌溉方式及施肥措施對棉田土壤呼吸及各組分貢獻的影響. 中國農(nóng)業(yè)科學, 2012, 45(12): 2420-2430. ZHANG Q B, YANG L, WANG J, LUO H H, ZHANG Y L, ZHANG W F. Effects of different irrigation methods and fertilization measures on soil respiration and its component contributions in cotton field in arid region. Scientia Agricultura Sinica, 2012, 45(12): 2420-2430. (in Chinese)

    [11] 張俊麗, 廖允成, 曾愛, 劉楊, 林杰, 王永平. 不同施氮水平下旱作玉米田土壤呼吸速率與土壤水熱關系. 農(nóng)業(yè)環(huán)境科學學報, 2013,32(7): 1382-1388. ZHANG J L, LIAO Y C, ZENG A, LIU Y, LIN J, WANG Y P. Effects of different levels of N fertilizer on soil respiration, and its relation to soil moisture and soil temperature under rainfed land of summer maize. Journal of Agro-Environment Science, 2013, 32(7): 1382-1388. (in Chinese)

    [12] DING W X, YU H Y, CAI Z C, HAN F X, XU Z H. Responses of soil respiration to N fertilization in a loamy soil under maize cultivation. Geoderma, 2010, 155(3): 381-389.

    [13] MORARU P I, RUSU T. Effect of tillage systems on soil moisture,soil temperature, soil respiration and production of wheat, maize and soybean crops. Journal of Food Agriculture & Environment, 2012,10(2): 445-448.

    [14] 涂純, 王俊, 官情, 劉文兆. 秸稈覆蓋對旱作冬小麥農(nóng)田土壤呼吸、作物產(chǎn)量及經(jīng)濟環(huán)境效益的影響. 中國生態(tài)農(nóng)業(yè)學報, 2013,21(8): 931-937. TU C, WANG J, GUAN Q, LIU W Z. Effect of straw mulching on soil respiration, crop yield, economy-environment benefit in rained winter wheat fields. Chinese Journal of Eco-Agriculture, 2013, 21(8):931-937. (in Chinese)

    [15] 周小平, 王效科, 張紅星, 龐軍柱, 任玉芬, 侯培強, 劉文兆. 黃土高原小麥田土壤呼吸季節(jié)和年際變化. 生態(tài)學報, 2013(23):7525-7536. ZHOU X P, WANG X K, ZHANG H X, PANG J Z, REN Y F, HOU P Q, LIU W Z. Seasonal and interannual variability soil respiration in wheat field of the Loess Plateau, China. Acta Ecologica Sinica,2013(23): 7525-7536. (in Chinese)

    [16] 孟凡喬, 關桂紅, 張慶忠, 史雅娟, 屈波, 況星. 華北高產(chǎn)農(nóng)田長期不同耕作方式下土壤呼吸及其季節(jié)變化規(guī)律. 環(huán)境科學學報,2006, 26(6): 992-999. MENG F Q, GUAN G H, ZHANG Q Z, SHI Y J, QU B, KUANG X. Seasonal variation in soil respiration under different long term cultivation practices on high yield farmland in the North China Plain. Acta Scientiae Circumstantiae, 2006, 26(6): 992-999. (in Chinese)

    [17] XU M, QI Y. Soil-surface CO2efflux and its spatial and temporal variations in a young ponderosa pine plantation in northern California. Global Change Biology, 2001, 7(6): 667-677.

    [18] DAVIDSON E A, VERCHOT L V, CATT?NIO J H, ACKERMAN I L, CARVALHO J E M. Effects of soil water content on soil respiration in forests and cattle pastures of eastern Amazonia. Biogeochemistry, 2000, 48(1): 53-69.

    [19] 金皖豫, 李銘, 何楊輝, 杜正剛, 邵鈞炯, 張國棟, 周靈燕, 周旭輝.不同施氮水平對冬小麥生長期土壤呼吸的影響. 植物生態(tài)學報,2015, 39(3): 249-257. JIN W Y, LI M, HE Y H, DU Z G, SHAO J J, ZHANG G D, ZHOU L Y, ZHOU X H. Effects of different levels of nitrogen fertilization on soil respiration during growing season in winter wheat (Triticum aestivum). Chinese Journal of Plant Ecology, 2015, 39(3): 249-257. (in Chinese)

    [20] 李銀坤, 陳敏鵬, 夏旭, 梅旭榮, 李昊儒, 郝衛(wèi)平. 不同氮水平下夏玉米農(nóng)田土壤呼吸動態(tài)變化及碳平衡研究. 生態(tài)環(huán)境學報,2013(1): 18-24. LI Y K, CHEN M P, XIA X, MEI X R, LI H R, HAO W P. Dynamic changes of soil respiration and carbon balance in summer maize field under different nitrogen levels. Ecology and Environmental Sciences,2013(1): 18-24. (in Chinese)

    [21] 齊玉春, 郭樹芳, 董云社, 彭琴, 賈軍強, 曹叢叢, 孫良杰, 閆鐘清,賀云龍. 灌溉對農(nóng)田溫室效應貢獻及土壤碳儲量影響研究進展.中國農(nóng)業(yè)科學, 2014, 47(9): 1764-1773. QI Y C, GUO S F, DONG Y S, PENG Q, JIA J Q, CAO C C, SUN L J,YAN Z Q, HE Y L. Advances in research on the effects of irrigation on the greenhouse gases emission and soil carbon sequestration in agro-ecosystem. Scientia Agricultura Sinica, 2014, 47(9): 1764-1773. (in Chinese)

    [22] YONEMURA S, NOUCHI I, NISHIMURA S, SAKURAI G,TOGAMI K, YAGI K. Soil respiration, N2O, and CH4emissions from an andisol under conventional-tillage and no-tillage cultivation for 4 years. Biology and Fertility of Soils, 2014, 50(1): 63-74.

    [23] SCHWEN A, JEITLER E, B?TTCHER J. Spatial and temporal variability of soil gas diffusivity, its scaling and relevance for soil respiration under different tillage. Geoderma, 2015, 259: 323-336.

    [24] LIU J M, HU L F, ZHANG A J. Research progress in greenhouse effect caused by conservation tillage. Chinese Agricultural Science Bulletin, 2006, 8(22): 246-249.

    [25] LAL R. Agricultural activities and the global carbon cycle. Nutrient Cycling in Agroecosystems. 2004, 70(2): 103-116.

    [26] 張宇, 張海林, 陳繼康, 陳阜. 耕作方式對冬小麥田土壤呼吸及各組分貢獻的影響. 中國農(nóng)業(yè)科學, 2009, 42(9): 3354-3360. ZHANG Y, ZHANG H L, CHEN J K, CHEN F. Tillage effects on soil respiration and contributions of its components in winter wheat field. Scientia Agricultura Sinica, 2009, 42(9): 3354-3360. (in Chinese)

    [27] 李昌珍, 張婷婷, 楊改河, 任廣鑫, 馮永忠. 秸稈覆蓋和施肥對關中灌區(qū)夏玉米生長后期土壤呼吸速率的影響. 生態(tài)環(huán)境學報, 2013,22(3): 411-416. LI C Z, ZHANG T T, YANG G H, REN G X, FENG Y Z. Effects of straw mulching and fertilization on the soil respiration rate in the late summer maize growth in Guanzhong irrigation district. Ecology and Environmental Sciences, 2013, 22(3): 411-416. (in Chinese)

    [28] SHI X, ZHANG X, YANG X, DRURY C F, MCLAUGHLIN N B,LIANG A Z, FAN R Q, JIA S X. Contribution of winter soil respiration to annual soil CO2emission in a mollisol under different tillage practices in northeast China. Global Biogeochemical Cycles,2012, 26(2).

    [29] JIANG J G, GUO S L, ZHANG Y J, LIU Q F, WANG R, WANG Z Q, LI N N, LI R J. Changes in temperature sensitivity of soil respiration in the phases of a three-year crop rotation system. Soil and Tillage Research, 2015(150): 139-146.

    [30] 駱亦其, 周旭輝. 土壤呼吸與環(huán)境. 北京: 高等教育出版社, 2007:92. LUO Y Q, ZHOU X H. Soil Respiration and the Envirnment. Beijing:Higher Education Press, 2007: 92. (in Chinese)

    [31] ZHAI L, LIU H, ZHANG J, HUANG J, WANG B W. Long-term application of organic manure and mineral fertilizer on N2O and CO2emissions in a red soil from cultivated maize-wheat rotation in China. Agricultural Sciences in China, 2011, 10(11): 1748-1757.

    [32] REICOSKY D C, LINDSTROM M J, SCHUMACHER T E, LOBB D E, MALO D D. Tillage-induced CO2loss across an eroded landscape. Soil and Tillage Research, 2005, 81(2): 183-194.

    [33] 李瑋, 張佳寶, 張叢志. 秸稈還田方式和氮肥類型對黃淮海平原夏玉米土壤呼吸的影響. 中國生態(tài)農(nóng)業(yè)學報, 2012, 20(7): 842-849. LI W, ZHANG J B, ZHANG C Z. Effects of straw incorporation and N fertilization on soil respiration during maize (Zea mays L.) growth in Huanghuaihai Plain. Chinese Journal of Eco-Agriculture, 2012,20(7): 842-849. (in Chinese)

    [34] 于愛忠, 黃高寶, 柴強. 不同耕作措施對西北綠洲灌區(qū)冬小麥農(nóng)田土壤呼吸的影響. 草業(yè)學報, 2012, 21(1): 273-278. YU A Z, HUANG G B, CHAI Q. Effects of different tillage measures on soil respiration of winter wheat in the irrigation area of Northwest China. Acta Prataculturae Sinica, 2012, 21(1): 273. (in Chinese)

    [35] 孫媛, 任廣鑫, 馮永忠, 張青, 李慧瑛, 楊改河. 秸稈還田和施氮對土壤水熱因子及呼吸速率的影響. 西北農(nóng)林科技大學學報(自然科學版), 2015, 43(3): 146-152. SUN Y, REN G X, FENG Y Z, ZHANG Q, LI H Y, YANG G H. Comprehensive influence of straw-returning and nitrogen fertilization on hydrothermal factors and soil respiration. Journal of Northwest A&F University (Natural Science Edition), 2015, 43(3): 146-152. (in Chinese)

    [36] 馬永良, 宇振榮, 江永紅, 羅維. 兩種還田模式下玉米秸稈分解速率的比較. 生態(tài)學雜志, 2002, 21(6): 68-70. MA Y L, YU Z R, JIANG Y H, LUO W. Comparison of decomposition rates of maize straw between two kinds of straw incorporation. Chinese Journal of Ecology, 2002, 21(6): 68-70. (in Chinese)

    [37] 付國占, 李潮海, 王俊忠, 王振林, 曹鴻鳴, 焦念元, 陳明燦. 殘茬覆蓋與耕作方式對土壤性狀及夏玉米水分利用效率的影響. 農(nóng)業(yè)工程學報, 2005, 21(1): 52-56. FU G Z, LI C H, WANG J Z, WANG Z L, CAO H M, JIAO N Y,CHEN M C. Effects of stubble mulch and tillage methods on soil properties and water use efficiency of summer maize. Transactions of the Chinese Society of Agricultural Engineering, 2005, 21(1): 52-56. (in Chinese)

    [38] 黃高寶, 李玲玲, 張仁陟, 蔡立群, KWONGYIN CHAN. 免耕秸稈覆蓋對旱作麥田土壤溫度的影響. 干旱地區(qū)農(nóng)業(yè)研究, 2006, 24(5):1-4. HUANG G B, LI L L, ZHANG R S, CAI L Q, KWONGYIN CHAN. No tillage effects on soil temperature of dryland wheat. Agricultural Research in the Arid Areas, 2006, 24(5): 1-4. (in Chinese)

    [39] 杜新艷, 楊路華, 脫云飛, 高惠嫣, 張振偉. 秸稈覆蓋對夏玉米農(nóng)田水分狀況、土壤溫度及生長發(fā)育的影響. 南水北調(diào)與水利科技,2006, 4(2): 24-26.DU X Y, YANG L H, TUO Y F, GAO H Y, ZHANG Z W. Moisture condition, soil temperature and growth condition in maize field mulched. South-to-North Water Transfers and Water Science & Technology, 2006, 4(2): 24-26. (in Chinese)

    [40] WANG X, JIA Z, LIANG L. Effect of straw incorporation on soil moisture, evapotranspiration, and rainfall-use efficiency of maize under dryland farming. Journal of Soil and Water Conservation, 2014,69(5): 449-455.

    [41] 王暉, 劉泉汝, 張圣勇, 沈加印, 趙丹丹, 于京平, 李全起. 秸稈覆蓋下超高產(chǎn)夏玉米農(nóng)田產(chǎn)量和土壤水分的動態(tài)變化. 水土保持學報, 2011, 25(5): 261-264. WANG H, LIU Q R, ZHANG S Y, SHEN J Y, ZHAO D D, YU J P, LI Q Q. Grain yield and soil water content of super-high-yield summer maize under straw mulching. Journal of Soil and Water Conservation,2011, 25(5): 261-264. (in Chinese)

    [42] 祿興麗, 廖允成. 保護性耕作對旱作夏玉米苗期土壤水熱及作物產(chǎn)量的影響. 土壤通報, 2014, 45(1): 147-150. LU X L, LIAO Y C. Conservation tillage in dryland soil water and heat of summer maize and crop yields. Chinese Journal of Soil Science, 2014, 45(1): 147-150. (in Chinese)

    [43] 劉爽, 嚴昌榮, 何文清, 劉勤. 不同耕作措施下旱作農(nóng)田土壤呼吸及其影響因素. 生態(tài)學報, 2010, 30(11): 2919-2924. LIU S, YAN C R, HE W Q, LIU Q. Soil respiration and it's affected factors under different tillage systems in dryland production systems. Acta Ecologica Sinica, 2010, 30(11): 2919-2924. (in Chinese)

    [44] WANG X G, ZHU B, WANG Y Q, ZHENG X H. Soil respiration and its sensitivity to temperature under different land use conditions. Acta Ecologica Sinica, 2007, 27(5): 1960-1968.

    [45] 張慶忠, 吳文良, 王明新, 周中仁, 陳淑峰. 秸稈還田和施氮對農(nóng)田土壤呼吸的影響. 生態(tài)學報, 2005, 25(11): 2883-2887. ZHANG Q Z, WU W L, WANG M X, ZHOU Z R, CHEN S F. The effects of crop residue amendment and N rate on soil respiration. Acta Ecologica Sinica, 2005, 25(11): 2883-2887. (in Chinese)

    [46] 孫小花, 張仁陟, 蔡立群, 陳強強. 不同耕作措施對黃土高原旱地土壤呼吸的影響. 應用生態(tài)學報, 2009(9): 2173-2180. SUN X H, ZHANG R S, CAI L Q, CHEN Q Q. Effects of different tillage measures on upland soil respiration in Loess Plateau. Chinese Journal of Applied Ecology, 2009(9): 2173-2180. (in Chinese)

    [47] ADHIKARI K R, DAHAL K R, CHEN Z S, TAN Y C, LAI J S. Rice-wheat cropping system: tillage, mulch, and nitrogen effects on soil carbon sequestration and crop productivity. Paddy and Water Environment, 2015. DOI 10.1007/s10333-015-0511-1.

    [48] 陳鮮妮, 岳西杰, 葛璽祖, 王旭東. 長期秸稈還田對塿土耕層土壤有機碳庫的影響. 自然資源學報, 2012, 27(1): 25-32. CHEN X N, YUE X J, GE X Z, WANG X D. Effect of long-term straw return on soil organic carbon pool in Lou soil. Journal of Resources, 2012, 27(1): 25-32. (in Chinese)

    [49] HAN H F, NING T Y, LI Z J, CAO H M. Soil respiration rate in summer maize field under different soil tillage and straw application. Maydica, 2014, 59(2): 185-190.

    [50] LENKA N K, LAL R. Soil aggregation and greenhouse gas flux after 15 years of wheat straw and fertilizer management in a no-till system. Soil and Tillage Research, 2013, 126: 78-89.

    [51] 武際, 郭熙盛, 魯劍巍, 王允清, 張曉玲, 許征宇. 連續(xù)秸稈覆蓋對土壤無機氮供應特征和作物產(chǎn)量的影響. 中國農(nóng)業(yè)科學, 2011,45(9): 1741-1749. WU J, GUO X S, LU J W, WANG Y Q, ZHANG X L, XU Z Y. Effects of continuous straw mulching on supply characteristics of soil inorganic nitrogen and crop yields. Scientia Agricultura Sinica, 2011,45(9): 1741-1749. (in Chinese)

    [52] 屈會娟, 李金才, 沈?qū)W善, 魏鳳珍, 武進東, 馬蓓. 秸稈全量還田對冬小麥不同小穗位和粒位結(jié)實粒數(shù)和粒重的影響. 中國農(nóng)業(yè)科學, 2010, 44(10): 2176-2183. QU H J, LI J C, CHEN X L, WEI F Z, WU J D, MA B. Effects of all straw returned to the field on grain number and grain weight at different spikelets and grain positions in winter wheat. Scientia Agricultura Sinica, 2010, 44(10): 2176-2183. (in Chinese)

    (責任編輯 楊鑫浩,李莉)

    Effects of Straw Returning to Field on Soil Respiration and Soil Water Heat in Winter Wheat - Summer Maize Rotation System Under No Tillage

    WANG Wei-yu1,2, QIAO Bo1,2, Kashif AKHTAR2, YUAN Shuai1,2, REN Guang-xin1,2, FENG Yong-zhong1,2
    (1College of Agronomy, Northwest A & F University, Yangling 712100, Shaanxi;2The Research Center of Recycle Agricultural Engineering and Technology of Shaanxi Province, Yangling 712100, Shaanxi)

    Abstract:【Objective】The purpose of the study was to investigate the effects of straw mulching on soil respiration, soiltemperature and soil moisture of winter wheat and summer maize rotation system in no tillage.【Method】From October 2011 to September 2014, a long-term location test of 6 different tillage treatments was set in Yangling, Shaanxi. The treatments used were total amount of straw returning to field+fertilization (S1F1), total amount of straw returning to field+no fertilization(S1F0), half the amount of straw to field+fertilization (S1/2F1), half the amount of straw to field+no fertilization(S1/2F0), without straw+fertilization (S0F1), without straw+no fertilization(S0F0), then the differences in soil respiration, soil temperature, soil moisture, crop yield and soil organic carbon of the soil plow layer under different treatments were measured and analyzed from 2011 to 2014. 【Result】During the growth period of winter wheat, soil respiration rate in each treatment showed an increase after the first downward trend and then a decrease. During the growth period of summer maize, the respiration rate in each treatment showed a trend of first increase and then decrease. During the same growth period, the following was S1F1>S1/2F1>S1/2F0>S0F1>S1F0>CK in the average rate and cumulative soil respiration. During different growth stages of the same crop, the cumulative soil respiration in different treatments showed a decreasing trend year by year. In the whole research term, the variation trend of soil temperature was similar with the trend of the monthly average temperature. The soil temperature variation trends of different treatments during the same growth period were similar. Simultaneously, the average soil temperatures during the same period of different treatments were decreased with the increase of soil depth. In winter, the soil temperature of different treatments was higher than that of CK, but the average soil temperature was lower than that of CK in a same growth period. The soil moisture decreased with the increase of soil depth, but the soil moisture was affected by rainfall, the change of soil moisture in different rotation cycles was larger, and the average water moisture of each treatment showed a trend of S1F0>S1F1>S1/2F0>S1/2F1>CK>S0F1 in a same growth period. Meanwhile, the soil moisture in different straw mulching treatments was significantly different from that of CK, respectively (P<0.05); 32.5%-60.4% variations in soil respiration rate could be recorded due to the soil temperature, and 38.4%-82.5% variations in it could be recorded due to the soil moisture. In different depths of soil, the correlation between the soil temperature and the soil respiration in 5 cm depth was the strongest, and the correlation between the soil moisture and the soil respiration at 10-20 cm was the highest. In the same year, the yields of winter wheat and summer maize in different treatments, respectively, showed a trend of S1F1 >S1F0>S1/2F1>S0F1>S1/2F0>CK. In the research cycle, the yield of wheat increased continuously in three years, and the yield of maize increased in the previous two years of three, but the yield in the third year significantly decreased because of the influence of extremely hot weather. After the harvest of single crop, soil organic carbon in different treatments at the same soil depth,respectively, demonstrated a trend of S1F1>S1/2F1>S1F0>S1/2F0>S0F1>CK. Meanwhile, the soil organic carbon in different treatments of straw mulching was increased continuously in the three years.【Conclusion】Long-term no tillage with straw returning can effectively reduce soil carbon emissions, improve soil water use efficiency and winter soil temperature, and increase crop yield and soil organic carbon. The effect of S1F0 treatment is the best among different treatments.

    Key words:no tillage; winter wheat-summer maize; rotation; straw returning; soil respiration; water and heat condition

    收稿日期:2015-10-16;接受日期:2016-03-06

    基金項目:陜西省科技統(tǒng)籌計劃(2015KTCL02-07,2016KTCL02-11)

    猜你喜歡
    免耕秸稈還田夏玉米
    小麥收割之后 如何種植夏玉米才能高產(chǎn)
    夏玉米高產(chǎn)的關鍵栽培技術措施
    一種能處理雜草的免耕玉米播種機
    小香蒜免耕種植技術
    玉米寬窄行免耕精量播種推廣技術
    玉米高產(chǎn)栽培秸稈還田保護性耕作技術措施
    赤壁市秸稈還田替代鉀肥在中稻上的應用效果研究
    論述秸稈還田對水稻病蟲害發(fā)生的影響
    大力推進聯(lián)耕聯(lián)種 構建新型農(nóng)業(yè)生產(chǎn)經(jīng)營模式
    美國土壤印痕裝置的耕地水土流失防治效果
    两人在一起打扑克的视频| 国产成人一区二区三区免费视频网站 | 久久热在线av| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区av电影网| 极品人妻少妇av视频| 又大又黄又爽视频免费| 黄色片一级片一级黄色片| 成人手机av| bbb黄色大片| 久久久精品国产亚洲av高清涩受| 波多野结衣一区麻豆| 日韩大码丰满熟妇| 亚洲熟女精品中文字幕| 别揉我奶头~嗯~啊~动态视频 | 麻豆av在线久日| 欧美精品一区二区大全| 日韩av不卡免费在线播放| 免费看十八禁软件| 99九九在线精品视频| 亚洲 欧美一区二区三区| 在线观看免费高清a一片| 99国产精品99久久久久| 中文字幕高清在线视频| 亚洲成人手机| 国产免费一区二区三区四区乱码| 91九色精品人成在线观看| 丰满少妇做爰视频| 我要看黄色一级片免费的| 欧美人与性动交α欧美软件| 一区二区三区四区激情视频| 久久久久久人人人人人| videosex国产| 婷婷丁香在线五月| 九色亚洲精品在线播放| videosex国产| 精品亚洲成国产av| 人妻 亚洲 视频| 人妻一区二区av| 亚洲精品美女久久av网站| 午夜免费观看性视频| 赤兔流量卡办理| 十八禁网站网址无遮挡| 人人妻人人澡人人爽人人夜夜| 人人妻人人澡人人爽人人夜夜| 国产精品av久久久久免费| 男女免费视频国产| 最新的欧美精品一区二区| 国产精品国产三级专区第一集| 在线观看免费日韩欧美大片| 一区在线观看完整版| 亚洲人成电影免费在线| 十分钟在线观看高清视频www| 国产一区二区在线观看av| 99热网站在线观看| 婷婷丁香在线五月| 国产高清videossex| 黄片播放在线免费| 亚洲av片天天在线观看| 下体分泌物呈黄色| 亚洲精品第二区| 久久久久久免费高清国产稀缺| 韩国高清视频一区二区三区| 国产爽快片一区二区三区| 午夜福利视频在线观看免费| 热99久久久久精品小说推荐| 19禁男女啪啪无遮挡网站| 波野结衣二区三区在线| 51午夜福利影视在线观看| 国产欧美日韩一区二区三 | 国产精品 欧美亚洲| 亚洲情色 制服丝袜| 亚洲色图综合在线观看| 午夜福利影视在线免费观看| 下体分泌物呈黄色| 婷婷成人精品国产| 日本五十路高清| 久久久精品免费免费高清| 亚洲成av片中文字幕在线观看| 大陆偷拍与自拍| 欧美激情高清一区二区三区| 婷婷色麻豆天堂久久| 国产精品成人在线| 日韩伦理黄色片| 18禁裸乳无遮挡动漫免费视频| av网站在线播放免费| 国产视频首页在线观看| 国精品久久久久久国模美| 欧美性长视频在线观看| 亚洲成人免费电影在线观看 | 国产色视频综合| bbb黄色大片| 一级片'在线观看视频| 丰满迷人的少妇在线观看| 看免费av毛片| av不卡在线播放| 国产精品一二三区在线看| 99久久人妻综合| 波多野结衣av一区二区av| 午夜福利一区二区在线看| 黄片小视频在线播放| 精品久久久久久久毛片微露脸 | 精品亚洲成国产av| 人妻 亚洲 视频| www.av在线官网国产| 男女下面插进去视频免费观看| 久久九九热精品免费| 亚洲av男天堂| 精品一区在线观看国产| 精品少妇久久久久久888优播| 中文字幕高清在线视频| 午夜激情av网站| 国产欧美日韩一区二区三区在线| 无遮挡黄片免费观看| 午夜免费男女啪啪视频观看| 美女脱内裤让男人舔精品视频| 美女大奶头黄色视频| 日韩一区二区三区影片| 老司机靠b影院| 啦啦啦啦在线视频资源| 女人被躁到高潮嗷嗷叫费观| 亚洲国产看品久久| www.999成人在线观看| 久久精品国产亚洲av高清一级| 一区在线观看完整版| 久久久久国产一级毛片高清牌| 黄片播放在线免费| 99久久99久久久精品蜜桃| 久热爱精品视频在线9| 国产在线一区二区三区精| 国产不卡av网站在线观看| 国产男女超爽视频在线观看| 18在线观看网站| 国产成人免费观看mmmm| 国产人伦9x9x在线观看| 亚洲欧美一区二区三区久久| 亚洲欧美中文字幕日韩二区| 午夜精品国产一区二区电影| 亚洲国产精品国产精品| 考比视频在线观看| 99精国产麻豆久久婷婷| 亚洲精品自拍成人| 亚洲国产av影院在线观看| 久久午夜综合久久蜜桃| 亚洲国产精品成人久久小说| 女人被躁到高潮嗷嗷叫费观| 中文字幕制服av| 欧美黑人精品巨大| 成人免费观看视频高清| 日韩av不卡免费在线播放| 18禁裸乳无遮挡动漫免费视频| 精品少妇久久久久久888优播| 国产精品99久久99久久久不卡| 黑人巨大精品欧美一区二区蜜桃| 尾随美女入室| 永久免费av网站大全| 国产极品粉嫩免费观看在线| 国产人伦9x9x在线观看| 99国产精品一区二区三区| 好男人电影高清在线观看| 国产不卡av网站在线观看| 中文字幕人妻丝袜制服| 久久久精品国产亚洲av高清涩受| 免费久久久久久久精品成人欧美视频| 黄片播放在线免费| 欧美日韩国产mv在线观看视频| 精品国产超薄肉色丝袜足j| 巨乳人妻的诱惑在线观看| 女性生殖器流出的白浆| av网站免费在线观看视频| 视频区图区小说| 在线观看国产h片| 五月天丁香电影| 久久精品国产a三级三级三级| 久久久欧美国产精品| 精品国产一区二区三区久久久樱花| 欧美亚洲 丝袜 人妻 在线| 爱豆传媒免费全集在线观看| 亚洲精品国产一区二区精华液| 欧美黄色淫秽网站| 中文字幕最新亚洲高清| 熟女少妇亚洲综合色aaa.| 热re99久久国产66热| 久久人人爽人人片av| 国产免费又黄又爽又色| 波野结衣二区三区在线| 最新的欧美精品一区二区| 亚洲精品乱久久久久久| 免费在线观看日本一区| 国产视频一区二区在线看| 免费黄频网站在线观看国产| 日韩,欧美,国产一区二区三区| 日韩人妻精品一区2区三区| 亚洲欧美成人综合另类久久久| 午夜久久久在线观看| 女人爽到高潮嗷嗷叫在线视频| 性色av一级| 国产淫语在线视频| 黑丝袜美女国产一区| 亚洲天堂av无毛| 国产视频首页在线观看| 亚洲情色 制服丝袜| 国产黄色视频一区二区在线观看| 精品亚洲乱码少妇综合久久| 在线精品无人区一区二区三| 国产成人精品无人区| 亚洲伊人久久精品综合| 亚洲精品av麻豆狂野| 国产精品国产三级国产专区5o| 亚洲国产看品久久| 天天添夜夜摸| 视频在线观看一区二区三区| 一级片'在线观看视频| 一级毛片我不卡| 日韩制服丝袜自拍偷拍| 久久久久精品人妻al黑| 男的添女的下面高潮视频| 成人国产av品久久久| 日韩中文字幕欧美一区二区 | www.精华液| 亚洲精品国产av成人精品| av有码第一页| 国产精品国产av在线观看| av在线播放精品| 一区二区三区四区激情视频| 欧美精品高潮呻吟av久久| 大话2 男鬼变身卡| 国产视频首页在线观看| 晚上一个人看的免费电影| 国产成人精品在线电影| 久久影院123| 妹子高潮喷水视频| 99国产综合亚洲精品| 国产精品亚洲av一区麻豆| 午夜久久久在线观看| 亚洲av综合色区一区| 黄频高清免费视频| 美女福利国产在线| 深夜精品福利| 免费av中文字幕在线| 伊人久久大香线蕉亚洲五| 男的添女的下面高潮视频| xxxhd国产人妻xxx| 免费看不卡的av| 久久久精品国产亚洲av高清涩受| 久久久久久免费高清国产稀缺| 肉色欧美久久久久久久蜜桃| 日本av手机在线免费观看| 波多野结衣av一区二区av| 国产色视频综合| 不卡av一区二区三区| 十八禁网站网址无遮挡| 久久久久精品国产欧美久久久 | 午夜福利视频在线观看免费| 久9热在线精品视频| 五月天丁香电影| 高清视频免费观看一区二区| 国产精品 欧美亚洲| 99久久人妻综合| 亚洲精品久久久久久婷婷小说| 亚洲欧美一区二区三区久久| 黄色一级大片看看| 青青草视频在线视频观看| 日本欧美视频一区| 国产精品秋霞免费鲁丝片| 侵犯人妻中文字幕一二三四区| 成人亚洲精品一区在线观看| 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| 中文字幕色久视频| 精品人妻1区二区| 视频在线观看一区二区三区| 美女国产高潮福利片在线看| 亚洲九九香蕉| 精品少妇内射三级| 国产一区二区三区综合在线观看| 亚洲av电影在线观看一区二区三区| 国产精品麻豆人妻色哟哟久久| 久9热在线精品视频| 97精品久久久久久久久久精品| 99精国产麻豆久久婷婷| 亚洲国产欧美网| 国产成人一区二区在线| 国产高清videossex| 国产成人系列免费观看| 视频区图区小说| 亚洲三区欧美一区| cao死你这个sao货| 不卡av一区二区三区| 久久女婷五月综合色啪小说| 高清不卡的av网站| 亚洲精品一二三| 亚洲av男天堂| 欧美+亚洲+日韩+国产| 一区二区三区精品91| 国产精品熟女久久久久浪| 岛国毛片在线播放| 男女高潮啪啪啪动态图| av欧美777| 99久久精品国产亚洲精品| 欧美成人午夜精品| 亚洲激情五月婷婷啪啪| 欧美日韩av久久| 亚洲九九香蕉| 日韩免费高清中文字幕av| 视频区欧美日本亚洲| av天堂久久9| 亚洲国产看品久久| 精品欧美一区二区三区在线| 免费黄频网站在线观看国产| 操美女的视频在线观看| 老熟女久久久| 人人妻人人澡人人看| 久久精品久久久久久久性| 1024视频免费在线观看| 色婷婷久久久亚洲欧美| 黄频高清免费视频| 免费日韩欧美在线观看| 亚洲欧美精品自产自拍| 亚洲自偷自拍图片 自拍| 国产高清不卡午夜福利| 高清不卡的av网站| 亚洲中文av在线| 久久久久久久久免费视频了| 欧美黑人欧美精品刺激| 国产一区二区 视频在线| 午夜免费观看性视频| 七月丁香在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av男天堂| 无遮挡黄片免费观看| 午夜老司机福利片| 欧美+亚洲+日韩+国产| 精品熟女少妇八av免费久了| 女人被躁到高潮嗷嗷叫费观| 欧美人与性动交α欧美精品济南到| 久久人人爽人人片av| 人成视频在线观看免费观看| 国产片特级美女逼逼视频| 免费女性裸体啪啪无遮挡网站| 色播在线永久视频| 精品国产一区二区三区久久久樱花| 中文字幕av电影在线播放| 亚洲国产精品999| www.熟女人妻精品国产| 少妇被粗大的猛进出69影院| 国产激情久久老熟女| 日韩视频在线欧美| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| 成年女人毛片免费观看观看9 | 免费观看a级毛片全部| 亚洲欧美激情在线| 国产一区二区三区av在线| 一边摸一边做爽爽视频免费| 日韩大码丰满熟妇| 亚洲成人免费电影在线观看 | 亚洲 国产 在线| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久免费高清国产稀缺| 亚洲五月婷婷丁香| 欧美日韩黄片免| 在线观看免费日韩欧美大片| 妹子高潮喷水视频| 欧美日韩一级在线毛片| 天天躁夜夜躁狠狠躁躁| 一区二区av电影网| 亚洲图色成人| 桃花免费在线播放| 亚洲国产毛片av蜜桃av| 免费久久久久久久精品成人欧美视频| 久久热在线av| 十八禁网站网址无遮挡| 国产成人a∨麻豆精品| 国产精品九九99| 亚洲色图综合在线观看| 国产亚洲欧美在线一区二区| 91九色精品人成在线观看| 老司机影院毛片| 热re99久久国产66热| 亚洲av日韩精品久久久久久密 | 亚洲av美国av| 人人妻人人添人人爽欧美一区卜| 亚洲国产日韩一区二区| 国产一区二区 视频在线| 宅男免费午夜| 亚洲色图综合在线观看| 久久久国产精品麻豆| 国产成人a∨麻豆精品| 午夜福利在线免费观看网站| 女人爽到高潮嗷嗷叫在线视频| 波野结衣二区三区在线| 一本久久精品| 最近最新中文字幕大全免费视频 | 波多野结衣一区麻豆| 久久精品久久久久久久性| 成人黄色视频免费在线看| 一边摸一边做爽爽视频免费| 欧美精品一区二区免费开放| 国产高清视频在线播放一区 | 中文字幕av电影在线播放| 女人久久www免费人成看片| av国产精品久久久久影院| 国产成人一区二区在线| 99国产精品99久久久久| 狠狠婷婷综合久久久久久88av| 高清欧美精品videossex| 天堂中文最新版在线下载| 老司机亚洲免费影院| 大型av网站在线播放| 9热在线视频观看99| 亚洲av国产av综合av卡| 欧美 日韩 精品 国产| 久久久精品国产亚洲av高清涩受| 天天躁狠狠躁夜夜躁狠狠躁| svipshipincom国产片| 黄片小视频在线播放| 欧美日韩视频高清一区二区三区二| 亚洲久久久国产精品| 啦啦啦中文免费视频观看日本| 亚洲av日韩精品久久久久久密 | 日韩,欧美,国产一区二区三区| 日日爽夜夜爽网站| 精品国产一区二区三区久久久樱花| 久久精品亚洲熟妇少妇任你| 久久免费观看电影| 精品国产超薄肉色丝袜足j| 波多野结衣一区麻豆| www.自偷自拍.com| 亚洲成国产人片在线观看| 日韩精品免费视频一区二区三区| 两个人看的免费小视频| 国产在线免费精品| av电影中文网址| 在线看a的网站| 黄色视频在线播放观看不卡| 大香蕉久久成人网| 狠狠精品人妻久久久久久综合| 精品国产一区二区久久| 亚洲国产毛片av蜜桃av| 国产人伦9x9x在线观看| 男的添女的下面高潮视频| 一二三四在线观看免费中文在| 国产成人影院久久av| 另类精品久久| 一区二区av电影网| 精品国产一区二区三区四区第35| 久久国产精品人妻蜜桃| 久久久久精品人妻al黑| 亚洲精品国产av成人精品| 久久久久久久久久久久大奶| 999精品在线视频| 日韩制服骚丝袜av| 欧美97在线视频| 青春草亚洲视频在线观看| 亚洲精品美女久久av网站| 亚洲成人手机| 亚洲av日韩在线播放| 久久人人爽人人片av| 欧美日韩一级在线毛片| 国产日韩欧美在线精品| 国产一区二区 视频在线| 国产亚洲欧美在线一区二区| 老司机影院成人| 亚洲一区二区三区欧美精品| 狠狠精品人妻久久久久久综合| 国产精品亚洲av一区麻豆| 赤兔流量卡办理| 亚洲色图综合在线观看| 欧美xxⅹ黑人| 啦啦啦啦在线视频资源| 日韩制服骚丝袜av| 午夜免费观看性视频| 黄色怎么调成土黄色| 久久女婷五月综合色啪小说| 欧美日韩精品网址| 99精品久久久久人妻精品| 久久久国产欧美日韩av| 国产男女内射视频| 在线精品无人区一区二区三| 国产在视频线精品| 1024视频免费在线观看| 欧美+亚洲+日韩+国产| av网站在线播放免费| 欧美黑人欧美精品刺激| 国产在视频线精品| 欧美亚洲 丝袜 人妻 在线| 久久久久久人人人人人| 爱豆传媒免费全集在线观看| kizo精华| 国产成人a∨麻豆精品| 国产成人精品无人区| 一区二区三区激情视频| 国产麻豆69| 精品一区二区三卡| 精品一区二区三区av网在线观看 | 极品少妇高潮喷水抽搐| 男女边吃奶边做爰视频| 在线天堂中文资源库| 女人爽到高潮嗷嗷叫在线视频| 一级,二级,三级黄色视频| 狠狠精品人妻久久久久久综合| 91麻豆av在线| a级片在线免费高清观看视频| 国产成人免费无遮挡视频| 性色av乱码一区二区三区2| 亚洲五月色婷婷综合| av福利片在线| 国产精品香港三级国产av潘金莲 | 你懂的网址亚洲精品在线观看| 亚洲欧美精品综合一区二区三区| 中文字幕人妻熟女乱码| 久久狼人影院| 亚洲精品久久午夜乱码| 男女免费视频国产| 欧美 日韩 精品 国产| av又黄又爽大尺度在线免费看| 蜜桃在线观看..| 美女福利国产在线| 男女边摸边吃奶| 纯流量卡能插随身wifi吗| 可以免费在线观看a视频的电影网站| 蜜桃在线观看..| 久久精品国产综合久久久| 老鸭窝网址在线观看| 晚上一个人看的免费电影| 色综合欧美亚洲国产小说| 精品人妻熟女毛片av久久网站| av不卡在线播放| 午夜av观看不卡| 亚洲精品一区蜜桃| 久久午夜综合久久蜜桃| 久久人妻福利社区极品人妻图片 | 日韩制服骚丝袜av| 久久天堂一区二区三区四区| 国产成人精品久久二区二区免费| 母亲3免费完整高清在线观看| 每晚都被弄得嗷嗷叫到高潮| xxx大片免费视频| 高清视频免费观看一区二区| 久久精品亚洲熟妇少妇任你| 69精品国产乱码久久久| 美女午夜性视频免费| 国产av国产精品国产| 在线看a的网站| 久久久久久人人人人人| 国产成人啪精品午夜网站| 黄色视频不卡| 曰老女人黄片| 亚洲人成电影免费在线| avwww免费| 丁香六月欧美| 久久久久久久久久久久大奶| 乱人伦中国视频| 激情五月婷婷亚洲| 国产高清不卡午夜福利| 国产亚洲一区二区精品| 亚洲成国产人片在线观看| 欧美激情 高清一区二区三区| 一本大道久久a久久精品| 中文欧美无线码| 好男人电影高清在线观看| 国产成人91sexporn| 国产精品麻豆人妻色哟哟久久| 午夜免费鲁丝| 男女之事视频高清在线观看 | 人妻一区二区av| 美国免费a级毛片| 国产精品一国产av| 国产xxxxx性猛交| 久久青草综合色| 天天操日日干夜夜撸| 久久久久久久国产电影| 国产精品三级大全| 99国产综合亚洲精品| 亚洲精品国产区一区二| 国产精品免费大片| videosex国产| 丰满饥渴人妻一区二区三| 中文精品一卡2卡3卡4更新| 国产一区二区三区综合在线观看| 97人妻天天添夜夜摸| 日本wwww免费看| 国产不卡av网站在线观看| 高清欧美精品videossex| 中文字幕人妻熟女乱码| 国产又色又爽无遮挡免| 国产在视频线精品| 中文字幕人妻熟女乱码| 国产欧美日韩一区二区三 | 侵犯人妻中文字幕一二三四区| 成人国产av品久久久| 亚洲熟女毛片儿| 久久久久国产一级毛片高清牌| 久久久久精品人妻al黑| 一二三四在线观看免费中文在| 成人国产av品久久久| 搡老岳熟女国产| 真人做人爱边吃奶动态| 一区二区三区精品91| 美女中出高潮动态图| 久久精品国产综合久久久| 后天国语完整版免费观看| 日韩电影二区| 一级毛片我不卡| 国产精品国产三级专区第一集| 巨乳人妻的诱惑在线观看| 国产1区2区3区精品| 亚洲精品日韩在线中文字幕| 国产亚洲av片在线观看秒播厂| 91老司机精品| 欧美在线一区亚洲| 男女午夜视频在线观看| xxxhd国产人妻xxx| 这个男人来自地球电影免费观看|