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

    兩種輪作模式下秸稈還田對土壤呼吸及其溫度敏感性的影響*

    2017-08-09 03:31:36丁瑞霞王維鈺
    關(guān)鍵詞:夏大豆土壤溫度夏玉米

    丁瑞霞, 王維鈺, 張 青

    ?

    兩種輪作模式下秸稈還田對土壤呼吸及其溫度敏感性的影響*

    丁瑞霞1, 王維鈺2,3, 張 青2

    (1. 西北農(nóng)林科技大學(xué)中國旱區(qū)節(jié)水農(nóng)業(yè)研究院/農(nóng)業(yè)部西北黃土高原作物生理生態(tài)與耕作重點實驗室 楊凌 712100; 2. 西北農(nóng)林科技大學(xué)農(nóng)學(xué)院 楊凌 712100; 3. 陜西省循環(huán)農(nóng)業(yè)工程技術(shù)研究中心 楊凌 712100)

    通過分析不同作物輪作模式下秸稈還田對土壤呼吸及其溫度敏感性的影響, 為深入探究關(guān)中地區(qū)農(nóng)田生態(tài)系統(tǒng)碳循環(huán)提供理論依據(jù)。試驗設(shè)置于陜西省楊凌地區(qū), 在2012年10月至2014年9月期間以冬小麥-夏玉米輪作模式和冬小麥-夏大豆輪作模式作為研究對象, 分別設(shè)置秸稈還田(SM)和秸稈不還田(NS)兩個處理, 測定分析不同處理下土壤呼吸、土壤溫度及土壤含水量的變化趨勢和差異, 并估算土壤呼吸的溫度敏感性(10)。結(jié)果表明: 土壤呼吸存在明顯的季節(jié)變化, 在作物生育期大部分時間內(nèi), SM處理的土壤呼吸速率均顯著高于NS處理(<0.05), 且SM處理的作物生育期土壤呼吸平均速率及土壤呼吸累計排放量也極顯著高于NS處理(<0.01); 不同作物生育期土壤呼吸平均速率依次為夏玉米>夏大豆>冬小麥, 土壤呼吸總量表現(xiàn)為冬小麥>夏玉米>夏大豆、冬小麥-夏玉米輪作>冬小麥-夏大豆輪作。冬小麥-夏玉米輪作與冬小麥-大豆輪作的土壤溫度間存在差異; 其中, 在冬小麥生育前期, 冬小麥-夏玉米輪作的土壤溫度顯著高于冬小麥-大豆輪作; 第2季夏玉米生育期內(nèi)5 cm深度的土壤溫度顯著低于同季的夏大豆; 相比NS處理, SM處理能提高冬季土壤的溫度, 并降低春季和夏季的土壤溫度; 在高溫少雨的時期內(nèi), SM處理能夠提高0~30 cm土壤的平均含水量, 不同的前茬作物引起兩種輪作模式中冬小麥耕作層土壤含水量間明顯的差異,夏玉米耕作層土壤含水量顯著高于夏大豆。相關(guān)分析表明, 土壤呼吸與5 cm和10 cm土壤溫度均存在極顯著的正相關(guān)性, 且與5 cm土壤溫度的相關(guān)性更好; 但土壤呼吸與0~30 cm的土壤平均含水量無顯著相關(guān)性。5 cm和10 cm土壤溫度變化能夠分別解釋土壤呼吸變化的64.6%~67.3%和51.5%~59.6%。整個研究周期內(nèi), 溫度敏感性(10)為1.70~2.01, 冬小麥-夏玉米輪作的溫度敏感性顯著高于冬小麥-大豆輪作, 且同一輪作模式下SM處理的溫度敏感性顯著低于NS處理。因此, 秸稈還田能夠提高農(nóng)田的土壤呼吸作用, 降低土壤呼吸的溫度敏感性, 同時能夠調(diào)節(jié)土壤的水熱狀況。

    輪作模式; 秸稈還田; 土壤呼吸; 溫度敏感性; 土壤水熱

    自工業(yè)革命以來, 人類社會快速發(fā)展導(dǎo)致自然環(huán)境受到前所未有的挑戰(zhàn), 其中氣候變化及其影響是當(dāng)前人類面臨的最大環(huán)境問題之一[1]。隨著環(huán)境的日益惡化, 與氣候變化密切相關(guān)的碳循環(huán)問題成為當(dāng)今國內(nèi)外研究熱點之一。土壤碳庫是陸地生態(tài)系統(tǒng)中最大的碳庫, 土壤呼吸作為土壤碳庫與大氣碳庫之間CO2交換的主要形式, 其排放量占整個陸地生態(tài)系統(tǒng)總排放量的2/3以上[2]。據(jù)統(tǒng)計, 每年全球超過75 Pg(1 Pg=1×1015g)的碳通過土壤呼吸的方式向大氣中釋放[3]。土壤呼吸微小的變化也將勢必導(dǎo)致大氣中碳濃度的變化[4], 進(jìn)而加劇全球氣候變暖, 危急人類今后的生存與環(huán)境。有研究認(rèn)為農(nóng)業(yè)生產(chǎn)過程排放的CO2占人為溫室氣體排放量的21%~25%[5]。我國作為一個農(nóng)業(yè)大國, 由農(nóng)田所引起的碳排放量在1999—2009這10年間的增幅達(dá)93.9%[6]。因此, 研究農(nóng)業(yè)生產(chǎn)過程的碳排放具有重要的現(xiàn)實意義。

    相比其他陸地生態(tài)系統(tǒng), 農(nóng)田生態(tài)系統(tǒng)受到人類活動的干擾最大, 耕作方式、施肥方式、灌溉方式的不同均會引起農(nóng)田生態(tài)系統(tǒng)土壤環(huán)境的差異, 進(jìn)而影響土壤呼吸速率[7-8]。一般認(rèn)為, 土壤表面所測定的土壤呼吸通量是根呼吸和微生物呼吸的總和[9]。土壤呼吸強度與土壤有機質(zhì)含量及礦化速率、土壤微生物種類及活性之間均存在密切的聯(lián)系[10]。秸稈還田作為目前農(nóng)業(yè)生產(chǎn)中廣泛應(yīng)用的耕作措施, 被認(rèn)為能夠改善土壤的水熱狀況, 提高土壤的透氣性。同時前人的研究也指出秸稈還田能夠顯著影響土壤的有機質(zhì)含量[11]及土壤的微生物量[12]。前人針對土壤呼吸對秸稈還田的響應(yīng)做出了大量的研究, 研究結(jié)果表明, 秸稈還田對于土壤呼吸具有顯著影響。其中于愛忠等[13]的研究表明免耕秸稈還田提高了土壤的有機質(zhì)含量, 促進(jìn)了有機質(zhì)分解釋放CO2, 因而提高了農(nóng)田土壤呼吸速率。也有研究者認(rèn)為免耕留茬能夠顯著降低土壤呼吸速率, 且土壤呼吸速率降低的幅度隨留茬量的增多而增大[14]。

    前人的研究表明, 土壤呼吸與溫度有密切的聯(lián)系[15], 這種關(guān)聯(lián)通常表現(xiàn)為: 溫度升高會導(dǎo)致土壤CO2通量增大[16]。而土壤呼吸對溫度的變化也具有適應(yīng)性, 這種適應(yīng)性會弱化及降低土壤呼吸和溫度變化之間的關(guān)系[17]。這種適應(yīng)性的強弱主要取決于土壤呼吸對溫度變化的敏感性。因此, 土壤呼吸的溫度敏感性與全球碳循環(huán)對氣候變化的響應(yīng)關(guān)系有非常緊密的聯(lián)系[18]。目前, 研究者已經(jīng)對此開展了大量的有關(guān)研究, 構(gòu)建了溫度響應(yīng)函數(shù), 而在這其中使用最廣泛的函數(shù)為10(土壤呼吸敏感性),10被定義為溫度增加10 ℃所引起的土壤呼吸的變化量[19]。在影響土壤呼吸溫度敏感性的眾多因子中, 溫度和水分對于土壤呼吸及其溫度敏感性的影響占主導(dǎo)地位。溫度和水分的差異不僅造成土壤生物種類和數(shù)量的差異, 同時也影響土壤微生物所需底物的供應(yīng)[20]。前人總結(jié)不同生態(tài)系統(tǒng)下的土壤呼吸溫度敏感性可知,10值介于1.28~4.75[21]。近年來, 有關(guān)農(nóng)田生態(tài)系統(tǒng)下土壤呼吸的溫度敏感性(10)的研究逐漸增多, 前人的研究發(fā)現(xiàn)耕作栽培措施的差異能夠引起土壤的水熱條件、養(yǎng)分含量及微生物狀況的變化, 最終導(dǎo)致土壤呼吸的溫度敏感性發(fā)生變化。但研究的主要側(cè)重點多集中在不同施肥條件或耕作模式下土壤呼吸溫度敏感性的變化[22-25]。而有關(guān)秸稈還田條件對土壤呼吸敏感性的研究仍不多見, 已有的研究也多將秸稈還田作為耕作措施的一種, 分析不同耕作措施對土壤呼吸及其溫度敏感性的影響[26-27]; 此外, 針對不同秸稈還田量對土壤呼吸及其溫度敏感性影響的研究也多為一種作物[28], 對于不同作物土壤呼吸溫度敏感性的研究鮮有報道。本研究以陜西省常見的兩種輪作模式為研究對象, 分析和比較不同輪作模式下秸稈還田對土壤呼吸溫度敏感性的影響, 探究秸稈還田對土壤呼吸與溫度變化之間反饋關(guān)系的影響, 為進(jìn)一步研究和分析關(guān)中地區(qū)農(nóng)田生態(tài)系統(tǒng)碳循環(huán)提供理論依據(jù)。

    1 材料與方法

    1.1 研究地基本信息

    試驗區(qū)域設(shè)于陜西省楊凌農(nóng)業(yè)示范區(qū)內(nèi)(108°07′E, 34°12′N), 該區(qū)域氣候類型為大陸性暖溫帶季風(fēng)區(qū)氣候, 年平均溫度12.9 ℃, 全年降水主要分布在每年的7—9月, 年平均降水量約為660 mm, 作物熟制為一年兩熟制。試驗土壤為土, 試驗地土壤基本理化性狀見表1。

    表1 試驗地土壤養(yǎng)分特性

    試驗地土壤養(yǎng)分特性于2011年試驗前測定。The soil chemical properties of the tested field were measured before experiment in 2011.

    1.2 試驗設(shè)計

    本試驗的周期為2012—2014年, 研究對象為冬小麥(L.)-夏玉米(L.)和冬小麥-夏大豆(L.)兩種輪作模式。其中, 冬小麥分別于2012年及2013年10月初播種, 并于2013年及2014年的6月收獲, 夏玉米及夏大豆在上季冬小麥?zhǔn)斋@后立即播種, 于2013年及2014年的10月初收獲。冬小麥品種為‘西農(nóng)889’, 夏玉米品種為‘漯單9號’, 夏大豆品種為‘東豆339’。由于灌溉刺激干旱土壤的土壤呼吸[29], 進(jìn)而影響土壤呼吸測定的準(zhǔn)確性, 在本研究周期內(nèi)盡可能避免灌溉, 每季冬小麥均只進(jìn)行一次冬灌, 時間約為每年的1月中旬(冬小麥越冬期), 其余時間均依靠降水補充作物生長所需; 每季夏玉米及夏大豆生育期內(nèi)均只灌溉一次, 以減緩作物的干旱脅迫, 時間為每年的7月中旬, 其余時間均依靠降水補充作物生長所需。所有灌溉均為漫灌。

    整個試驗地均實行免耕播種, 每個輪作模式下均設(shè)置2個處理: 秸稈還田(SM)和秸稈不還田(NS, 對照)。其中, 冬小麥-夏玉米輪作模式下的處理為SM1和NS1, 冬小麥-夏大豆輪作模式下的處理為SM2和NS2。所有試驗處理均設(shè)置為8.6 m×8 m的小區(qū), 小區(qū)相互間隔0.5 m, 相同處理均設(shè)置3個重復(fù)。秸稈還田處理在前茬作物收獲后, 直接將作物殘茬還田; 秸稈不還田處理在前茬作物收獲后, 清理所有根茬。冬小麥?zhǔn)┯媚蛩丶傲姿岫@(尿素375 kg×hm-2+磷酸二銨375 kg×hm-2), 夏玉米僅施用尿素(375 kg×hm-2), 夏大豆僅施用磷酸二銨(150 kg×hm-2), 所有肥料均作為基肥施用。不同作物秸稈還田情況及施氮量見表2。各處理其余的田間管理均保持一致。

    表2 秸稈還田處理的不同作物秸稈還田量及秸稈還田和不還田處理的施氮量

    施氮量根據(jù)所購買的尿素及磷酸二銨的總氮含量估算, 其中尿素的總氮含量≥46%, 磷酸二銨總氮含量≥17.4%。The amount of nitrogen was calculated based on nitrogen content in urea and diammonium phosphate, respectively. Urea and diammonium phosphate contain 46% and 17.4% nitrogen, respectively.

    1.3 測定項目

    1.3.1 氣體測定

    土壤呼吸速率采用GXH-3010E1型便攜式紅外線分析器(華云分析儀器研究所有限公司)測定。在作物播種后, 在每個小區(qū)中部分別放置PVC腔室(腔室高10 cm, 底部直徑為16 cm), 埋置深度為5 cm, 每個小區(qū)重復(fù)3次, 腔室呈“三角形”排列, 相互間隔1 m。在測量一定時間范圍內(nèi)的土壤CO2通量的過程中, 在腔室的頂部放置1個小風(fēng)扇, 持續(xù)開啟3 min, 使腔室內(nèi)的氣體混合均勻, 然后測量。每次測量均在每日的9:00—11:00完成。

    土壤呼吸速率計算公式:

    (2-1)Δ(1)

    式中:為土壤呼吸速率(μmol·m-2·s-1);為換算系數(shù),=1.80(25 ℃, 1個標(biāo)準(zhǔn)大氣壓);1、2分別為測定時CO2初始質(zhì)量分?jǐn)?shù)和測定結(jié)束后即時質(zhì)量分?jǐn)?shù)(%);為容器高(m); Δ為測定時間變化(s)。

    土壤呼吸速率自第1季冬小麥播種后(2012年10月)開始測量, 由于冬小麥在返青前的生長較慢, 故分別于冬小麥的苗期(播種后15 d)、分蘗期(播種后50 d)、越冬期(播種后95 d)、返青期(播種后155 d)各測定一次, 返青期后每15 d測定一次, 全生育期內(nèi)共測定9次; 下茬作物播種后, 每15 d測定一次, 全生育期內(nèi)共測定7次。在降雨天氣出現(xiàn)后, 可適當(dāng)后延3~5 d測定土壤呼吸速率, 以降低降水引起的干擾。

    1.3.2 土壤溫度測定

    作物播種后, 將一組不同深度的直角地溫計埋至距PVC腔室10 cm處, 每個處理3組。測定土壤呼吸速率的當(dāng)天內(nèi), 讀取土壤CO2質(zhì)量分?jǐn)?shù)的同時讀取5 cm、10 cm、15 cm 層土壤溫度值。

    1.3.3 土壤含水量測定

    采用烘干稱量法測定土壤含水量, 每次土壤呼吸速率測定的當(dāng)天, 用土鉆采集0~30 cm土壤樣品(每10 cm為一層), 3次重復(fù), 帶回實驗室測定。

    1.3.4 土壤呼吸累計釋放量測定

    土壤呼吸累計釋放量參考Zhai等[29]的計算公式:

    R=(R+R+1)/2×3 600×24×44×10-8(2)

    式中:R為土壤呼吸累計釋放量[t(CO2)·hm-2],R為第次測量的土壤呼吸速率,R+1為間隔天測量的土壤呼吸速率,為相鄰兩次測量間隔天數(shù)。

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

    有關(guān)土壤呼吸速率對溫度響應(yīng)的研究多使用指數(shù)函數(shù)[30], 函數(shù)表達(dá)式為:

    =′e(3)

    式中:為測定的土壤呼吸速率(μmol·m-2·s-1),為測定的土壤溫度(℃),和為回歸系數(shù)。

    土壤呼吸對溫度的敏感性通常用10來描述,10指溫度升高10 ℃所引起的土壤呼吸的變化值:

    10=e10b(4)

    式中:通過(2)式獲得。

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

    2 結(jié)果與分析

    2.1 不同輪作模式下各處理土壤溫度的動態(tài)變化

    如圖1所示, 2012年10月—2013年9月期間, 月平均氣溫最低的月份為2013年1月, 為-1.6 ℃; 月平均氣溫最高的月份為2013年8月, 為23.1 ℃。2013年10月—2014年9月期間, 月平均氣溫最低的月份為2013年12月, 為-1.3 ℃; 月平均氣溫最高的月份為2014年6月, 為28.2 ℃。2013—2014連續(xù)兩年內(nèi), 6—8月的月平均溫度均高于21 ℃。在整個研究周期內(nèi), 大部分時間的月平均降雨量均低于230 mm, 僅2013年4月、2014年2月、2014年4月、2014年9月的降雨量超過300 mm。

    在冬小麥-夏玉米(麥玉)和冬小麥-夏大豆(麥豆)兩種輪作模式下, 冬小麥生育期的土壤溫度以播種后95 d為節(jié)點, 呈現(xiàn)先降低后升高的變化趨勢(圖2)。在連續(xù)兩季冬小麥播種后95 d時, 麥玉輪作SM處理5 cm深度的土壤溫度比NS處理高0.1 ℃和0.8 ℃, 麥豆輪作SM處理5 cm深度的土壤溫度比NS處理高0.7 ℃和0.2 ℃; 同理, 麥玉輪作SM處理10 cm深度的土壤溫度比NS處理高0.6 ℃和0.2 ℃, 麥豆輪作SM處理10 cm深度的土壤溫度比NS處理高0.8 ℃和0.4 ℃。在冬小麥播種95~230 d, SM處理在5 cm及10 cm深度的土壤溫度均低于NS處理, 其中麥玉輪作2012—2013年冬小麥生育期分別下降4.4%~10.7%和1.3%~5.3%, 2013—2014年分別下降2.5%~10%和0.5%~4.2%; 同理, 麥豆輪作2012—2013年冬小麥生育期SM處理較NS處理分別下降1.9%~5.4%和1.1%~4.6%, 而在2013—2014年冬小麥生育期內(nèi)分別下降2.1%~8.1%和1.6%~4.4%。

    第一季后茬作物生育期內(nèi)(2013年6—9月)2013年夏大豆生育期內(nèi)土壤溫度呈降低趨勢; 夏玉米播種后15~60 d土壤溫度變化幅度不大, 60~105 d土壤溫度逐漸下降; 在第二季后茬作物生育期內(nèi)(2014年6—9月)夏玉米和夏大豆生育期土壤溫度以播種后45 d為節(jié)點均呈現(xiàn)出先升高后降低的變化趨勢。第一季夏玉米生育期內(nèi)(2013年), 5 cm深度和10 cm深度, SM處理土壤溫度分別比NS低0.3%~7.7%和0.3%~7.4%, 在第二季生育期內(nèi)(2014年)分別低1.7%~9.8%和0.35%~0.4%。同理, 在2013年夏大豆生育期內(nèi), SM處理土壤溫度分別比NS處理低0.3%~4.4%和0.7%~6.1%, 而在2014年夏大豆生育期內(nèi)則分別低0.7%~4.6%和0.7%~3.7%。

    在整個研究周期內(nèi), 麥玉輪作在冬小麥生育前期的土壤溫度顯著高于麥豆輪作; 而在第2茬作物生育期內(nèi), 夏玉米5 cm深度的土壤溫度也顯著低于夏大豆。同一輪作模式下, SM處理與NS處理間的差異存在顯著差異; 此外, 除冬小麥播種后50 d至95 d外, 其余時段內(nèi)各處理土壤溫度均隨土壤深度的增加而顯著降低。

    2.2 不同輪作模式下各處理土壤含水量變化

    不同輪作模式下各處理的0~30 cm平均含水量的變化如圖3所示。不同輪作模式間差異顯著(<0.05)在2012—2013年冬小麥播種后185 d、200 d和230 d, 麥玉輪作SM處理的土壤平均含水量比NS處理分別高為9.58%、1.77%和7.46%; 麥豆輪作高0.08%、6.21%和3.53%。而在2013—2014年冬小麥生育期內(nèi), 播種后215 d、230 d麥玉和麥豆輪作模式下SM處理土壤平均含水量均高于NS處理, 增幅分別為6.62%、2.92%和3.08%、5.39%。

    兩季夏玉米生育期內(nèi), 播種后30~90 d, SM處理的土壤平均含水量均高于NS處理, 增幅依次為0.5%~3.18%和0.06%~4.17%。2012—2013年夏大豆整個生育期內(nèi), SM處理的土壤平均含水量均高于NS處理, 增幅為0.78%~3.50%, 而2013—2014年夏大豆生育期內(nèi), 除播種后90 d外, 其余時間內(nèi)SM處理的土壤平均含水量均高于NS處理, 增幅為0.44%~13.85%。

    2.3 不同輪作模式下各處理土壤水熱因子對土壤呼吸的影響

    通過Pearson相關(guān)分析法分析不同處理土壤呼吸與土壤水熱因子間的相關(guān)性。分析結(jié)果表明, 兩種輪作模式下各處理土壤呼吸與土壤溫度極顯著正相關(guān)(表3)。麥玉輪作下, 土壤呼吸與土壤溫度間的相關(guān)系數(shù)為0.611~0.687, 麥豆輪作為0.227~0.397。此外, 兩種輪作模式下, SM處理土壤呼吸與不同深度土壤溫度間的相關(guān)系數(shù)大于NS處理, 且所有處理土壤呼吸與地下5 cm的土壤溫度間的相關(guān)系數(shù)均大于地下10 cm。但不同處理土壤呼吸與0~30 cm深度下的土壤平均含水量間的相關(guān)性不顯著。

    通過擬合不同處理土壤溫度對土壤呼吸的影響可知, 土壤呼吸與土壤溫度間的關(guān)系可用指數(shù)函數(shù)解釋(表4)。從土壤深度來看, 地下5 cm土壤的溫度能夠解釋土壤呼吸變化的64.6%~67.3%, 而地下10 cm土壤溫度能夠解釋土壤呼吸變化的51.5%~59.6%。

    實線為冬小麥-夏玉米輪作, 虛線為冬小麥-夏大豆輪作; NS為秸稈不還田處理, SM為秸稈還田處理。The solid line indicates winter wheat-summer maize rotation system and the dashed line indicates winter wheat-summer soybean rotation system. NS: no straw incorporation treatment; SM: straw incorporation treatment

    2.4 不同輪作模式下各處理土壤呼吸速率的動態(tài)變化

    由圖4所示, 冬小麥播種后, 不同處理土壤呼吸速率均逐漸下降, 至播種后95 d, 土壤呼吸速率降至最低。但不同處理的土壤呼吸速速率在播種后155~230 d的變化規(guī)律有所差異。2012—2013年冬小麥播種后95~185 d, 僅麥豆輪作的SM處理的土壤呼吸速率持續(xù)增大, 其余處理均表現(xiàn)為先升高后降低的變化趨勢。2013—2014年冬小麥播種后, 所有處理的土壤呼吸速率95~200 d持續(xù)升高, 但215~230 d均變?yōu)榻档?。在后茬作物生育期?nèi), 不同處理的土壤呼吸速率均基本呈現(xiàn)先升高后降低的變化趨勢。2012—2013年后茬作物生育期內(nèi), 夏玉米和夏大豆土壤呼吸速率最大值出現(xiàn)的時期有所差異, 夏大豆出現(xiàn)在播種后45 d, 而夏玉米出現(xiàn)在播種后60 d。2013—2014年后茬作物生育期內(nèi), 麥玉輪作的播種后30~105 d土壤呼吸率始終高于麥豆輪作, 且所有處理土壤呼吸速率的最大值均出現(xiàn)在播種后60 d。

    麥玉輪作, 除冬小麥2012—2013年播種后230 d以及2013—2014年播種后200 d外, SM處理的土壤呼吸速率在其余時期均顯著高于NS處理, 增幅為1.08%~36.53%; 同理, 麥豆輪作下, SM處理在2012—2014年播種后95~155 d之外的其余時期內(nèi)的土壤呼吸速率均顯著高于NS處理, 增幅為1.51%~99.43%。在連續(xù)兩季夏玉米生育期內(nèi), 除2013—2014年播后30 d及75 d外, 其余時期內(nèi)SM處理的土壤呼吸速率均顯著高于NS處理, 增幅為0.70%~56.24%; 而在連續(xù)兩季夏大豆生育期內(nèi), SM處理的土壤呼吸速率均顯著高于NS處理, 增幅為0.73%~11.64%。

    實線為冬小麥-夏玉米輪作, 虛線為冬小麥-夏大豆輪作; NS為秸稈不還田處理, SM為秸稈還田處理。The solid line indicates winter wheat-summer maize rotation system and the dashed line indicates winter wheat-summer soybean rotation system. NS: no straw incorporation treatment; SM: straw incorporation treatment.

    表3 兩種輪作模式下不同處理土壤呼吸與土壤水熱因子的相關(guān)性

    NS為秸稈不還田處理, SM為秸稈還田處理。**和*分別表示極顯著(<0.01)和顯著(<0.05)相關(guān)。NS: no straw incorporation treatment; SM: straw incorporation treatment. ** and * indicate significant correlation at 0.01 and 0.05 levels, respectively.

    表4 兩種輪作模式下不同處理土壤呼吸與土壤溫度擬合方程

    NS為秸稈不還田處理, SM為秸稈還田處理。NS: no straw incorporation treatment; SM: straw incorporation treatment.

    冬小麥生育期土壤呼吸平均速率為2.673~3.141 μmol·m-2·s-1, 夏玉米生育期土壤呼吸平均速率為3.401~4.810 μmol·m-2·s-1, 夏大豆生育期土壤呼吸平均速率為3.390~3.762 μmol·m-2·s-1(表5)。整個研究周期內(nèi)SM處理的作物生育期土壤呼吸平均速率均顯著高于NS處理(<0.01); 兩季麥玉輪作模式中, 冬小麥SM處理的作物生育期土壤呼吸平均速率比NS處理分別高14.13%和9.35%, 而麥豆輪作模式中冬小麥SM處理比NS處理分別高12.61%和7.29%, 兩季夏玉米中SM處理比NS處理分別高16.55%和16.70%, 兩季夏大豆中SM處理分別比NS處理高4.01%和2.33%。

    實線為冬小麥-夏玉米輪作, 虛線為冬小麥-夏大豆輪作; NS為秸稈不還田處理, SM為秸稈還田處理。The solid line indicates winter wheat-summer maize rotation system and the dashed line indicates winter wheat-summer soybean rotation system. NS: no straw incorporation treatment; SM: straw incorporation treatment.

    2.5 不同輪作模式下各處理土壤呼吸累計釋放量的動態(tài)變化

    土壤呼吸總量能夠描述整個作物生育期內(nèi)土壤釋放的CO2總量。如表6所示, 不同作物生育期內(nèi)土壤呼吸總量有很大的差異, 表現(xiàn)為冬小麥>夏玉米>夏大豆, 受此影響, 麥玉輪作的土壤呼吸總量[34.68~40.81 t(CO2)·hm-2]高于麥豆輪作[30.04~33.86 t(CO2)·hm-2]。不同作物生育期內(nèi)土壤呼吸總量均為SM>NS, 且處理間差異極顯著(<0.01), 這與作物生育期土壤呼吸平均速率的變化規(guī)律相同。

    2.6 不同輪作模式下各處理土壤呼吸溫度敏感性(10)的變化

    利用指數(shù)函數(shù)=e對連續(xù)兩個輪作周期內(nèi)不同處理土壤呼吸速率對地下5 cm土壤溫度的響應(yīng)關(guān)系進(jìn)行模擬, 通過模擬, 進(jìn)而估算不同處理土壤呼吸的溫度敏感性(10), 不同處理的10值如圖5所示。整個研究周期內(nèi), 不同處理的10值介于1.70~2.01, 相同處理的10表現(xiàn)為麥玉輪作>麥豆輪作(<0.05); 兩種輪作模式中, SM處理的10值均極顯著低于NS處理。在麥玉輪作中, SM的10值比NS處理的低0.20, 在麥豆輪作中, SM的10值比NS處理的低0.16。

    表5 2012年10月(2012-10)—2014年6月(2014-06)兩種輪作模式下作物生育期不同處理的土壤呼吸平均速率

    NS為秸稈不還田處理, SM為秸稈還田處理。不同字母表示同一生育期內(nèi)不同處理間差異顯著(<0.01)。NS: no straw incorporation treatment; SM: straw incorporation treatment. Different letters indicate significant difference among treatments in the same growth period at< 0.01 level.

    表6 2012年10月(2012-10)—2014年6月(2014-06)兩種輪作模式下作物生育期不同處理的土壤呼吸總量

    NS為秸稈不還田處理, SM為秸稈還田處理。不同字母表示同一生育期內(nèi)不同處理間差異顯著(<0.01)。NS: no straw incorporation treatment; SM: straw incorporation treatment. Different letters indicate significant difference among treatments in the same growth period at< 0.01 level.

    NS為秸稈不還田處理, SM為秸稈還田處理。不同大寫字母表示同一輪作模式不同處理間差異顯著(<0.01)。NS: no straw incorporation treatment; SM: straw incorporation treatment. Different letters indicate significant difference between treatments for the same rotation system at< 0.01 level.

    3 討論與結(jié)論

    土壤溫度的動態(tài)變化對氣溫的變化有明顯的響應(yīng), 同時也與地表接收的太陽輻射量有關(guān)。在本研究中, 不同處理土壤溫度的季節(jié)變化趨勢與氣溫變化趨勢相同。但土壤溫度隨土層深度變化的規(guī)律并不一致, 具體表現(xiàn)為: 冬季地下10 cm土壤溫度高于5 cm, 而春季和夏季則正好相反??赡艿脑蚴钱?dāng)外界氣溫較低時, 土壤熱量由深層向淺層傳導(dǎo), 而當(dāng)外界氣溫升高時, 土壤中的熱量則由淺層向深層傳導(dǎo)[31]。前人的研究結(jié)果證實, 秸稈還田能夠調(diào)節(jié)耕層土壤的溫度[32], 本研究所得結(jié)果也驗證了這一結(jié)論, 具體表現(xiàn)為: 當(dāng)外界溫度較低時, 秸稈還田處理的土壤溫度高于不還田處理, 當(dāng)外界溫度逐漸升高時, 秸稈還田處理的土壤溫度均低于不還田處理。在兩種輪作模式下的冬小麥生育期內(nèi), 由于前茬作物(夏玉米和夏大豆)的秸稈還田量不同, 導(dǎo)致土壤向大氣散失的熱量出現(xiàn)差異, 最終引起兩種輪作模式下冬小麥生育前期土壤溫度出現(xiàn)顯著差異, 對于后茬作物而言, 由于夏玉米和夏大豆的作物群體不同, 進(jìn)入作物的生育后期, 植株能夠在土壤的上方形成一個覆蓋層[33], 在外界氣溫較高時, 夏玉米較大的葉片, 減少了地表接收的太陽輻射, 進(jìn)而降低了土壤的溫度。

    前人的研究結(jié)果表明, 秸稈還田能夠提高作物生育期內(nèi)土壤的含水量[34-35]。但本研究結(jié)果顯示, 在作物生育期內(nèi), 秸稈還田的土壤含水量并非在所有時間均高于秸稈不還田處理, 通過分析可知, 在外界氣溫較高且降雨量較低的月份里, 秸稈還田處理耕作層土壤平均含水量高于秸稈不還田處理。這可能是由于在這些月份中, 較高的氣溫加劇了土壤水分的蒸散作用, 單次降雨量較少以及降雨的間隔時間較長, 導(dǎo)致土壤的整體濕潤程度較低, 秸稈還田能夠減緩株間無效蒸散, 進(jìn)而提高土壤的含水量。此外, 研究還發(fā)現(xiàn), 夏玉米耕作層土壤含水量高于夏大豆, 可能的原因是由于夏玉米的葉片較大, 一定程度上減少了太陽輻射與土壤的直接接觸, 從而減緩了土壤的蒸散作用; 對于兩種輪作模式內(nèi)冬小麥耕作層土壤含水量的差異, 可能的原因是由于輪作體系下, 作物根系的數(shù)量及分布均有所差異, 一定程度上改變了土壤的入滲速率[36]。

    有研究者認(rèn)為, 免耕秸稈還田能夠降低土壤呼吸[37-38]。而本試驗基于免耕條件下實行秸稈還田, 研究結(jié)果表明, 秸稈還田處理顯著提高了土壤呼吸作用, 這與Blanco-Canqui等[39-40]研究結(jié)果一致。原因可能是秸稈的腐解是一個緩慢的過程, 秸稈中的碳不能完全轉(zhuǎn)化為土壤碳, 未轉(zhuǎn)化的部分通常以CO2或CH4的形式向大氣中釋放。同時, 由于研究中在秸稈還田的基礎(chǔ)上額外施用氮肥的量較大, 引起土壤C/N發(fā)生變化, 更有利于微生物分解有機質(zhì), 進(jìn)而增加了土壤有機質(zhì)中碳的分解與釋放[41]。此外, 秸稈還田處理的作物生育期土壤呼吸平均速率也高于秸稈不還田處理, 這一研究結(jié)果也與前人相近[42]。研究發(fā)現(xiàn), 在不同輪作模式下土壤呼吸速率有明顯的差異, 這種差異主要來源于不同生長特性的作物交替種植改變了土壤微生物的群體分布和多樣性, 以及土壤養(yǎng)分的循環(huán)[43]。同時, 在輪作模式下, 不同作物生長過程中, 土壤溫度及土壤濕度也存在差異[44]。此外, 前人的研究還指出施氮量、秸稈還田量的不同能夠引起土壤呼吸的差異[45-46], 這些都是引起不同輪作模式間土壤呼吸速率差異的潛在因素。土壤溫度和土壤含水量被認(rèn)為是影響土壤呼吸的重要環(huán)境因子, 在本研究中, 土壤呼吸與地下5 cm及地下10 cm土壤溫度均呈顯著正相關(guān), 且與地下5 cm土壤溫度的相關(guān)性更好, 這與李昌珍等[47]的研究結(jié)果一致。通過回歸分析可知, 地下5 cm土壤的溫度能夠解釋土壤呼吸變化的64.6%~67.3%, 略低于張俊麗等[48]研究結(jié)果中的62.31%~78.66%。此外, 本研究顯示土壤呼吸與0~30 cm土壤平均含水量的相關(guān)性不顯著, 這一結(jié)果也與前人的研究結(jié)果一致[49-50]。

    在整個研究周期內(nèi), 不同處理的10值介于1.70~2.01。數(shù)值略低于前人研究結(jié)果中我國10值范圍(2.25±0.28)[51]。本研究結(jié)果顯示, 麥玉輪作的10略高于麥豆輪作, 可能的原因是由于第2輪作周期內(nèi), 冬小麥生育后期以及后茬作物生育期內(nèi), 麥豆輪作的土壤溫度均顯著高于麥玉輪作, 較低的溫度能夠豐富土壤微生物的種類, 進(jìn)而潛在地提高10[52]。本研究結(jié)果還顯示, 秸稈還田處理的10低于秸稈不還田處理, 這與前人的研究結(jié)果一致[53-54]。導(dǎo)致這一現(xiàn)象的原因可能是由于外界氣候條件發(fā)生變化所致, 前人在研究中指出溫度升高及干旱都會引起土壤呼吸溫度敏感性降低[55-56]。在本研究周期內(nèi), 由于試驗區(qū)域的夏季溫度較高, 且夏玉米及夏大豆生育期內(nèi)根據(jù)當(dāng)年的天氣情況僅灌溉一次, 可能引發(fā)潛在的干旱脅迫, 進(jìn)而降低土壤呼吸的溫度敏感性。由于本研究屬于長期定位試驗的一個部分, 研究周期也僅為兩年, 所得結(jié)論還需進(jìn)一步的驗證, 今后有關(guān)秸稈還田對土壤呼吸的溫度敏感性長期影響的試驗將進(jìn)一步完善和開展。

    2012—2014連續(xù)兩年內(nèi)不同輪作模式下秸稈還田對土壤呼吸及其溫度敏感性的研究結(jié)果顯示, 不同作物生育期內(nèi), 土壤呼吸均存在明顯的季節(jié)變化趨勢, 且在生育期的大部分時間內(nèi), SM處理的土壤呼吸速率均高于NS處理, 相比NS處理, SM顯著提高了作物生育期土壤呼吸平均速率和生育期土壤呼吸累計釋放量。不同輪作模式下, 麥玉輪作的生育期土壤呼吸累計釋放量高于麥豆輪作。秸稈還田具有調(diào)節(jié)土壤溫度的作用, 在冬季能夠提高土壤的溫度, 但在春季和夏季表現(xiàn)為降低土壤溫度。對于土壤含水量而言, 在高溫少雨的時期內(nèi), 秸稈還田能夠有效提高土壤含水量。土壤溫度與土壤呼吸間存在顯著的正相關(guān)關(guān)系, 地下5 cm和地下10 cm的土壤溫度變化能夠分別解釋土壤呼吸變化的64.6%~67.3%和51.5%~59.6%, 但土壤呼吸與地下0~30 cm土壤平均含水量間無顯著相關(guān)性。此外, 秸稈還田還能夠降低土壤呼吸的溫度敏感性, 對農(nóng)田生態(tài)系統(tǒng)碳循環(huán)對全球氣候變化響應(yīng)的研究具有一定實際意義。

    [1] Keeling C D. Climate change and carbon dioxide: An introduce[J]. Proceedings of the National Academy of Sciences of the United States of America, 1997, 94(16): 8273–8274

    [2] Bond-Lamberty B, Thomson A. A global database of soil respiration data[J]. Biogeosciences, 2010, 7(6): 1915–1926

    [3] Hibbard K A, Law B E, Reichstein M, et al. An analysis of soil respiration across northern hemisphere temperate ecosystems[J]. Biogeochemistry, 2005, 73(1): 29–70

    [4] Schlesinger W H, Andrews J A. Soil respiration and the global carbon cycle[J]. Biogeochemistry, 2000, 48(1): 7–20

    [5] Li H, Qiu J J, Wang L G, et al. Modelling impacts of alternative farming management practices on greenhouse gas emissions from a winter wheat–maize rotation system in China[J]. Agriculture, Ecosystems & Environment, 2010, 135(1/2): 24–33

    [6] 段華平, 張悅, 趙建波, 等. 中國農(nóng)田生態(tài)系統(tǒng)的碳足跡分析[J]. 水土保持學(xué)報, 2011, 25(5): 203–208 Duan H P, Zhang Y, Zhao J B, et al. Carbon footprint analysis of farmland ecosystem in China[J]. Journal of Soil and Water Conservation, 2011, 25(5): 203–208

    [7] Raich J W, Tufekciogul A. Vegetation and soil respiration: Correlations and controls[J]. Biogeochemistry, 2000, 48(1): 71–90

    [8] Lal R. Soil carbon sequestration impacts on global climate change and food security[J]. Science, 2004, 304(5677): 1623–1627

    [9] Luo Y Q, Zhou X H. 土壤呼吸與環(huán)境[M]. 姜麗芬, 曲來葉, 周玉梅, 等, 譯. 北京: 高等教育出版社, 2007 Luo Y Q, Zhou X H. Soil Respiration and Environment[M]. Jiang L F, Qu L Y, Zhou Y M, et al., trans. Beijing: Higher Education Press, 2007

    [10] 李海防, 夏漢平, 熊燕梅, 等. 土壤溫室氣體產(chǎn)生與排放影響因素研究進(jìn)展[J]. 生態(tài)環(huán)境, 2007, 16(6): 1781–1788 Li H F, Xia H P, Xiong Y M, et al. Mechanism of greenhouse gases fluxes from soil and its controlling factors: A review[J]. Ecology and Environment, 2007, 16(6): 1781–1788

    [11] 南雄雄, 游東海, 田霄鴻, 等. 關(guān)中平原農(nóng)田作物秸稈還田對土壤有機碳和作物產(chǎn)量的影響[J]. 華北農(nóng)學(xué)報, 2011, 26(5): 222–229 Nan X X, You D H, Tian X H, et al. Effect of returning of cropland straw to field on soil organic carbon and grain yield in Guanzhong plain[J]. Acta Agriculturae Boreali-Sinica, 2011, 26(5): 222–229

    [12] Buman R A, Alesii B A, Hatfield J L, et al. Profit, yield, and soil quality effects of tillage systems in corn-soybean rotations[J]. Journal of Soil and Water Conservation, 2004, 59(6): 260–270

    [13] 于愛忠, 黃高寶, 柴強. 不同耕作措施對西北綠洲灌區(qū)冬小麥農(nóng)田土壤呼吸的影響[J]. 草業(yè)學(xué)報, 2012, 21(1): 273–278 Yu A Z, Huang G B, Chai Q. Effect of different tillage treatments on soil respiration of winter-wheat farmland in oasis irrigated area Northwest China[J]. Acta Prataculturae Sinica, 2012, 21(1): 273–278

    [14] 崔鳳娟, 李立軍, 劉景輝, 等. 免耕留茬覆蓋對土壤呼吸和土壤酶活性及養(yǎng)分的影響[J]. 中國農(nóng)學(xué)通報, 2011, 27(21): 147–153 Cui F J, Li L J, Liu J H, et al. No-tillage stubble with residues on soil respiration and the soil enzyme activity and nutrient influence[J]. Chinese Agricultural Science Bulletin, 2011, 27(21): 147–153

    [15] Irvine J, Law B E. Contrasting soil respiration in young and old-growth ponderosas pine forests[J]. Global Change Biology, 2002, 8(12): 1183–1194

    [16] Zhou X H, Sherry R A, An Y, et al. Main and interactive effects of warming, clipping, and doubled precipitation on soil CO2efflux in a grassland ecosystem[J]. Global Biogeochemical Cycles, 2006, 20(1): GB1003

    [17] 楊毅, 黃玫, 劉洪升, 等. 土壤呼吸的溫度敏感性和適應(yīng)性研究進(jìn)展[J]. 自然資源學(xué)報, 2011, 26(10): 1811–1820 Yang Y, Huang M, Liu H S, et al. A The interrelation between temperature sensitivity and adaptability of soil respiration[J]. Journal of Natural Resources, 2011, 26(10): 1811–1820

    [18] 楊慶朋, 徐明, 劉洪升, 等. 土壤呼吸溫度敏感性的影響因素和不確定性[J]. 生態(tài)學(xué)報, 2011, 31(8): 2301–2311 Yang Q P, Xu M, Liu H S, et al. Impact factors and uncertainties of the temperature sensitivity of soil respiration[J]. Acta Ecologica Sinica, 2011, 31(8): 2301–2311

    [19] Raich J W, Schlesinger W H. The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate[J]. Tellus B: Chemical and Physical Meteorology, 1992, 44(2): 81–99

    [20] Davidson E A, Janssens I A, Luo Y Q. On the variability of respiration in terrestrial ecosystems: Moving beyond10[J]. Global Change Biology, 2006, 12(2): 154–164

    [21] Zheng Z M, Yu G R, Fu Y L, et al. Temperature sensitivity of soil respiration is affected by prevailing climatic conditions and soil organic carbon content: A trans-China based case study[J]. Soil Biology and Biochemistry, 2009, 41(7): 1531–1540

    [22] Peng Q, Dong Y S, Qi Y C, et al. Effects of nitrogen fertilization on soil respiration in temperate grassland in Inner Mongolia, China[J]. Environmental Earth Sciences, 2011, 62(6): 1163–1171

    [23] Wang R, Sun Q Q, Wang Y, et al. Temperature sensitivity of soil respiration: Synthetic effects of nitrogen and phosphorus fertilization on Chinese Loess Plateau[J]. Science of the Total Environment, 2017, 574: 1665–1673

    [24] 劉爽, 嚴(yán)昌榮, 何文清, 等. 不同耕作措施下旱地農(nóng)田土壤呼吸及其影響因素[J]. 生態(tài)學(xué)報, 2010, 30(11): 2919–2924 Liu S, Yan C R, He W Q, et al. Soil respiration and it’s affected factors under different tillage systems in dryland production systems[J]. Acta Ecologica Sinica, 2010, 30(11): 2919–2924

    [25] 姜繼韶, 郭勝利, 王蕊, 等. 施氮對黃土旱塬區(qū)春玉米土壤呼吸和溫度敏感性的影響[J]. 環(huán)境科學(xué), 2015, 36(5): 1802–1809 Jiang J S, Guo S L, Wang R, et al. Effects of nitrogen fertilization on soil respiration and temperature sensitivity in spring maize field in semi-arid regions on Loess Plateau[J]. Environmental Science, 2015, 36(5): 1802–1809

    [26] 代快, 蔡典雄, 王燕, 等. 不同耕作措施對旱作春玉米農(nóng)田土壤呼吸影響的研究——土壤溫度對土壤呼吸速率的影響[J]. 中國土壤與肥料, 2010, (6): 64–69 Dai K, Cai D X, Wang Y, et al. Effects of different tillage on soil respiration rate of spring corn in the dryland of North China — Effects of soil temperature on soil respiration rate[J]. Soils and Fertilizer Sciences in China, 2010, (6): 64–69

    [27] 張賽, 羅海秀, 王龍昌, 等. 保護(hù)性耕作下大豆農(nóng)田土壤呼吸及影響因素分析[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2013, 21(8): 913–920 Zhang S, Luo H X, Wang L C, et al. Analysis of soil respiration and the influencing factors in soybean fields under conservation tillage[J]. Chinese Journal of Eco-Agriculture, 2013, 21(8): 913–920

    [28] 涂純, 王俊, 劉文兆. 不同覆蓋條件下旱作農(nóng)田土壤呼吸及其影響因素[J]. 植物營養(yǎng)與肥料學(xué)報, 2012, 18(5): 1107–1114 Tu C, Wang J, Liu W Z. Variation in soil respiration and its driving factors in rainfed winter wheat fields with different mulching measures[J]. Plant Nutrition and Fertilizer Science, 2012, 18(5): 1107–1114

    [29] Zhai L M, Liu H B, Zhang J Z, et al. Long-term application of organic manure and mineral fertilizer on N2O and CO2emissions in a red soil from cultivated maize-wheat rotation in China[J]. Agricultural Sciences in China, 2011, 10(11): 1748–1757

    [30] Lloyd J, Taylor J A. On the temperature dependence of soil respiration[J]. Functional Ecology, 1994, 8(3): 315–323

    [31] 陳繼康, 李素娟, 張宇, 等. 不同耕作方式麥田土壤溫度及其對氣溫的響應(yīng)特征——土壤溫度日變化及其對氣溫的響應(yīng)[J]. 中國農(nóng)業(yè)科學(xué), 2009, 42(7): 2592–2600Chen J K, Li S J, Zhang Y, et al. Characteristics of soil temperature and response to air temperature under different tillage systems — Diurnal dynamic of soil temperature and its response to air temperature[J]. Scientia Agricultura Sinica, 2009, 42(7): 2592–2600

    [32] 付國占, 李潮海, 王俊忠, 等. 殘茬覆蓋與耕作方式對土壤性狀及夏玉米水分利用效率的影響[J]. 農(nóng)業(yè)工程學(xué)報, 2005, 21(1): 52–56 Fu G Z, Li C H, Wang J Z, et al. Effects of stubble mulch and tillage managements on soil physical properties and water use efficiency of summer maize[J]. Transactions of the CSAE, 2005, 21(1): 52–56

    [33] 戴云仙, 張向前, 郭曉霞, 等. 免耕輪作對旱作區(qū)土壤水熱動態(tài)及速效養(yǎng)分含量的影響[J]. 中國土壤與肥料, 2016, (6): 14–20 Dai Y X, Zhang X Q, Guo X X, et al. Effects of no tillage with rotation on soil hydrothermal dynamics and available nutrients in dry farming area[J]. Soil and Fertilizer Sciences in China, 2016, (6): 14–20

    [34] Wang X, Jia Z, Liang L. Effect of straw incorporation on soil moisture, evapotranspiration, and rainfall-use efficiency of maize under dryland farming[J]. Journal of Soil and Water Conservation, 2014, 69(5): 449–455

    [35] 王暉, 劉泉汝, 張圣勇, 等. 秸稈覆蓋下超高產(chǎn)夏玉米農(nóng)田產(chǎn)量和土壤水分的動態(tài)變化[J]. 水土保持學(xué)報, 2011, 25(5): 261–264 Wang H, Liu Q R, Zhang S Y, et al. Grain yield and soil water content of super-high-yield summer maize under straw mulching[J]. Journal of Soil and Water Conservation, 2011, 25(5): 261–264

    [36] Kar G, Verma H N, Singh R. Effects of winter crop and supplemental irrigation on crop yield, water use efficiency and profitability in rainfed rice based cropping system of eastern India[J]. Agricultural Water Management, 2006, 79(3): 280–292

    [37] 李曉莎, 武寧, 劉玲, 等. 不同秸稈還田和耕作方式對夏玉米農(nóng)田土壤呼吸及微生物活性的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2015, 26(6): 1765–1771 Li X S, Wu N, Liu L, et al. Effects of different straw recycling and tillage methods on soil respiration and microbial activity[J]. Chinese Journal of Applied Ecology, 2015, 26(6): 1765–1771

    [38] 祿興麗, 廖允成. 不同耕作措施對旱作夏玉米田土壤呼吸及根呼吸的影響[J]. 環(huán)境科學(xué), 2015, 36(6): 2266–2273 Lu X L, Liao Y C. Effects of tillage on soil respiration and root respiration under rain-fed summer corn field[J]. Environmental Science, 2015, 36(6): 2266–2273

    [39] Blanco-Canqui H, Lal R. Soil structure and organic carbon relationships following 10 years of wheat straw management in no-till[J]. Soil and Tillage Research, 2007, 95(1/2): 240–254

    [40] 劉曉芳, 李立軍, 閆睿卿, 等. 內(nèi)蒙古旱作農(nóng)田土壤有機碳對免耕留茬的響應(yīng)[J]. 中國農(nóng)學(xué)通報, 2015, 31(7): 218–222 Liu X F, Li L J, Yan R Q, et al. Response of soil organic carbon of rain fed cropland to stubble tillage in Inner Mongolia[J]. Chinese Agricultural Science Bulletin, 2015, 31(7): 218–222

    [41] 孫媛, 任廣鑫, 馮永忠, 等. 秸稈還田和施氮對土壤水熱因子及呼吸速率的影響[J]. 西北農(nóng)林科技大學(xué)學(xué)報: 自然科學(xué)版, 2015, 43(3): 146–152Sun Y, Ren G X, Feng Y Z, et al. Comprehensive influence of straw-returning and nitrogen fertilization on hydrothermal factors and soil respiration[J]. Journal of Northwest A&F University: Natural Science Edition, 2015, 43(3): 146–152

    [42] 王維鈺, 喬博, Akhtar K, 等. 免耕條件下秸稈還田對冬小麥-夏玉米輪作系統(tǒng)土壤呼吸及土壤水熱狀況的影響[J]. 中國農(nóng)業(yè)科學(xué), 2016, 49(11): 2136–2152 Wang W Y, Qiao B, Akhtar K, et al. Effects of straw returning to field on soil respiration and soil water heat in winter wheat-summer maize rotation system under no tillage[J]. Scientia Agricultura Sinica, 2016, 49(11): 2136–2152

    [43] McDaniel M D, Grandy A S, Tiemann L K, et al. Crop rotation complexity regulates the decomposition of high and low quality residues[J]. Soil Biology and Biochemistry, 2014, 78: 243–254

    [44] Jiang J S, Guo S L, Zhang Y J, et al. Changes in temperature sensitivity of soil respiration in the phases of a three-year crop rotation system[J]. Soil and Tillage Research, 2015, 150: 139–146

    [45] 張慶忠, 吳文良, 王明新, 等. 秸稈還田和施氮對農(nóng)田土壤呼吸的影響[J]. 生態(tài)學(xué)報, 2005, 25(11): 2883–2887 Zhang Q Z, Wu W L, Wang M X, et al. The effects of crop residue amendment and N rate on soil respiration[J]. Acta Ecologica Sinica, 2005, 25(11): 2883–2887

    [46] 強學(xué)彩, 袁紅莉, 高旺盛. 秸稈還田量對土壤CO2釋放和土壤微生物量的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2004, 15(3): 469–472 Qiang X C, Yuan H L, Gao W S. Effect of crop-residue incorporation on soil CO2emission and soil microbial biomass[J]. Chinese Journal of Applied Ecology, 2004, 15(3): 469–472

    [47] 李昌珍, 張婷婷, 楊改河, 等. 秸稈覆蓋和施肥對關(guān)中灌區(qū)夏玉米生長后期土壤呼吸速率的影響[J]. 生態(tài)環(huán)境學(xué)報, 2013, 22(3): 411–416 Li C Z, Zhang T T, Yang G H, et al. Effect of straw mulching and fertilization on soil respiration from summer-maize farmland during later growth period in Guanzhong irrigation district[J]. Ecology and Environmental Sciences, 2013, 22(3): 411–416

    [48] 張俊麗, 廖允成, 曾愛, 等. 不同施氮水平下旱作玉米田土壤呼吸速率與土壤水熱關(guān)系[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2013, 32(7): 1382–1388 Zhang J L, Liao Y C, Zeng A, et al. 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[J]. Journal of Agro-Environment Science, 2013, 32(7): 1382–1388

    [49] 楊玲, 廖允成, 高茂盛, 等. 不同耕作措施下旱作麥田CO2排放速率與土壤水熱關(guān)系分析[J]. 西北農(nóng)業(yè)學(xué)報, 2011, 20(1): 70–75 Yang L, Liao Y C, Gao M S, et al. Relationship analysis of CO2emission fluxes and soil temperature and moisture in rainfed wheat field with different tillage[J]. Acta Agriculturae Boreali-Occidentalis Sinica, 2011, 20(1): 70–75

    [50] 張宇, 張海林, 陳繼康, 等. 耕作措施對華北農(nóng)田CO2排放影響及水熱關(guān)系分析[J]. 農(nóng)業(yè)工程學(xué)報, 2009, 25(4): 47–53 Zhang Y, Zhang H L, Chen J K, et al. Effects of different tillage practices on CO2emission fluxes from farmland in North China Plain and the analysis of soil temperature and moisture[J]. Transactions of the CSAE, 2009, 25(4): 47–53

    [51] Peng S S, Piao S L, Wang T, et al. Temperature sensitivity of soil respiration in different ecosystems in China[J]. Soil Biology and Biochemistry, 2009, 41(5): 1008–1014

    [52] Chen B Y, Liu S R, Ge J P, et al. Annual and seasonal variations of10soil respiration in the sub-alpine forests of the Eastern Qinghai-Tibet Plateau, China[J]. Soil Biology and Biochemistry, 2010, 42(10): 1735–1742

    [53] He X H, Du Z L, Wang Y D, et al. Sensitivity of soil respiration to soil temperature decreased under deep biochar amended soils in temperate croplands[J]. Applied Soil Ecology, 2016, 108: 204–210

    [54] Hu Z H, Ling H, Chen S T, et al. Soil respiration, nitrification, and denitrification in a wheat farmland soil under different managements[J]. Communications in Soil Science and Plant Analysis, 2013, 44(21): 3092–3102

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

    [56] Almagro M, López J, Querejeta J I, et al. Temperature dependence of soil CO2efflux is strongly modulated by seasonal patterns of moisture availability in a Mediterranean ecosystem[J]. Soil Biology and Biochemistry, 2009, 41(3): 594–605

    Effect of straw mulching on soil respiration and its’ temperature sensitivity under different crop rotation systems*

    DING Ruixia1, WANG Weiyu2,3, ZHANG Qing2

    (1. Institute of Water-saving Agriculture in Arid Areas of China, Northwest A&F University / Key Laboratory of Crop Physiology and Tillage Sciences in Northwestern Loess Plateau, Ministry of Agriculture, Yangling 712100, China; 2. College of Agronomy, Northwest A & F University, Yangling 712100, China; 3. Research Center of Recycle Agricultural Engineering and Technology of Shaanxi Province, Yangling 712100, China)

    Soil respiration is the second largest source of terrestrial carbon (C) flux between the atmosphere and the terrestrial ecosystems. It is critical for regulating global soil C dynamics. As soil temperature and soil moisture would exert stronger effects on soil respiration in the future, a thorough understanding of the response of soil microbes to temperature change can provide a novel method of studying the effects of drought on soil respiration and of predicting drought-induced changes in future terrestrial C cycle. Temperature sensitivity of soil respiration can explain the relationship between soil respiration and soil temperature. The objective of this study was to explore the effects of straw mulch on the linkages between the changes in soil respiration and temperature. The study aimed to lay the basis of C cycle process in agro-ecosystems in Yangling, Shaanxi province. To that end, a 2-year field experiment (October 2012 to September 2014) was conducted to study the linkages under different crop rotation systems. It included two treatments of no straw (NS) and straw mulch (SM) in winter wheat-summer maize rotation and winter wheat-summer soybean rotation systems. Soil respiration rate, temperature, and moisture were analyzed under different crop rotation systems. In addition, the10(with10value as the multiplier in determining soil respiration rate after temperature increase of 10 ℃) was used to determine the effect of soil temperature change on soil respiration. It was noted that SM significantly (< 0.05) increased soil respiration rate during crop growth period. Mean soil respiration rate and cumulative soil respiration during crop growth period significantly increased under straw mulch (< 0.01). The order of mean soil respiration rate under various crops was as follows: maize (3.401–4.810 μmol×m-2×s-1) > soybean (3.390–3.762 μmol×m-2×s-1) > wheat (2.673–3.141 μmol×m-2×s-1). Then the order of cumulative soil respiration among different rotations was as follows: wheat-maize [34.68–40.81 t(CO2)×hm-2] > wheat-soybean [30.04–33.86 t(CO2)×hm-2]. In addition, soil temperature varied significantly (< 0.05) among different crops. Particularly, soil temperature under wheat-maize rotation system was higher than that under wheat-soybean rotation system during the growth stage of wheat. Soil temperature at 5 cm soil depth in maize field was higher than that in soybean field during the summer of 2014. It was noted that SM treatment was a major regulator of soil temperature — significantly increasing it in winter and then significantly decreasing it in spring and summer. Moreover, mean soil moisture content in the 0–30 cm soil layer was significantly higher under SM treatment than under NS treatment during the dry season. Further, mean soil moisture content in the 0–30 cm layer during wheat growth period varied significantly among different crop rotation systems, which was associated with root characteristics under crop rotation. The mean soil moisture content for the 0–30 cm soil layer during maize growth period was significantly higher than that during soybean growth period. Soil temperature at the 5–10 cm soil layer was positively correlated with soil respiration. However, soil moisture at the 0–30 cm soil layer was not significantly correlated with soil respiration. Changes in soil temperature at 5 and 10 cm soil depths were respectively 64.6%–67.3% and 51.5%–59.6% explained by the variance in soil respiration. In this study,10varied within 1.70?2.01 across different crop rotation systems and was significantly higher under wheat-maize than wheat-soybean rotation system. In addition,10was significantly higher under SM treatment than under NS treatment. Therefore, SM treatment was more advantageous in terms of the ability to effectively reduce temperature sensitivity of soil respiration and to accurately predict soil moisture and soil heat conditions.

    Crop rotation system; Straw mulch; Soil respiration; Temperature sensitivity; Soil moisture and heat condition

    10.13930/j.cnki.cjea.170078

    S344.13

    A

    1671-3990(2017)08-1106-13

    *陜西省科技統(tǒng)籌計劃項目(2015KTCL02-07)資助

    丁瑞霞, 主要從事旱區(qū)農(nóng)業(yè)研究。E-mail: dingruixia@sohu.com

    2017-01-21 接受日期: 2017-04-05

    *This work was funded by the Science and Technology Overall Planning Project in Shaanxi Province of China (2015KTCL02-07).

    , DING Ruixia, E-mail: dingruixia@sohu.com

    Jan. 21, 2017; accepted Apr. 5, 2017

    丁瑞霞, 王維鈺, 張青. 兩種輪作模式下秸稈還田對土壤呼吸及其溫度敏感性的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2017, 25(8): 1106-1118

    Ding R X, Wang W Y, Zhang Q. Effect of straw mulching on soil respiration and its’ temperature sensitivity under different crop rotation systems[J]. Chinese Journal of Eco-Agriculture, 2017, 25(8): 1106-1118

    猜你喜歡
    夏大豆土壤溫度夏玉米
    不同年份氣象因子對夏大豆生長及產(chǎn)量的影響
    種子科技(2024年4期)2024-03-28 12:00:58
    菏澤市夏大豆高產(chǎn)穩(wěn)產(chǎn)栽培技術(shù)
    優(yōu)質(zhì)夏大豆綠色高效種植模式思考研究
    西部論叢(2020年6期)2020-08-16 10:41:32
    小麥?zhǔn)崭钪?如何種植夏玉米才能高產(chǎn)
    夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
    遼東山區(qū)3種人工林土壤呼吸對土壤溫度和土壤水分的響應(yīng)
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    黃淮海地區(qū)高蛋白夏大豆高產(chǎn)栽培技術(shù)
    大豆科技(2016年6期)2016-03-11 15:30:22
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    稻草還田對油菜生長、土壤溫度及濕度的影響
    一边摸一边抽搐一进一小说| 欧美日本视频| 成人永久免费在线观看视频| 麻豆久久精品国产亚洲av| 国产在线观看jvid| 国产91精品成人一区二区三区| 听说在线观看完整版免费高清| 毛片女人毛片| 精品久久久久久久久久免费视频| 欧美绝顶高潮抽搐喷水| 精华霜和精华液先用哪个| 久久久久久久午夜电影| 真人一进一出gif抽搐免费| 熟女少妇亚洲综合色aaa.| 美女扒开内裤让男人捅视频| 成人国产综合亚洲| 精华霜和精华液先用哪个| 久久中文字幕一级| 首页视频小说图片口味搜索| 超碰成人久久| 男女下面进入的视频免费午夜| 舔av片在线| 久久午夜亚洲精品久久| 久久伊人香网站| 欧美 亚洲 国产 日韩一| 国产av在哪里看| av有码第一页| 免费在线观看亚洲国产| 午夜日韩欧美国产| 免费看a级黄色片| 亚洲精品美女久久av网站| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 亚洲人成网站在线播放欧美日韩| 久久天堂一区二区三区四区| 国产亚洲av嫩草精品影院| 国内精品久久久久精免费| 久久午夜综合久久蜜桃| 全区人妻精品视频| 亚洲国产欧洲综合997久久,| 国产欧美日韩一区二区精品| 高潮久久久久久久久久久不卡| 久久国产精品影院| 亚洲成人国产一区在线观看| 午夜精品在线福利| 啦啦啦韩国在线观看视频| 精品一区二区三区视频在线观看免费| 亚洲av成人av| 精品一区二区三区四区五区乱码| 搞女人的毛片| 亚洲av成人不卡在线观看播放网| 中文字幕熟女人妻在线| 国产精品亚洲av一区麻豆| 麻豆一二三区av精品| 午夜福利高清视频| 国产一区二区三区在线臀色熟女| 国产日本99.免费观看| 波多野结衣巨乳人妻| ponron亚洲| 亚洲精品粉嫩美女一区| 一进一出抽搐动态| 真人做人爱边吃奶动态| 又粗又爽又猛毛片免费看| 日韩三级视频一区二区三区| 国产视频内射| 搞女人的毛片| 草草在线视频免费看| 亚洲一码二码三码区别大吗| 国产成人啪精品午夜网站| 亚洲午夜精品一区,二区,三区| 五月伊人婷婷丁香| 久久精品91无色码中文字幕| 欧美久久黑人一区二区| 黄色 视频免费看| 美女高潮喷水抽搐中文字幕| 免费在线观看黄色视频的| 精品第一国产精品| 午夜福利高清视频| 久久久国产精品麻豆| 国产麻豆成人av免费视频| 一级毛片精品| 成人三级黄色视频| 亚洲五月天丁香| 日本在线视频免费播放| 精品一区二区三区视频在线观看免费| 天堂√8在线中文| 日韩有码中文字幕| 好看av亚洲va欧美ⅴa在| 国产男靠女视频免费网站| 日韩欧美三级三区| 五月玫瑰六月丁香| 免费观看人在逋| 国产精品久久视频播放| 男插女下体视频免费在线播放| 人妻丰满熟妇av一区二区三区| 久久人妻福利社区极品人妻图片| 亚洲国产看品久久| 日韩欧美 国产精品| 18禁观看日本| av天堂在线播放| 午夜福利18| 亚洲中文字幕一区二区三区有码在线看 | 男女做爰动态图高潮gif福利片| 波多野结衣高清无吗| 亚洲av成人精品一区久久| 99久久综合精品五月天人人| 老司机午夜福利在线观看视频| 免费在线观看亚洲国产| 久久久久精品国产欧美久久久| 淫妇啪啪啪对白视频| 怎么达到女性高潮| 国产高清videossex| 欧美一级a爱片免费观看看 | 欧美日韩亚洲综合一区二区三区_| 视频区欧美日本亚洲| 手机成人av网站| 欧美日韩精品网址| 精品福利观看| 久久婷婷成人综合色麻豆| 国产男靠女视频免费网站| 精品久久久久久,| 99国产综合亚洲精品| 国产日本99.免费观看| 看黄色毛片网站| 国产一级毛片七仙女欲春2| 91字幕亚洲| 少妇粗大呻吟视频| 久久欧美精品欧美久久欧美| 欧美最黄视频在线播放免费| 岛国在线观看网站| 久久国产精品影院| www日本在线高清视频| 午夜视频精品福利| 悠悠久久av| 午夜福利欧美成人| 在线视频色国产色| 无遮挡黄片免费观看| 国产爱豆传媒在线观看 | 亚洲欧美一区二区三区黑人| 午夜亚洲福利在线播放| 操出白浆在线播放| 在线观看日韩欧美| 曰老女人黄片| 精品乱码久久久久久99久播| 免费在线观看黄色视频的| 欧美日韩国产亚洲二区| 久久天堂一区二区三区四区| 桃红色精品国产亚洲av| 高潮久久久久久久久久久不卡| 大型黄色视频在线免费观看| 久久久久久国产a免费观看| 一级作爱视频免费观看| 欧美日韩国产亚洲二区| 在线永久观看黄色视频| 国产麻豆成人av免费视频| 超碰成人久久| 宅男免费午夜| 一进一出抽搐动态| 亚洲自偷自拍图片 自拍| 国产又色又爽无遮挡免费看| 久久九九热精品免费| 色尼玛亚洲综合影院| 久久久国产成人精品二区| 色综合亚洲欧美另类图片| 18禁国产床啪视频网站| 日本a在线网址| 成年女人毛片免费观看观看9| 免费一级毛片在线播放高清视频| 一区二区三区国产精品乱码| 小说图片视频综合网站| 国内少妇人妻偷人精品xxx网站 | 亚洲色图av天堂| 日韩中文字幕欧美一区二区| 女同久久另类99精品国产91| 欧美一区二区精品小视频在线| 一进一出好大好爽视频| 琪琪午夜伦伦电影理论片6080| 亚洲国产欧美一区二区综合| 可以在线观看的亚洲视频| 校园春色视频在线观看| 一级毛片高清免费大全| 精品少妇一区二区三区视频日本电影| 天天躁夜夜躁狠狠躁躁| 一a级毛片在线观看| 露出奶头的视频| 夜夜夜夜夜久久久久| 天天一区二区日本电影三级| 国产亚洲欧美在线一区二区| 国内久久婷婷六月综合欲色啪| 999精品在线视频| 露出奶头的视频| 成熟少妇高潮喷水视频| 久久久久亚洲av毛片大全| 51午夜福利影视在线观看| 亚洲天堂国产精品一区在线| 亚洲最大成人中文| 婷婷亚洲欧美| 亚洲国产欧美人成| 老鸭窝网址在线观看| 黄片大片在线免费观看| 性欧美人与动物交配| 黑人欧美特级aaaaaa片| 欧美精品啪啪一区二区三区| 悠悠久久av| 又粗又爽又猛毛片免费看| 日韩 欧美 亚洲 中文字幕| 国产成年人精品一区二区| 一边摸一边抽搐一进一小说| 欧美极品一区二区三区四区| 国产成人av激情在线播放| 亚洲欧美精品综合久久99| 黄色毛片三级朝国网站| 在线国产一区二区在线| 男人舔女人的私密视频| or卡值多少钱| 日本三级黄在线观看| 欧美日韩黄片免| 免费看美女性在线毛片视频| 免费在线观看日本一区| 久久精品aⅴ一区二区三区四区| 女警被强在线播放| 亚洲欧美精品综合一区二区三区| 露出奶头的视频| 级片在线观看| 日韩成人在线观看一区二区三区| 黄色片一级片一级黄色片| 欧美不卡视频在线免费观看 | 正在播放国产对白刺激| 国产野战对白在线观看| 禁无遮挡网站| 国产乱人伦免费视频| 久久精品影院6| 91麻豆av在线| 亚洲精品一区av在线观看| 丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲人成伊人成综合网2020| 国产精品98久久久久久宅男小说| 我要搜黄色片| 国产99久久九九免费精品| 一级a爱片免费观看的视频| 男插女下体视频免费在线播放| 无限看片的www在线观看| 国产成年人精品一区二区| 欧美3d第一页| 日韩欧美 国产精品| 国产成人精品无人区| 日韩大尺度精品在线看网址| 宅男免费午夜| 欧美丝袜亚洲另类 | 国产三级在线视频| 成年人黄色毛片网站| www日本在线高清视频| 日本精品一区二区三区蜜桃| 欧美极品一区二区三区四区| 2021天堂中文幕一二区在线观| 午夜免费激情av| 久久精品亚洲精品国产色婷小说| or卡值多少钱| 国产99久久九九免费精品| av在线播放免费不卡| 精品无人区乱码1区二区| 在线a可以看的网站| 91av网站免费观看| 中文字幕久久专区| 婷婷六月久久综合丁香| 亚洲人与动物交配视频| 91九色精品人成在线观看| 精品久久久久久,| 午夜免费观看网址| 丝袜人妻中文字幕| 亚洲第一电影网av| 日本熟妇午夜| 精品久久久久久,| а√天堂www在线а√下载| 国产v大片淫在线免费观看| 久久久久免费精品人妻一区二区| 高潮久久久久久久久久久不卡| 两个人的视频大全免费| 床上黄色一级片| 热99re8久久精品国产| 日韩欧美在线乱码| 国产一区在线观看成人免费| 叶爱在线成人免费视频播放| 丝袜人妻中文字幕| 特大巨黑吊av在线直播| 成人欧美大片| 窝窝影院91人妻| 国产亚洲精品久久久久久毛片| 欧美黑人精品巨大| 91av网站免费观看| a级毛片在线看网站| 精品欧美一区二区三区在线| 亚洲午夜理论影院| 亚洲精品一区av在线观看| 久久中文字幕一级| 99热6这里只有精品| 国产精品野战在线观看| 丁香欧美五月| 欧美黑人巨大hd| 久久国产精品人妻蜜桃| 色综合亚洲欧美另类图片| 欧美丝袜亚洲另类 | 男人舔奶头视频| 波多野结衣高清作品| 午夜福利在线在线| 久久亚洲精品不卡| 欧美性猛交黑人性爽| 亚洲欧美激情综合另类| 欧美日本视频| 国产一区二区在线观看日韩 | 欧美在线一区亚洲| 国产激情久久老熟女| 亚洲 欧美 日韩 在线 免费| 亚洲精品久久成人aⅴ小说| 天堂√8在线中文| 可以在线观看的亚洲视频| 亚洲avbb在线观看| 亚洲欧美日韩东京热| 亚洲 欧美 日韩 在线 免费| 国产三级在线视频| 久久久精品欧美日韩精品| 亚洲一卡2卡3卡4卡5卡精品中文| 日日摸夜夜添夜夜添小说| 精品福利观看| 国产在线观看jvid| 亚洲九九香蕉| 亚洲成a人片在线一区二区| 国产精品免费视频内射| 色综合站精品国产| 国产97色在线日韩免费| 天天一区二区日本电影三级| 男女床上黄色一级片免费看| 午夜日韩欧美国产| 好男人电影高清在线观看| 国产一区二区激情短视频| 国产亚洲av嫩草精品影院| 在线播放国产精品三级| 久久国产精品影院| 国产精品一及| 日韩av在线大香蕉| 久久久久国产一级毛片高清牌| 人人妻人人看人人澡| 夜夜看夜夜爽夜夜摸| 日本 av在线| 欧美成人一区二区免费高清观看 | 男女视频在线观看网站免费 | 亚洲 欧美 日韩 在线 免费| 黄片小视频在线播放| 真人一进一出gif抽搐免费| 久久热在线av| 亚洲 欧美 日韩 在线 免费| 亚洲国产欧美人成| 可以免费在线观看a视频的电影网站| 熟妇人妻久久中文字幕3abv| 少妇被粗大的猛进出69影院| 精品一区二区三区四区五区乱码| 欧美不卡视频在线免费观看 | 校园春色视频在线观看| 国产午夜福利久久久久久| 国产成年人精品一区二区| 免费在线观看成人毛片| 动漫黄色视频在线观看| 麻豆国产97在线/欧美 | 高清在线国产一区| 国产野战对白在线观看| 在线观看午夜福利视频| 欧美3d第一页| 老司机午夜十八禁免费视频| 91九色精品人成在线观看| 丁香欧美五月| 丰满人妻一区二区三区视频av | 亚洲狠狠婷婷综合久久图片| av在线天堂中文字幕| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站高清观看| 国产私拍福利视频在线观看| 亚洲精品国产一区二区精华液| 午夜福利18| 亚洲人成电影免费在线| 久久精品人妻少妇| 国产成人aa在线观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲国产欧美一区二区综合| 黑人操中国人逼视频| 久久婷婷人人爽人人干人人爱| 毛片女人毛片| 日韩三级视频一区二区三区| 999久久久精品免费观看国产| 美女免费视频网站| 一夜夜www| 久久人妻福利社区极品人妻图片| 成人亚洲精品av一区二区| 女同久久另类99精品国产91| 啦啦啦免费观看视频1| 久9热在线精品视频| 亚洲欧洲精品一区二区精品久久久| 最新美女视频免费是黄的| 亚洲国产高清在线一区二区三| 亚洲激情在线av| 91国产中文字幕| 在线观看舔阴道视频| 亚洲成人精品中文字幕电影| 午夜福利视频1000在线观看| 变态另类成人亚洲欧美熟女| 午夜精品在线福利| 午夜福利高清视频| 99riav亚洲国产免费| 又黄又粗又硬又大视频| 美女免费视频网站| 18禁黄网站禁片免费观看直播| 久久中文字幕人妻熟女| 桃色一区二区三区在线观看| 欧美日韩福利视频一区二区| 777久久人妻少妇嫩草av网站| 国产精品,欧美在线| 国产精品久久久久久人妻精品电影| 麻豆国产av国片精品| 亚洲天堂国产精品一区在线| 熟女电影av网| 舔av片在线| 亚洲av成人精品一区久久| 欧美在线黄色| 成人三级做爰电影| 久久精品夜夜夜夜夜久久蜜豆 | 欧美成人午夜精品| 国产av又大| 999精品在线视频| 后天国语完整版免费观看| 日本免费一区二区三区高清不卡| 亚洲人成伊人成综合网2020| 男女午夜视频在线观看| 成年女人毛片免费观看观看9| www日本黄色视频网| 国内久久婷婷六月综合欲色啪| 国产精品一及| 黄色毛片三级朝国网站| 99国产综合亚洲精品| 不卡av一区二区三区| 久久久久久久精品吃奶| 在线观看免费视频日本深夜| 在线观看免费午夜福利视频| 欧美日本亚洲视频在线播放| 午夜成年电影在线免费观看| 亚洲乱码一区二区免费版| 人妻久久中文字幕网| 欧美黑人精品巨大| 正在播放国产对白刺激| 国产精品九九99| 久久伊人香网站| 亚洲美女黄片视频| 亚洲精品色激情综合| 国产午夜精品久久久久久| a级毛片在线看网站| 国产成+人综合+亚洲专区| 欧美性猛交╳xxx乱大交人| 成人永久免费在线观看视频| 麻豆av在线久日| 女同久久另类99精品国产91| 操出白浆在线播放| 亚洲自拍偷在线| 国产激情久久老熟女| 午夜免费观看网址| 可以在线观看毛片的网站| 亚洲无线在线观看| 亚洲最大成人中文| 在线观看免费日韩欧美大片| 又粗又爽又猛毛片免费看| 狂野欧美激情性xxxx| 草草在线视频免费看| 国产又黄又爽又无遮挡在线| 麻豆国产97在线/欧美 | 五月伊人婷婷丁香| 正在播放国产对白刺激| 欧美一区二区国产精品久久精品 | 变态另类丝袜制服| 中文资源天堂在线| 男女视频在线观看网站免费 | 最新在线观看一区二区三区| 欧美黑人精品巨大| 国内精品久久久久精免费| 91成年电影在线观看| 国产伦一二天堂av在线观看| 男女视频在线观看网站免费 | 99久久精品国产亚洲精品| 日韩欧美精品v在线| 国产片内射在线| 日本精品一区二区三区蜜桃| netflix在线观看网站| 岛国在线免费视频观看| 午夜两性在线视频| 无人区码免费观看不卡| 日韩欧美精品v在线| 久久伊人香网站| 亚洲国产欧洲综合997久久,| 国产伦一二天堂av在线观看| 1024视频免费在线观看| 老熟妇仑乱视频hdxx| 欧美3d第一页| 亚洲,欧美精品.| 99国产精品99久久久久| 免费人成视频x8x8入口观看| 亚洲电影在线观看av| 欧美人与性动交α欧美精品济南到| 国内少妇人妻偷人精品xxx网站 | 搡老熟女国产l中国老女人| 国产高清有码在线观看视频 | 午夜免费观看网址| 一级毛片高清免费大全| 成人午夜高清在线视频| 999久久久精品免费观看国产| x7x7x7水蜜桃| 欧美三级亚洲精品| 免费在线观看成人毛片| 18禁观看日本| av片东京热男人的天堂| 国产三级在线视频| 少妇的丰满在线观看| 最近最新中文字幕大全电影3| 在线视频色国产色| 别揉我奶头~嗯~啊~动态视频| 一个人免费在线观看的高清视频| 看免费av毛片| av欧美777| 日韩欧美在线二视频| 五月伊人婷婷丁香| 久久伊人香网站| 不卡一级毛片| www.www免费av| 国产99白浆流出| 久久九九热精品免费| 国产精品一区二区三区四区免费观看 | 国产熟女xx| 亚洲性夜色夜夜综合| 亚洲av片天天在线观看| 国产精品自产拍在线观看55亚洲| 9191精品国产免费久久| 在线观看美女被高潮喷水网站 | 99精品在免费线老司机午夜| 美女高潮喷水抽搐中文字幕| 伦理电影免费视频| 色综合亚洲欧美另类图片| 国产av一区在线观看免费| 久久这里只有精品中国| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人中文字幕在线播放| 在线看三级毛片| 亚洲国产精品999在线| 久99久视频精品免费| 午夜激情av网站| 午夜福利免费观看在线| 免费一级毛片在线播放高清视频| 麻豆国产97在线/欧美 | 小说图片视频综合网站| 欧美zozozo另类| 9191精品国产免费久久| 亚洲无线在线观看| 亚洲成av人片在线播放无| a级毛片a级免费在线| 黄色毛片三级朝国网站| 亚洲av日韩精品久久久久久密| 可以在线观看毛片的网站| 久久精品综合一区二区三区| 99热这里只有是精品50| 91麻豆精品激情在线观看国产| 91麻豆av在线| 国产不卡一卡二| 日本在线视频免费播放| 久久99热这里只有精品18| 国产久久久一区二区三区| 国产91精品成人一区二区三区| 天堂影院成人在线观看| 男女下面进入的视频免费午夜| svipshipincom国产片| 99久久精品热视频| 色av中文字幕| 国产免费av片在线观看野外av| 人妻久久中文字幕网| 亚洲自拍偷在线| 国产精品美女特级片免费视频播放器 | 男女那种视频在线观看| 88av欧美| 亚洲欧洲精品一区二区精品久久久| 美女高潮喷水抽搐中文字幕| 国产精品一区二区免费欧美| 长腿黑丝高跟| 成人一区二区视频在线观看| 男人舔女人下体高潮全视频| 国产精品 欧美亚洲| 日本 av在线| 精品久久久久久久久久免费视频| 麻豆国产97在线/欧美 | 国产亚洲av高清不卡| 国产成人影院久久av| 国产高清视频在线播放一区| 久久这里只有精品19| 亚洲国产精品999在线| 成人18禁在线播放| 99精品欧美一区二区三区四区| 中文字幕高清在线视频| 香蕉丝袜av| 精品福利观看| 制服诱惑二区| 老司机靠b影院| 成人国语在线视频| 亚洲国产看品久久| 一本精品99久久精品77| 在线免费观看的www视频| 99国产精品99久久久久| 久久香蕉国产精品| 黑人操中国人逼视频| 天天添夜夜摸| 757午夜福利合集在线观看| 免费看十八禁软件| 久久久久久免费高清国产稀缺|