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

    不同秸稈還田方式對(duì)玉米農(nóng)田土壤CO2排放量和碳平衡的影響

    2023-08-15 01:22:10李金任立軍李曉宇畢潤(rùn)學(xué)金鑫鑫虞娜張玉玲鄒洪濤張玉龍
    關(guān)鍵詞:影響研究

    李金,任立軍,李曉宇,畢潤(rùn)學(xué),金鑫鑫,虞娜,張玉玲,鄒洪濤,張玉龍

    不同秸稈還田方式對(duì)玉米農(nóng)田土壤CO2排放量和碳平衡的影響

    李金,任立軍,李曉宇,畢潤(rùn)學(xué),金鑫鑫,虞娜,張玉玲,鄒洪濤,張玉龍

    沈陽農(nóng)業(yè)大學(xué)土地與環(huán)境學(xué)院/農(nóng)業(yè)農(nóng)村部東北耕地保育重點(diǎn)實(shí)驗(yàn)室/土肥高效利用國(guó)家工程研究中心,沈陽 110866

    【目的】探究秸稈還田方式對(duì)土壤CO2排放特征及碳平衡的影響,為東北地區(qū)農(nóng)田土壤固碳減排和秸稈還田方式的選擇提供科學(xué)依據(jù)?!痉椒ā坎捎锰镩g微區(qū)試驗(yàn),以玉米為供試作物,設(shè)置3種秸稈還田方式:秸稈淺層還田(QH)、秸稈深層還田(SH)和秸稈覆蓋還田(FG),無秸稈還田(CK)處理為對(duì)照。利用LI-8100A土壤碳通量自動(dòng)測(cè)定儀監(jiān)測(cè)玉米生長(zhǎng)季不同秸稈還田方式下土壤CO2的排放特征,探討土壤溫度、含水量、pH、微生物量碳及氮磷鉀速效養(yǎng)分和全量養(yǎng)分對(duì)土壤CO2排放的影響,并分析不同還田方式下的土壤碳平衡。【結(jié)果】在玉米生長(zhǎng)季,各處理土壤CO2排放速率均表現(xiàn)為先升高后降低的趨勢(shì)。土壤CO2累積排放量表現(xiàn)為FG>QH>SH>CK處理,相較于SH處理,F(xiàn)G和QH處理土壤CO2累積排放量分別增加了14.0%和6.4%,各處理間差異顯著(<0.05)。不同還田方式下土壤CO2排放速率與土壤溫度、土壤含水量進(jìn)行單因素模型擬合,均呈二次函數(shù)相關(guān)關(guān)系,且達(dá)到了顯著水平(<0.05),土壤溫度和土壤含水量分別解釋68.2%—73.7%和21.3%—82.8%的土壤CO2排放速率變化,但土壤溫度和土壤含水量的雙因素復(fù)合模型能更好地解釋土壤CO2排放速率的變化,解釋度達(dá)到78.5%—82.8%。相關(guān)性分析表明,土壤CO2累積排放量與速效鉀、微生物量碳呈極顯著相關(guān)關(guān)系(<0.01),與土壤有機(jī)質(zhì)、堿解氮、全氮和pH呈顯著的相關(guān)關(guān)系(<0.05)。秸稈還田處理下土壤碳平衡均為正值,為大氣CO2碳匯。SH處理下土壤碳平衡和固碳潛力顯著高于QH、FG處理,提高幅度分別為23.4%、475.7%和7.1%、30.7%(<0.05),表現(xiàn)出較強(qiáng)的碳匯功能。在兩年收獲期,秸稈還田顯著提高了玉米產(chǎn)量,其中SH處理最高,但與QH和FG處理間無顯著差異。【結(jié)論】本試驗(yàn)條件下,綜合考慮固碳減排效應(yīng)和產(chǎn)量,3種秸稈還田方式相比,秸稈深層還田(SH)是一種較好的還田方式。

    秸稈還田方式;土壤CO2排放;土壤溫度;土壤含水量;玉米產(chǎn)量;碳平衡

    0 引言

    【研究意義】溫室氣體排放導(dǎo)致的全球氣候變化和糧食安全已成為人類面臨的重大挑戰(zhàn)[1]。CO2作為一種主要溫室氣體成分,20%左右來源于土壤[2],農(nóng)田土壤是陸地生態(tài)系統(tǒng)向大氣排放CO2的重要源頭,其微小的波動(dòng)會(huì)對(duì)大氣中CO2濃度產(chǎn)生較大影響[3],并影響著全球碳循環(huán)和碳平衡過程。我國(guó)每年作物秸稈產(chǎn)量為7億—9億t[4],其中大部分秸稈被焚燒浪費(fèi),導(dǎo)致嚴(yán)重的環(huán)境污染并加速了全球溫室氣體的排放[5]。近年來,有研究學(xué)者發(fā)現(xiàn)秸稈是一種便于利用的自然資源,其還田后在提高土壤有機(jī)碳含量、改善土壤質(zhì)量和增加作物產(chǎn)量等方面有著重要的作用[6],但在秸稈還田的同時(shí)其對(duì)農(nóng)田土壤CO2排放的影響也是不容忽視的。因此,探究不同秸稈還田方式對(duì)土壤CO2排放特征及碳平衡的影響,對(duì)農(nóng)田土壤固碳減排,作物秸稈合理利用具有重大現(xiàn)實(shí)意義?!厩叭搜芯窟M(jìn)展】目前研究學(xué)者們關(guān)于秸稈還田與土壤CO2排放的關(guān)系仍存在分歧。YANG等[7]研究認(rèn)為,秸稈還田可以促進(jìn)有機(jī)物分解和養(yǎng)分循環(huán),但增加了土壤CO2排放。白銀萍等[8]研究表明,與不還田相比,秸稈覆蓋還田使土壤CO2排放增加了12.7%;而秸稈翻埋還田使土壤CO2排放降低了2.5%。WANG等[9]發(fā)現(xiàn),在冬小麥-春玉米輪作系統(tǒng)中,無論是翻耕還是免耕,秸稈覆蓋還田較不還田相比,土壤CO2排放分別降低了8.5%和14.5%。土壤碳平衡可以通過一定時(shí)間內(nèi)碳輸入和輸出的差異來計(jì)算[10],當(dāng)碳平衡值為正值時(shí),表現(xiàn)為大氣CO2的吸收“匯”,反之則為排放“源”[11]。DOSSOU-YOVO等[12]研究發(fā)現(xiàn),在貝寧北部地區(qū),秸稈覆蓋還田顯著降低了土壤碳平衡,但劉平奇等[13]研究表明,在我國(guó)東北黑土地區(qū),單獨(dú)進(jìn)行深松處理的土壤碳平衡表現(xiàn)為負(fù)值,而深松結(jié)合秸稈還田處理的土壤碳平衡卻表現(xiàn)為正值,表明秸稈還田能夠增加土壤碳平衡值。【本研究切入點(diǎn)】目前,相關(guān)學(xué)者對(duì)秸稈還田下土壤CO2排放的研究主要集中在覆蓋還田和表層淺混還田[14],而深還田下土壤CO2排放特征及其影響因素和土壤碳平衡的研究有待深入?!緮M解決的關(guān)鍵問題】本研究在田間試驗(yàn)條件下,探究秸稈還田方式對(duì)土壤CO2排放特征及碳平衡的影響,明確最佳的秸稈還田方式,以期為東北地區(qū)農(nóng)田土壤固碳減排和秸稈還田方式的選擇提供理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)地概況

    試驗(yàn)在沈陽農(nóng)業(yè)大學(xué)科研試驗(yàn)基地進(jìn)行(41°82′ N,123°56′ E),該區(qū)域?qū)儆跍貛О霛駶?rùn)大陸性氣候,試驗(yàn)地為雨養(yǎng)農(nóng)業(yè),年均溫度為9.00 ℃,年均降雨量為881.07 mm。供試土壤為棕壤,其基本理化性質(zhì)為有機(jī)質(zhì)含量12.21 g·kg-1、全氮0.67 g·kg-1、pH 6.70、土壤容重1.39 g·cm-3,土壤速效鉀、速效磷和堿解氮含量分別為132.52、62.19、64.45 mg·kg-1。

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

    該試驗(yàn)始于2020年5月,已連續(xù)兩年進(jìn)行秸稈還田試驗(yàn)。采用隨機(jī)區(qū)組設(shè)計(jì),設(shè)3種秸稈還田方式,分別為秸稈淺層還田(QH)、秸稈深層還田(SH)、秸稈覆蓋還田(FG),無秸稈還田(CK)為對(duì)照,共4個(gè)處理,每個(gè)處理設(shè)3次重復(fù),各試驗(yàn)微區(qū)面積為2.08 m2(1.6 m×1.3 m),每個(gè)微區(qū)四周用混凝土澆筑并做防水處理。2020年秋季玉米收獲后將秸稈粉碎至3 cm左右還入田中,各處理秸稈還田量相同,為10 500 kg·hm-2,具體操作見表1。

    表1 試驗(yàn)處理

    該區(qū)域種植模式為一年一熟制度,供試作物為玉米,品種為鄭單958。于2020年和2021年5月下旬播種,9月中旬收獲。每個(gè)小區(qū)種植2行玉米,每行4株,行距為0.4 m,株距為0.3 m。各處理施肥量相同,施用過磷酸鈣(P2O5含量為12%)75 kg·hm-2,硫酸鉀(K2O含量為50%)105 kg·hm-2,尿素(N含量為46.4%)240 kg·hm-2,作為基肥在播種前一次性施入,田間管理與當(dāng)?shù)爻R?guī)玉米栽培一致。

    1.3 試驗(yàn)方法與數(shù)據(jù)處理

    1.3.1 樣品采集與測(cè)定 采用LI-8100A土壤碳通量自動(dòng)測(cè)定儀(Li-Cor,Lincon,NE,USA)及其自帶呼吸室測(cè)定2021年玉米季土壤CO2排放速率。為平衡土壤中氣體,降低土壤結(jié)構(gòu)擾動(dòng)帶來的影響,在測(cè)定前2—3 d,先將PVC基座(直徑20 cm,高10 cm)均勻插入土壤中[15],露出土面2 cm左右。同時(shí)為減小試驗(yàn)測(cè)量誤差,盡量將PVC基座安放在各試驗(yàn)區(qū)同一位置。測(cè)定時(shí),將呼吸室放在PVC基座上,收集從土壤釋放的CO2,每個(gè)處理布置3次重復(fù),每個(gè)試驗(yàn)區(qū)重復(fù)測(cè)定3次,若3次測(cè)定值偏差大于15%,則需繼續(xù)進(jìn)行測(cè)定直至偏差小于15%。自施肥日起,開始測(cè)定土壤CO2排放速率,第一周每2 d測(cè)定一次,第二周每2—3 d測(cè)定一次,之后平均每7 d測(cè)定一次(具體根據(jù)天氣降雨等條件變化,測(cè)定時(shí)間會(huì)有微小調(diào)整),測(cè)定時(shí)間為上午8:00—12:00,在測(cè)定氣體的同時(shí),監(jiān)測(cè)土壤溫度和土壤含水量。試驗(yàn)期間,土壤CO2累積排放量()計(jì)算公式[16]為:

    式中,為土壤CO2累積排放量(kg·hm-2);為第次測(cè)定土壤CO2排放速率(μmol·m-2·s-1);(t-t)為連續(xù)兩次測(cè)定間隔時(shí)間(d);為測(cè)定的總次數(shù)。

    土壤呼吸碳累積排放量():

    =×0.27 (2)

    式中,為土壤呼吸碳累積排放量(kg·hm-2);為土壤CO2累積排放量(kg·hm-2);0.27為C占CO2分子量的比例。

    碳平衡()[10-11]:

    =-(3)

    =(+)×0.45 (4)

    =/2.1 (5)

    ×0.865 (6)

    =/(7)

    式中,為生態(tài)系統(tǒng)碳平衡(kg·hm-2);為凈初級(jí)生產(chǎn)力固碳量(kg·hm-2);為土壤微生物異養(yǎng)呼吸的碳釋放量(kg·hm-2);為地上生物量(kg·hm-2);為根生物量(kg·hm-2);0.45為作物地上部與根部含碳量;2.1為作物地上生物量與根生物量的比例;0.865為土壤微生物異養(yǎng)呼吸轉(zhuǎn)化系數(shù);為生態(tài)系統(tǒng)土壤固碳潛力。

    玉米收獲后,于2021年用五點(diǎn)法采集0—20 cm土層土樣,剔除根系和秸稈等雜物,混勻后風(fēng)干、過篩測(cè)定土壤理化性質(zhì)。土壤容重采用環(huán)刀法測(cè)定;土壤pH采用pH計(jì)測(cè)定(土水比=1﹕2.5);土壤有機(jī)質(zhì)和全氮采用元素分析儀(Vario EL Ⅲ,Elementar,Germany)測(cè)定。土壤微生物碳采用氯仿熏蒸法測(cè)定。土壤速效和全量氮磷鉀養(yǎng)分采用《土壤農(nóng)化分析》中方法測(cè)定[17]。土壤溫度采用長(zhǎng)度為10 cm的地溫計(jì)測(cè)定;10 cm土層體積含水量(全文簡(jiǎn)稱“含水量”)采用土壤水分測(cè)定儀(EC-5)測(cè)定;在2020年和2021年玉米收獲期,測(cè)量其穗行數(shù)、穗粒數(shù)、百粒重和籽粒產(chǎn)量,同時(shí)采集玉米植株地上部分,105 ℃殺青0.5 h后,在60 ℃下烘干至恒重,稱重。

    1.3.2 數(shù)據(jù)處理 采用Microsoft Excel 2019和SPSS 25.0進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析;采用單因素(one-way ANOVA)和Duncan法進(jìn)行方差分析和多重比較(α=0.05)。利用Origin 2017軟件繪圖并進(jìn)行單因素非線性曲線擬合和雙因素非線性曲面模型擬合。

    2 結(jié)果

    2.1 秸稈還田方式對(duì)土壤CO2排放速率動(dòng)態(tài)變化的影響

    由圖1可知,玉米生育期內(nèi)各處理土壤CO2排放速率變化趨勢(shì)基本一致,呈先升高后降低的波動(dòng)趨勢(shì)。各處理土壤CO2排放速率均在6月9日達(dá)到峰值,QH、FG、SH和CK處理土壤CO2排放速率分別為7.65、7.40、6.31、5.64 μmol·m-2·s-1。隨后土壤CO2排放速率逐漸下降,在8月2日出現(xiàn)降雨后,土壤CO2排放速率呈上升趨勢(shì),并在8月4日再次達(dá)到峰值,此時(shí)FG處理下土壤CO2排放速率最高,為6.70 μmol·m-2·s-1,其次是SH和QH處理,分別為5.84和5.54 μmol·m-2·s-1,CK處理最低,為4.50 μmol·m-2·s-1。

    圖1 不同秸稈還田方式下土壤CO2排放速率的動(dòng)態(tài)變化

    2.2 秸稈還田方式對(duì)土壤CO2累積排放量的影響

    由圖2可知,在玉米生長(zhǎng)季,不同處理間土壤CO2累積排放量均達(dá)到顯著差異。FG處理土壤CO2累積排放量最高,為23 897.17 kg·hm-2;其次是QH、SH處理,分別為22 298.22和20 962.72 kg·hm-2;CK處理土壤CO2累積排放量最少,為16 602.24 kg·hm-2。

    不同字母表示差異顯著(P<0.05)

    3種還田處理中,相較于SH處理,F(xiàn)G和QH處理土壤CO2累積排放量分別增加了14.0%和6.4%。由此可知,秸稈還田增加了土壤CO2排放,但秸稈深還相較于其他還田方式能夠減少土壤CO2排放。

    2.3 土壤CO2排放影響因素

    2.3.1 土壤溫度 由圖3可知,不同還田方式下土壤溫度的變化趨勢(shì)基本一致,呈先升高后降低的趨勢(shì)。在玉米生育期內(nèi),CK處理土壤溫度變化范圍最大,為17.3—33.3 ℃;其次為QH和FG處理,分別為16.0—31.0和15.0—30.0 ℃;SH處理的變化范圍最小,為15.8—30.5 ℃。各處理土壤溫度均在7月28日達(dá)到最大值,CK處理土壤溫度最高,為32 ℃,其次是QH、SH和FG處理,分別為30.6、30.2和29.9 ℃。

    圖3 不同秸稈還田方式下土壤溫度的動(dòng)態(tài)變化

    土壤溫度與土壤CO2排放速率關(guān)系采用非線性擬合的方法進(jìn)行擬合,結(jié)果如圖4所示,不同處理下土壤溫度與土壤CO2排放速率的擬合方程均達(dá)到了顯著水平。不同秸稈還田方式下,土壤溫度可以解釋68.2%—73.7%的土壤CO2排放速率變化,其中SH處理2最大,可以解釋73.7%土壤CO2排放速率變化。此外,CK處理下,當(dāng)土壤溫度<27.2 ℃時(shí),土壤CO2排放速率隨著土壤溫度的升高而增加,當(dāng)土壤溫度>27.2 ℃時(shí),表現(xiàn)出相反的變化趨勢(shì);其他3個(gè)秸稈還田處理表現(xiàn)出相似的變化規(guī)律,但土壤CO2排放速率出現(xiàn)的最高點(diǎn)不同,QH、FG和SH處理分別為24.7、23.9和24.8 ℃。

    2.3.2 土壤含水量 不同還田方式土壤含水量的動(dòng)態(tài)變化趨勢(shì)基本一致,呈現(xiàn)先升高后降低,隨后再升高達(dá)到峰值,之后下降的趨勢(shì)(圖5)。土壤含水量的變化受到了自然降雨的影響,整體隨著降雨的增加,土壤含水量增加,變化范圍在13.9%—39.6%。在大多數(shù)情況下,秸稈還田處理的土壤含水量相較于CK處理而言,均有所提高,說明秸稈還田后土壤的保水能力有所增加。土壤含水量平均值表現(xiàn)為FG>SH>QH>CK處理。

    *表示在0.05水平上差異顯著(n=60)。下同

    圖5 不同秸稈還田方式下土壤含水量的動(dòng)態(tài)變化

    對(duì)不同秸稈還田方式下土壤含水量與土壤CO2排放速率進(jìn)行非線性擬合發(fā)現(xiàn),土壤含水量解釋了21.3%—37.5%土壤CO2排放速率變化,低于土壤溫度對(duì)土壤CO2排放速率的解釋能力(圖6)。在土壤含水量較低時(shí),土壤CO2排放速率隨著土壤含水量的增加而增加,當(dāng)土壤含水量高于一定的數(shù)值時(shí),土壤CO2排放速率呈下降的趨勢(shì),但FG處理的轉(zhuǎn)折點(diǎn)要晚于其他處理。

    圖6 不同秸稈還田方式下土壤含水量與土壤CO2排放速率的模型擬合分析

    2.3.3 水熱因子的綜合作用 由圖7可知,在不同秸稈還田方式下,以土壤含水量()和土壤溫度()為自變量,土壤CO2排放速率()為因變量,建立雙因素復(fù)合模型(Y=z+aW+bT+cW+dT,<0.05),CK、QH、FG、SH處理的擬合方程見表2。各處理雙因素復(fù)合模型的2在0.785—0.828之間,相較于單因素(2為0.213—0.737)模型擬合較好,說明土壤含水量和土壤溫度的雙因素復(fù)合模型可以更好地解釋土壤CO2排放速率的變化。

    2.3.4 土壤理化性狀 將玉米收獲期0—20 cm土層理化性狀指標(biāo)與土壤CO2累積排放量進(jìn)行相關(guān)性分析,結(jié)果如圖8所示。土壤CO2累積排放量與速效鉀、微生物量碳呈極顯著正相關(guān)關(guān)系(<0.01),相關(guān)系數(shù)為0.83、0.74;與土壤有機(jī)質(zhì)、堿解氮和全氮呈顯著的正相關(guān)關(guān)系(<0.05),而與pH呈顯著的負(fù)相關(guān)關(guān)系,表明土壤CO2排放與土壤養(yǎng)分間有密切的聯(lián)系。

    2.4 秸稈還田方式對(duì)玉米生長(zhǎng)季碳平衡的影響

    由表3可知,SH處理的凈初級(jí)生產(chǎn)力固碳量最高,為6 775.72 kg·hm-2,比QH、FG和CK處理分別提高了0.7%、14.7%和24.0%。不同還田方式對(duì)土壤微生物異養(yǎng)呼吸的碳釋放量有顯著影響,整體表現(xiàn)為FG>QH>SH>CK處理。不同還田方式各處理的碳平衡值均為正值,說明該小區(qū)生態(tài)系統(tǒng)是大氣CO2的吸收“匯”,整體表現(xiàn)為SH>CK>QH>FG處理,SH與QH、FG處理相比,土壤碳平衡值分別提高了23.37%和475.71%,且均達(dá)到了顯著性差異。土壤的固碳潛力表現(xiàn)為CK處理最高,其次是SH處理,兩個(gè)處理間,差異沒有達(dá)到顯著水平。而SH處理土壤固碳潛力顯著高于QH和FG處理,分別提高了7.1%和30.7%,表明不同秸稈還田方式中SH處理具有較強(qiáng)的固碳潛力。

    表2 不同秸稈還田方式下水熱雙因素復(fù)合模型的擬合參數(shù)

    擬合方程為2+2;:土壤CO2排放速率;:土壤含水量;:土壤溫度

    The fitting equation is2+2;: Soil CO2emission rate;: Soil moisture content;: Soil temperature

    圖7 不同秸稈還田方式下土壤溫度和含水量對(duì)土壤CO2排放速率的影響

    表3 不同還田方式下玉米生長(zhǎng)季的碳平衡

    :地上生物量;:凈初級(jí)生產(chǎn)力固碳量;:土壤微生物異養(yǎng)呼吸的碳釋放量;:生態(tài)系統(tǒng)碳平衡;當(dāng)是正值時(shí),表明該小區(qū)能夠作為大氣CO2的吸收“匯”,反之則是排放“源”;:土壤固碳潛力。同列不同字母表示差異顯著(<0.05)。下同

    : Aboveground biomass;: Net primary productivity carbon sequestration;: Carbon emission by soil microbial heterotrophic respiration;: Ecosystem carbon balance; Whenwas positive, it indicated that the plot could be used as a "sink" of atmospheric CO2absorption; otherwise, it was the “source” of emissions;: Soil carbon sequestration potential. Different letters in the same column represent significant difference (<0.05). The same as below

    圖8 土壤CO2累積排放量與土壤理化性狀相關(guān)分析

    2.5 秸稈還田方式對(duì)玉米產(chǎn)量及其構(gòu)成因素的影響

    由表4可知,不同秸稈還田方式下玉米的穗粒數(shù)、百粒重和產(chǎn)量存在一定差異,但穗行數(shù)沒有達(dá)到差異顯著水平。2020年,各處理間玉米產(chǎn)量表現(xiàn)為SH>FG>QH>CK處理,在8 201.12—10 117.79 kg·hm-2間變化;2021年則表現(xiàn)為SH>QH>FG>CK處理,變化范圍為11 466.67—13 533.33 kg·hm-2。各處理玉米穗粒數(shù)與產(chǎn)量的變化趨勢(shì)相似。不同還田處理中SH處理顯著提高了玉米百粒重,與QH和FG處理相比,2年分別增加了4.63%、9.61%和4.11%、6.55%。

    表4 不同秸稈還田方式下的玉米產(chǎn)量及其構(gòu)成因素

    3 討論

    3.1 不同秸稈還田方式下土壤CO2的排放

    本研究中,在玉米生長(zhǎng)季前期QH處理與FG處理土壤CO2排放速率差異不大,但從6月中旬開始,F(xiàn)G處理土壤CO2排放速率明顯高于QH處理,可能是因?yàn)榍捌诮斩捙c土壤充分混合,使秸稈與土壤微生物直接接觸,刺激了微生物活性[18],到6月中旬,隨著易降解有機(jī)物質(zhì)的不斷消耗[19],QH處理下土壤CO2排放速率變小。在整個(gè)玉米生長(zhǎng)季,秸稈還田均增加了土壤CO2排放,可能有以下原因:(1)秸稈添加到土壤中,豐富了活性有機(jī)碳含量[20],進(jìn)而促進(jìn)土壤微生物代謝活動(dòng),增加了CO2排放,本研究中微生物量碳含量與CO2排放顯著的正相關(guān)關(guān)系也證實(shí)了這一點(diǎn)。(2)秸稈還田可以改善土壤肥力,促進(jìn)作物根系生長(zhǎng)和代謝[21],從而使來自根系的CO2排放量增加。(3)秸稈進(jìn)入土壤中加速了原有土壤有機(jī)質(zhì)的分解并引起正激發(fā)效應(yīng)[22],這都會(huì)導(dǎo)致土壤CO2排放增加。本研究結(jié)果表明,QH處理和SH處理下土壤CO2累積排放量顯著低于FG處理,這與閆翠萍等[23]研究發(fā)現(xiàn),秸稈還田至20 cm左右時(shí),土壤CO2累積排放量顯著高于覆蓋還田的結(jié)果存在差異。這可能與進(jìn)行秸稈還田前對(duì)土壤擾動(dòng)不同有關(guān),本文覆蓋還田和其他兩個(gè)還田方式均對(duì)不同土層土壤進(jìn)行翻動(dòng),對(duì)土壤擾動(dòng)一致,但其研究中只有秸稈還田至20 cm時(shí)對(duì)土壤進(jìn)行了擾動(dòng),而覆蓋還田為免耕播種,沒有擾動(dòng)土壤,導(dǎo)致郁閉在土壤中的氣體難以排放,故土壤CO2排放較低。同時(shí),本研究結(jié)果表明,不同秸稈還田方式中SH處理下土壤CO2累積排放量最低,一方面是因?yàn)楸韺油寥琅c底層土壤相比更易接觸氧氣,含有較高的C、N基質(zhì)供微生物活動(dòng),促進(jìn)了有機(jī)物的分解[24],另一方面,SH處理下土壤含水量明顯高于其余兩個(gè)處理,導(dǎo)致CO2在水中的溶解度變大,向地面擴(kuò)散的CO2減少[25]。此外,不同秸稈還田方式對(duì)土壤有機(jī)碳固定不同,有研究表明[26]秸稈深翻還田可以通過增加土壤有機(jī)碳固存來減少土壤CO2排放,于建光等[27]研究也發(fā)現(xiàn)土壤有機(jī)碳的穩(wěn)定性隨著土層深度的增加而增強(qiáng),深層土壤有機(jī)碳不易被微生物分解利用,所引起的激發(fā)效應(yīng)要小于表層土壤[28],因此,SH處理還田至20—40 cm土層的土壤呼吸較弱,CO2排放較低。

    3.2 影響土壤CO2排放的因素

    土壤溫度和土壤含水量是影響土壤CO2排放的重要環(huán)境因子。有研究指出[29],土壤溫度與土壤CO2排放速率呈指數(shù)關(guān)系,這與本文研究結(jié)果不一致。本研究中,土壤溫度與土壤CO2排放速率呈現(xiàn)較好的二次函數(shù)相關(guān)關(guān)系,當(dāng)土壤溫度達(dá)到一定值的時(shí)候,土壤CO2排放速率下降,并沒有隨著溫度的升高而增大,該結(jié)果與陳全勝等[30]研究發(fā)現(xiàn)土壤呼吸對(duì)土壤溫度的響應(yīng)存在一個(gè)適宜的范圍結(jié)論相一致,當(dāng)土壤溫度較高時(shí),一方面,土壤呼吸基質(zhì)數(shù)量和質(zhì)量發(fā)生變化,一些與呼吸作用相關(guān)的酶活性降低;另一方面,土壤中氮素的礦化速率增加,提高了可利用氮素含量,植物獲取氮素相對(duì)容易,使光合產(chǎn)物向地下根系分配的比例降低,影響了植物根系的代謝活動(dòng),進(jìn)而抑制了土壤呼吸。土壤含水量對(duì)土壤CO2排放的影響較為復(fù)雜,XU等[31]報(bào)道稱土壤含水量在20%左右有一個(gè)臨界值,超過這個(gè)臨界值,土壤含水量對(duì)土壤CO2排放的影響由正相關(guān)關(guān)系變?yōu)樨?fù)相關(guān)關(guān)系,與本文的研究結(jié)果一致,但也有研究表明[32]土壤CO2排放和土壤含水量之間沒有顯著的相關(guān)關(guān)系,原因可能是其研究區(qū)域土壤含水量變化范圍較窄,觀察到的值對(duì)土壤CO2排放的影響較弱[33];此外,土壤溫度是調(diào)節(jié)土壤含水量對(duì)CO2排放影響的重要因素,DING等[34]發(fā)現(xiàn)當(dāng)排除土壤溫度的混雜影響時(shí),極大地改善了土壤含水量與土壤CO2排放的關(guān)系,這是因?yàn)閮蓚€(gè)環(huán)境因素是同時(shí)變化的,土壤含水量對(duì)土壤CO2排放的影響被土壤溫度所掩蓋。有研究表明[35],當(dāng)土壤溫度或土壤含水量中的某一個(gè)因素處于較高或較低的范圍時(shí),另一個(gè)因素可能成為影響土壤CO2排放的主要因素。當(dāng)土壤含水量高于萎蔫系數(shù)時(shí),土壤CO2排放與土壤溫度密切相關(guān)[36];而當(dāng)土壤溫度較高時(shí),土壤CO2排放則隨著土壤含水量的變化而改變[37]。因此,單因素模型忽略了各因素間相互作用的影響,不能很好地描述土壤CO2排放。GOU等[38]認(rèn)為土壤溫度和土壤含水量的相互作用與土壤CO2排放之間具有更明顯的相關(guān)性,這一結(jié)果在本研究中也得到了證實(shí)。本研究發(fā)現(xiàn),土壤含水量和土壤溫度雙因素復(fù)合模型(2為0.785— 0.828)較單因素模型(2為0.213—0.737)相比,更好地解釋了土壤CO2排放速率變化。

    除了受土壤溫度和含水量的影響外,土壤CO2排放還受到土壤理化性狀的影響。葛高飛等[39]研究表明,微生物量碳與土壤CO2排放有顯著的正相關(guān)關(guān)系,相關(guān)系數(shù)為0.80,本文也得到了相似的研究結(jié)果,這是因?yàn)橥寥牢⑸锖粑菢?gòu)成土壤CO2排放的重要部分,秸稈還田為微生物提供了豐富的碳源和能源,增加土壤微生物數(shù)量的同時(shí)也增加了土壤微生物量碳,提高微生物新陳代謝強(qiáng)度,加速土壤養(yǎng)分礦化進(jìn)而釋放更多的CO2[40]。有研究通過主成分分析表明[41],土壤理化性質(zhì)中土壤有機(jī)質(zhì)和碳氮鉀含量是土壤CO2排放的主要影響因子,與本文研究結(jié)果一致,而土壤磷含量對(duì)土壤CO2排放沒有影響,可能是因?yàn)榱撞⒉皇窃搮^(qū)域農(nóng)田的限制營(yíng)養(yǎng)元素。本研究中,pH與土壤CO2排放呈顯著的負(fù)相關(guān)關(guān)系,與WANG等[42]研究結(jié)果不一致,主要?dú)w因于土壤本身pH不同,本試驗(yàn)土壤呈弱酸性,而WANG等試驗(yàn)土壤環(huán)境是堿性的,堿性環(huán)境能夠更好的溶解和分散土壤有機(jī)物質(zhì),使微生物充分接觸有機(jī)質(zhì)并進(jìn)行分解,增加了微生物呼吸[43],進(jìn)而增加土壤CO2排放。

    3.3 不同秸稈還田方式下玉米生長(zhǎng)季的碳平衡

    本研究中,各還田處理的碳平衡值為326.53— 1 879.87 kg·hm-2,說明進(jìn)行秸稈還田的各小區(qū)均是大氣CO2的吸收“匯”,這與在黃土高原地區(qū),秸稈還田后旱作春玉米田可作為一個(gè)碳匯的研究結(jié)果相符[44]。FG和QH較SH處理顯著降低了農(nóng)田碳匯效應(yīng),主要是因?yàn)镕G處理和QH處理兩種還田方式具有較高的土壤碳排放量,抵消了部分系統(tǒng)凈初級(jí)生產(chǎn)力固碳量,從而導(dǎo)致生態(tài)系統(tǒng)碳平衡值減小,弱化了農(nóng)田的碳匯作用。然而,也有一些研究表明[45],秸稈還田對(duì)農(nóng)田碳匯的影響微乎其微,甚至是負(fù)面的,可能是因?yàn)檫@些研究都是長(zhǎng)期試驗(yàn),長(zhǎng)期的秸稈還田使其土壤有機(jī)碳含量已經(jīng)達(dá)到了非常高的水平,進(jìn)而無法達(dá)到更高的碳平衡[46]。此外,LIU等[47]對(duì)已發(fā)表的176項(xiàng)研究進(jìn)行了meta分析,結(jié)果表明當(dāng)秸稈還田12年后,土壤碳接近飽和狀態(tài),農(nóng)田將達(dá)到新的碳平衡。由于本試驗(yàn)土壤有機(jī)質(zhì)相對(duì)較低,且秸稈還田只進(jìn)行了兩年,因此,我們認(rèn)為較低的土壤有機(jī)碳含量可能是該區(qū)域農(nóng)田表現(xiàn)為碳匯的主要原因,并且在有機(jī)碳含量未達(dá)到飽和之前,進(jìn)行秸稈還田作為農(nóng)田碳匯的效應(yīng)還將持續(xù)。WANG等[48]對(duì)中國(guó)農(nóng)田不同管理措施下土壤固碳潛力進(jìn)行預(yù)測(cè),發(fā)現(xiàn)秸稈還田措施具有最大的固碳潛力,可達(dá)34.4 Tg C·a-1。本研究發(fā)現(xiàn),SH處理土壤固碳潛力略低于CK處理,主要是因?yàn)榻斩掃€田向土壤中帶入了大量的有機(jī)質(zhì)[49],導(dǎo)致玉米生長(zhǎng)季內(nèi)有較高的土壤呼吸,增大了碳輸出量,進(jìn)而降低了土壤的固碳潛力。但劉平齊等[13]認(rèn)為,秸稈還田后,由于秸稈碳輸入導(dǎo)致的土壤呼吸增加,這并不能簡(jiǎn)單看作是土壤有機(jī)碳的損失或增加CO2排放的過程,此時(shí)土壤呼吸更多表示的是土壤動(dòng)物和微生物的活性強(qiáng)弱。若將本研究中不還田處理的秸稈全部進(jìn)行燃燒,根據(jù)玉米秸稈燃燒排放因子進(jìn)行計(jì)算[50],可產(chǎn)生19 833.33—21 262.50 kg·hm-2的CO2,加上不還田處理排放的CO2,其土壤CO2釋放總量達(dá)到36 435.57— 36 786.74 kg·hm-2,接近秸稈還田的2倍,因此綜合來看,SH處理的固碳潛力較優(yōu)。

    4 結(jié)論

    在玉米生長(zhǎng)季,秸稈還田處理顯著增加了土壤CO2排放,不同秸稈還田方式下土壤CO2累積排放量表現(xiàn)為秸稈覆蓋>秸稈淺層還田>秸稈深層還田>無秸稈還田處理。土壤溫度和土壤含水量雙因素復(fù)合模型可以解釋78.5%—82.8%的土壤CO2排放速率變化,土壤pH和土壤養(yǎng)分(碳氮鉀)含量對(duì)土壤CO2排放也存在顯著影響。不同秸稈還田方式下土壤碳平衡均為正值,其中秸稈深還處理有較高的土壤碳平衡值,達(dá)到1 879.87 kg·hm-2,表現(xiàn)出較強(qiáng)的碳“匯”功能,同時(shí)能夠增加土壤固碳潛力和提高作物產(chǎn)量。因此,本試驗(yàn)條件下,綜合考慮固碳減排效應(yīng)和產(chǎn)量,3種秸稈還田方式相比,秸稈深層還田(SH)是一種較好的還田方式。

    [1] PERKINS K M, MUNGUIA N, ELLENBECKER M, MOURE- ERASO R, VELAZQUEZ L. COVID-19 pandemic lessons to facilitate future engagement in the global climate crisis. Journal of Cleaner Production, 2021, 290: 125178.

    [2] LINQUIST B, VAN GROENIGEN K J, ADVIENTO-BORBE M A, PITTELKOW C, VAN KESSEL C. An agronomic assessment of greenhouse gas emissions from major cereal crops. Global Change Biology, 2012, 18(1): 194-209.

    [3] 武開闊, 張哲, 武志杰, 馮良山, 宮平, 白偉, 馮晨, 張麗莉. 不同秸稈還田量和氮肥配施對(duì)玉米田土壤CO2排放的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2022, 33(3): 664-670.

    WU K K, ZHANG Z, WU Z J, FENG L S, GONG P, BAI W, FENG C, ZHANG L L. Effects of different amounts of straw return and nitrogen fertilizer application on soil CO2emission from maize fields. Chinese Journal of Applied Ecology, 2022, 33(3): 664-670. (in Chinese)

    [4] 靳玉婷, 李先藩, 蔡影, 胡宏祥, 劉運(yùn)峰, 付思偉, 張博睿. 秸稈還田配施化肥對(duì)稻-油輪作土壤酶活性及微生物群落結(jié)構(gòu)的影響. 環(huán)境科學(xué), 2021, 42(8): 3985-3996.

    JIN Y T, LI X F, CAI Y, HU H X, LIU Y F, FU S W, ZHANG B R. Effects of straw returning with chemical fertilizer on soil enzyme activities and microbial community structure in rice-rape rotation. Environmental Science, 2021, 42(8): 3985-3996. (in Chinese)

    [5] ZHANG W S, LI H, LIANG L Y, WANG S X, LAKSHMANAN P, JIANG Z C, LIU C Y, YANG H, ZHOU M L, CHEN X P. An integrated straw-tillage management increases maize crop productivity, soil organic carbon, and net ecosystem carbon budget. Agriculture, Ecosystems & Environment, 2022, 340: 108175.

    [6] ZUO S S, WU D, DU Z L, XU C C, WU W L. Effects of white-rot fungal pretreatment of corn straw return on greenhouse gas emissions from the North China Plain soil. Science of the Total Environment, 2022, 807: 150837.

    [7] YANG X, MENG J, LAN Y, CHEN W F, YANG T X, YUAN J, LIU S N, HAN J. Effects of maize stover and its biochar on soil CO2emissions and labile organic carbon fractions in Northeast China. Agriculture, Ecosystems & Environment, 2017, 240: 24-31.

    [8] 白銀萍, 海江波, 楊剛, 黃晶, 董發(fā)勤, 劉明學(xué), 聶小琴. 稻田土壤呼吸及酶活性對(duì)不同秸稈還田方式的響應(yīng). 應(yīng)用與環(huán)境生物學(xué)報(bào), 2017, 23(1): 28-32.

    BAI Y P, HAI J B, YANG G, HUANG J, DONG F Q, LIU M X, NIE X Q. Effect of the straw returning pattern on soil respiration and enzyme activities. Chinese Journal of Applied and Environmental Biology, 2017, 23(1): 28-32. (in Chinese)

    [9] WANG H, WANG S L, YU Q, ZHANG Y J, WANG R, LI J, WANG X L. No tillage increases soil organic carbon storage and decreases carbon dioxide emission in the crop residue-returned farming system. Journal of Environmental Management, 2020, 261: 110261.

    [10] WU L F, LI B B, QIN Y, GREGORICH E. Soil CO2emission and carbon budget of a wheat/maize annual double-cropped system in response to tillage and residue management in the North China Plain. International Journal of Agricultural Sustainability, 2017, 15(3): 253-263.

    [11] 殷文, 史倩倩, 郭瑤, 馮福學(xué), 趙財(cái), 于愛忠, 柴強(qiáng). 秸稈還田、一膜兩年用及間作對(duì)農(nóng)田碳排放的短期效應(yīng). 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2016, 24(6): 716-724.

    YIN W, SHI Q Q, GUO Y, FENG F X, ZHAO C, YU A Z, CHAI Q. Short-term response of farmland carbon emission to straw return, two-year plastic film mulching and intercropping. Chinese Journal of Eco-Agriculture, 2016, 24(6): 716-724. (in Chinese)

    [12] DOSSOU-YOVO E R, BRüGGEMANN N, AMPOFO E, IGUE A M, JESSE N, HUAT J, AGBOSSOU E K. Combining no-tillage, rice straw mulch and nitrogen fertilizer application to increase the soil carbon balance of upland rice field in northern Benin. Soil and Tillage Research, 2016, 163: 152-159.

    [13] 劉平奇, 張夢(mèng)璇, 王立剛, 王迎春. 深松秸稈還田措施對(duì)東北黑土土壤呼吸及有機(jī)碳平衡的影響. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2020, 39(5): 1150-1160.

    LIU P Q, ZHANG M X, WANG L G, WANG Y C. Effects of subsoiling and straw return on soil respiration and soil organic carbon balance in black soil of northeast China. Journal of Agro-Environment Science, 2020, 39(5): 1150-1160. (in Chinese)

    [14] 馬玲, 王丹蕾, 韓昌東, 范麗娟, 金鑫鑫, 葉旭紅, 鄒洪濤. 秸稈還田方式對(duì)東北農(nóng)田土壤NH3揮發(fā)和N2O排放的影響. 環(huán)境科學(xué)研究, 2020, 33(10): 2351-2360.

    MA L, WANG D L, HAN C D, FAN L J, JIN X X, YE X H, ZOU H T. Effects of straw returning mode on NH3volatilization and N2O emission of farmland in northeast China. Research of Environmental Sciences, 2020, 33(10): 2351-2360. (in Chinese)

    [15] 張曉龍, 沈冰, 權(quán)全, 董樑, 田開迪. 渭河平原農(nóng)田冬小麥土壤呼吸及其影響因素. 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27(8): 2551-2560.

    ZHANG X L, SHEN B, QUAN Q, DONG L, TIAN K D. Soil respiration rates and its affecting factors in winter wheat land in the Weihe Plain, Northwest China. Chinese Journal of Applied Ecology, 2016, 27(8): 2551-2560. (in Chinese)

    [16] LI J H, LI H, ZHANG Q, SHAO H B, GAO C H, ZHANG X Z. Effects of fertilization and straw return methods on the soil carbon pool and CO2emission in a reclaimed mine spoil in Shanxi Province, China. Soil and Tillage Research, 2019, 195: 104361.

    [17] 鮑士旦. 土壤農(nóng)化分析. 3版. 北京: 中國(guó)農(nóng)業(yè)出版社, 2000.

    BAO S D. Soil and Agricultural Chemistry Analysis. 3rd ed. Beijing: China Agriculture Press, 2000. (in Chinese)

    [18] ZHAO X M, HE L, ZHANG Z D, WANG H B, ZHAO L P. Simulation of accumulation and mineralization (CO2release) of organic carbon in chernozem under different straw return ways after corn harvesting. Soil and Tillage Research, 2016, 156: 148-154.

    [19] TIAN P, SUI P X, LIAN H L, WANG Z Y, MENG G X, SUN Y, WANG Y Y, SU Y H, MA Z Q, QI H, JIANG Y. Maize straw returning approaches affected straw decomposition and soil carbon and nitrogen storage in northeast China. Agronomy, 2019, 9(12): 818.

    [20] XIA L L, WANG S W, YAN X Y. Effects of long-term straw incorporation on the net global warming potential and the net economic benefit in a rice–wheat cropping system in China. Agriculture, Ecosystems & Environment, 2014, 197: 118-127.

    [21] HU Z H, CUI H L, CHEN S T, SHEN S H, LI H M, YANG Y P, LI C Z. Soil respiration and N2O flux response to UV-B radiation and straw incorporation in a soybean–winter wheat rotation system. Water, Air, & Soil Pollution, 2013, 224(1): 1394.

    [22] WU L, ZHANG W J, WEI W J, HE Z L, KUZYAKOV Y, BOL R, HU R G. Soil organic matter priming and carbon balance after straw addition is regulated by long-term fertilization. Soil Biology and Biochemistry, 2019, 135: 383-391.

    [23] 閆翠萍, 張玉銘, 胡春勝, 董文旭, 王玉英, 李曉欣, 秦樹平. 不同耕作措施下小麥–玉米輪作農(nóng)田溫室氣體交換及其綜合增溫潛勢(shì). 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2016, 24(6): 704-715.

    YAN C P, ZHANG Y M, HU C S, DONG W X, WANG Y Y, LI X X, QIN S P. Greenhouse gas exchange and comprehensive global warming potential under different wheat-maize rotation patterns. Chinese Journal of Eco-Agriculture, 2016, 24(6): 704-715. (in Chinese)

    [24] HEITK?TTER J, HEINZE S, MARSCHNER B. Relevance of substrate quality and nutrients for microbial C-turnover in top- and subsoil of a Dystric Cambisol. Geoderma, 2017, 302: 89-99.

    [25] 朱曉晴, 安晶, 馬玲, 陳松嶺, 李嘉琦, 鄒洪濤, 張玉龍. 秸稈還田深度對(duì)土壤溫室氣體排放及玉米產(chǎn)量的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2020, 53(5): 977-989. doi: 10.3864/j.issn.0578-1752.2020.05.010.

    ZHU X Q, AN J, MA L, CHEN S L, LI J Q, ZOU H T, ZHANG Y L. Effects of different straw returning depths on soil greenhouse gas emission and maize yield. Scientia Agricultura Sinica, 2020, 53(5): 977-989. doi: 10.3864/j.issn.0578-1752.2020.05.010. (in Chinese)

    [26] 郭海斌, 張軍剛, 王文文, 馮曉曦, 許波, 王成業(yè), 白冬, 楊永乾. 耕作方式與秸稈還田對(duì)土壤呼吸及其影響因素的影響. 陜西農(nóng)業(yè)科學(xué), 2018, 64(6): 1-3.

    GUO H B, ZHANG J G, WANG W W, FENG X X, XU B, WANG C Y, BAI D, YANG Y Q. Effects of tillage methods and straw returning on soil respiration and its influencing factors. Shaanxi Journal of Agricultural Sciences, 2018, 64(6): 1-3. (in Chinese)

    [27] 于建光, 李輝信, 胡鋒, 沈英. 施用秸稈及接種蚯蚓對(duì)土壤顆粒有機(jī)碳及礦物結(jié)合有機(jī)碳的影響. 生態(tài)環(huán)境, 2006, 15(3): 606-610.

    YU J G, LI H X, HU F, SHEN Y. Effects of straw application and earthworm inoculation on soil particulate organic carbon and mineral-associated organic carbon. Ecology and Environment, 2006, 15(3): 606-610. (in Chinese)

    [28] ZHANG H, LIANG S, WANG Y H, LIU S W, SUN H D. Greenhouse gas emissions of rice straw return varies with return depth and soil type in paddy systems of Northeast China. Archives of Agronomy and Soil Science, 2021, 67(12): 1591-1602.

    [29] 王維鈺, 喬博, Kashif AKHTAR, 袁率, 任廣鑫, 馮永忠. 免耕條件下秸稈還田對(duì)冬小麥-夏玉米輪作系統(tǒng)土壤呼吸及土壤水熱狀況的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2016, 49(11): 2136-2152. doi: 10.3864/ j.issn.0578-1752.2016.11.010.

    WANG W Y, QIAO B, AKHTAR K, YUAN S, REN G X, FENG Y Z. Effects of straw returning to field on soil respiration and soil water heat in winter wheat-summer maize rotation system under no tillage. Scientia Agricultura Sinica, 2016, 49(11): 2136-2152. doi: 10.3864/ j.issn.0578-1752.2016.11.010. (in Chinese)

    [30] 陳全勝, 李凌浩, 韓興國(guó), 董云社, 王智平, 熊小剛, 閻志丹. 土壤呼吸對(duì)溫度升高的適應(yīng). 生態(tài)學(xué)報(bào), 2004, 24(11): 2649-2655.

    CHEN Q S, LI L H, HAN X G, DONG Y S, WANG Z P, XIONG X G, YAN Z D. Acclimatization of soil respiration to warming. Acta Ecologica Sinica, 2004, 24(11): 2649-2655. (in Chinese)

    [31] 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.

    [32] DONG W Y, LIU E K, WANG J B, YAN C R, LI J, ZHANG Y Q. Impact of tillage management on the short- and long-term soil carbon dioxide emissions in the dryland of Loess Plateau in China. Geoderma, 2017, 307: 38-45.

    [33] LI L J, YOU M Y, SHI H A, DING X L, QIAO Y F, HAN X Z. Soil CO2emissions from a cultivated Mollisol: effects of organic amendments, soil temperature, and moisture. European Journal of Soil Biology, 2013, 55: 83-90.

    [34] 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/4): 381-389.

    [35] LELLEI-KOVáCS E, KOVáCS-LáNG E, BOTTA-DUKáT Z, KALAPOS T, EMMETT B, BEIER C. Thresholds and interactive effects of soil moisture on the temperature response of soil respiration. European Journal of Soil Biology, 2011, 47(4): 247-255.

    [36] DILUSTRO J J, COLLINS B, DUNCAN L S, CRAWFORD C. Moisture and soil texture effects on soil CO2efflux components in southeastern mixed pine forests. Forest Ecology and Management, 2005, 204(1): 87-97.

    [37] MAESTRE F T, CORTINA J. Small-scale spatial variation in soil CO2efflux in a Mediterranean semiarid steppe. Applied Soil Ecology, 2003, 23(3): 199-209.

    [38] GOU Z W, YIN W, CHAI Q. Straw and residual film management enhances crop yield and weakens CO2emissions in wheat–maize intercropping system. Scientific Reports, 2021, 11: 14077.

    [39] 葛高飛, 梁永超. 玉米生長(zhǎng)過程中施肥對(duì)土壤呼吸和微生物量碳的影響. 中國(guó)農(nóng)學(xué)通報(bào), 2011, 27(18): 73-78.

    GE G F, LIANG Y C. Effect of fertilization on soil respiration and microbial biomass carbon during maize growth period. Chinese Agricultural Science Bulletin, 2011, 27(18): 73-78. (in Chinese)

    [40] 孫美佳, 周志勇, 王勇強(qiáng), 沈穎, 夏威. 有機(jī)物添加對(duì)山西太岳山油松林土壤呼吸及碳組分的影響. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2023, 47(1): 67-75.

    SUN M J, ZHOU Z Y, WANG Y Q, SHEN Y, XIA W. The effect of organic matter addition on soil respiration and carbon component in Pinus tabuliformis forests in Taiyue Mountain, Shanxi Province, China. Journal of Nanjing Forestry University (Natural Sciences Edition), 2023, 47(1): 67-75. (in Chinese)

    [41] 陳仕陽, 楊榮, 蘇永中, 杜澤玉. 河西綠洲農(nóng)田土壤呼吸特征及其對(duì)長(zhǎng)期施肥的響應(yīng). 中國(guó)沙漠, 2022, 42(3): 178-186.

    CHEN S Y, YANG R, SU Y Z, DU Z Y. Characteristics of soil respiration in farmland of Hexi oasis and its response to long-term fertilization. Journal of Desert Research, 2022, 42(3): 178-186. (in Chinese)

    [42] WANG J L, TENG D X, HE X M, QIN L, YANG X D, LV G H. Spatial non-stationarity effects of driving factors on soil respiration in an arid desert region. Catena, 2021, 207: 105617.

    [43] 張俊麗, Sikander Khan Tanveer, 溫曉霞, 陳月星, 高明博, 劉楊, 廖允成. 不同耕作方式下旱作玉米田土壤呼吸及其影響因素. 農(nóng)業(yè)工程學(xué)報(bào), 2012, 28(18): 192-199.

    ZHANG J L, TANVEER S K, WEN X X, CHEN Y X, GAO M B, LIU Y, LIAO Y C. Soil respiration and its affecting factors in dry-land maize field under different tillage systems. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(18): 192-199. (in Chinese)

    [44] GAO X, GU F X, HAO W P, MEI X R, LI H R, GONG D Z, MAO L L, ZHANG Z G. Carbon budget of a rainfed spring maize cropland with straw returning on the Loess Plateau, China. Science of the Total Environment, 2017, 586: 1193-1203.

    [45] WANG Y Y, HU C S, DONG W X, LI X X, ZHANG Y M, QIN S P, OENEMA O. Carbon budget of a winter-wheat and summer-maize rotation cropland in the North China Plain. Agriculture, Ecosystems & Environment, 2015, 206: 33-45.

    [46] BUYSSE P, ROISIN C, AUBINET M. Fifty years of contrasted residue management of an agricultural crop: impacts on the soil carbon budget and on soil heterotrophic respiration. Agriculture, Ecosystems & Environment, 2013, 167: 52-59.

    [47] LIU C, LU M, CUI J, LI B, FANG C M. Effects of straw carbon input on carbon dynamics in agricultural soils: a meta-analysis. Global Change Biology, 2014, 20(5): 1366-1381.

    [48] WANG X K, LU F, HAN B, OUYANG Z. Carbon sequestration by cropland soil in China: potential and feasibilty. IOP Conference Series: Earth and Environmental Science, 2009, 6(24): 242041.

    [49] 張莉, 王婧, 逄煥成, 張珺穜, 郭建軍, 董國(guó)豪. 短期秸稈顆粒還田對(duì)小麥-玉米系統(tǒng)作物產(chǎn)量與土壤呼吸的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2018, 29(2): 565-572.

    ZHANG L, WANG J, PANG H C, ZHANG J T, GUO J J, DONG G H. Effects of short-term granulated straw incorporation on grain yield and soil respiration in a winter wheat-summer maize cropping system. Chinese Journal of Applied Ecology, 2018, 29(2): 565-572. (in Chinese)

    [50] LI X H, WANG S X, DUAN L, HAO J M, LI C, CHEN Y S, YANG L. Particulate and trace gas emissions from open burning of wheat straw and corn stover in China. Environmental Science & Technology, 2007, 41(17): 6052-6058.

    Effects of Different Straw Returning Patterns on Soil CO2Emission and Carbon Balance in Maize Field

    LI Jin, REN LiJun, LI XiaoYu, BI RunXue, JIN XinXin, YU Na, ZHANG YuLing, ZOU HongTao, ZHANG YuLong

    College of Land and Environment, Shenyang Agricultural University/Key Laboratory of Arable Land Conservation (Northeast China), Ministry of Agriculture and Rural Affairs/National Engineering Research Center for Efficient Utilization of Soil and Fertilizer Resources, Shenyang 110866

    【Objective】The effects of different straw returning patterns on soil carbon dioxide (CO2) emission characteristics and carbon balance were discussed, which provided a scientific basis for carbon (C) sequestration and emission reduction as well as the selection of straw returning patterns in Northeast China. 【Method】A field micro-plot experiment were conducted with maize as the experimental crop, and three straw returning patterns were set up, including straw shallow returning (QH), straw deep returning (SH), and straw mulching (FG). No straw returning (CK) was used as the control treatment. The LI-8100A automatic soil C flux tester was used to monitor soil CO2emission characteristics under different straw returning patterns during the maize growth period. Effects of soil temperature, soil moisture content, pH, MBC, available nutrients and total nutrients of nitrogen, phosphorus and potassium on soil CO2emissions were analyzed, and soil carbon balance was investigated too. 【Result】During the maize season, soil CO2emission rates showed a trend of first increasing and then decreasing under different straw returning patterns. The cumulative soil CO2emissions were as follows: FG>QH>SH>CK treatment. Compared with SH treatment, the cumulative soil CO2emissions under FG and QH treatments increased by 14.0% and 6.4%, respectively. There was a significant difference between the treatments (<0.05). The single factor model fitting of soil CO2emission rates, soil temperature and soil moisture content under different straw returning patterns showed a quadratic function correlation, and reached a significant level (<0.05), soil temperature could explain the variation of soil CO2emission rate of 68.2%-73.7%, and soil moisture content could explain 21.3%-37.5%. However, the two-factor composite model of soil temperature and soil moisture content could better explain the variation of soil CO2emission rate, with an explanation of 78.5%-82.8%. Correlation analysis showed that cumulative CO2emissions were significantly correlated with available potassium and MBC (<0.01), and significantly correlated with soil organic matter, available nitrogen, total nitrogen, and pH (<0.05). The soil carbon balance was positive under different straw returning patterns, which were the "sink" of atmospheric carbon dioxide. The soil carbon balance and carbon sequestration potential under the SH treatment were significantly higher than the QH and FG treatments by increased of 23.4%, 475.7% and 7.1%, 30.7% (<0.05), respectively. Compared with other treatments, the SH treatment showed a strong carbon "sink" function. In the two-year harvest period, straw returning treatments significantly increased maize yield, SH treatment had the highest maize yield, but there was no significant difference with QH and FG treatments. 【Conclusion】Therefore, under the conditions of this experiment, taking into account the carbon sequestration and emission reduction effect and yield, SH was a better straw returning pattern compared with the three patterns.

    straw returning patterns; soil CO2emission; soil temperature; soil moisture content; maize yield; carbon balance

    10.3864/j.issn.0578-1752.2023.14.009

    2022-07-08;

    2022-11-14

    遼寧省教育廳重點(diǎn)項(xiàng)目(LSNZD202001)、遼寧省“興遼英才計(jì)劃”領(lǐng)軍人才項(xiàng)目(XLYC1905010)

    李金,Tel:15524365461;E-mail:lijin149917@163.com。通信作者鄒洪濤,E-mail:zht@syau.edu.cn

    (責(zé)任編輯 李云霞)

    猜你喜歡
    影響研究
    FMS與YBT相關(guān)性的實(shí)證研究
    是什么影響了滑動(dòng)摩擦力的大小
    2020年國(guó)內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    哪些顧慮影響擔(dān)當(dāng)?
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    新版C-NCAP側(cè)面碰撞假人損傷研究
    沒錯(cuò),痛經(jīng)有時(shí)也會(huì)影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    擴(kuò)鏈劑聯(lián)用對(duì)PETG擴(kuò)鏈反應(yīng)與流變性能的影響
    一二三四社区在线视频社区8| 日韩欧美一区二区三区在线观看| 国产精品av视频在线免费观看| 国产欧美日韩一区二区三| 亚洲中文日韩欧美视频| 婷婷六月久久综合丁香| 午夜精品久久久久久毛片777| 国产爱豆传媒在线观看 | 免费在线观看日本一区| 日韩高清综合在线| 黄色视频不卡| 亚洲精品一区av在线观看| 看黄色毛片网站| 少妇的丰满在线观看| 色精品久久人妻99蜜桃| 国产成年人精品一区二区| 亚洲精品一区av在线观看| 女同久久另类99精品国产91| 日韩欧美在线乱码| 成人18禁在线播放| 国产精华一区二区三区| 黄色视频不卡| 国产精品av视频在线免费观看| 国产乱人伦免费视频| 国产激情偷乱视频一区二区| 校园春色视频在线观看| 日韩有码中文字幕| 中国美女看黄片| 国产精品av久久久久免费| 美女高潮喷水抽搐中文字幕| 久久久国产欧美日韩av| 欧美三级亚洲精品| 日本熟妇午夜| 9191精品国产免费久久| www.熟女人妻精品国产| 国产精品亚洲av一区麻豆| 精品欧美国产一区二区三| 午夜福利免费观看在线| 精品国产乱码久久久久久男人| 午夜亚洲福利在线播放| 国产一区二区在线av高清观看| av欧美777| 亚洲精品在线观看二区| 12—13女人毛片做爰片一| 国产乱人伦免费视频| 午夜日韩欧美国产| 精品一区二区三区av网在线观看| 成人亚洲精品av一区二区| 午夜激情福利司机影院| 桃红色精品国产亚洲av| 欧美在线一区亚洲| 欧美性猛交╳xxx乱大交人| 九色成人免费人妻av| 草草在线视频免费看| 国产主播在线观看一区二区| 在线永久观看黄色视频| 国产av在哪里看| 久99久视频精品免费| av福利片在线| 成人国产综合亚洲| 18禁裸乳无遮挡免费网站照片| 97碰自拍视频| 精品一区二区三区视频在线观看免费| 黑人操中国人逼视频| 国产精品爽爽va在线观看网站| 男女之事视频高清在线观看| 三级男女做爰猛烈吃奶摸视频| 免费在线观看成人毛片| 91成年电影在线观看| 国产野战对白在线观看| 中文亚洲av片在线观看爽| 亚洲欧美精品综合久久99| 久久午夜亚洲精品久久| 99国产精品一区二区蜜桃av| 亚洲18禁久久av| 又粗又爽又猛毛片免费看| 欧美日韩瑟瑟在线播放| 婷婷丁香在线五月| 亚洲av日韩精品久久久久久密| 99国产精品99久久久久| 亚洲欧美一区二区三区黑人| 亚洲专区中文字幕在线| 国产成人av激情在线播放| 亚洲一码二码三码区别大吗| 女人高潮潮喷娇喘18禁视频| 免费一级毛片在线播放高清视频| 俺也久久电影网| 在线观看免费视频日本深夜| 夜夜夜夜夜久久久久| 淫秽高清视频在线观看| www.999成人在线观看| 无遮挡黄片免费观看| 天天一区二区日本电影三级| 欧美最黄视频在线播放免费| 亚洲人成网站高清观看| 亚洲avbb在线观看| 18禁黄网站禁片午夜丰满| 午夜免费激情av| 麻豆av在线久日| 两个人视频免费观看高清| 90打野战视频偷拍视频| 桃色一区二区三区在线观看| 日本精品一区二区三区蜜桃| 麻豆成人av在线观看| 观看免费一级毛片| 亚洲人成伊人成综合网2020| 亚洲av熟女| 久久精品aⅴ一区二区三区四区| 村上凉子中文字幕在线| 亚洲成人久久性| 999久久久精品免费观看国产| 女同久久另类99精品国产91| 老熟妇乱子伦视频在线观看| 色老头精品视频在线观看| 欧美日韩国产亚洲二区| 久久精品国产综合久久久| 在线国产一区二区在线| 亚洲一区二区三区不卡视频| 中出人妻视频一区二区| 亚洲精品美女久久av网站| 99久久久亚洲精品蜜臀av| 特大巨黑吊av在线直播| 一边摸一边抽搐一进一小说| 免费搜索国产男女视频| 国产高清有码在线观看视频 | a级毛片在线看网站| 国产精品爽爽va在线观看网站| 久久久水蜜桃国产精品网| 两个人免费观看高清视频| 一级a爱片免费观看的视频| 午夜视频精品福利| 激情在线观看视频在线高清| 日本一本二区三区精品| 欧美大码av| 免费看日本二区| 18美女黄网站色大片免费观看| 国产精品自产拍在线观看55亚洲| 亚洲熟妇中文字幕五十中出| 色综合站精品国产| 色综合站精品国产| 岛国在线免费视频观看| 中亚洲国语对白在线视频| 黄片大片在线免费观看| 国产成人av教育| 母亲3免费完整高清在线观看| 免费在线观看完整版高清| 日韩欧美精品v在线| 999久久久国产精品视频| 亚洲 国产 在线| 神马国产精品三级电影在线观看 | 免费电影在线观看免费观看| 久久精品夜夜夜夜夜久久蜜豆 | 黄片小视频在线播放| 亚洲精品美女久久av网站| 深夜精品福利| 国产亚洲欧美在线一区二区| 国产黄a三级三级三级人| АⅤ资源中文在线天堂| 国产精品爽爽va在线观看网站| 制服丝袜大香蕉在线| 国产成人系列免费观看| 一二三四在线观看免费中文在| 在线观看舔阴道视频| √禁漫天堂资源中文www| 中文字幕精品亚洲无线码一区| 久久午夜亚洲精品久久| 一个人免费在线观看电影 | 五月玫瑰六月丁香| 亚洲 欧美 日韩 在线 免费| 久久香蕉精品热| 久久精品人妻少妇| 美女免费视频网站| 丰满人妻一区二区三区视频av | 极品教师在线免费播放| 一个人免费在线观看的高清视频| 国产精品久久久av美女十八| 国产精品国产高清国产av| 欧美高清成人免费视频www| 日日摸夜夜添夜夜添小说| 最新美女视频免费是黄的| 欧美一级a爱片免费观看看 | 亚洲精品久久成人aⅴ小说| 18禁国产床啪视频网站| 久久久水蜜桃国产精品网| 亚洲九九香蕉| 我要搜黄色片| 国产在线精品亚洲第一网站| 久久 成人 亚洲| 精品久久久久久久人妻蜜臀av| 天堂av国产一区二区熟女人妻 | 老司机深夜福利视频在线观看| 在线观看www视频免费| 午夜视频精品福利| 欧美又色又爽又黄视频| 午夜视频精品福利| 国产一区在线观看成人免费| 国产黄a三级三级三级人| 国产成人精品久久二区二区91| www.www免费av| e午夜精品久久久久久久| 欧美久久黑人一区二区| 天堂av国产一区二区熟女人妻 | 免费一级毛片在线播放高清视频| 丁香六月欧美| 久热爱精品视频在线9| 欧美+亚洲+日韩+国产| av在线播放免费不卡| av欧美777| 变态另类丝袜制服| 99久久无色码亚洲精品果冻| 一本一本综合久久| 老司机福利观看| 欧美一级a爱片免费观看看 | 国产又色又爽无遮挡免费看| 国产精品,欧美在线| 天天添夜夜摸| 91国产中文字幕| 少妇的丰满在线观看| 国产探花在线观看一区二区| 老司机在亚洲福利影院| 一区二区三区国产精品乱码| 日韩大尺度精品在线看网址| 国产精品永久免费网站| 亚洲国产精品久久男人天堂| 欧美中文日本在线观看视频| 日韩有码中文字幕| 亚洲欧美激情综合另类| 好男人在线观看高清免费视频| av超薄肉色丝袜交足视频| 久久久精品国产亚洲av高清涩受| 国产成人aa在线观看| 亚洲一区高清亚洲精品| 91大片在线观看| 成人高潮视频无遮挡免费网站| 老司机靠b影院| 亚洲精品久久成人aⅴ小说| 99在线人妻在线中文字幕| 久久亚洲真实| 国内精品久久久久久久电影| 日韩大尺度精品在线看网址| www.精华液| 亚洲欧洲精品一区二区精品久久久| 国内精品久久久久久久电影| 国产视频内射| 精品第一国产精品| 欧美久久黑人一区二区| 又黄又粗又硬又大视频| 欧美成人午夜精品| 午夜福利18| 天堂动漫精品| 久久久精品欧美日韩精品| 18美女黄网站色大片免费观看| 99热6这里只有精品| 亚洲国产精品久久男人天堂| 国产激情偷乱视频一区二区| 国产精品一区二区三区四区免费观看 | 夜夜看夜夜爽夜夜摸| 夜夜看夜夜爽夜夜摸| 久久婷婷人人爽人人干人人爱| 在线观看美女被高潮喷水网站 | 天堂动漫精品| 亚洲精品中文字幕一二三四区| 国产成人精品久久二区二区91| 午夜福利视频1000在线观看| 欧美成人免费av一区二区三区| 91在线观看av| 精品一区二区三区视频在线观看免费| 国产一区在线观看成人免费| 极品教师在线免费播放| 嫁个100分男人电影在线观看| 亚洲第一电影网av| 日韩精品中文字幕看吧| 久久久国产欧美日韩av| 国产三级黄色录像| 成年人黄色毛片网站| 国产av一区在线观看免费| 亚洲va日本ⅴa欧美va伊人久久| www.自偷自拍.com| 99热只有精品国产| 色精品久久人妻99蜜桃| 亚洲avbb在线观看| 亚洲va日本ⅴa欧美va伊人久久| 欧美乱色亚洲激情| 亚洲午夜精品一区,二区,三区| 欧美乱码精品一区二区三区| 久久精品91蜜桃| 国产一区二区三区视频了| 午夜福利在线观看吧| 国产亚洲av嫩草精品影院| 麻豆国产av国片精品| 免费看日本二区| 色av中文字幕| 免费无遮挡裸体视频| 男人舔奶头视频| 听说在线观看完整版免费高清| 欧美日韩福利视频一区二区| 脱女人内裤的视频| 国产一区二区三区视频了| 美女高潮喷水抽搐中文字幕| 18美女黄网站色大片免费观看| 日韩欧美在线二视频| 欧美日韩瑟瑟在线播放| 五月伊人婷婷丁香| 亚洲欧美一区二区三区黑人| 可以免费在线观看a视频的电影网站| 天堂影院成人在线观看| 脱女人内裤的视频| 老司机在亚洲福利影院| 午夜免费成人在线视频| 久久香蕉国产精品| 精品欧美国产一区二区三| 最新在线观看一区二区三区| 久久热在线av| 午夜激情av网站| 亚洲自拍偷在线| 国产三级黄色录像| 国产精华一区二区三区| 国产激情久久老熟女| 欧美乱色亚洲激情| 国产主播在线观看一区二区| 日本一二三区视频观看| 日本一本二区三区精品| 日韩精品中文字幕看吧| 亚洲欧美日韩无卡精品| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| 欧美最黄视频在线播放免费| 一级黄色大片毛片| 一区福利在线观看| 午夜福利在线观看吧| 色综合欧美亚洲国产小说| 色老头精品视频在线观看| 欧美+亚洲+日韩+国产| 变态另类丝袜制服| 在线观看免费午夜福利视频| 国产1区2区3区精品| 少妇粗大呻吟视频| av有码第一页| aaaaa片日本免费| 国产成+人综合+亚洲专区| 丁香欧美五月| 日韩av在线大香蕉| 亚洲天堂国产精品一区在线| 黄色 视频免费看| 国产午夜精品久久久久久| 日韩大码丰满熟妇| 高清在线国产一区| 国产精品亚洲一级av第二区| 免费电影在线观看免费观看| 亚洲五月婷婷丁香| 一级毛片女人18水好多| 神马国产精品三级电影在线观看 | 日韩有码中文字幕| 色综合婷婷激情| 亚洲国产精品999在线| 亚洲五月天丁香| 精品第一国产精品| 亚洲全国av大片| 日本免费a在线| www.精华液| 国产午夜福利久久久久久| a级毛片a级免费在线| www国产在线视频色| 一二三四在线观看免费中文在| 岛国视频午夜一区免费看| 成人av在线播放网站| 香蕉久久夜色| 男人舔女人下体高潮全视频| 又大又爽又粗| 亚洲国产精品久久男人天堂| 日韩欧美在线乱码| 老汉色∧v一级毛片| 欧美黑人精品巨大| 一个人免费在线观看电影 | 18禁美女被吸乳视频| 最新在线观看一区二区三区| 国产成人av激情在线播放| 欧美日韩黄片免| 久久精品aⅴ一区二区三区四区| 大型av网站在线播放| 特大巨黑吊av在线直播| 777久久人妻少妇嫩草av网站| 欧美午夜高清在线| 久久久国产成人免费| 亚洲一区二区三区不卡视频| 草草在线视频免费看| 久久久久九九精品影院| 国产亚洲欧美98| 成人高潮视频无遮挡免费网站| 国产三级在线视频| 国产一区二区激情短视频| 少妇的丰满在线观看| 国产精品免费视频内射| 十八禁网站免费在线| x7x7x7水蜜桃| 色综合婷婷激情| 99国产极品粉嫩在线观看| 香蕉丝袜av| 欧美成人免费av一区二区三区| 变态另类成人亚洲欧美熟女| 久久国产乱子伦精品免费另类| 精品久久久久久久毛片微露脸| 麻豆成人午夜福利视频| 天堂av国产一区二区熟女人妻 | 免费av毛片视频| 亚洲欧美日韩高清专用| 亚洲午夜理论影院| 美女午夜性视频免费| a在线观看视频网站| 亚洲欧美日韩高清在线视频| 神马国产精品三级电影在线观看 | 91九色精品人成在线观看| 97人妻精品一区二区三区麻豆| 18禁美女被吸乳视频| 欧美另类亚洲清纯唯美| 亚洲色图 男人天堂 中文字幕| 人人妻人人看人人澡| 亚洲自拍偷在线| 中亚洲国语对白在线视频| 九色国产91popny在线| 国产三级中文精品| cao死你这个sao货| 99热这里只有是精品50| 午夜影院日韩av| 午夜福利在线观看吧| 一区二区三区国产精品乱码| 亚洲国产看品久久| 一本一本综合久久| 欧美乱色亚洲激情| 欧美午夜高清在线| 在线观看日韩欧美| 精品久久久久久久人妻蜜臀av| 国产亚洲精品第一综合不卡| 黄频高清免费视频| 制服人妻中文乱码| 波多野结衣高清无吗| 最近最新免费中文字幕在线| 国产激情偷乱视频一区二区| 日韩精品免费视频一区二区三区| 欧美色视频一区免费| 老司机深夜福利视频在线观看| 国产av一区在线观看免费| 欧美日韩福利视频一区二区| 日韩精品中文字幕看吧| 亚洲国产看品久久| 啪啪无遮挡十八禁网站| 人妻丰满熟妇av一区二区三区| 麻豆成人av在线观看| 国产午夜精品久久久久久| 亚洲人成电影免费在线| netflix在线观看网站| 99国产极品粉嫩在线观看| 久久中文字幕一级| 欧美一级毛片孕妇| 午夜激情av网站| 五月玫瑰六月丁香| 99精品久久久久人妻精品| 又黄又粗又硬又大视频| 亚洲精品色激情综合| 人妻久久中文字幕网| 好男人电影高清在线观看| 国模一区二区三区四区视频 | 人妻久久中文字幕网| 999久久久精品免费观看国产| 亚洲av日韩精品久久久久久密| 性欧美人与动物交配| 国产乱人伦免费视频| 亚洲熟妇中文字幕五十中出| 99久久国产精品久久久| www日本在线高清视频| 免费电影在线观看免费观看| 免费在线观看黄色视频的| 亚洲一区中文字幕在线| 波多野结衣高清作品| 精品无人区乱码1区二区| 不卡av一区二区三区| 18禁国产床啪视频网站| 国产精品美女特级片免费视频播放器 | 身体一侧抽搐| 国产精品精品国产色婷婷| 国产激情久久老熟女| 国产午夜福利久久久久久| 久久久久久大精品| 校园春色视频在线观看| 午夜激情av网站| 亚洲av成人av| 久久精品国产99精品国产亚洲性色| av免费在线观看网站| 19禁男女啪啪无遮挡网站| 91麻豆av在线| 老司机靠b影院| 一个人免费在线观看电影 | 亚洲狠狠婷婷综合久久图片| 欧美另类亚洲清纯唯美| 夜夜看夜夜爽夜夜摸| 国产午夜精品久久久久久| 99热6这里只有精品| 久久久久国产精品人妻aⅴ院| 一进一出好大好爽视频| 国产成人aa在线观看| 久久久久国产精品人妻aⅴ院| 国产黄a三级三级三级人| 国产69精品久久久久777片 | 人妻久久中文字幕网| 久久久久久人人人人人| 两个人视频免费观看高清| 国产1区2区3区精品| 久久这里只有精品中国| 欧美色视频一区免费| 日韩国内少妇激情av| 草草在线视频免费看| 操出白浆在线播放| 亚洲电影在线观看av| 欧美成人性av电影在线观看| 韩国av一区二区三区四区| 天堂动漫精品| 国产高清videossex| 亚洲美女黄片视频| 在线观看免费午夜福利视频| 国产精品亚洲一级av第二区| 欧美zozozo另类| 熟妇人妻久久中文字幕3abv| 90打野战视频偷拍视频| 国内精品久久久久精免费| 久久天堂一区二区三区四区| 成人亚洲精品av一区二区| 久久国产乱子伦精品免费另类| 在线观看免费午夜福利视频| 久久香蕉精品热| 久久亚洲真实| 亚洲最大成人中文| 黄色视频,在线免费观看| 免费一级毛片在线播放高清视频| 美女高潮喷水抽搐中文字幕| 法律面前人人平等表现在哪些方面| 日韩欧美精品v在线| 欧美日韩一级在线毛片| 国内精品一区二区在线观看| 一二三四社区在线视频社区8| 亚洲男人的天堂狠狠| 成人国产一区最新在线观看| 久久精品成人免费网站| 狂野欧美激情性xxxx| 在线观看66精品国产| videosex国产| 久久人妻av系列| 久久婷婷人人爽人人干人人爱| 久久这里只有精品19| 全区人妻精品视频| 中国美女看黄片| 成年版毛片免费区| а√天堂www在线а√下载| 香蕉久久夜色| 叶爱在线成人免费视频播放| 真人做人爱边吃奶动态| 国产黄a三级三级三级人| 欧美一级a爱片免费观看看 | 桃色一区二区三区在线观看| 日本五十路高清| 老司机午夜十八禁免费视频| 国产三级中文精品| 免费在线观看亚洲国产| 久久久久久九九精品二区国产 | 在线观看免费视频日本深夜| 久久亚洲真实| 国产精品98久久久久久宅男小说| 亚洲av第一区精品v没综合| 欧美精品啪啪一区二区三区| 日本a在线网址| 亚洲午夜精品一区,二区,三区| 亚洲aⅴ乱码一区二区在线播放 | 国产精品,欧美在线| 首页视频小说图片口味搜索| 亚洲成人久久性| 给我免费播放毛片高清在线观看| 久久国产精品影院| 69av精品久久久久久| 久久久久国产精品人妻aⅴ院| 一夜夜www| 亚洲自偷自拍图片 自拍| 一个人免费在线观看的高清视频| 精品一区二区三区av网在线观看| 欧美黑人欧美精品刺激| 欧美大码av| 国产精品免费一区二区三区在线| 天堂影院成人在线观看| 他把我摸到了高潮在线观看| 中文字幕最新亚洲高清| 国产精品一区二区三区四区久久| 脱女人内裤的视频| 亚洲精品美女久久av网站| 波多野结衣巨乳人妻| 九九热线精品视视频播放| 亚洲av中文字字幕乱码综合| 淫秽高清视频在线观看| 在线观看舔阴道视频| 国内毛片毛片毛片毛片毛片| 老司机深夜福利视频在线观看| 搡老妇女老女人老熟妇| www日本在线高清视频| 久久久久久久精品吃奶| 亚洲av五月六月丁香网| 一个人观看的视频www高清免费观看 | 两个人免费观看高清视频| 亚洲第一欧美日韩一区二区三区| 男人舔女人下体高潮全视频| 色av中文字幕| 色噜噜av男人的天堂激情| 精品久久久久久成人av| 两个人免费观看高清视频| 中文字幕最新亚洲高清| АⅤ资源中文在线天堂| 男男h啪啪无遮挡| 在线观看美女被高潮喷水网站 | 91麻豆av在线|