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

    長江中下游稻區(qū)不同水旱輪作模式和氮肥水平對稻田CH4排放的影響

    2019-08-08 08:06:24劉少文殷敏褚光徐春梅王丹英章秀福陳松
    中國農(nóng)業(yè)科學(xué) 2019年14期
    關(guān)鍵詞:產(chǎn)甲烷菌輪作通量

    劉少文,殷敏,褚光,徐春梅,王丹英,章秀福,陳松

    (中國水稻研究所/水稻生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室,杭州 311400)

    0 引言

    【研究意義】全球氣候變暖正日益影響著人類的生存環(huán)境和質(zhì)量[1],而其發(fā)生與溫室氣體濃度的增加有關(guān)[2]。CH4作為主要的溫室氣體之一,其增溫效應(yīng)約占全部溫室氣體總效應(yīng)的 15%—20%,且每年以1.1%的速度增長[3]。研究發(fā)現(xiàn)農(nóng)田是大氣CH4的重要排放源,約17.9%的CH4排放來源于稻田[4]。水稻是中國重要的糧食作物,其種植面積約占世界水稻種植總面積的 23%[5]。深入研究稻田 CH4排放特點(diǎn)及其影響因素,并在此基礎(chǔ)上制定相應(yīng)的減排措施,對緩解全球增溫效應(yīng)意義重大?!厩叭搜芯窟M(jìn)展】稻田CH4的產(chǎn)生是由產(chǎn)甲烷菌將土壤中有機(jī)物質(zhì)通過厭氧發(fā)酵轉(zhuǎn)化而成[6],CH4從厭氧沉積層或淹水土壤向大氣傳輸主要經(jīng)由分子擴(kuò)散、氣泡迸發(fā)和通過植株排放[7],其中超過90%的稻田CH4通過水稻植株排放到大氣[8]。對稻田CH4排放及相應(yīng)的減排措施已有大量研究。氣候條件、土壤特性、耕作方式以及田間管理是影響稻田CH4排放的重要因素[9-11]。一般認(rèn)為水旱輪作能夠改善土壤的理化性質(zhì),降低土壤還原物質(zhì)含量,減少 CH4的產(chǎn)生[12];而稻田持續(xù)淹水,會阻礙氧氣的擴(kuò)散,降低土壤氧化還原電位,易于產(chǎn)甲烷菌的快速繁殖,導(dǎo)致 CH4排放的增加[13];中期烤田則能減少36%—65%以上的CH4排放[14-15]。此外,研究表明,秸稈添加能顯著增加稻田 CH4排放量,且排放總量和秸稈還田量呈正相關(guān)[16]。在生產(chǎn)上,水旱輪作過程中往往伴隨著冬季作物秸稈全量還田,長江中下游地區(qū)主要的水旱輪作包括稻麥、稻油、稻綠、稻閑等模式[17-18],而近年來,秸稈覆蓋種植馬鈴薯也被廣泛推廣[19]。合理的秸稈還田能增加土壤有機(jī)質(zhì)含量,培肥地力、改善土壤結(jié)構(gòu)[20],但不同屬性的秸稈還田對溫室氣體排放的影響,目前尚未見系統(tǒng)報(bào)道。相關(guān)研究表明,田間氮肥施用對稻田 CH4排放亦有顯著影響,但目前尚未有一致的結(jié)論,有研究認(rèn)為氮肥會增加稻田 CH4排放[21-22],也有研究發(fā)現(xiàn)隨著施氮量的增加,稻田 CH4排放減少[23]。由于氮肥會影響還田秸稈的分解及其 CH4轉(zhuǎn)化過程[24],因此氮肥對 CH4排放的影響可能還與稻田的秸稈類型有關(guān)?!颈狙芯壳腥朦c(diǎn)】目前關(guān)于水旱輪作模式(不同屬性秸稈還田處理)與氮肥水平對稻田 CH4的排放特性的互作效應(yīng)研究報(bào)道較少,而對其相關(guān)排放機(jī)理的研究則尚未見報(bào)道?!緮M解決的關(guān)鍵問題】鑒此,本研究以多年長期定位試驗(yàn)(2003至今)為基礎(chǔ),選擇 4種不同水旱輪作模式(不同屬性秸稈還田處理),結(jié)合3種不同氮肥水平,通過田間氣體采集與土壤樣品分析,主要目的在于:(1)評估還田秸稈屬性與氮肥水平互作對稻季 CH4的排放規(guī)律及累積排放量的影響;(2)結(jié)合土壤甲烷代謝微生物、土壤理化性質(zhì)與秸稈屬性綜合分析,明確不同秸稈還田與氮肥處理下稻季CH4排放差異機(jī)理。以期為我國稻田溫室氣體排放以及合理減排措施的制定提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)區(qū)概況

    試驗(yàn)于2016—2017年,在浙江省富陽市中國水稻研究所農(nóng)場進(jìn)行(120.2°E,30.3°N,海拔11 m)。該農(nóng)場位于長江中下游平原,屬亞熱帶季風(fēng)氣候,年平均氣溫為13—20℃,其中7月份為35℃;年平均降水量為1200—1600 mm,其中4月份和9月份降雨量約占80%。本研究所需的稻田種植模式定位試驗(yàn)于2003年開始實(shí)施,于 2016—2017年測定稻田CH4排放,兩年土壤基本性質(zhì)年度間差異不顯著,具體指標(biāo)見表1。

    表1 不同輪作及氮肥水平下土壤基本性質(zhì)(兩年平均)Table1 Basic properties of soil under different rotation systems and different N rates(the two years’ average)

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

    本研究以不同水旱輪作模式長期定位試驗(yàn)為基礎(chǔ)(2003至今)。于2016—2017年測定稻田CH4排放。定位試驗(yàn)采用裂區(qū)隨機(jī)區(qū)組設(shè)計(jì),4種不同輪作模式為主區(qū),分別為水稻-休閑(RF),水稻-紫云英(RC-G),水稻-小麥(RW)和水稻-稻草覆蓋種植馬鈴薯(RP);水稻季氮肥水平為副區(qū),分別為N0(0)、N1(142.5 kg N·hm-2)和 N2(202.5 kg N·hm-2),小區(qū)位置和面積長期固定,每個小區(qū)長8 m、寬8 m,3次田間區(qū)組重復(fù)。

    冬季作物栽培:RF,田間不作任何處理,空閑放置,次年水稻移栽前,翻耕耘田;RC-G,紫云英于10月中下旬播種,出苗后施入尿素15 kg·hm-2作為基肥,次年盛花期(3月底4月初)將其割翻,翻耕還田;RW,小麥于10月下旬11月初撒直播,苗期和抽穗期分別施入復(fù)合肥75 kg·hm-2(基肥)和尿素15 kg·hm-2(追肥),次年5月中旬收獲,除籽粒外,其余秸稈翻耕還田;RP,稻草全程覆蓋種植馬鈴薯,于1月下旬在田面上直接擺放種薯,并覆蓋8—10 cm厚稻草,施入復(fù)合肥75 kg·hm-2作為基肥,5月中旬收馬鈴薯塊莖,其余部分翻耕還田,還田的有機(jī)物為腐爛的稻草和馬鈴薯莖葉等。

    圖1 不同輪作模式下秸稈還田量(A)和還田秸稈C/N比(B)Fig.1 Amount (A) and C/N ratio (B) of straw incorporated under different rotations

    單季稻種植:供試品種為嘉58(常規(guī)粳稻)。水稻田于6月上旬翻耕。水稻浸種催芽后,5月中上旬播于秧田,至3葉1心左右移栽,行株距26.7 cm×13.3 cm。水稻季施氮肥,氮肥總量根據(jù)氮肥處理設(shè)置,其中基肥∶分蘗肥∶穗肥按照 40%∶30%∶30%施入;磷肥(過磷酸鈣,60 kg·hm-2)作為底肥一次性施入;鉀肥(氯化鉀,135 kg hm-2)按基、蘗肥各50%分次施入。

    秸稈還田處理:水稻人工收獲,其秸稈全部移除用于冬季馬鈴薯覆蓋栽培。冬季作物秸稈全量還田,不同輪作模式下秸稈還田量及還田秸稈C/N比如圖1所示。不同輪作模式下秸稈還田量差異顯著(圖1-A,除2017年RW和RC-G差異不顯著),與RP相比,RW和RC-G的秸稈還田量分別減少45.7%和78.3%(2016),51.6%和62.5%(2017)。長期還田秸稈的C/N比(圖1-B)表現(xiàn)為RW顯著大于RP和RC-G(RP和RC-G差異不顯著);相較于RW,RP和RC-G輪作下C/N分別降低58.0%和71.3%(2016),58.9%和74.1%(2017)。因此,根據(jù)秸稈屬性,可以將 4種輪作模式劃分為RP(高秸稈量、低C/N比)、RW(高秸稈量、高C/N比)、RC-G(低秸稈量、低C/N比)和RF(對照-無秸稈)。

    1.3 樣品采集與分析

    CH4氣體采集于2016和2017年6—10月間開展,采用密閉式靜態(tài)(透明)箱-氣相色譜法測定田間CH4排放通量。氣體采樣箱由底座和采樣箱組成,靜態(tài)箱底橫截面積為0.5 m×0.6 m,采樣箱由PVC材質(zhì)制成,箱體高度隨水稻高度而增加(拔節(jié)前 0.5 m、拔節(jié)后1.2 m),中部固定風(fēng)扇,箱內(nèi)嵌有溫度計(jì),方便記錄溫度變化,底座為正方形U槽,嵌入地表,放置箱體后加水使之密封。

    采集頻率:水稻自移栽后第3天起,根據(jù)實(shí)際情況每5—10天采氣1次。采樣時間為8:00—11:00。氣體采集前打開箱內(nèi)頂部 12 V小風(fēng)扇以充分混勻箱內(nèi)氣體,采樣時,將采樣箱垂直安放在底座5 cm深的凹槽內(nèi)并加水密封,樣品采集分別在關(guān)箱后的0、6、12、24 min進(jìn)行。氣體樣中 CH4濃度由氣象色譜(Agilent 7890A,美國安捷倫科技有限公司)分析測定。CH4檢測器是 FID(氫火焰離子化檢測器),檢測溫度300℃,柱溫 60℃,載氣為 99.999%高純 N2,流速 30 mL·min-1。

    土壤及秸稈屬性測定:有機(jī)質(zhì)含量采用水合熱重鉻酸鉀氧化-比色法測定[67];全氮用全自動凱氏定氮儀(海能 K9840,濟(jì)南海能儀器股份有限公司)測定;堿解氮采用堿解擴(kuò)散法進(jìn)行測定[67];土壤pH用pH計(jì)(F2,梅特勒-托利多儀器(上海)有限公司)進(jìn)行測定。

    鑒于田間試驗(yàn)處理,年度間基本一致,且經(jīng)過多年定位試驗(yàn),土壤性質(zhì)基本穩(wěn)定,因此于2018年分蘗期測定土壤氧化還原電位與土壤m(xù)crA和pmoA基因豐度以探索CH4排放機(jī)理差異。

    土壤氧化還原電位(Eh)測定(2018年):于水稻移栽后每5—7天測定一次。利用土壤氧化還原電位儀(FJA-6型,南京傳滴儀器設(shè)備有限公司)測定根際土壤中的Eh值。

    唐玉煙不好意思地笑了笑,道:“實(shí)在抱歉,那天,我蒙蒙眬眬地感覺四周不停地顛簸搖晃,就像地震了一般,后來又感覺一陣窒息,我在窒息中醒來,一睜眼,就發(fā)現(xiàn)自己正被一只怪物攜著向前跑,我很害怕,于是拔出刀,割斷了它的脖子?!?/p>

    土壤m(xù)crA和pmoA基因豐度測定:分別于水稻移栽后 3 d(2018/6/28)、14 d(2018/7/9)、28 d(2018/7/23)利用取土器采集每個小區(qū)距稻樁20 cm,耕層(0—20 cm)土壤樣品(5點(diǎn)取樣法取樣),去雜質(zhì)后混合待用。土壤m(xù)crA和pmoA基因豐度采用實(shí)時熒光定量PCR分析。即用土壤DNA快速提取試劑盒提取土壤 DNA,用0.8%的瓊脂糖凝膠電泳檢測所提 DNA片段大小,并用超微量紫外分光光度計(jì)測定DNA濃度及質(zhì)量。設(shè)計(jì)產(chǎn)甲烷菌的mcrA基因、甲烷氧化菌的pmoA基因引物,進(jìn)行實(shí)時熒光定量PCR分析,反應(yīng)體系為20 μL[25]。各基因定量PCR分析所用的引物和反應(yīng)條件見表2,其中采集信號溫度為72℃。

    表2 光實(shí)時定量PCR擴(kuò)增引物及反應(yīng)條件Table2 Aplification primer and reaction condition of quantitative PCR

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

    溫室氣體排放通量由該氣體在箱中濃度隨時間的變化率計(jì)算得出,計(jì)算方程為氣體排放通量計(jì)算公式如下:

    式中,F(xiàn)為排放通量(mg·m-2·h-1);dc/dt為采樣過程中箱內(nèi)氣體濃度隨時間的變化率(mL·m-3·h-1));h為采樣箱體高度(m);ρ為標(biāo)準(zhǔn)狀態(tài)下氣體的密度(0.714 kg·m-3);T為采樣時箱內(nèi)溫度(℃),273為氣態(tài)方程常數(shù)。每次觀測時各處理氣體排放通量用3次重復(fù)的平均值表示;作物生長季節(jié)氣體平均排放通量的計(jì)算則以相鄰兩次采樣時間間隔為權(quán)重,將各次排放通量進(jìn)行加權(quán)平均;氣體累積排放量用平均排放通量和全生長季時間的乘積表示。

    本文采用Excel 2010和SAS 9.2軟件對數(shù)據(jù)結(jié)果進(jìn)行計(jì)算和分析,采用 Duncan法分析數(shù)據(jù)間的差異顯著性(P=0.05)。

    2 結(jié)果

    2.1 稻季CH4排放通量

    水稻生長期稻田CH4排放通量的變化受水分管理影響顯著,呈現(xiàn)明顯的階段性變化(圖2),可觀測到兩個排放高峰期—分蘗期和孕穗期CH4排放階段。其中分蘗期CH4排放階段(約移栽后7—30 d)持續(xù)時間占全生育期的19.0%—19.7%。分蘗中后期(移栽后 37 d左右),稻田逐漸進(jìn)入擱田期,各處理 CH4排放通量都表現(xiàn)不同程度下降。擱田復(fù)水后,CH4排放小幅上升,進(jìn)入孕穗期 CH4排放階段(約移栽后50—70 d),該階段持續(xù)時間約占生育期的 16.7%—17.2%。之后進(jìn)入低排放通量周期。

    輪作模式對稻田CH4排放的影響因氮肥水平差異而有所不同。N0水平下,分蘗階段CH4平均排放通量,RP、RW和RC-G無顯著差異,但均顯著高于RF,分別提高 72.1%—168.7%、70.8%—145.2%和 56.5%—145.2%(圖2- A 和2-B)。孕穗階段CH4排放通量,各模式間無顯著差異,但2016年RP顯著高于其他。此外,N0水平下分蘗期CH4排放表現(xiàn)為雙峰特點(diǎn)(約為移栽后第12—15 d和25—30 d)。各輪作在移栽后25—30 d出現(xiàn)最大排放峰,RP、RW、RC-G的最高排放通量(兩年平均)分別為55.0、42.9和38.6 mg·m-2·h-1。相比RF,RP、RW、RC-G的最高排放通量分別提高182.1%、120.0%和97.9%。

    隨著施氮量的增加,CH4排放通量對輪作的響應(yīng)出現(xiàn)分化。RF、RP和RW的分蘗期和孕穗期CH4平均排放通量隨著施氮量的增加而顯著提高。N2水平下,RP 分蘗期 CH4平均排放通量為 36.2 mg·m-2·h-1,要比N1和N0分別高14.8%和40.4%;與RP相似,RW則要高13.9%和43.9%。RF則表現(xiàn)為分蘗期CH4排放通量,N1和N2顯著高于N0,尤其在移栽后25—30 d,其平均排放通量達(dá)到23.8和27.5 mg·m-2·h-1(兩年平均),相比N0,分別提高29.9%和50.3%。而RC-G分蘗期CH4排放通量保持相對穩(wěn)定(2017)或略有下降(2016)。與 N0相比,增施氮肥下 RC-G分蘗期CH4平均排放通量降低20.0%—29.8%(N1)和11.9%—20.1%(N2)。

    圖2 2016(A、C和E)和2017(B、D和F)年度不同輪作模式及氮肥梯度下稻田CH4排放通量的季節(jié)變化Fig.2 Seasonal variations of CH4 fluxes in rice season under different rotation systems and N fertilizer rates in 2016 (A, C, E), 2017 (B, D, F)

    2.2 稻季CH4累積排放量

    由圖3中可以看出,不同輪作模式及氮肥水平對稻季CH4排放量的影響存在顯著互作。N0條件下,休閑(RF)模式的稻季CH4累積排放量要顯著低于冬季作物種植(RP、RW 和 RC-G)。RF稻季年均 CH4累積排放僅 18.4—21.2 g·m-2,分別比 RP、RW 和 RC-G降低42.6%—59.4%、40.7%—52.1%和36.0%— 52.1%。

    隨著氮肥的增加,RP、RW和RF處理下稻季CH4累積排放量均有不同程度的增加。N2水平下,RP、RW 和 RF的 CH4累積排放量分別為 51.2—55.8、45.3—51.5 和 25.0—30.5 g·m-2,比 N0水平提高23.0%—38.4%、26.7%—33.7%和 35.3%—43.5%;而與 N1相比,則提高 9.9%—19.7%、20.8%—23.1%和17.4%—18.8%(2016,RF模式下,N1和N2無顯著差異)。而RC-G稻季CH4累積排放量則隨著施氮量的增加而顯著降低。與N0相比,N1和N2下稻季CH4累積排放量分別降低20.7%—42.4%和10.6%—16.6%。不同輪作間比較,N1和N2水平下,RF和RC-G的稻季CH4累積排放量無顯著差異,但顯著低于RW和RP。

    圖3 2016(A)和2017年(B)不同輪作模式及肥料梯度下稻季CH4累積排放量Fig.3 CH4 accumulation emissions in rice season under different rotations and N fertilizer rates in 2016 (A) and 2017 (B)

    稻季CH4排放通量還與水稻生育階段相關(guān)。為了進(jìn)一步明確不同輪作模式與氮肥水平下稻季CH4的排放規(guī)律,將水稻全生育期分為分蘗期、孕穗期和灌漿期3個主要生育階段??傮w而言,稻季CH4的累積排放量分蘗期>孕穗期>灌漿期,其中分蘗期、孕穗期和灌漿期占比分別為51.85%—72.31%、22.9%—31.3%和6.0%—20.7%,這表明CH4排放主要集中在生育前中期,而后期占比較小(圖4)。輪作模式和氮肥水平對水稻不同生育階段CH4累積排放量影響顯著,差異主要集中體現(xiàn)在分蘗期和孕穗期,而灌漿期則無顯著差異。N0水平下,與RP、RW和RC-G相比,RF的分蘗期 CH4累積排放量降低 95.9%—171.1%、70.7%—144.1%和74.2%—165.8%;隨著氮肥的施用,各輪作分蘗期CH4累積排放量占比均顯著提高(2016年RC-G除外),增幅為7.0%—48.3%。孕穗期CH4累積排放量,RP和RW模式顯著高于RC-G和RF,而RC-G和RF無顯著差異。灌漿期,CH4排放量較小,各處理之間無顯著差異。

    2.3 產(chǎn)甲烷菌(mcrA)和甲烷氧化菌(pmoA)基因豐度

    土壤中產(chǎn)甲烷菌(mcrA)和甲烷氧化菌(pmoA)是影響稻季CH4排放的重要微生物。為進(jìn)一步明確不同輪作模式下稻季CH4排放微生物差異,于2018年水稻移栽后3、14和28 d采集并測定土壤中mcrA和pmoA基因豐度(圖5)。不同輪作模式與氮肥水平下,mcrA和pmoA基因豐度差異顯著。N0水平下,移栽后3 d,mcrA基因豐度量較少,而移栽14 d和28 d后豐度顯著增加。移栽后28 d mcrA基因豐度大小順序?yàn)?RP>RW>RC-G>RF,相較于RF,RP、RW 和 RC-G模式下mcrA基因豐度增加195.7%、121.3%和45.7%,這說明秸稈還田有利于提高土壤中產(chǎn)甲烷菌豐度。施氮后,各輪作模式mcrA基因豐度均顯著增加,相較于N0,N1水平下RP、RW、RC-G和RF的mcrA基因豐度分別增加191.4%、135.6%、43.1%和105.7%(移栽后28 d),其中以RC-G增加最少。

    圖4 2016(A、C和E)和2017(B、D和F)年度不同輪作模式及氮肥梯度下水稻不同生育階段CH4累積排放量Fig.4 CH4 accumulation emissions at different rice growth periods under different rotation systems and N fertilizer rates in 2016(A, C and E) and 2017 (B, D and F)

    與mcrA不同,N0水平下,移栽后14和28 d,RP和RW的pmoA基因豐度顯著高于RF和RC-G。增氮后,RF、RW和RP的pmoA基因豐度顯著減少降幅分別為42.7%、14.1%和13.3%(移栽后28 d),說明氮肥一定程度上抑制了甲烷氧化菌的活性,而RC-G則剛好相反,增施氮肥后pmoA基因豐度增加72.8%(移栽后28 d)。

    以上結(jié)果表明,隨著氮肥的增施,RP、RW和RF的產(chǎn)甲烷菌顯著提高,但RC-G增加較少;而甲烷氧化菌的豐度變化則與CH4排放剛好相反,表現(xiàn)為RP、RW和RF pmoA基因豐度減少,而RC-G增加。進(jìn)一步分析mcrA/pmoA比值(圖6-E和6-F),發(fā)現(xiàn)N0水平下,不同輪作模式下 mcrA/pmoA比值大小順序?yàn)镽P>RW、RC-G>RF(移栽后14 d和28 d),其變化規(guī)律與稻季不同輪作方式下CH4排放變化基本一致。增施氮肥后RF、RP和RW mcrA/pmoA基因豐度比值顯著提高,其增加比例分別為 71.4%—141.1%、197.1%—258.2%和84.6%—165.5%;而RC-G則下降26.8%—42.3%。這與田間CH4的排放規(guī)律相似。

    圖5 2018年不同輪作模式及氮肥梯度下產(chǎn)甲烷菌基因(mcrA;A和B)和甲烷氧化菌基因(pmoA;C和D)拷貝數(shù)和mcrA/pmoA比值(E和F)Fig.5 Abundance of methanogenic gene mcrA (A and B) and methanotrophic gene pmoA (C and D) and mcrA/pmoA ratio (E and F) under different rotation systems and N rates in 2018

    2.4 氧化還原電位差異性

    不同輪作模式下氧化還原電位(Eh)變化規(guī)律在不同氮肥處理下基本一致(圖6),擱田前(稻田一直處于淹水狀態(tài)),各輪作模式 Eh持續(xù)下降,在移栽后18—25 d降到最低,然后開始回升,移栽30—35 d以后Eh>0,并呈波浪式變化。N0水平下,水稻移栽前期(3—7 d),此時各輪作模式 Eh大小順序?yàn)镽P>RW>RC-G>RF(輪作間 Eh差異不顯著),可能是由于秸稈的添加與機(jī)械翻耕增加了土壤的通氣性。移栽后22 d,RP、RW和RC-G模式之間無顯著差異,但均顯著低于 RF模式,降幅為 118.4%、102.6%和80.5%,說明隨著淹水時間的延長,秸稈的降解消耗了土壤中的O2,加速了土壤中氧化還原電位的降低。增施氮肥后,Eh變化幅度因輪作模式而有所差異(移栽后22 d),N2水平下,RP、RW、RC-G和RF的Eh分別為-198.9、-188.5、-156.3和-80.3 mV(RP、RW 和 RC-G模式之間無顯著差異,但均顯著低于RF模式),與N0相比,分別降低33.1%、34.5%、29.7%和24.2%。而與N1相比,則降低7.2%、20.2%、38.4%和 15.9%。無機(jī)氮肥的投入,一定程度上促進(jìn)秸稈的腐解速率,并進(jìn)一步加速土壤中氧化還原電位的降低。

    3 討論

    3.1 輪作模式對稻田CH4排放的影響

    圖6 不同輪作模式及氮肥梯度下稻季土壤氧化還原電位變化規(guī)律Fig.6 Variations of soil redox potential in rice season under different rotation systems and N fertilizer rates

    輪作處理與氮肥水平并未改變稻季CH4總的排放規(guī)律。受田間水分管理和生長發(fā)育影響,稻季CH4排放主要集中在水稻分蘗期,其它時期(如拔節(jié)抽穗期、灌漿結(jié)實(shí)期)相對較少[28-29]。在水稻生長季,CH4排放表現(xiàn)為分蘗期前期升高,中后期特別是擱田開始后降低,擱田結(jié)束稻田復(fù)水后略有升高,并保持低排放量直至水稻收獲,呈現(xiàn)波浪式變化[5]。

    本研究發(fā)現(xiàn)輪作模式與氮肥水平對稻季CH4不同階段排放通量及其累積排放量均有顯著影響,其影響主要集中在分蘗期,這與張?jiān)婪糩30]和黃太慶[31]等的部分研究結(jié)論一致。另外,張嘯林等[32]研究認(rèn)為添加秸稈和綠肥作物均增加稻季CH4排放,而添加高生物量秸稈比添加低生物量綠肥作物有更多的CH4排放。這與本研究施氮(N1和 N2)處理下不同輪作模式 CH4排放規(guī)律基本一致;但在N0下,RP、RW和RC-G無顯著差異。其變化規(guī)律還與產(chǎn)甲烷菌基因豐度密切聯(lián)系,即添加秸稈生物量與CH4排放之間的關(guān)系受氮肥施用調(diào)控,并受到產(chǎn)甲烷菌活性的制約。微生物增殖與激活需要底物和能量,而土壤有效碳量是影響土壤微生物活性的限制因子[33]。稻田CH4排放主要受限于產(chǎn)甲烷菌和甲烷氧化菌的豐度及其活性。產(chǎn)甲烷菌是稻田土壤中甲烷產(chǎn)生的主要生物源[34]。而甲烷氧化菌,是以甲烷作為唯一碳源和能源進(jìn)行同化和異化代謝的微生物[35]。李大明等[36]研究指出產(chǎn)甲烷菌豐度與稻田CH4排放量呈正相關(guān),這與本研究結(jié)果一致。秸稈還田與氮肥處理影響了有機(jī)碳的轉(zhuǎn)化和產(chǎn)甲烷菌與甲烷氧化菌可利用的有機(jī)碳及能量,進(jìn)而影響到CH4排放[37]。在種植系統(tǒng)中,秸稈全量還田能夠顯著增加產(chǎn)甲烷菌豐度,進(jìn)而導(dǎo)致稻田CH4排放量提高[38]。本研究發(fā)現(xiàn)N0下,產(chǎn)甲烷菌豐度變化規(guī)律與秸稈還田量有關(guān),這與BAO[38]和ARIF[39]等關(guān)于大量秸稈碳的投入增加稻田產(chǎn)甲烷菌數(shù)量的結(jié)果一致。另一方面,稻田CH4排放量也因秸稈還田的有機(jī)物料種類而異。前人研究認(rèn)為添加的有機(jī)物料C/N越高,能源物質(zhì)越豐富,土壤有關(guān)微生物活性就越高,在環(huán)境適宜的條件下,CH4排放就越多[40]。本研究中RW、RP和RC-G還田秸稈的C/N比分別為70.39、29.25和19.26(圖1 B),其中RP處理C/N較低,但其CH4排放通量顯著高于RC-G,這與其秸稈還田量顯著高于RC-G相符(圖1-A)。顯然,以還田秸稈 C/N比作為判斷還田秸稈對 CH4排放的影響亦不全面。筆者認(rèn)為與C/N相比,還田秸稈中有機(jī)碳含量可能與 CH4排放關(guān)系更加密切,但其效應(yīng)受氮肥施用的影響。此外,長期定位輪作后,RP土壤中較高的有機(jī)質(zhì)含量可能也是導(dǎo)致 RP輪作模式 CH4排放量最大的重要因素(表1)。

    3.2 輪作模式和氮肥水平對稻田CH4排放的互作效應(yīng)

    氮肥能夠影響稻田CH4排放,其效應(yīng)受種類、施用量、施用方式和施用時間等影響[23]。長期還田秸稈的C/N比能干擾增施氮肥帶來的CH4排放效應(yīng)。在本研究的RP、RW和RF輪作模式中,我們發(fā)現(xiàn),增施氮肥能夠顯著增加稻田CH4排放(圖3),與前人研究印證[41]。但RC-G的表現(xiàn)則剛好相反。這與ZHU[42]和胡安永等[5]關(guān)于紫云英-水稻輪作模式不同施氮量下CH4排放研究結(jié)論一致。氮肥與有機(jī)物料的互作效應(yīng)對稻季CH4排放的影響主要與有機(jī)物料分解及其后續(xù)有機(jī)物種類有關(guān)[43]。一般而言,具有較高C/N比的有機(jī)底物,其分解速率隨著氮輸入量的增加而增加[44]。另外秸稈的分解速率還與相關(guān)微生物豐度密切相關(guān),后者則受到土壤中有效氮含量的制約[45],研究發(fā)現(xiàn),在氮受限的環(huán)境中,增加氮的輸入提高土壤中氮的有效性,這相應(yīng)提高了土壤中產(chǎn)甲烷菌所需有機(jī)底物的有效性,促進(jìn)了產(chǎn)甲烷微生物活性,增加CH4產(chǎn)生[46]。本研究進(jìn)一步發(fā)現(xiàn)微生物對氮肥的響應(yīng)機(jī)制因輪作模式而有所差異,增施氮肥促進(jìn)產(chǎn)甲烷菌的增殖,卻抑制了甲烷氧化菌的增殖(圖5-A,5-B),但其變化幅度因輪作處理而有所不同。這與許欣等[25]關(guān)于土壤中可利用氮含量(銨態(tài)氮與硝態(tài)氮)與mcrA基因豐度正相關(guān)而與 pmoA基因豐度呈顯著負(fù)相關(guān)的結(jié)論一致。筆者認(rèn)為CH4排放對氮肥的響應(yīng),原因可能是:(1)產(chǎn)甲烷菌主要以NH4+-N為氮源,外源無機(jī)氮的輸入提高了產(chǎn)甲烷菌繁殖所需底物的有效性[46],會加快產(chǎn)甲烷菌的生理活動[47],從而促進(jìn)其群體的快速增長[48];(2)甲烷單氧酶(MMO)是甲烷氧化菌的關(guān)鍵功能酶,可以在氧氣的作用下催化甲烷等低碳烷烴或烯烴羥基化或環(huán)氧化[49],NH4+是MMO的競爭性抑制劑,通過與MMO的競爭,減少了甲烷氧化,促進(jìn)了甲烷排放[50];(3)氮素增加還會引起植物根系分泌量增加,進(jìn)而為產(chǎn)甲烷菌帶來更多生活底物[51];(4)由于產(chǎn)甲烷菌代謝過程嚴(yán)格厭氧,只能生活在氧化還原電位很低的環(huán)境中[52]。增加無機(jī)氮的投入,能夠加速有機(jī)物料的腐解速率,腐解的同時消耗了水中的O2,降低了土壤中的氧化還原電位(圖6-B和6-C),進(jìn)一步增加了產(chǎn)甲烷菌的數(shù)量,進(jìn)而促進(jìn)CH4的排放。另外,WATANABE等[53]研究發(fā)現(xiàn),稻田CH4排放與秸稈還田量顯著正相關(guān),當(dāng)秸稈還田量為 2、4和 6 t·hm-2時,CH4排放量分別增加 19%、97%和 228%,其中秸稈中的碳對CH4排放的貢獻(xiàn)分別為19%、51%和60%。本研究中,RW和RP中秸稈還田量達(dá)到15.15和 7.85 t·hm-2,屬于高還田量栽培。微生物利用碳源進(jìn)行繁殖增長的同時,需要外界提供較高的氮源[46]。而RC-G處理本身帶入秸稈量較少(4.30 t·hm-2),含氮量較高。鑒此,筆者認(rèn)為稻田系統(tǒng)中秸稈還田C/N的相對含量可能是干擾氮肥水平對稻田CH4排放作用的關(guān)鍵。另外,本研究發(fā)現(xiàn),稻田土壤CH4排放量并不是取決于產(chǎn)甲烷菌或甲烷氧化菌中的某一個,而產(chǎn)甲烷菌與甲烷氧化菌基因豐度比值(mcrA/pmoA)能更客觀的反映CH4產(chǎn)生、氧化及排放之間的數(shù)量關(guān)系,這與LIU[52]等研究結(jié)果一致。當(dāng)土壤中生成的CH4通過擴(kuò)散進(jìn)入氧化區(qū)域時,大量CH4被甲烷氧化菌氧化。因此稻田土壤的CH4排放量不僅取決于CH4的生成能力,還決定于土壤對內(nèi)源CH4的氧化能力[54]。本研究表明,相較N0水平下,增施氮肥(N1)顯著增加RF、RW和RP輪作模式mcrA/pmoA比值;而RC-G輪作模式卻正好相反。其變化規(guī)律與CH4排放量基本一致。表明紫云英綠肥還田配施氮肥(N1),對土壤中甲烷氧化菌的促進(jìn)作用強(qiáng)于產(chǎn)甲烷菌,有助于減少稻田CH4排放。

    3.3 水稻植株及冬季作物種植對稻季CH4排放的影響

    稻田甲烷排放是稻田甲烷產(chǎn)生、氧化和傳輸綜合作用的結(jié)果,而水稻的生長強(qiáng)烈影響這3個過程[55]。水稻根系作為重要的中間載體,其根系形態(tài)和生理特征嚴(yán)重影響CH4排放[56],而投入的有機(jī)物料及其分解后有機(jī)物種類對根系的影響程度不同,最終會導(dǎo)致甲烷排放的差異。王紅妮等[57]發(fā)現(xiàn)水稻根系為了維持正常的呼吸作用,大量消耗根系附近的氧氣,降低了根系土壤的氧化還原電位,促進(jìn)了稻田土壤厭氧環(huán)境的形成,這也可能是導(dǎo)致不同輪作模式下稻田 Eh差異性的原因之一(圖6);同時,水稻根系通過分泌物和凋落物的形式為產(chǎn)甲烷菌提供底物,促進(jìn)產(chǎn)甲烷的增殖[58]。另外,水稻植株也會影響稻田CH4氧化作用,AULAKH等[59]研究指出大氣中的氧氣可以通過水稻植株的通氣組織進(jìn)入根系,形成一個好氧微區(qū)域,促進(jìn)甲烷氧化菌的生長,氧化稻田土壤中產(chǎn)生的甲烷;稻田土壤中產(chǎn)生的甲烷大約有36%[55]到80%甚至更多[60]在根區(qū)微氧區(qū)域被甲烷氧化細(xì)菌消耗。本研究中紫云英綠肥還田配施氮肥后CH4排放相對較低,是否與水稻植株根系發(fā)育程度及其通氣組織有關(guān)還有待更深入的研究。筆者認(rèn)為大量有機(jī)物料及無機(jī)氮肥的投入對水稻植株及其根系影響程度存在差異,進(jìn)而造成稻田甲烷產(chǎn)生、氧化和傳輸過程的不同,最終導(dǎo)致CH4排放存在差異。

    本研究發(fā)現(xiàn)系統(tǒng)中C/N相對含量與CH4排放密切聯(lián)系。冬季輪作及冬作后殘留氮會影響土壤中有效氮含量[61-62],同時也會影響稻季秸稈還田和氮肥施用,這也可能是導(dǎo)致溫室氣體差異的原因之一。土壤中有效氮源(碳源)是微生物繁殖的限制因素,而不同冬季輪作種植能夠影響土壤有機(jī)質(zhì)含量、培肥地力及改善土壤結(jié)構(gòu)[20],進(jìn)而改變了微生物代謝的有效有機(jī)碳/氮的組成與數(shù)量[63],從而顯著影響相關(guān)微生物的生存能力和活性[33]。本研究中,相較于 RF,RP、RW 和RC-G輪作系統(tǒng)下土壤全氮分別提高24.3%、14.0%和7.9%,有機(jī)碳分別提高24.2%、7.2%和6.3%(表1)。這與粱愛珍[64]和汪娟[65]等的部分結(jié)論一致。土壤全氮是能夠反映土壤供氮能力的重要指標(biāo),而有機(jī)碳在一定程度上影響著土壤結(jié)構(gòu)及其微生物多樣性,同時又是土壤肥力和基礎(chǔ)地力的最重要的物質(zhì)基礎(chǔ)[66]。盡管不同處理下,土壤屬性有所差異,但輪作研究中土壤屬性的變化緩慢而微量。而與秸稈還田和氮肥施用帶來迅速且高通量的CH4排放相比,土壤地力的變化是否是影響當(dāng)季CH4排放的重要因子,有待進(jìn)一步研究。

    4 結(jié)論

    N0處理不同輪作模式下稻田CH4排放差異顯著,具體表現(xiàn)為RP>RW>RC-G>RF;增施氮肥后,輪作模式與氮肥水平存在顯著互作,RW、RP和RF的CH4排放顯著增加,而RC-G則略有下降;產(chǎn)甲烷菌與甲烷氧化菌基因豐度比值(mcrA/pmoA)與CH4排放規(guī)律基本一致。鑒此,筆者認(rèn)為稻田系統(tǒng)中秸稈還田C/N的相對含量可能是干擾氮肥水平對稻田CH4排放作用的關(guān)鍵,當(dāng)系統(tǒng)中碳冗余時,相關(guān)微生物活性受到土壤中有效氮制約,投入無機(jī)氮可以減輕氮的限制作用從而顯著提高CH4排放;而碳不足時,繼續(xù)投入無機(jī)氮,相關(guān)微生物繁殖由于受到土壤中有限碳源的限制其活性也會受到抑制,CH4排放相對減少。

    猜你喜歡
    產(chǎn)甲烷菌輪作通量
    冬小麥田N2O通量研究
    農(nóng)民瓜菜“接茬輪作”增效益
    玉米淀粉水解液的制備及對油藏中產(chǎn)甲烷菌的激活
    智富時代(2019年5期)2019-07-05 18:43:32
    輪作休耕,讓土地歇一會兒
    耕地輪作效益調(diào)查報(bào)告
    基于GenBank中16S rRNA分析的產(chǎn)甲烷菌的地域分布特征
    農(nóng)藝措施對稻田土壤產(chǎn)甲烷菌的影響研究進(jìn)展
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    未培養(yǎng)技術(shù)在瘤胃產(chǎn)甲烷菌群研究中的應(yīng)用
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    国产免费一区二区三区四区乱码| 久久久久久久久久久久大奶| 亚洲av免费高清在线观看| 亚洲av男天堂| 大话2 男鬼变身卡| 久久狼人影院| 一区二区三区四区激情视频| 精品人妻偷拍中文字幕| 欧美日韩亚洲高清精品| 亚洲成av片中文字幕在线观看 | 国产精品一区二区在线观看99| 欧美亚洲日本最大视频资源| 青青草视频在线视频观看| 男女啪啪激烈高潮av片| 亚洲精品一二三| 啦啦啦在线观看免费高清www| 日日啪夜夜爽| 欧美日韩视频精品一区| 秋霞在线观看毛片| 多毛熟女@视频| 免费女性裸体啪啪无遮挡网站| 午夜激情av网站| 久久免费观看电影| 看免费成人av毛片| 亚洲精品美女久久av网站| 高清在线视频一区二区三区| 亚洲情色 制服丝袜| 欧美日韩视频精品一区| 人人妻人人澡人人爽人人夜夜| 少妇精品久久久久久久| 丝袜脚勾引网站| 亚洲人成电影观看| 亚洲在久久综合| 成年av动漫网址| 久久99热这里只频精品6学生| 女人久久www免费人成看片| 中文字幕精品免费在线观看视频| 只有这里有精品99| 国产伦理片在线播放av一区| 日本-黄色视频高清免费观看| 亚洲av成人精品一二三区| 91久久精品国产一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产亚洲欧美精品永久| 日本av免费视频播放| 伊人亚洲综合成人网| 人人妻人人澡人人爽人人夜夜| 少妇被粗大的猛进出69影院| 一区二区日韩欧美中文字幕| 国产不卡av网站在线观看| 水蜜桃什么品种好| 在现免费观看毛片| 热re99久久国产66热| 日日撸夜夜添| 一级,二级,三级黄色视频| 亚洲男人天堂网一区| 日韩精品免费视频一区二区三区| 免费日韩欧美在线观看| 国产精品 欧美亚洲| 国产日韩欧美亚洲二区| 国产男女内射视频| 最新的欧美精品一区二区| 中文欧美无线码| 久久免费观看电影| 亚洲人成77777在线视频| 亚洲国产精品成人久久小说| 亚洲成色77777| 久久青草综合色| 一个人免费看片子| 久久99精品国语久久久| 久久精品国产亚洲av涩爱| 精品人妻熟女毛片av久久网站| 亚洲成人av在线免费| 国产亚洲精品第一综合不卡| 美女国产高潮福利片在线看| 久久青草综合色| 女人精品久久久久毛片| 一级片免费观看大全| 欧美人与善性xxx| 国产日韩欧美在线精品| 秋霞在线观看毛片| 一本久久精品| 国产精品麻豆人妻色哟哟久久| 亚洲av日韩在线播放| 一边亲一边摸免费视频| 久久久久久免费高清国产稀缺| 欧美激情 高清一区二区三区| 精品人妻偷拍中文字幕| 波多野结衣一区麻豆| 五月伊人婷婷丁香| 中文欧美无线码| 国产国语露脸激情在线看| 永久网站在线| 亚洲av日韩在线播放| 亚洲精品aⅴ在线观看| 亚洲成色77777| av线在线观看网站| 不卡av一区二区三区| 亚洲成人手机| 久久久精品94久久精品| 黄色怎么调成土黄色| 女性被躁到高潮视频| 九九爱精品视频在线观看| 亚洲欧美色中文字幕在线| 国产又色又爽无遮挡免| 黄色一级大片看看| 中国三级夫妇交换| 国产精品久久久久久精品电影小说| 亚洲国产看品久久| 只有这里有精品99| 国产一区二区 视频在线| 最近的中文字幕免费完整| 春色校园在线视频观看| 亚洲精品国产一区二区精华液| 狠狠婷婷综合久久久久久88av| 亚洲精品美女久久av网站| 欧美老熟妇乱子伦牲交| 久久av网站| 18禁国产床啪视频网站| 日韩精品免费视频一区二区三区| 国产精品久久久久久精品电影小说| 亚洲国产欧美在线一区| 夫妻性生交免费视频一级片| 蜜桃国产av成人99| 国产一区二区 视频在线| 麻豆av在线久日| 中文欧美无线码| 色婷婷av一区二区三区视频| 黑人欧美特级aaaaaa片| 亚洲欧美中文字幕日韩二区| 97在线人人人人妻| 伦精品一区二区三区| 免费观看av网站的网址| 少妇人妻 视频| 亚洲人成电影观看| 九草在线视频观看| 国产乱来视频区| 午夜91福利影院| 中文字幕精品免费在线观看视频| 国产成人免费无遮挡视频| 日本爱情动作片www.在线观看| 电影成人av| 国产 一区精品| 制服丝袜香蕉在线| 天天操日日干夜夜撸| 美女大奶头黄色视频| 肉色欧美久久久久久久蜜桃| 人体艺术视频欧美日本| 在线观看免费日韩欧美大片| a级毛片在线看网站| 欧美日韩成人在线一区二区| av网站免费在线观看视频| 人成视频在线观看免费观看| 在线观看免费日韩欧美大片| 中文欧美无线码| 国产成人91sexporn| 中文字幕另类日韩欧美亚洲嫩草| 捣出白浆h1v1| 赤兔流量卡办理| 久久亚洲国产成人精品v| 亚洲第一青青草原| 热99久久久久精品小说推荐| 久久久久精品性色| 亚洲av电影在线观看一区二区三区| 18在线观看网站| 久久这里只有精品19| 一区二区三区四区激情视频| 婷婷成人精品国产| 在线精品无人区一区二区三| 免费高清在线观看视频在线观看| 久久久精品94久久精品| 亚洲欧美中文字幕日韩二区| 久久久久久久亚洲中文字幕| 国产成人一区二区在线| 精品人妻在线不人妻| 日韩精品免费视频一区二区三区| 久久精品aⅴ一区二区三区四区 | 天堂俺去俺来也www色官网| 午夜激情av网站| 一边摸一边做爽爽视频免费| 国产激情久久老熟女| 精品一区在线观看国产| 久久影院123| 国产1区2区3区精品| 国产精品久久久久久精品电影小说| 下体分泌物呈黄色| 亚洲精华国产精华液的使用体验| 寂寞人妻少妇视频99o| 观看av在线不卡| 欧美日韩一级在线毛片| 久久精品久久久久久噜噜老黄| 五月伊人婷婷丁香| 在线观看三级黄色| 爱豆传媒免费全集在线观看| av免费在线看不卡| 9色porny在线观看| 黄片播放在线免费| 亚洲一区中文字幕在线| 免费在线观看完整版高清| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产色片| 美女高潮到喷水免费观看| 一本色道久久久久久精品综合| 一级毛片电影观看| 9色porny在线观看| 性色avwww在线观看| 精品99又大又爽又粗少妇毛片| 永久网站在线| 91在线精品国自产拍蜜月| 乱人伦中国视频| 最近中文字幕2019免费版| 日韩一卡2卡3卡4卡2021年| 精品福利永久在线观看| 美女脱内裤让男人舔精品视频| 满18在线观看网站| 五月伊人婷婷丁香| 国产欧美日韩综合在线一区二区| 亚洲综合色网址| 国产亚洲午夜精品一区二区久久| 欧美成人午夜免费资源| 久久精品亚洲av国产电影网| av女优亚洲男人天堂| 国产成人一区二区在线| 欧美精品一区二区免费开放| 汤姆久久久久久久影院中文字幕| av在线观看视频网站免费| 中国三级夫妇交换| 久久韩国三级中文字幕| 久久久久久免费高清国产稀缺| 欧美精品高潮呻吟av久久| 寂寞人妻少妇视频99o| 久久久精品国产亚洲av高清涩受| 日日撸夜夜添| www.精华液| 热re99久久国产66热| 啦啦啦视频在线资源免费观看| 91精品伊人久久大香线蕉| 18禁动态无遮挡网站| 天天躁狠狠躁夜夜躁狠狠躁| 一区二区日韩欧美中文字幕| 欧美日本中文国产一区发布| 男女边摸边吃奶| 国产综合精华液| 国产又爽黄色视频| 色婷婷av一区二区三区视频| 午夜精品国产一区二区电影| 欧美日韩成人在线一区二区| 国产成人aa在线观看| 波野结衣二区三区在线| 午夜免费鲁丝| 宅男免费午夜| 免费久久久久久久精品成人欧美视频| 精品亚洲成a人片在线观看| 国产精品 国内视频| 久久午夜综合久久蜜桃| 一本大道久久a久久精品| 考比视频在线观看| 男人添女人高潮全过程视频| 80岁老熟妇乱子伦牲交| 你懂的网址亚洲精品在线观看| 久久免费观看电影| 国产精品人妻久久久影院| 久久午夜综合久久蜜桃| 性色avwww在线观看| 18禁观看日本| 亚洲成av片中文字幕在线观看 | 国产不卡av网站在线观看| 高清黄色对白视频在线免费看| 亚洲天堂av无毛| 久久影院123| av又黄又爽大尺度在线免费看| 男女国产视频网站| 久久久久久久久免费视频了| 91精品伊人久久大香线蕉| 麻豆精品久久久久久蜜桃| 一区福利在线观看| 女人久久www免费人成看片| 婷婷色麻豆天堂久久| 午夜av观看不卡| 一本—道久久a久久精品蜜桃钙片| 亚洲综合色惰| 97在线视频观看| 一区二区三区乱码不卡18| 99热全是精品| 精品亚洲乱码少妇综合久久| 十八禁网站网址无遮挡| 肉色欧美久久久久久久蜜桃| 成人国产麻豆网| 激情五月婷婷亚洲| 国产1区2区3区精品| 亚洲四区av| 日本91视频免费播放| 亚洲欧美一区二区三区黑人 | 欧美老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 国产日韩欧美亚洲二区| 精品国产乱码久久久久久男人| 亚洲综合色惰| 91在线精品国自产拍蜜月| 国产无遮挡羞羞视频在线观看| 国产午夜精品一二区理论片| 黄片小视频在线播放| 中文字幕人妻丝袜制服| 国产精品人妻久久久影院| 色视频在线一区二区三区| 中文乱码字字幕精品一区二区三区| 国产免费又黄又爽又色| a级片在线免费高清观看视频| 国产高清不卡午夜福利| 91午夜精品亚洲一区二区三区| 国产亚洲午夜精品一区二区久久| 亚洲色图综合在线观看| 女人被躁到高潮嗷嗷叫费观| 欧美日本中文国产一区发布| 亚洲av日韩在线播放| 中文精品一卡2卡3卡4更新| 国产男女内射视频| 最近中文字幕高清免费大全6| 日韩免费高清中文字幕av| 久久人人97超碰香蕉20202| 少妇的逼水好多| 国产黄色免费在线视频| 美女主播在线视频| 免费av中文字幕在线| 欧美人与性动交α欧美软件| 三上悠亚av全集在线观看| h视频一区二区三区| 91精品三级在线观看| freevideosex欧美| 成人国产av品久久久| 一级片'在线观看视频| 中文字幕制服av| 两性夫妻黄色片| 国产黄色视频一区二区在线观看| 人成视频在线观看免费观看| 超碰成人久久| 在线看a的网站| 午夜老司机福利剧场| 色94色欧美一区二区| 一区二区三区四区激情视频| 捣出白浆h1v1| 亚洲,一卡二卡三卡| 日韩电影二区| 婷婷色综合大香蕉| 春色校园在线视频观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲第一av免费看| 黄色一级大片看看| av免费在线看不卡| 午夜老司机福利剧场| 在线亚洲精品国产二区图片欧美| 如日韩欧美国产精品一区二区三区| 嫩草影院入口| 久久婷婷青草| 一级片免费观看大全| 日韩三级伦理在线观看| 不卡av一区二区三区| 午夜福利乱码中文字幕| 青春草亚洲视频在线观看| 亚洲伊人久久精品综合| 毛片一级片免费看久久久久| 男人添女人高潮全过程视频| 汤姆久久久久久久影院中文字幕| 777米奇影视久久| 热99国产精品久久久久久7| 国产精品二区激情视频| 这个男人来自地球电影免费观看 | 国产精品秋霞免费鲁丝片| 最新中文字幕久久久久| 午夜福利视频精品| 欧美亚洲日本最大视频资源| 欧美人与善性xxx| 香蕉国产在线看| 亚洲国产最新在线播放| 一区福利在线观看| 日韩av在线免费看完整版不卡| 亚洲情色 制服丝袜| 十八禁高潮呻吟视频| 寂寞人妻少妇视频99o| 最近最新中文字幕免费大全7| 夜夜骑夜夜射夜夜干| 国产午夜精品一二区理论片| 这个男人来自地球电影免费观看 | 日韩大片免费观看网站| 亚洲五月色婷婷综合| 亚洲男人天堂网一区| 免费在线观看黄色视频的| 国产 精品1| 久久久久精品人妻al黑| 午夜影院在线不卡| 2022亚洲国产成人精品| 国精品久久久久久国模美| 久久久久久久大尺度免费视频| 十八禁网站网址无遮挡| 日韩大片免费观看网站| 午夜激情av网站| 精品少妇黑人巨大在线播放| 波多野结衣一区麻豆| 久久国产亚洲av麻豆专区| 天天影视国产精品| 欧美精品人与动牲交sv欧美| 黄色怎么调成土黄色| 亚洲少妇的诱惑av| 中文字幕人妻丝袜制服| 久久久国产一区二区| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久久久久| 多毛熟女@视频| 涩涩av久久男人的天堂| 午夜福利网站1000一区二区三区| 女人精品久久久久毛片| 国产精品国产三级国产专区5o| 欧美变态另类bdsm刘玥| 久久精品亚洲av国产电影网| 亚洲综合精品二区| 秋霞伦理黄片| 婷婷色综合大香蕉| 国产精品熟女久久久久浪| 99久久综合免费| 少妇精品久久久久久久| 欧美精品人与动牲交sv欧美| 亚洲欧美一区二区三区久久| 日韩精品免费视频一区二区三区| 在线观看免费日韩欧美大片| 久久 成人 亚洲| h视频一区二区三区| 最近最新中文字幕免费大全7| 国产乱来视频区| 亚洲精品美女久久av网站| 极品人妻少妇av视频| 一级毛片我不卡| 亚洲三级黄色毛片| 91国产中文字幕| 亚洲av免费高清在线观看| 一区二区日韩欧美中文字幕| freevideosex欧美| 日韩成人av中文字幕在线观看| 女人高潮潮喷娇喘18禁视频| 午夜福利网站1000一区二区三区| 午夜免费鲁丝| 久久久欧美国产精品| 少妇被粗大猛烈的视频| 亚洲国产精品成人久久小说| 成人亚洲欧美一区二区av| 国产视频首页在线观看| tube8黄色片| 国产极品天堂在线| 好男人视频免费观看在线| 国产伦理片在线播放av一区| 国产精品 欧美亚洲| 国产欧美日韩综合在线一区二区| 在线 av 中文字幕| 久久鲁丝午夜福利片| 国产亚洲午夜精品一区二区久久| 久久精品国产a三级三级三级| 两性夫妻黄色片| av国产精品久久久久影院| 午夜福利乱码中文字幕| 99热国产这里只有精品6| 一级黄片播放器| 香蕉丝袜av| 自拍欧美九色日韩亚洲蝌蚪91| 天天躁夜夜躁狠狠躁躁| 麻豆av在线久日| 9191精品国产免费久久| 亚洲男人天堂网一区| 亚洲av福利一区| 亚洲美女搞黄在线观看| 一级片免费观看大全| 91精品国产国语对白视频| 18+在线观看网站| 中文字幕人妻丝袜制服| 高清不卡的av网站| 欧美精品高潮呻吟av久久| av免费观看日本| a 毛片基地| 久久久久久久久久人人人人人人| 欧美日韩亚洲国产一区二区在线观看 | 自线自在国产av| 日产精品乱码卡一卡2卡三| 国产免费一区二区三区四区乱码| 亚洲精品美女久久av网站| 十分钟在线观看高清视频www| 激情视频va一区二区三区| 国产成人精品无人区| xxxhd国产人妻xxx| 午夜免费观看性视频| av免费在线看不卡| 18禁观看日本| 国产 一区精品| 9热在线视频观看99| 久久女婷五月综合色啪小说| 搡老乐熟女国产| 精品人妻一区二区三区麻豆| 少妇被粗大猛烈的视频| 热99国产精品久久久久久7| 欧美人与善性xxx| 亚洲一级一片aⅴ在线观看| 欧美亚洲日本最大视频资源| 亚洲av日韩在线播放| 亚洲国产精品一区二区三区在线| 国产高清国产精品国产三级| h视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 狠狠精品人妻久久久久久综合| 欧美精品一区二区大全| 国产1区2区3区精品| 欧美人与善性xxx| 国产亚洲欧美精品永久| 午夜影院在线不卡| 一本—道久久a久久精品蜜桃钙片| 国产深夜福利视频在线观看| 日韩av不卡免费在线播放| 一级毛片我不卡| 高清视频免费观看一区二区| 爱豆传媒免费全集在线观看| 18禁观看日本| 久久99蜜桃精品久久| 国产精品久久久久久精品电影小说| 永久网站在线| 亚洲国产精品一区三区| 精品亚洲成a人片在线观看| 26uuu在线亚洲综合色| 人妻系列 视频| 大话2 男鬼变身卡| 99九九在线精品视频| 亚洲成色77777| 王馨瑶露胸无遮挡在线观看| 美女高潮到喷水免费观看| 老司机亚洲免费影院| 国产一区二区三区av在线| 免费黄网站久久成人精品| 搡女人真爽免费视频火全软件| 亚洲精品国产av蜜桃| 波多野结衣一区麻豆| 亚洲精品国产av成人精品| 色94色欧美一区二区| 看免费av毛片| kizo精华| 中国国产av一级| 成年女人在线观看亚洲视频| 又粗又硬又长又爽又黄的视频| 久久久久久久国产电影| 中文字幕人妻熟女乱码| 欧美亚洲日本最大视频资源| 日本午夜av视频| 国产精品国产三级国产专区5o| 久久久久久久久久久久大奶| a级毛片黄视频| 熟女电影av网| av卡一久久| 国产在线视频一区二区| 中文字幕av电影在线播放| 精品酒店卫生间| 男男h啪啪无遮挡| 免费大片黄手机在线观看| 午夜免费鲁丝| 国产成人精品福利久久| 欧美中文综合在线视频| 国产免费又黄又爽又色| 777米奇影视久久| 亚洲国产看品久久| 777米奇影视久久| 妹子高潮喷水视频| 日韩伦理黄色片| 精品99又大又爽又粗少妇毛片| 欧美日韩视频精品一区| 国产一级毛片在线| 免费高清在线观看视频在线观看| 777久久人妻少妇嫩草av网站| 成人毛片60女人毛片免费| 美女国产高潮福利片在线看| 一级片免费观看大全| 一二三四在线观看免费中文在| 日韩av不卡免费在线播放| 丰满少妇做爰视频| 日本免费在线观看一区| 亚洲国产精品成人久久小说| 2022亚洲国产成人精品| 婷婷成人精品国产| 美女脱内裤让男人舔精品视频| 精品卡一卡二卡四卡免费| 欧美最新免费一区二区三区| 精品少妇久久久久久888优播| 少妇的逼水好多| 国产精品免费大片| 国产又爽黄色视频| 日韩av免费高清视频| 国产又色又爽无遮挡免| 中国三级夫妇交换| 免费少妇av软件| www.精华液| 亚洲精品久久成人aⅴ小说| av网站在线播放免费| 国产不卡av网站在线观看| 欧美日韩视频高清一区二区三区二| 1024香蕉在线观看| 久久综合国产亚洲精品| 深夜精品福利| 欧美日韩视频高清一区二区三区二| 亚洲成av片中文字幕在线观看 | 国产亚洲最大av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 美女福利国产在线| 大陆偷拍与自拍| 咕卡用的链子| 久久99一区二区三区| 777久久人妻少妇嫩草av网站| 一级,二级,三级黄色视频| 免费看av在线观看网站| 国产精品一国产av| 精品少妇一区二区三区视频日本电影 | 大码成人一级视频| 亚洲国产看品久久| 日日撸夜夜添|