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

    秸稈還田年限及還田量對稻田凈溫室效應(yīng)的影響①

    2022-12-14 02:52:10朱小莉沈皖豫張廣斌
    土壤 2022年5期
    關(guān)鍵詞:田量溫室效應(yīng)通量

    黃 瓊,朱小莉,沈皖豫,樊 迪,張廣斌,馬 靜,徐 華*

    秸稈還田年限及還田量對稻田凈溫室效應(yīng)的影響①

    黃 瓊1,2,朱小莉1,2,沈皖豫1,2,樊 迪1,2,張廣斌1,馬 靜1,徐 華1*

    (1 土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點實驗室(中國科學(xué)院南京土壤研究所),南京 210008;2 中國科學(xué)院大學(xué),北京 100049)

    本文研究了不同秸稈還田年限(2 a和13 a)和還田量(0、1.6、3.2和4.8 t/hm2)對稻田CH4和N2O排放、土壤固碳速率和凈溫室效應(yīng)的影響。結(jié)果表明:稻田CH4排放量在第2年和第13年均隨秸稈還田量的增加呈線性增加,但第13年單位秸稈還田量對CH4排放的平均促進(jìn)效應(yīng)比第2年降低3%;各處理N2O排放量在第2年無顯著差異,第13年顯著促進(jìn)N2O排放55% ~ 171%,且隨秸稈還田量的增加呈線性增加;土壤固碳速率在第13年的變幅為0.40 ~ 0.50 t/(hm2·a),各處理間無顯著差異;凈溫室效應(yīng)在第13年隨秸稈還田量的增加呈線性增加,表現(xiàn)為源,這說明秸稈還田量對凈溫室效應(yīng)的影響受還田年限調(diào)控,其增加的溫室效應(yīng)逐步完全抵消了土壤固碳的減排效益,加劇了全球氣候變暖。因此,稻田在長期秸稈還田條件下是一個重要的碳源,亟待優(yōu)化秸稈還田技術(shù),進(jìn)一步固碳減排,以實現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展。

    還田年限;還田量;溫室氣體排放;土壤固碳;凈溫室效應(yīng)

    甲烷(CH4)和氧化亞氮(N2O)是大氣中僅次于二氧化碳(CO2)的重要溫室氣體,其百年尺度的全球增溫潛勢(global warming potential,GWP)分別為CO2的28倍和265倍[1]。2019年,大氣中CH4和N2O濃度分別達(dá)到1 877和332 ppb,是工業(yè)化前的260%和123%[2]。稻田是這兩種溫室氣體的重要排放源,據(jù)統(tǒng)計,中國稻田CH4和N2O排放量分別為5.3 Tg/a和N 31.1 Gg/a[3-4]。因此,加強稻田溫室氣體研究對減緩全球氣候變化有重要意義。

    秸稈是重要的可再生能源之一。我國的秸稈資源豐富,年產(chǎn)量約為8.3 × 108t[5]。秸稈還田是一項大力推行的保護性耕作措施,它能夠增加土壤大團聚體,改善土壤物理結(jié)構(gòu),增加作物產(chǎn)量,減少焚燒帶來的空氣污染,提高土壤有機質(zhì)含量[6-8]。我國秸稈肥料化利用率達(dá)47.3%[9],其中江蘇省秸稈還田率達(dá)64%,夏季麥秸還田率更高達(dá)85%[10]。秸稈還田一方面促進(jìn)稻田土壤有機碳固定,增加碳匯,緩解全球增溫效應(yīng)[11-12];另一方面顯著增加稻田甲烷排放,加劇全球增溫效應(yīng)[13-14]。所以,在發(fā)展循環(huán)農(nóng)業(yè)的過程中,綜合評估秸稈還田對稻田生態(tài)系統(tǒng)的凈溫室效應(yīng)十分必要。

    凈溫室效應(yīng)/凈全球增溫潛勢(net global warming potential,Net GWP)是給定生態(tài)系統(tǒng)中土壤固碳和溫室氣體排放的綜合評估指標(biāo),該指標(biāo)的正值和負(fù)值分別表示該生態(tài)系統(tǒng)溫室氣體的源和匯[15-16],可用于定量評估秸稈還田對稻田生態(tài)系統(tǒng)的全球增溫效應(yīng)的綜合影響。有關(guān)秸稈還田影響稻田溫室氣體排放和土壤固碳的研究報道眾多[17-21],而綜合評價其對稻田凈溫室效應(yīng)影響的報道較少,且研究結(jié)果不盡相同。湯宏等[22]報道秸稈還田1 a后顯著降低稻田凈溫室效應(yīng)295% ~ 1 128%;而張岳芳等[15]研究發(fā)現(xiàn)秸稈還田2 a后顯著增加凈溫室效應(yīng)57%;劉曉雨等[23]報道稻田的凈溫室效應(yīng)在秸稈還田21 a后顯著增加48%;Xia等[24]也報道秸稈還田20 ~ 21 a后顯著增加稻田凈溫室效應(yīng)25% ~ 52%。

    本研究以太湖流域稻麥輪作稻田為研究對象,以當(dāng)?shù)爻R?guī)模式為對照,通過田間原位觀測不同年限不同秸稈還田量稻田的CH4和N2O排放規(guī)律,估算土壤的固碳速率,綜合評估秸稈還田對稻田凈溫室效應(yīng)的影響,探明其隨不同年限變化的規(guī)律,分析其變化可能的原因,以期為實現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展及“雙碳”目標(biāo)提供參考數(shù)據(jù)。

    1 材料與方法

    1.1 試驗區(qū)概況

    試驗區(qū)位于江蘇省句容市白兔鎮(zhèn)行香村(31°58′N,119°18′E),該地區(qū)主要種植制度為稻–麥輪作,長期定位試驗開始于2006年。該地區(qū)屬于亞熱帶季風(fēng)氣候區(qū),2007年和2018年的年平均氣溫分別為16.9 ℃和16.7 ℃,年總降水量分別為1 162.5 mm和1 374.0 mm。供試土壤是發(fā)育于下蜀黃土的爽水性水稻土,耕層土壤(0 ~ 20 cm)初始基本理化性質(zhì)為:pH 5.9,容重1.2 g/cm3,有機碳含量8.2 g/kg,全氮含量1.1 g/kg。

    1.2 試驗設(shè)計

    本試驗選取2007年(2 a)和2018年(13 a)兩個處理年限為觀測期。共設(shè)置4個秸稈還田水平處理:S0(0 t/hm2)、S1(1.6 t/hm2)、S2(3.2 t/hm2)、S3(4.8 t/hm2),這4個秸稈還田量分別是當(dāng)?shù)爻R?guī)秸稈還田量的0、1/3、2/3和全量。各處理還田的秸稈切短至10 cm 左右,每年在稻麥季前施入,秸稈與農(nóng)田耕作層土壤均勻混合。每個處理設(shè)置3次重復(fù),隨機區(qū)組排列,各小區(qū)面積為3 m × 2.5 m。

    1.3 田間管理

    田間管理措施與當(dāng)?shù)爻R?guī)管理模式一致。水稻生長季尿素施用量為N 300 kg/hm2,施用比例為:基肥∶分蘗肥∶穗肥 = 2∶1∶1;磷肥(過磷酸鈣,Ca(H2PO4)2)和鉀肥(氯化鉀,KCl)作為基肥,一次性施入,施用量分別為450 kg/hm2和225 kg/hm2。田間水分管理方式為間歇灌溉(前期淹水–中期烤田–復(fù)水–后期干濕交替),觀測期內(nèi)田間管理具體時間見表1。

    表1 2007年和2018年水稻生長季作物管理和水分管理時間

    1.4 樣品采集與測定

    1.4.1 氣體樣品 采用手動密封箱法[25]采集稻田溫室氣體(CH4和N2O),水稻移栽后第2天開始采樣,3 ~ 7 d采樣1次,同時記錄相關(guān)環(huán)境因素數(shù)據(jù),包括土溫、氣溫、采樣箱內(nèi)溫度、水層深度和天氣情況等,采樣時間為上午9:00—11:30。采樣時,將采樣箱(箱體由透明有機玻璃制成,尺寸為0.5 m × 0.5 m × 0.5/1 m)放置于采樣底座(尺寸為0.5 m × 0.5 m × 0.15 m)凹槽(深度為4 cm)內(nèi),并用水密封,底座內(nèi)作物種植密度(6穴,每穴3株)與底座外保持一致。用雙通針將氣樣從采樣箱導(dǎo)入真空氣瓶待測,每次采樣時間間隔為10 min,共采樣4次。采集的氣樣用氣相色譜儀進(jìn)行測定(2007年:Shimadzu GC-14B,日本;2018年:Agilent 7890B,美國),其中CH4濃度用氫火焰離子化檢測器(FID)測定,N2O濃度用63Ni電子捕獲檢測器(ECD)測定,標(biāo)準(zhǔn)氣體由中國計量科學(xué)研究院提供。

    1.4.2 土壤樣品 土壤樣品在水稻收獲后使用內(nèi)徑2.5 cm、長20 cm的取土器進(jìn)行采集。在田間按“之”字形采樣,土壤充分混勻后揀去植物殘根和石礫等,經(jīng)風(fēng)干、研磨和過篩(100目),用碳氮元素分析儀進(jìn)行SOC含量測定(Elementar Vario MAX CN,德國)。環(huán)刀法測定土壤容重[26]。

    1.5 數(shù)據(jù)處理與統(tǒng)計分析

    1.5.1 稻田溫室氣體排放量計算 CH4和N2O的排放通量根據(jù)下式計算[25]:

    =×/×/× 273/(1)

    式中:是CH4或N2O排放通量,其中CH4排放通量單位為mg/(m2·h),N2O排放通量單位為μg/(m2·h)(以N計);是標(biāo)準(zhǔn)狀態(tài)下CH4或N2O的密度,其值是CH4為0.714 kg/m3,N2O 為1.25 kg/m3(以N計);是采樣箱內(nèi)有效體積,m3;是采樣箱所覆蓋的土壤面積,m2;/是單位時間內(nèi)采樣箱內(nèi)CH4或N2O濃度的變化,CH4單位為μl/(L·h),N2O單位為 nl/(L·h)(以N計);是采樣過程中采樣箱內(nèi)的溫度,K。

    CH4和N2O的排放量根據(jù)下式計算:

    E= ∑[(F+1+F)/2] × (D+1–D) × 24/1 000 (2)

    式中:E是CH4或N2O排放量,CH4排放量單位為g/m2,N2O排放量單位為mg/m2(以N計);FF+1分別是第和+1次采樣時CH4或N2O排放通量,CH4排放通量單位為mg/(m2·h),N2O排放通量單位為μg/(m2·h)(以N計);DD+1分別是第和+1次的采樣時間,d。

    1.5.2 土壤固碳速率計算 基于有機碳儲量的增量計算每年有機碳固定速率(SOCsequestrate rate)。

    有機碳儲量(SOCstock)計算公式如下[27]:

    SOCstock= SOCcontent××× 100 (3)

    式中:SOCstock是土壤有機碳儲存量,t/hm2(以C計);SOCcontent是土壤有機碳含量,g/kg(以C計);是土壤密度,g/cm3;是土壤樣品取樣深度(20 cm)。

    有機碳固定速率(SOCsequestrate rate)計算公式如下:

    SOCsequestrate rate= (SOCstock n–SOCstock i) /(4)

    式中:SOCsequestrate rate是土壤有機碳固定速率,t/(hm2·a) (以C計);SOCstock n和SOCstock i分別是觀測年和初始年的土壤有機碳儲存量值,t/hm2(以C計);是觀測年處理年限值,a。

    1.5.3 稻田凈溫室效應(yīng)計算 凈溫室效應(yīng)(Net GWP)是稻田排放CH4和N2O的綜合增溫潛勢(GWP)與土壤固碳減緩全球變暖貢獻(xiàn)的差值,計算公式如下[28]:

    GWP = CH4× 28 + N2O × 265 (5)

    Net GWP = GWP – SOCsequestrate rate× 44/12 (6)

    式中:GWP是稻田排放CH4和N2O的綜合增溫潛勢,t/hm2(以CO2-eq計);CH4是CH4排放量,t/hm2;N2O是N2O排放量,t/hm2;28和265分別是百年尺度上CH4和N2O的GWP值;Net GWP是凈溫室效應(yīng),t/(hm2·a) (以CO2-eq計);44/12是將土壤有機碳換算為CO2當(dāng)量的系數(shù)。

    1.5.4 統(tǒng)計分析 試驗結(jié)果均以每個處理3次重復(fù)測量值的平均值與標(biāo)準(zhǔn)差來表示,數(shù)據(jù)處理用Microsoft Excel for Windows 2010完成,用SPSS 25進(jìn)行統(tǒng)計分析,采用Origin lab 9.1進(jìn)行作圖。

    2 結(jié)果

    2.1 CH4排放

    2007年和2018年水稻生長季,稻田CH4排放通量的季節(jié)變化如圖1A和1B所示。各處理CH4排放通量的變化趨勢基本一致,即在水稻移栽初期緩慢上升,并在持續(xù)淹水期的第11 ~ 39 天達(dá)到季節(jié)排放峰值(6.0 ~ 58.9 mg/(m2·h)),烤田時降低到最低值,烤田后復(fù)水及后期干濕交替期間排放通量均較低。在第2年,與S0相比,S1、S2和S3的最大峰值分別增加126%、378% 和720%;第13年,與S0相比,S1、S2和S3的最大峰值分別增加167%、538% 和886%,且其最大峰值的到達(dá)時間分別前移8、16和16 d。

    如表2所示,各處理稻田CH4季節(jié)排放總量為27.7 ~ 301.0 kg/hm2,隨著秸稈還田量的增加,稻田CH4排放量呈線性增加(圖2A,<0.01)。在第2年,S1、S2和S3的CH4排放量比S0分別高2.1倍、4.8倍和9.4倍(<0.05);第13年,S1、S2和S3的CH4排放量比S0分別高2.2倍、5.3倍和7.9倍(<0.05),但其單位秸稈還田量對CH4排放的平均促進(jìn)效應(yīng)比第2年降低3%(表2)。

    2.2 N2O排放

    2007年和2018年水稻生長季,稻田N2O排放通量的季節(jié)變化如圖1C和1D所示。各處理N2O排放通量的變化趨勢基本一致,即在水稻移栽初期和施分蘗肥時均有少量N2O排放,N2O排放通量在烤田期間迅速上升并達(dá)到季節(jié)排放峰值(470 ~ 1 219 μg/(m2·h)),復(fù)水后迅速下降,施穗肥和干濕交替期間有少量N2O排放。隨著秸稈還田量的增加,各處理的N2O排放通量在第2年和第13年的變化趨勢不同。在第2年,各處理的N2O排放通量隨著秸稈還田量的增加呈現(xiàn)下降趨勢,其中,與S0相比,S1、S2和S3的最大峰值分別降低了19%、21% 和47%;而第13年,各處理的N2O排放通量隨著秸稈還田量的增加呈現(xiàn)先上升后下降的趨勢,其中,與S0相比,S1、S2和S3的最大峰值分別增加了63%、237% 和159%。

    圖1 2007年和2018年水稻生長季CH4和N2O排放通量的季節(jié)變化

    表2 2007年和2018年稻田CH4排放量、N2O排放量、土壤有機碳含量、土壤固碳速率和凈溫室效應(yīng)

    注:同列數(shù)據(jù)小寫字母不同表示處理間差異顯著(<0.05)。

    圖2 2007年和2018年水稻生長季秸稈還田量與CH4排放量(A)、N2O排放量(B)、土壤固碳速率(C)和凈溫室效應(yīng)(D)的關(guān)系

    如表2所示,各處理稻田N2O排放量的變幅是N 0.60 ~ 2.22 kg/hm2,隨著秸稈還田量的增加,各處理的N2O排放量在第2年和第13年的變化趨勢不同。在第2年,各處理的N2O排放量隨秸稈還田量的增加呈現(xiàn)下降趨勢,均未達(dá)到顯著水平,與秸稈還田量無顯著線性關(guān)系(圖2B);而第13年,各處理的N2O排放量隨秸稈還田量的增加呈線性增加(圖2B,<0.01),與S0相比,S1、S2和S3的N2O排放量分別增加了55%、171% 和122%,S2和S3達(dá)到顯著水平(<0.05)。

    2.3 土壤有機碳含量和固碳速率

    2007年和2018年稻田土壤的有機碳含量如表2所示,各處理的有機碳含量變幅為9.97 ~ 12.23 g/kg,在第2年和第13年,各處理的有機碳含量均隨秸稈還田量的增加而增加。在第2年,與S0相比,S1、S2和S3的SOC含量增加了1% ~ 5%,均未達(dá)到顯著水平;第13年,與S0相比,S1、S2和S3的SOC含量分別增加了6%、17% 和22%,S2和S3達(dá)到顯著水平(<0.05)。

    如表2所示,在第13年,各處理的土壤固碳速率變幅為0.40 ~ 0.50 t/(hm2·a),各處理間無顯著差異,與秸稈還田量無線性關(guān)系(圖2C)。

    2.4 稻田凈溫室效應(yīng)

    2018年的稻田凈溫室效應(yīng)如表2所示,各處理的稻田凈溫室效應(yīng)變幅為 –0.35 ~ 6.00 t/(hm2·a),隨秸稈還田量的增加呈線性增加(圖2D,<0.01)。在第13年,秸稈還田處理的凈溫室效應(yīng)為正值,表現(xiàn)為稻田溫室系統(tǒng)溫室氣體的源,與S0相比,S1、S2和S3的凈溫室效應(yīng)分別顯著增加了1.61、4.31和6.35 t/(hm2·a)(<0.05)。

    3 討論

    研究顯示,秸稈無論是短期還田還是長期還田,均顯著促進(jìn)稻田CH4排放(圖1和表2),且稻田CH4排放量與秸稈還田量之間呈線性正相關(guān)(圖2A),這與以往大量研究結(jié)果一致[13-14, 29]。一方面,還田的秸稈可以為產(chǎn)甲烷菌提供充足的產(chǎn)甲烷基質(zhì);另一方面,秸稈在淹水條件下分解,加速降低土壤的氧化還原電位(Eh),為產(chǎn)甲烷菌提供適宜的生長代謝環(huán)境,最終促進(jìn)稻田的CH4排放[25]。隨著秸稈還田量的增加,提供的產(chǎn)甲烷基質(zhì)隨之增加,秸稈分解消耗更多的O2,進(jìn)一步降低土壤Eh,稻田CH4排放量也隨之增加[14]。然而,隨著秸稈還田年限的增加,其單位秸稈還田量對稻田CH4排放的平均促進(jìn)效應(yīng)呈下降趨勢(表 2),這說明增排CH4的作用在減弱。其可能原因是:長期秸稈還田更有利于稻田土壤大團聚體形成,土壤空隙增多(本研究中,土壤容重由2007年的1.2 g/cm3下降至2018年的1.0 g/cm3),土壤含氧量增加,有利于稻田土壤中甲烷氧化菌的生長代謝,促進(jìn)甲烷氧化過程,從而抵消部分秸稈還田對稻田CH4排放的促進(jìn)效應(yīng)[30-32]。Jiang 等[30]的田間觀測也發(fā)現(xiàn),在稻–麥輪作系統(tǒng),秸稈還田13 a對CH4排放的促進(jìn)效應(yīng)比秸稈還田8 a降低73%。

    稻田排放的N2O主要由土壤微生物主導(dǎo)的硝化和反硝化作用過程產(chǎn)生[25]。秸稈分解可以為硝化和反硝化作用提供能量和基質(zhì)[25],進(jìn)而影響稻田N2O排放。目前,秸稈還田對稻田N2O排放的影響還存在爭議[10, 29, 33-35]。本研究結(jié)果顯示,在第2年,秸稈還田減少N2O排放,但影響并不顯著(圖1和表2),其原因可能是:高C/N的秸稈還田后,微生物在分解秸稈過程中需要利用土壤中的有效氮,使硝化和反硝化作用的底物減少,從而降低N2O排放[14, 29, 36]。鄒建文等[35]在短期秸稈還田后,也觀測到N2O排放總量減少18%。柴凱斌[37]的研究進(jìn)一步發(fā)現(xiàn),與不還田相比,1/3秸稈還田、2/3秸稈還田和全量秸稈還田后,N2O排放量分別降低了18%、39% 和38%。第13年,秸稈還田顯著增加N2O排放,且排放量隨秸稈還田量的增加呈線性增加趨勢(圖2B),其原因可能是:秸稈富含氮、磷、鉀和有機質(zhì)等,長期還田后,土壤全氮含量增加[38](本研究中,土壤全氮含量由2007年的1.1 g/kg增至2018年的1.5 g/kg),可為稻田土壤中氮轉(zhuǎn)化提供更多的基質(zhì),促進(jìn)稻田土壤的硝化和反硝化作用,從而增加稻田N2O排放[33, 39]。Wang 等[39]在稻麥輪作系統(tǒng)中也觀測到秸稈還田7 ~ 9 a顯著增加土壤全氮含量(從初始的1.4 g/kg增至1.6 ~ 1.8 g/kg),從而促進(jìn)稻田N2O排放58% ~ 155%。

    秸稈還田促進(jìn)稻田土壤固碳(表 2),這與前人的研究結(jié)果一致[11, 19]。秸稈還田后,在土壤微生物的作用下分解,部分有機碳進(jìn)入土壤,增加土壤碳庫,從而促進(jìn)稻田土壤固碳[17]。隨著秸稈還田量的增加,進(jìn)入土壤碳庫的有機碳隨之增加(表2)。李成芳等[17]也觀測到隨著秸稈還田量的增加,土壤的固碳量也隨之增加。但秸稈還田對土壤固碳的促進(jìn)作用受還田年限的影響,隨著秸稈還田年限的增加,一方面土壤固碳量增加,與初始狀態(tài)相比,其固碳潛力有所下降;另一方面長期秸稈還田會增強對碳的分解潛力,大量的碳被分解,只有少部分作為穩(wěn)定性碳保留在土壤中,因此長期秸稈還田對土壤固碳的促進(jìn)作用減弱[40-41]。Tian等[42]通過Meta分析也發(fā)現(xiàn),土壤固碳速率隨著秸稈還田年限的增加有相似的下降。

    本研究中各處理的稻田凈溫室效應(yīng)范圍是–0.35 ~ 6.00 t /(hm2·a)(表2),在前期研究結(jié)果的變化范圍內(nèi)(–11.7 ~ 9.9 t/(hm2·a))[15, 19, 22-23, 43]。在第13年,秸稈還田顯著增加稻田GWP(CH4+ N2O)1.7 ~ 5.8倍(表2),而土壤固碳速率的范圍是0.40 ~ 0.50 t/(hm2·a) (表 2),此時秸稈還田處理的土壤固碳對減緩溫室效應(yīng)的貢獻(xiàn)無法抵消秸稈還田增排溫室氣體所引起的溫室效應(yīng),因此,凈溫室效應(yīng)為正值,表現(xiàn)為碳源。Xia等[24]在太湖流域稻–麥輪作系統(tǒng)中進(jìn)行長期定位試驗,也發(fā)現(xiàn)長期秸稈還田對稻田生態(tài)系統(tǒng)溫室效應(yīng)的加劇作用大于減緩作用,其增排溫室氣體造成的溫室效應(yīng)比固碳所減少的溫室效應(yīng)高3.2倍~ 3.9倍。

    4 結(jié)論

    秸稈還田第2年顯著促進(jìn)稻田CH4排放和土壤固碳,N2O排放無顯著差異;但第13年顯著促進(jìn)稻田CH4和N2O排放,其增加的溫室效應(yīng)完全抵消了土壤固碳的減排效益,稻田生態(tài)系統(tǒng)的凈溫室效應(yīng)為碳源。這意味著長期秸稈還田下,稻田是一個很重要的碳源。在當(dāng)前水稻輕簡化和機械化栽培情況下,秸稈還田幾乎是全量還田,這與本研究中的全量處理情況很相似。也就是說,長期秸稈全量還田大背景下,稻田的碳源強度會很大。因此未來亟待研究適于輕簡化和機械化栽培的秸稈還田技術(shù)措施,可以通過秸稈多種綜合利用方式(如生物質(zhì)炭和秸稈快速好氧分解等),以及農(nóng)業(yè)區(qū)域共同體內(nèi)的秸稈循環(huán)利用(如將生產(chǎn)的秸稈施用到茶園、果園或菜地中),以實現(xiàn)稻田的固碳減排和農(nóng)業(yè)可持續(xù)發(fā)展。

    [1] Landman W. Climate change 2013: The physical cience basis[Z]. IPCC, Geneva. 2013.

    [2] WMO. The state of greenhouse gases in the atmosphere based on global observations through 2019[Z]. 2020. https://library.wmo.int/doc_num.php?explnum_id=10437

    [3] FAO. FAO statistical databases for crops[Z]. Roma. 2019. http://www.fao.org/faostat/en/#data/QC

    [4] Yue Q, Ledo A, Cheng K, et al. re-assessing nitrous oxide emissions from croplands across Mainland China[J]. Agriculture, Ecosystems & Environment, 2018, 268: 70–78.

    [5] 張曉慶, 王梓凡, 參木友, 等. 中國農(nóng)作物秸稈產(chǎn)量及綜合利用現(xiàn)狀分析[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報, 2021, 26(9): 30–41.

    [6] 董建新, 叢萍, 劉娜, 等. 秸稈深還對黑土亞耕層土壤物理性狀及團聚體分布特征的影響[J]. 土壤學(xué)報, 2021, 58(4): 921–934.

    [7] 趙穎, 周楓, 羅佳琳, 等. 水稻秸稈還田配施肥料對小麥產(chǎn)量和氮素利用的影響[J]. 土壤, 2021, 53(5): 937–944.

    [8] 吳鵬年, 王艷麗, 侯賢清, 等. 秸稈還田配施氮肥對寧夏揚黃灌區(qū)滴灌玉米產(chǎn)量及土壤物理性狀的影響[J]. 土壤, 2020, 52(3): 470–475.

    [9] 石祖梁. 中國秸稈資源化利用現(xiàn)狀及對策建議[J]. 世界環(huán)境, 2018(5): 16–18.

    [10] 張剛. 太湖地區(qū)稻麥兩熟制農(nóng)田秸稈還田綜合效應(yīng)研究[D]. 南京: 南京林業(yè)大學(xué), 2020.

    [11] 逯非, 王效科, 韓冰, 等. 稻田秸稈還田: 土壤固碳與甲烷增排[J]. 應(yīng)用生態(tài)學(xué)報, 2010, 21(1): 99–108.

    [12] 史然, 陳曉娟, 沈建林, 等. 稻田秸稈還田的土壤增碳及溫室氣體排放效應(yīng)和機理研究進(jìn)展[J]. 土壤, 2013, 45(2): 1193–1198.

    [13] 馬靜, 徐華, 蔡祖聰. 施肥對稻田甲烷排放的影響[J]. 土壤, 2010, 42(2): 153–163.

    [14] 蔣靜艷, 黃耀, 宗良綱. 水分管理與秸稈施用對稻田CH4和N2O排放的影響[J]. 中國環(huán)境科學(xué), 2003, 23(5): 552–556.

    [15] 張岳芳, 陳留根, 朱普平, 等. 秸稈還田對稻麥兩熟高產(chǎn)農(nóng)田凈增溫潛勢影響的初步研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2012, 31(8): 1647–1653.

    [16] Zhang M, Li B, Xiong Z Q. Effects of organic fertilizer on net global warming potential under an intensively managed vegetable field in southeastern China: A three-year field study[J]. Atmospheric Environment, 2016, 145: 92–103.

    [17] 李成芳, 寇志奎, 張枝盛, 等. 秸稈還田對免耕稻田溫室氣體排放及土壤有機碳固定的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2011, 30(11): 2362–2367.

    [18] 張翰林, 呂衛(wèi)光, 鄭憲清, 等. 不同秸稈還田年限對稻麥輪作系統(tǒng)溫室氣體排放的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2015, 23(3): 302–308.

    [19] 馮曉赟, 萬鵬, 李潔, 等. 秸稈還田與氮肥配施對中南地區(qū)稻田土壤固碳和溫室氣體排放的影響[J]. 農(nóng)業(yè)資源與環(huán)境學(xué)報, 2016, 33(6): 508–517.

    [20] 李靜, 馮淑怡, 陳利根, 等. 秸稈還田對稻田溫室氣體排放的影響: Meta分析——以長江中下游地區(qū)為例[J]. 中國人口·資源與環(huán)境, 2016, 26(5): 91–100.

    [21] 張麗. 南方典型稻作系統(tǒng)綠肥和秸稈還田對溫室氣體排放特征和環(huán)境效應(yīng)的影響[D]. 南京: 南京農(nóng)業(yè)大學(xué), 2018.

    [22] 湯宏, 曾掌權(quán), 沈健林, 等. 秸稈與水分管理稻田的溫室氣體排放和碳固定[J]. 環(huán)境科學(xué)與技術(shù), 2021, 44(1): 41–48.

    [23] 劉曉雨, 李志鵬, 潘根興, 等. 長期不同施肥下太湖地區(qū)稻田凈溫室效應(yīng)和溫室氣體排放強度的變化[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2011, 30(9): 1783–1790.

    [24] 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[J]. Agriculture, Ecosystems & Environment, 2014, 197: 118–127.

    [25] 蔡祖聰, 徐華, 馬靜. 稻田生態(tài)系統(tǒng)CH4和N2O排放[M]. 合肥: 中國科學(xué)技術(shù)大學(xué)出版社, 2009.

    [26] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社, 2000.

    [27] Lu F, Wang X K, Han B, et al. Soil carbon sequestrations by nitrogen fertilizer application, straw return and no-tillage in China’s cropland[J]. Global Change Biology, 2009, 15(2): 281–305.

    [28] Shang Q Y, Yang X X, Gao C M, et al. Net annual global warming potential and greenhouse gas intensity in Chinese double rice-cropping systems: A 3-year field measurement in long-term fertilizer experiments[J]. Global Change Biology, 2011, 17(6): 2196–2210.

    [29] 張岳芳, 鄭建初, 陳留根, 等. 麥秸還田與土壤耕作對稻季CH4和N2O排放的影響[J]. 生態(tài)環(huán)境學(xué)報, 2009, 18(6): 2334–2338.

    [30] Jiang Y, Qian H Y, Huang S, et al. Acclimation of methane emissions from rice paddy fields to straw addition[J]. Science Advances, 2019, 5(1): eaau9038.

    [31] Fonte S J, Yeboah E, Ofori P, et al. Fertilizer and residue quality effects on organic matter stabilization in soil aggregates[J]. Soil Science Society of America Journal, 2009, 73(3): 961–966.

    [32] Hütsch B W. Methane oxidation in non-flooded soils as affected by crop production—invited paper[J]. European Journal of Agronomy, 2001, 14(4): 237–260.

    [33] 王海候, 沈明星, 陸長嬰, 等. 不同秸稈還田模式對稻麥兩熟農(nóng)田稻季甲烷和氧化亞氮排放的影響[J]. 江蘇農(nóng)業(yè)學(xué)報, 2014, 30(4): 758–763.

    [34] 張艷波, 彭其安. 水稻秸稈還田對稻田土壤N2O排放的影響[J]. 湖北農(nóng)業(yè)科學(xué), 2016, 55(10): 2539–2543, 2554.

    [35] 鄒建文, 黃耀, 宗良綱, 等. 不同種類有機肥施用對稻田CH4和N2O排放的綜合影響[J]. 環(huán)境科學(xué), 2003, 24(4): 7–12.

    [36] 馬靜, 徐華, 蔡祖聰, 等. 墑溝埋草對稻田CH4和N2O排放的影響[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2008, 24(4): 27–31.

    [37] 柴凱斌. 秸稈還田對稻麥系統(tǒng)作物產(chǎn)量及溫室氣體排放的影響[D]. 武漢: 華中農(nóng)業(yè)大學(xué), 2018.

    [38] 周子軍, 郭松, 陳琨, 等. 長期秸稈覆蓋對免耕稻-麥產(chǎn)量、土壤氮組分及微生物群落的影響[J]. 土壤學(xué)報, 2022, 59(4): 1148–1159.

    [39] Wang H H, Shen M X, Hui D F, et al. Straw incorporation influences soil organic carbon sequestration, greenhouse gas emission, and crop yields in a Chinese rice (L.)-wheat (L.) cropping system[J]. Soil and Tillage Research, 2019, 195: 104377.

    [40] Zhou P, Sheng H, Li Y, et al. Lower C sequestration and N use efficiency by straw incorporation than manure amendment on paddy soils[J]. Agriculture, Ecosystems & Environment, 2016, 219: 93–100.

    [41] Cheng W G, Padre A T, Sato C, et al. Changes in the soil C and N contents, C decomposition and N mineralization potentials in a rice paddy after long-term application of inorganic fertilizers and organic matter[J]. Soil Science and Plant Nutrition, 2016, 62(2): 212–219.

    [42] Tian K, Zhao Y C, Xu X H, et al. Effects of long-term fertilization and residue management on soil organic carbon changes in paddy soils of China: A meta- analysis[J]. Agriculture, Ecosystems & Environment, 2015, 204: 40–50.

    [43] 李桂花, 周吉祥, 張建峰, 等. 有機肥和緩控肥替代部分化肥降低雙季稻田綜合凈溫室效應(yīng)[J]. 植物營養(yǎng)與肥料學(xué)報, 2020, 26(6): 1017–1024.

    Effects of Straw Incorporation Years and Rates on Net Global Warming Potential in Paddy Fields

    HUANG Qiong1,2, ZHU Xiaoli1,2, SHEN Wanyu1,2, FAN Di1,2, ZHANG Guangbin1, MA Jing1, XU Hua1*

    (1 State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China; 2 University of Chinese Academy of Sciences, Beijing 100049, China)

    In this paper, the effects of different straw incorporation years (2 and 13 years) and rates (0, 1.6, 3.2 and 4.8 t/hm2) on CH4and N2O emissions from paddy fields, soil organic carbon (SOC) sequestration rate and net global warming potential (Net GWP) were studied. The results show that CH4emissions is increased linearly with the increase of straw incorporation rate both in 2 and 13 years, but the average promoting effect of unit straw incorporation amount on CH4emissions in 13 years is 3% lower than that in 2 years. For N2O emissions, no significant difference is found in 2 years, but the emissions is significantly increased by 55% ~ 171% and shows a linear increase with the increase of straw incorporation rate in 13 years. Meanwhile, SOC sequestration rate varies from 0.40 to 0.50 t/(hm2·a) in 13 years without significant difference among the treatments. Therefore, Net GWP shows a linear increase with the increase of straw incorporation rate, which is carbon source in 13 years. The findings suggest that the effect of straw incorporation amount on Net GWP is regulated by the duration, and the increased GWP gradually completely offsets the emission reduction benefit of soil carbon sequestration and then exacerbates global warming.Therefore, the paddy field is an important carbon source after long-term straw incorporation, and it is urgent to optimize straw incorporation technology to further sequester carbon and reduce emissions to achieve sustainable agricultural development.

    Duration of straw incorporation; Straw incorporation rate; Greenhouse gas emissions; Soil organic carbon sequestration; Net global warming potential

    S131;X511

    A

    10.13758/j.cnki.tr.2022.05.006

    黃瓊,朱小莉,沈皖豫, 等. 秸稈還田年限及還田量對稻田凈溫室效應(yīng)的影響. 土壤, 2022, 54(5): 912–919.

    國家自然科學(xué)基金項目(41877325,42177233)和中國科學(xué)院青年創(chuàng)新促進(jìn)會項目(2018349)資助。

    (hxu@issas.ac.cn)

    黃瓊(1987—),女,江蘇邳州人,博士研究生,主要從事農(nóng)田溫室氣體排放及土壤碳氮循環(huán)研究。E-mail: qhuang@issas.ac.cn

    猜你喜歡
    田量溫室效應(yīng)通量
    模擬地球的氣候變化—溫室效應(yīng)小實驗
    不同秸稈還田量對設(shè)施番茄生長的影響
    冬小麥田N2O通量研究
    不同農(nóng)作物秸稈直接還田及其還田量對大豆生長的影響
    秸稈還田下灌水量對土壤水分運移特性研究
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    淺議城市進(jìn)化程度與溫室效應(yīng)
    河南科技(2014年18期)2014-02-27 14:15:09
    溫室效應(yīng)、低碳經(jīng)濟與我國減排對策分析
    江西省碳通量時空演變與模型構(gòu)建
    涩涩av久久男人的天堂| 丰满迷人的少妇在线观看| 有码 亚洲区| 国产男女内射视频| 国产综合精华液| 在线观看免费高清a一片| 夜夜骑夜夜射夜夜干| 在现免费观看毛片| 9热在线视频观看99| 人妻人人澡人人爽人人| 久久人人爽人人片av| 久久久久久人人人人人| 一区福利在线观看| 老司机影院毛片| 男女下面插进去视频免费观看| 国产福利在线免费观看视频| 亚洲综合色惰| 天美传媒精品一区二区| 2021少妇久久久久久久久久久| 亚洲国产精品成人久久小说| 国产成人精品一,二区| 久久久精品免费免费高清| 亚洲男人天堂网一区| 欧美激情 高清一区二区三区| 丝袜美腿诱惑在线| 香蕉国产在线看| 国产成人91sexporn| 国产成人精品婷婷| 99久国产av精品国产电影| 婷婷成人精品国产| 亚洲国产毛片av蜜桃av| a级片在线免费高清观看视频| 尾随美女入室| 免费观看在线日韩| 亚洲av日韩在线播放| 高清黄色对白视频在线免费看| 天天操日日干夜夜撸| 在线观看免费日韩欧美大片| 亚洲精品乱久久久久久| 一区二区三区四区激情视频| 黄色视频在线播放观看不卡| 伊人久久国产一区二区| 欧美变态另类bdsm刘玥| 国产欧美亚洲国产| 大话2 男鬼变身卡| 国产黄频视频在线观看| 国产成人精品婷婷| 一本大道久久a久久精品| 少妇熟女欧美另类| 亚洲第一av免费看| 日韩熟女老妇一区二区性免费视频| 欧美xxⅹ黑人| 欧美黄色片欧美黄色片| 蜜桃国产av成人99| 人妻少妇偷人精品九色| 国产在视频线精品| 蜜桃在线观看..| 免费看av在线观看网站| 日韩 亚洲 欧美在线| 少妇人妻精品综合一区二区| 国产精品熟女久久久久浪| 亚洲av免费高清在线观看| www日本在线高清视频| 日本欧美国产在线视频| www.熟女人妻精品国产| 欧美激情高清一区二区三区 | 免费观看a级毛片全部| 亚洲精品久久成人aⅴ小说| 18禁观看日本| 国产精品蜜桃在线观看| 韩国av在线不卡| 欧美精品亚洲一区二区| 王馨瑶露胸无遮挡在线观看| 日日摸夜夜添夜夜爱| 丰满迷人的少妇在线观看| 久久久久人妻精品一区果冻| 成人影院久久| 国产成人精品婷婷| 最新中文字幕久久久久| 国产男女内射视频| 热99国产精品久久久久久7| 韩国高清视频一区二区三区| 午夜免费观看性视频| 亚洲精华国产精华液的使用体验| 久久女婷五月综合色啪小说| 美女大奶头黄色视频| 男人操女人黄网站| 国产精品欧美亚洲77777| 人人澡人人妻人| 午夜日韩欧美国产| 久久久久久久国产电影| 老汉色av国产亚洲站长工具| 国产av精品麻豆| 精品酒店卫生间| 日韩一区二区视频免费看| 在线亚洲精品国产二区图片欧美| 蜜桃国产av成人99| 在线观看www视频免费| 国产精品99久久99久久久不卡 | 免费看不卡的av| 亚洲精品一二三| 精品亚洲成国产av| 久久久国产欧美日韩av| av在线app专区| 亚洲第一av免费看| 看十八女毛片水多多多| 亚洲人成电影观看| 男女边吃奶边做爰视频| 亚洲国产精品成人久久小说| 亚洲少妇的诱惑av| 久久女婷五月综合色啪小说| 一区二区av电影网| 久久久久久久精品精品| 久久久精品区二区三区| videos熟女内射| 日韩在线高清观看一区二区三区| 两性夫妻黄色片| 男女国产视频网站| 97人妻天天添夜夜摸| 国产乱人偷精品视频| 亚洲男人天堂网一区| 在线精品无人区一区二区三| 久久人妻熟女aⅴ| 亚洲精品自拍成人| 国产熟女欧美一区二区| 久久久久精品人妻al黑| 日韩在线高清观看一区二区三区| 免费少妇av软件| 国产精品久久久av美女十八| 午夜福利视频在线观看免费| 亚洲国产精品国产精品| 亚洲av在线观看美女高潮| 久久精品国产亚洲av天美| 丝瓜视频免费看黄片| 婷婷成人精品国产| www.自偷自拍.com| 欧美97在线视频| 亚洲成色77777| 国产精品成人在线| 亚洲av电影在线观看一区二区三区| 欧美人与性动交α欧美软件| 午夜福利视频在线观看免费| 欧美日韩综合久久久久久| 午夜av观看不卡| 多毛熟女@视频| 一本大道久久a久久精品| av线在线观看网站| 国产又爽黄色视频| 亚洲欧美成人精品一区二区| 午夜福利一区二区在线看| 18禁国产床啪视频网站| 99国产精品免费福利视频| 亚洲精品美女久久久久99蜜臀 | 色网站视频免费| 久久久久久久久久人人人人人人| av.在线天堂| 天天操日日干夜夜撸| 91aial.com中文字幕在线观看| 欧美精品一区二区免费开放| 99热国产这里只有精品6| 日韩一卡2卡3卡4卡2021年| 日韩三级伦理在线观看| av福利片在线| 亚洲美女视频黄频| 久久久久久伊人网av| 有码 亚洲区| 国产一区二区 视频在线| 美女xxoo啪啪120秒动态图| 亚洲av.av天堂| 女人久久www免费人成看片| 国产精品不卡视频一区二区| 欧美日本中文国产一区发布| 97人妻天天添夜夜摸| 精品国产乱码久久久久久男人| 电影成人av| 国产男女内射视频| 日韩欧美精品免费久久| 精品99又大又爽又粗少妇毛片| 久久久久精品人妻al黑| 99热国产这里只有精品6| 一二三四中文在线观看免费高清| 中文字幕制服av| 亚洲精品美女久久久久99蜜臀 | 99久久综合免费| 久久久久久久大尺度免费视频| 亚洲国产精品999| 精品人妻偷拍中文字幕| 亚洲综合色网址| 久久久久人妻精品一区果冻| 黄色毛片三级朝国网站| 亚洲av免费高清在线观看| 韩国高清视频一区二区三区| 亚洲精品一区蜜桃| 国产精品香港三级国产av潘金莲 | 午夜福利乱码中文字幕| 高清视频免费观看一区二区| 国产成人精品在线电影| 午夜福利在线观看免费完整高清在| 亚洲精品,欧美精品| 最黄视频免费看| 久久久亚洲精品成人影院| 久久精品久久久久久噜噜老黄| 中文字幕色久视频| 精品第一国产精品| 人人妻人人添人人爽欧美一区卜| 精品少妇内射三级| 免费日韩欧美在线观看| 国产 精品1| 麻豆乱淫一区二区| 18禁观看日本| av在线老鸭窝| 欧美日韩亚洲高清精品| 伊人久久国产一区二区| h视频一区二区三区| 国产成人一区二区在线| 国产白丝娇喘喷水9色精品| 亚洲国产精品一区二区三区在线| 只有这里有精品99| a级毛片在线看网站| av福利片在线| 色哟哟·www| 毛片一级片免费看久久久久| 久久人人97超碰香蕉20202| 国产一区二区三区综合在线观看| 欧美日韩精品成人综合77777| 久久久久久久久免费视频了| 男男h啪啪无遮挡| 不卡av一区二区三区| 一边摸一边做爽爽视频免费| 午夜精品国产一区二区电影| 亚洲av在线观看美女高潮| 80岁老熟妇乱子伦牲交| av有码第一页| 精品人妻偷拍中文字幕| 久久这里只有精品19| 一边亲一边摸免费视频| 成人18禁高潮啪啪吃奶动态图| 毛片一级片免费看久久久久| 国产精品久久久久成人av| 9191精品国产免费久久| 狂野欧美激情性bbbbbb| 18禁动态无遮挡网站| 老熟女久久久| 久久精品熟女亚洲av麻豆精品| 国产成人一区二区在线| 国产有黄有色有爽视频| 熟女电影av网| 电影成人av| 国产欧美日韩综合在线一区二区| 最近最新中文字幕免费大全7| 中文精品一卡2卡3卡4更新| 久久久久久久久免费视频了| 亚洲四区av| 亚洲精品在线美女| 男女无遮挡免费网站观看| 九草在线视频观看| av.在线天堂| 欧美人与善性xxx| 精品酒店卫生间| 国产成人精品福利久久| 国产成人精品在线电影| 国产成人精品久久久久久| 青草久久国产| 久久久久视频综合| 我的亚洲天堂| 婷婷色av中文字幕| 巨乳人妻的诱惑在线观看| 99国产精品免费福利视频| 成人亚洲精品一区在线观看| 一级,二级,三级黄色视频| 亚洲国产av新网站| 丰满迷人的少妇在线观看| 欧美国产精品va在线观看不卡| 欧美成人精品欧美一级黄| 久久毛片免费看一区二区三区| 国产成人精品久久二区二区91 | 精品少妇一区二区三区视频日本电影 | av在线app专区| 亚洲欧美一区二区三区黑人 | 中文字幕制服av| 大香蕉久久网| 亚洲精品乱久久久久久| 麻豆精品久久久久久蜜桃| 国产欧美亚洲国产| 一级毛片 在线播放| 国产免费一区二区三区四区乱码| 国产无遮挡羞羞视频在线观看| 69精品国产乱码久久久| www.熟女人妻精品国产| 日韩欧美一区视频在线观看| 亚洲精品国产一区二区精华液| 亚洲国产精品国产精品| 成人国产麻豆网| 超碰97精品在线观看| 成人毛片a级毛片在线播放| 一区二区三区乱码不卡18| 国产成人精品福利久久| 亚洲精品av麻豆狂野| 欧美黄色片欧美黄色片| 亚洲一级一片aⅴ在线观看| 日韩 亚洲 欧美在线| 欧美bdsm另类| 日韩精品免费视频一区二区三区| 天天躁夜夜躁狠狠躁躁| av在线观看视频网站免费| 日韩在线高清观看一区二区三区| 熟妇人妻不卡中文字幕| 最近中文字幕高清免费大全6| xxxhd国产人妻xxx| 日本欧美视频一区| 亚洲美女视频黄频| 亚洲欧美精品综合一区二区三区 | 欧美变态另类bdsm刘玥| 999久久久国产精品视频| 成人影院久久| 日韩制服骚丝袜av| 久久精品aⅴ一区二区三区四区 | 国产精品免费大片| 观看av在线不卡| 91精品伊人久久大香线蕉| 亚洲精华国产精华液的使用体验| 日韩三级伦理在线观看| 熟妇人妻不卡中文字幕| 国产 一区精品| av女优亚洲男人天堂| 少妇的逼水好多| 一级毛片 在线播放| 超碰97精品在线观看| 欧美xxⅹ黑人| 一本久久精品| 久久久久久久久久久免费av| 观看美女的网站| 国产av一区二区精品久久| 精品久久久久久电影网| 久热久热在线精品观看| 少妇人妻久久综合中文| 在线免费观看不下载黄p国产| www.av在线官网国产| av一本久久久久| 丝袜喷水一区| 日韩一卡2卡3卡4卡2021年| 亚洲五月色婷婷综合| 免费观看av网站的网址| 午夜日韩欧美国产| 在线精品无人区一区二区三| 另类亚洲欧美激情| 最近的中文字幕免费完整| 精品人妻熟女毛片av久久网站| 韩国av在线不卡| 一级毛片电影观看| 99热网站在线观看| 亚洲美女黄色视频免费看| 亚洲精品日韩在线中文字幕| 色吧在线观看| 日韩制服丝袜自拍偷拍| 哪个播放器可以免费观看大片| 亚洲欧美精品综合一区二区三区 | 一本色道久久久久久精品综合| 中文字幕亚洲精品专区| 亚洲色图 男人天堂 中文字幕| 国产片特级美女逼逼视频| 精品国产一区二区三区久久久樱花| 日本爱情动作片www.在线观看| 国产不卡av网站在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产午夜精品一二区理论片| av网站在线播放免费| 91成人精品电影| 最新中文字幕久久久久| 亚洲国产欧美日韩在线播放| 91aial.com中文字幕在线观看| 一级a爱视频在线免费观看| 寂寞人妻少妇视频99o| 好男人视频免费观看在线| 国产一区二区三区综合在线观看| 国产一区二区 视频在线| 日韩欧美一区视频在线观看| 一级片'在线观看视频| 国产亚洲av片在线观看秒播厂| 免费播放大片免费观看视频在线观看| 黄色 视频免费看| 日产精品乱码卡一卡2卡三| av女优亚洲男人天堂| 寂寞人妻少妇视频99o| 啦啦啦中文免费视频观看日本| 女人高潮潮喷娇喘18禁视频| 婷婷色麻豆天堂久久| 国产成人免费无遮挡视频| 91成人精品电影| 免费高清在线观看视频在线观看| 天天躁日日躁夜夜躁夜夜| 在线观看免费视频网站a站| 国产一区二区三区av在线| 亚洲欧美一区二区三区久久| 久久久欧美国产精品| 免费在线观看黄色视频的| 国产成人免费无遮挡视频| 青春草视频在线免费观看| 婷婷色综合大香蕉| 亚洲一区中文字幕在线| 精品酒店卫生间| av网站在线播放免费| 欧美 日韩 精品 国产| 久久精品国产亚洲av高清一级| 9色porny在线观看| 亚洲国产成人一精品久久久| 欧美精品一区二区大全| 国精品久久久久久国模美| 日本91视频免费播放| 1024视频免费在线观看| 99久久综合免费| 99re6热这里在线精品视频| 国产一区二区 视频在线| 人人妻人人爽人人添夜夜欢视频| 啦啦啦中文免费视频观看日本| 国产精品.久久久| 嫩草影院入口| 日本免费在线观看一区| 亚洲欧美中文字幕日韩二区| 亚洲精品视频女| 欧美成人午夜精品| 亚洲色图 男人天堂 中文字幕| 国产欧美日韩一区二区三区在线| 亚洲综合色惰| 欧美在线黄色| 精品一区在线观看国产| 欧美黄色片欧美黄色片| 美女福利国产在线| 国产熟女欧美一区二区| 曰老女人黄片| 亚洲精品自拍成人| 啦啦啦啦在线视频资源| 99香蕉大伊视频| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| av在线老鸭窝| 色哟哟·www| 国产伦理片在线播放av一区| 在线观看免费高清a一片| 免费高清在线观看日韩| 欧美国产精品va在线观看不卡| 永久免费av网站大全| 老司机影院成人| 国产精品 国内视频| 最近2019中文字幕mv第一页| 久久精品亚洲av国产电影网| 国产精品蜜桃在线观看| 国产麻豆69| 日日啪夜夜爽| 大片电影免费在线观看免费| 成人午夜精彩视频在线观看| 青春草亚洲视频在线观看| 涩涩av久久男人的天堂| 国产精品人妻久久久影院| 一区二区日韩欧美中文字幕| 校园人妻丝袜中文字幕| 国产精品一区二区在线观看99| 又黄又粗又硬又大视频| 纵有疾风起免费观看全集完整版| 国产精品一二三区在线看| 免费在线观看黄色视频的| 少妇 在线观看| 美女大奶头黄色视频| 青春草国产在线视频| 波多野结衣一区麻豆| 亚洲美女搞黄在线观看| 一本—道久久a久久精品蜜桃钙片| 青春草视频在线免费观看| 在线观看国产h片| 亚洲精品美女久久av网站| 亚洲国产最新在线播放| 久久久久久久久久久久大奶| 性色av一级| 国产国语露脸激情在线看| 一级爰片在线观看| 日韩av免费高清视频| 少妇的逼水好多| 久久97久久精品| 午夜福利在线免费观看网站| 亚洲国产精品成人久久小说| 国产精品av久久久久免费| 欧美国产精品va在线观看不卡| 日韩一本色道免费dvd| 日韩中文字幕欧美一区二区 | 人妻 亚洲 视频| 国产精品女同一区二区软件| 日本黄色日本黄色录像| 少妇猛男粗大的猛烈进出视频| 香蕉丝袜av| 欧美精品亚洲一区二区| 波多野结衣一区麻豆| 精品亚洲成国产av| 欧美另类一区| 国产欧美亚洲国产| 最近中文字幕高清免费大全6| av视频免费观看在线观看| 国产麻豆69| 蜜桃国产av成人99| 中文字幕人妻丝袜制服| 欧美日韩视频高清一区二区三区二| 亚洲国产精品国产精品| 高清黄色对白视频在线免费看| 久久99热这里只频精品6学生| 精品视频人人做人人爽| 精品久久久久久电影网| 国产无遮挡羞羞视频在线观看| 中文字幕人妻丝袜制服| 亚洲色图综合在线观看| 狠狠婷婷综合久久久久久88av| 亚洲欧美一区二区三区国产| av免费观看日本| 国产一区有黄有色的免费视频| 制服人妻中文乱码| 久久女婷五月综合色啪小说| 亚洲欧美精品综合一区二区三区 | 韩国av在线不卡| 国产麻豆69| 国产成人免费观看mmmm| 少妇的逼水好多| 一二三四中文在线观看免费高清| 在线看a的网站| 亚洲精品中文字幕在线视频| 日韩大片免费观看网站| 久久精品aⅴ一区二区三区四区 | 80岁老熟妇乱子伦牲交| 少妇猛男粗大的猛烈进出视频| 纵有疾风起免费观看全集完整版| 女的被弄到高潮叫床怎么办| 久久久久国产一级毛片高清牌| 午夜福利,免费看| 久久精品国产综合久久久| 国产福利在线免费观看视频| 国产亚洲一区二区精品| 久久久久久伊人网av| 两个人免费观看高清视频| 一区二区三区四区激情视频| 国产免费福利视频在线观看| 另类精品久久| 成年女人毛片免费观看观看9 | 久久精品久久久久久久性| 久久这里只有精品19| 欧美成人精品欧美一级黄| 久久精品亚洲av国产电影网| 久久精品久久久久久久性| 亚洲av在线观看美女高潮| 亚洲av中文av极速乱| 精品少妇一区二区三区视频日本电影 | 侵犯人妻中文字幕一二三四区| 亚洲av中文av极速乱| 天堂俺去俺来也www色官网| 久久99蜜桃精品久久| 最近2019中文字幕mv第一页| 久久久久久久久久久久大奶| 免费黄色在线免费观看| 在线看a的网站| 日韩 亚洲 欧美在线| 国产精品成人在线| 久久久久精品性色| 久久av网站| 欧美日韩一级在线毛片| 国产精品久久久久久精品电影小说| 免费看不卡的av| 最近最新中文字幕免费大全7| 中文字幕制服av| 看非洲黑人一级黄片| 国产精品久久久久久精品古装| 国产精品国产三级专区第一集| 曰老女人黄片| 又大又黄又爽视频免费| 精品人妻在线不人妻| 男男h啪啪无遮挡| 91午夜精品亚洲一区二区三区| 亚洲av在线观看美女高潮| 成人漫画全彩无遮挡| 午夜福利在线观看免费完整高清在| 国产无遮挡羞羞视频在线观看| 国产成人精品婷婷| 成人免费观看视频高清| 少妇熟女欧美另类| 制服诱惑二区| 国产成人精品一,二区| 黄色毛片三级朝国网站| 这个男人来自地球电影免费观看 | 国产av码专区亚洲av| 激情视频va一区二区三区| 一本色道久久久久久精品综合| 欧美日本中文国产一区发布| 亚洲欧美中文字幕日韩二区| 国产一区二区在线观看av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美黄色片欧美黄色片| 老司机影院成人| 精品国产一区二区久久| 亚洲精品一二三| 久久国内精品自在自线图片| 国产熟女午夜一区二区三区| 精品人妻一区二区三区麻豆| 高清av免费在线| 欧美人与善性xxx| 香蕉丝袜av| 久久午夜综合久久蜜桃| 国产欧美亚洲国产| 婷婷色综合www| 亚洲欧美精品自产自拍| 国产日韩欧美视频二区| 中文字幕人妻熟女乱码| videos熟女内射| 欧美精品国产亚洲| 亚洲,欧美,日韩| 黄色毛片三级朝国网站| 精品国产露脸久久av麻豆| 曰老女人黄片| 日本色播在线视频| 老女人水多毛片| 国产激情久久老熟女|