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

    不同灌溉施肥措施對(duì)夏玉米-冬小麥農(nóng)田N2O排放和產(chǎn)量的影響

    2018-08-22 03:18:14李昊儒郝衛(wèi)平梅旭榮

    李昊儒,郝衛(wèi)平,梅旭榮,郭 瑞

    (中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所,農(nóng)業(yè)農(nóng)村部旱作節(jié)水農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室,北京 100081)

    0 引 言

    氧化亞氮(N2O)是大氣中主要的溫室氣體之一,對(duì)全球氣候變化起到重要作用,也是導(dǎo)致臭氧層破壞的光化學(xué)反應(yīng)主要參與者[1],其體積分?jǐn)?shù)已由 1750年的270×10–9上升到 2005 年的 319×10–9,相比工業(yè)化前提高了18%[2]。農(nóng)業(yè)活動(dòng)是N2O濃度增加的主要原因之一,2005年中國農(nóng)業(yè)活動(dòng)導(dǎo)致的N2O排放占N2O總排放量的73.79%[3]。華北平原是中國重要的糧食生產(chǎn)基地,該地區(qū)耕地面積占全國的 27%,小麥和玉米產(chǎn)量分別占全國的50%和 35%[4]。夏玉米-冬小麥一年兩熟是華北平原主要的種植制度,大水漫灌和撒施氮肥是該地區(qū)普遍應(yīng)用的灌溉施肥方式,常規(guī)的周年施氮量約為550~600 kg/hm2[5],嚴(yán)重高于作物需求氮量[6]。這種“大水高肥”的粗放型管理措施不僅導(dǎo)致水肥利用效率低,而且還會(huì)增加N2O等溫室氣體排放,引起一系列環(huán)境問題。李虎等[7]利用DNDC(DeNitrification-DeComposition)模型模擬在傳統(tǒng)農(nóng)業(yè)管理措施下,該區(qū)域麥–玉輪作農(nóng)田每年氮素通過N2O排放損失的量達(dá)到17.71 kg/hm2。因此,在保證糧食產(chǎn)量的同時(shí),減少氮肥用量、降低N2O排放成為華北平原農(nóng)業(yè)科學(xué)研究的重點(diǎn)內(nèi)容之一。

    不同的灌溉施肥方式能改變土壤結(jié)構(gòu)和理化性質(zhì),造成土壤通氣性、有效氮分布和微生物有效性發(fā)生改變,從而對(duì)土壤硝化、反硝化以及 N2O排放產(chǎn)生影響[8]。與傳統(tǒng)的管理方式相比,滴灌水肥一體化技術(shù)可以根據(jù)作物需求通過微噴灌系統(tǒng)將水肥施入作物根區(qū),使水分和養(yǎng)分在土壤中均勻分布,達(dá)到農(nóng)田局部集中施肥和灌水的效果,以保證養(yǎng)分被根系快速吸收,大幅度地提高了肥料的利用率,降低了因過量施肥而造成的環(huán)境污染問題[9]。2012年,東北三省區(qū)實(shí)施“節(jié)水增糧”行動(dòng),投資380億元推廣噴滴灌施肥技術(shù)(253萬hm2)。2014年,啟動(dòng)華北地下水超采區(qū)綜合治理試點(diǎn),將滴灌施肥列為關(guān)鍵技術(shù),大力推廣小麥和玉米滴灌技術(shù)應(yīng)用面積,在河北省進(jìn)行大面積推廣。農(nóng)業(yè)部也視滴灌施肥為現(xiàn)代農(nóng)業(yè)“一號(hào)技術(shù)”,2016年4月,農(nóng)業(yè)部辦公廳印發(fā)《推進(jìn)水肥一體化實(shí)施方案(2016-2020年)》,進(jìn)一步明確了新形勢(shì)下推進(jìn)滴灌施肥技術(shù)發(fā)展的總體思路和目標(biāo)任務(wù)。目前,關(guān)于滴灌和滴灌水肥一體化對(duì)華北平原農(nóng)田 N2O排放的影響研究[10-11]主要分別針對(duì)玉米田和麥田,對(duì)玉米-小麥輪作農(nóng)田土壤N2O周年排放特征的研究少有報(bào)道,并且對(duì)滴灌水肥一體化下不同施氮量農(nóng)田土壤N2O排放規(guī)律和特征的研究較少。針對(duì)中國目前的資源環(huán)境現(xiàn)狀,在國家各項(xiàng)政策的扶持推動(dòng)下,滴灌水肥一體化技術(shù)在大田作物生產(chǎn)的應(yīng)用日趨廣泛,本文研究對(duì)等氮量下不同灌溉施肥技術(shù)以及在滴灌水肥一體化下減氮的產(chǎn)量效果和環(huán)境效應(yīng)進(jìn)行了比較,具有一定的生產(chǎn)指導(dǎo)意義。本研究以華北平原夏玉米-冬小麥輪作系統(tǒng)為研究對(duì)象,設(shè)置了常規(guī)施氮量下傳統(tǒng)灌溉施肥、滴灌+傳統(tǒng)施肥、滴灌水肥一體化以及滴灌水肥一體化下不同施氮量處理,以滴灌+不施氮肥為對(duì)照,研究了不同施灌溉施肥模式對(duì)土壤N2O排放特征及作物產(chǎn)量的影響;明確了既保證糧食產(chǎn)量又減少土壤N2O排放的合理氮用量及灌溉施肥模式,為滴灌水肥一體化技術(shù)的應(yīng)用提供理論依據(jù)和技術(shù)指導(dǎo)。

    1 材料與方法

    1.1 試驗(yàn)地概況

    試驗(yàn)于2015年6月至2016年6月在中國農(nóng)業(yè)科學(xué)院河北廊坊試驗(yàn)基地進(jìn)行。該試驗(yàn)站位于廊坊市萬莊鎮(zhèn),屬暖溫帶大陸性季風(fēng)氣候,四季分明,冷熱干濕差異明顯。廊坊市年平均氣溫(1971-2004年)為11.9 ℃,年平均降水量(1971-2004年)為554.9 mm,降水季節(jié)分布不均,多集中在 6-8月份。該地區(qū)土壤類型為砂壤質(zhì)潮土,基本性質(zhì)見表1。距試驗(yàn)田塊30 m處安裝有自動(dòng)氣象站,可以連續(xù)觀測(cè)氣溫、風(fēng)速、風(fēng)向、輻射、降雨等氣象數(shù)據(jù),觀測(cè)期內(nèi)降雨量和氣溫見圖1。

    表1 廊坊試驗(yàn)基地土壤理化性質(zhì)Table 1 Soil physicochemical properties of in Langfang experimental station

    圖1 觀測(cè)期降雨量和氣溫的變化Fig.1 Precipitation and air temperature during whole observed period

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

    供試夏玉米品種為廊坊地區(qū)主栽品種鄭單958,播種日期為2015年6月15日,收獲日期為2015年9月30日。供試冬小麥品種為廊坊地區(qū)主栽品種廊研 43,播種日期為2015年10月3日,收獲日期為2016年6月12日。試驗(yàn)設(shè)置了常規(guī)施氮量(夏玉米:205.5 kg/hm2;冬小麥:250 kg/hm2)下傳統(tǒng)灌溉施肥(FP100%)、滴灌+傳統(tǒng)施肥(DN100%)、滴灌水肥一體化(FN100%)以及滴灌水肥一體化下不同施氮量(減氮60%(FN40%)、減氮 30%(FN70%)、常規(guī)氮量(FN100%)和增氮 30%(FN130%))和滴灌+不施氮肥(CK)共 7個(gè)處理。每個(gè)處理3次重復(fù),試驗(yàn)小區(qū)面積為42 m2(6 m′7 m),每個(gè)小區(qū)之間設(shè)2m寬的隔離帶,冬小麥平均行距為30 cm,平均播種密度為4.2′106株/hm2,夏玉米平均行距為60 cm,株距為30 cm,播種密度為5.6′104株/hm2。

    FP100%處理播種前施用復(fù)合肥(N∶P∶K=15∶15∶15)作為基肥,拔節(jié)期施用尿素追肥,施氮肥方式為撒施,灌水方式為溝灌;DN100%處理施氮肥種類為尿素,分別在播種前施基肥,拔節(jié)期追肥,施氮肥方式為撒施;滴灌水肥一體化處理施氮肥種類均為尿素,施氮肥方式為隨灌水施入氮肥;各處理均施用相同的磷、鉀肥做底肥,夏玉米季用量為67.5 kg/hm2;冬小麥季用量為165 kg/hm2。滴灌和滴灌水肥一體化處理在作物關(guān)鍵需水需肥期,根據(jù)測(cè)定的土壤實(shí)際含水量設(shè)置灌溉量,夏玉米整個(gè)生育期共進(jìn)行 4次灌水,冬小麥季整個(gè)生育期共進(jìn)行6次灌水,具體施肥灌溉時(shí)間和量詳見表2。

    試驗(yàn)采用以色列Netafim公司同一型號(hào)的滴灌帶,灌水過程中的滴頭工作壓力為 15~20 kPa,滴頭流量為0.6~0.65 L/h,額定流量為1.9 L/h。為了解試驗(yàn)田塊灌水的均勻性,在灌水季節(jié)內(nèi),對(duì)滴頭流量在試驗(yàn)田塊內(nèi)的變化進(jìn)行了多次測(cè)試,結(jié)果表明,滴頭流量均勻系數(shù)在95%以上。滴頭間距為30 cm,夏玉米滴灌帶布置間距為60 cm(每行玉米布設(shè)一條毛管),冬小麥滴灌帶布置間距為30 cm(每?jī)尚行←湶荚O(shè)1條毛管)。試驗(yàn)水源采用地下井水,經(jīng)滴灌首部過濾后進(jìn)入比例式施肥泵(Mis Rite Model 2504,Tefen,以色列)進(jìn)水管道,每個(gè)小區(qū)各接一個(gè)獨(dú)立的施肥泵,并聯(lián)接一個(gè)儲(chǔ)液罐,施肥開始前按各小區(qū)所需氮肥分別加入儲(chǔ)液罐,將儲(chǔ)液罐充滿水,充分?jǐn)嚢?,使其完全溶解?/p>

    1.3 N2O的采集和分析方法

    采用密閉式靜態(tài)箱法[12]測(cè)定N2O,采集N2O的裝置為下端開口上部密封,直徑25 cm×高25 cm的柱形圓筒,由PVC材料制成,圓柱體頂部安裝有一個(gè)帶有螺旋開關(guān)的氣密性氣體抽氣閥門,在試驗(yàn)小區(qū)中部,以滴灌帶為直徑,滴頭為圓心,安裝一個(gè)插入土壤 5 cm深的 PVC塑料圓形底座,底座內(nèi)徑為21.5 cm,底座外徑為31.5 cm。采樣時(shí),將箱體置于底座上,向底座中倒入少量水,進(jìn)行水封確保密封性良好。密封后立即用100 mL注射器抽取氣體并注入帶有聚乙烯涂層的鋁箔真空氣袋內(nèi)(大連普萊特公司)內(nèi),每隔10 min取樣1次,共取樣5次,并準(zhǔn)確記錄采樣時(shí)間和箱內(nèi)氣溫。取樣結(jié)束后,立即將采樣箱從底座移開。每次滴灌施肥后逐日采樣7 d;每次10 mm以上日降雨后,逐日采樣3 d;每次播種后連續(xù)取樣3 d;其他情況下,3-11月每周采樣1次,12月、1月和2月份每2周采樣一次。為了便于比較和消除日變化所導(dǎo)致的通量差異,取樣監(jiān)測(cè)時(shí)間維持在當(dāng)?shù)貢r(shí)間的9:00-11:00之間進(jìn)行。每次取氣樣的同時(shí),測(cè)定氣溫、土壤5 cm深度溫度及15 cm深度土壤含水量。N2O氣體樣品在中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所采用改進(jìn)的氣相色譜儀(Agilent 7890A)分析測(cè)定。

    表2 不同灌溉施肥處理下夏玉米-冬小麥?zhǔn)┑獣r(shí)期和施氮量、灌溉時(shí)期和灌溉量Table 2 Nitrogen (N) application date, rates and irrigation date, amounts under different treatments during summer maize- winter wheat growing stage

    1.4 計(jì)算公式及數(shù)據(jù)分析

    土壤N2O排放通量計(jì)算公式為[13]

    式中F為N2O的排放通量,mg/m2·h,正值為排放,負(fù)值為吸收;r為標(biāo)準(zhǔn)大氣壓下N2O的密度,g/L;H為采樣箱氣室高度,cm;T為采樣箱內(nèi)氣溫,℃;P為采樣時(shí)氣壓,kPa;P0為標(biāo)準(zhǔn)大氣壓,kPa,P/P0≈1;/ctDD為采樣箱內(nèi) N2O 濃度的變化速率,uL/L·min。用內(nèi)插法計(jì)算相鄰兩次監(jiān)測(cè)之間未觀測(cè)日期的排放通量,然后將觀測(cè)值和未觀測(cè)日計(jì)算值逐日累加得到N2O 排放總量。

    土壤孔隙含水率(water-filled pore spaces, WFPS)計(jì)算公式為[14]

    式中vq為土壤體積含水量,%;BD為土壤容重,g/cm3;2.65為土壤密度,g/cm3。

    采用SAS9.2統(tǒng)計(jì)分析軟件處理試驗(yàn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析和差異顯著性檢驗(yàn),采用LSD方法對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行方差分析和顯著性測(cè)驗(yàn)。作圖采用Origin 8.0軟件完成。

    2 結(jié)果與分析

    2.1 夏玉米季農(nóng)田土壤N2O排放日通量的變化特征

    由圖2可以看出,在夏玉米整個(gè)觀測(cè)期內(nèi),N2O排放通量的增加多發(fā)生在降水或灌溉施肥后的1~4 d,4 d后N2O排放通量逐漸減弱,最后直至消失,這說明玉米季土壤硝化、反硝化作用基本在降水或灌溉施肥后4 d內(nèi)完成。CK處理除了7月16日-18日集中降雨后,N2O排放出現(xiàn)明顯的排放峰外,其余時(shí)間N2O排放波動(dòng)較小,排放通量最高值發(fā)生在7月19日,達(dá)到0.17 mg/m2·h。滴灌水肥一體化不同施氮量處理N2O排放動(dòng)態(tài)變化趨勢(shì)基本一致,且施氮量越低,N2O排放通量越小,F(xiàn)N40%處理N2O排放通量最高值僅為0.10 mg/m2·h,這說明N2O排放通量與施氮量密切相關(guān)。DN100%處理的 N2O排放通量具有較大的波動(dòng),變化范圍在0.00~1.64 mg/m2·h之間,最大排放峰出現(xiàn)在7月18日降雨后,持續(xù)5 d左右,峰值出現(xiàn)在7月19日。FP100%處理的第一次排放峰出現(xiàn)在基肥+第一次灌溉后,持續(xù)5 d左右,峰值出現(xiàn)在6月 20日,達(dá)到 0.27 mg/m2·h;與 DN100%處理一致,F(xiàn)P100%處理的最大排放峰出現(xiàn)在7月18日降雨后,峰值出現(xiàn)在7月19日,達(dá)到0.34 mg/m2·h,這說明在傳統(tǒng)撒施氮肥下,降雨對(duì)玉米田N2O排放影響較大。FN100%處理的N2O排放通量最高值(0.27 mg/m2·h)出現(xiàn)在抽雄期滴灌施肥后。夏玉米季 FN40%處理(DN100%除外)的 N2O平均排放通量較 FP100%處理顯著降低 70.8%(P<0.05)(圖 2)。

    圖2 不同灌溉施肥方式下2015年夏玉米季土壤N2O排放通量Fig.2 N2O emission flux under different management in summer maize fields during 2015 growing season

    2.2 冬小麥季農(nóng)田土壤N2O排放日通量的變化特征

    由圖 3可以看出,在冬小麥整個(gè)觀測(cè)期內(nèi),CK和FP100%處理的N2O最大排放通量出現(xiàn)在撒施基肥后,分別達(dá)到0.31和0.21 mg/m2·h,其余時(shí)間波動(dòng)較小,這說明傳統(tǒng)灌溉施肥下,施肥對(duì)麥田 N2O影響較大。DN100%和FN100%處理排放規(guī)律基本一致,第一個(gè)排放峰出現(xiàn)在第一次灌溉后,持續(xù)6 d左右,峰值出現(xiàn)在10月10日,分別達(dá)0.20和0.10 mg/m2·h;第二個(gè)排放峰出現(xiàn)在返青期灌水后,峰值出現(xiàn)在灌溉后第三天(3月10日),分別達(dá)到0.88和0.56 mg/m2·h。與玉米季相同,滴灌水肥一體化不同施氮量處理N2O排放動(dòng)態(tài)變化趨勢(shì)基本一致,且施氮量越低,N2O排放通量越小,F(xiàn)N40%處理的最大排放通量發(fā)生在第一次滴灌施肥以及4月2日降雨后的第二天,均為0.05 mg/m2·h。與FP100%處理相比,F(xiàn)N40%處理在冬小麥季的 N2O平均排放通量顯著降低 66.7%(P<0.05)(圖 3)。

    圖3 不同灌溉施肥方式下2015-2016年冬小麥季土壤N2O排放通量Fig.3 N2O emission flux under different management in winter wheat fields during 2015-2016 growing season

    2.3 農(nóng)田N2O排放的影響因素

    土壤溫濕度影響微生物活性,進(jìn)而影響土壤硝化反硝化作用,最終影響農(nóng)田土壤N2O的產(chǎn)生和向大氣中的擴(kuò)散,為明確不同灌溉施肥方式下農(nóng)田N2O排放的影響因素,需對(duì)各處理土壤溫濕度變化進(jìn)行進(jìn)一步分析。研究結(jié)果顯示,在滴灌水肥一體化處理中,減氮60%的N2O排放通量最小,顯著改善環(huán)境效應(yīng),因此,在農(nóng)田土壤N2O排放的影響因素和作物產(chǎn)量構(gòu)成分析中只采用FN40%處理試驗(yàn)結(jié)果與其他處理進(jìn)行分析比較。

    2.3.1 土壤溫度

    各處理0~5 cm土層平均溫度均表現(xiàn)出相似的變化趨勢(shì),夏玉米季先上升后下降,冬小麥季先下降后上升,7月10日土壤溫度達(dá)到最高值(28.4~30.9 ℃),冬小麥越冬期至返青期,土壤溫度低于0 ℃,1月16日土壤溫度達(dá)到最低值(–10.0 ℃)(圖4)。通過剔除灌溉施肥引起的N2O排放峰值,對(duì)不同灌溉施肥措施N2O排放通量與5 cm深度土壤溫度的相關(guān)分析表明,整個(gè)輪作周期內(nèi)CK和FN40%處理5 cm深度土壤溫度與N2O排放通量呈現(xiàn)極顯著的正相關(guān)關(guān)系(P<0.01),其他處理5 cm深度土壤溫度與N2O排放通量之間無顯著相關(guān)性(表3),這說明對(duì)整個(gè)輪作周期而言,低施氮量處理土壤溫度升高可顯著促進(jìn)農(nóng)田 N2O排放。而對(duì)玉米和小麥生長季分別進(jìn)行分析發(fā)現(xiàn),在夏玉米季,DN100%和FN100%處理N2O排放通量與5 cm土壤溫度呈現(xiàn)顯著正相關(guān)關(guān)系,F(xiàn)N40%處理呈顯著負(fù)相關(guān)關(guān)系(P<0.05),CK和FP100%處理N2O排放通量與5 cm深度土壤溫度的關(guān)系均無統(tǒng)計(jì)顯著性。在冬小麥季,DN100%和FP100%處理N2O排放通量與5 cm土壤溫度呈現(xiàn)顯著負(fù)相關(guān)關(guān)系,F(xiàn)N40%呈顯著正相關(guān)關(guān)系(P<0.05),CK和FN100%處理無顯著相關(guān)性。

    表3 不同灌溉施肥措施下夏玉米和冬小麥農(nóng)田N2O排放通量與0~5 cm土層地溫相關(guān)性分析Table 3 Correlation analysis between N2O emission fluxes and soil temperature in 0-5 cm soil depth under different management of summer maize and winter wheat fields

    圖4 不同灌溉施肥措施下夏玉米-冬小麥田(0~5 cm)土壤平均溫度變化特征Fig.4 Dynamics of average soil temperature with different irrigation and fertilization methods in layer (0–5 cm) of summer maize-winter wheat field

    2.3.2 土壤濕度

    夏玉米季整個(gè)觀測(cè)期土壤孔隙含水率在 19.4%~42.2%范圍內(nèi)變化,冬小麥返青期后土壤孔隙含水率在3.7%~58.1%范圍內(nèi)變化,每次灌溉降雨后土壤水分含量略有上升(圖5)。通過對(duì)各處理N2O排放通量與土壤孔隙含水率進(jìn)行相關(guān)分析發(fā)現(xiàn),周年內(nèi)兩者之間存在顯著的正相關(guān)關(guān)系(P<0.05),其中,F(xiàn)N100%和FP100%達(dá)極顯著水平(P<0.01)(表4),這說明對(duì)整個(gè)輪作周期而言,增加土壤濕度可顯著提高農(nóng)田 N2O排放。就 FP100%處理而言,高于整個(gè)觀測(cè)期平均值(0.05 mg/m2·h)的N2O排放通量中有57.1%出現(xiàn)在WFPS在29%~35.5%的條件下,其余42.9%出現(xiàn)在WFPS<29%和WFPS>23.5%的水分條件下,從該結(jié)果可以看出,29%~35.5%的濕度條件有利于N2O的產(chǎn)生和排放。對(duì)各處理玉米季和小麥季分別進(jìn)行分析發(fā)現(xiàn),夏玉米季,CK、FP100%和 FN40%處理N2O排放通量與WFPS呈現(xiàn)極顯著正相關(guān)關(guān)系(P<0.01);冬小麥季,DN100%、FN100%和FP100%處理N2O排放通量與WFPS呈現(xiàn)極顯著正相關(guān)關(guān)系(P<0.01)(表4)。

    圖5 不同耕作措施下夏玉米和冬小麥田(0~20 cm)土壤充水孔隙度變化特征Fig.5 Dynamics of water filled porosity space under irrigation and fertilization managements in layer (0–20 cm) of summer maize and winter wheat field

    表4 不同灌溉施肥措施下夏玉米和冬小麥農(nóng)田N2O排放通量與土壤孔隙含水率相關(guān)性分析Table 4 Correlation analysis between N2O emission fluxes and WFPS under different management of summer maize and winter wheat fields

    2.4 不同灌溉施肥方式下夏玉米-冬小麥產(chǎn)量與排放系數(shù)

    與CK處理相比,夏玉米季各施肥處理產(chǎn)量顯著增加了21.8%~31.4%;冬小麥季產(chǎn)量增加了42.3%~123.3%;周年產(chǎn)量顯著增加了 28.8%~53.4%。與 FP100%處理相比,F(xiàn)N40%處理夏玉米季產(chǎn)量顯著降低了7.4%,冬小麥產(chǎn)量和周年產(chǎn)量并沒有顯著性差異(P<0.05)。從產(chǎn)量構(gòu)成要素來看,夏玉米FN40%的穗粒數(shù)和有效穗數(shù)沒有顯著減少,千粒重顯著低于其他灌溉施肥處理(DN100%除外),與CK相比,F(xiàn)N40%的有效穗數(shù)和千粒質(zhì)量均顯著增加(表5);冬小麥FN40%的每穗粒數(shù)、有效穗數(shù)和千粒質(zhì)量均沒有顯著低于其他灌溉施肥處理,且千粒質(zhì)量顯著高于DN100%,與CK相比,F(xiàn)N40%的有效穗數(shù)顯著增加(表 5)。綜上所述,與其他灌溉施肥處理相比,夏玉米FN40%處理產(chǎn)量的減少主要是千粒質(zhì)量降低的結(jié)果,而冬小麥FN40%處理產(chǎn)量沒有顯著性降低是由于每穗粒數(shù)、有效穗數(shù)和千粒質(zhì)量均沒有顯著性減少。

    表5 不同灌溉施肥處理夏玉米和冬小麥產(chǎn)量及其構(gòu)成要素Table 5 Maize and wheat grain yield and yield components under different irrigation and fertilization management

    從夏玉米-冬小麥周年輪作來看,不同處理夏玉米-冬小麥輪作農(nóng)田土壤N2O排放總量為0.98~8.84 kg/hm2,排放次序?yàn)椋篋N100%>FN100%>FP100%>FN40%>CK,CK處理N2O年排放總量最低,為(0.98±0.28)kg/hm2;與FP100%相比,F(xiàn)N40%處理N2O年排放總量顯著減少了1.80 kg/hm2,降幅為62.9%;DN100%處理N2O年排放總量顯著高于其他處理。對(duì)各處理玉米和小麥生長季的N2O排放總量分別進(jìn)行分析發(fā)現(xiàn),玉米季N2O排放總量占全年 N2O排放總量的 25.8%~49.2%,小麥季占50.8%~74.2%。在夏玉米季,與 FP100%相比,CK和FN40%處理分別減少62.3%和58.4%(P<0.05),CK處理排放N2O最少,F(xiàn)N40%次之;而DN100%處理N2O排放總量增加208.0%,這與N2O排放通量規(guī)律一致。在冬小麥季,與FP100%相比,CK和FN40%處理的N2O排放總量分別降低68.4%和66.4%,CK處理排放N2O最少,F(xiàn)N40%排放N2O次之,但二者差異不顯著(表6)。

    夏玉米-冬小麥輪作農(nóng)田土壤 N2O的排放系數(shù)介于0.04%~1.72%之間,其次序是:DN100%>FN100%>FP100%>FN40%,其中,F(xiàn)N100%、FP100%和FN40%處理排放系數(shù)均低于Bouwman[15]提供的糧田土壤N2O排放系數(shù)1.25%。不同灌溉施肥方式下冬小麥生長季N2O排放系數(shù)為0.01%~1.79%,其中,F(xiàn)P100%和FN40%處理低于IPCC建議氮素肥料N2O-N排放系數(shù)的1%,夏玉米生長季N2O排放系數(shù)在0.06%~1.65%之間。綜合考慮各處理N2O的排放量、玉米-小麥產(chǎn)量和施肥量等因素,本研究認(rèn)為,F(xiàn)N40%玉米-小麥輪作周期內(nèi)農(nóng)田土壤N2O排放總量為1.06 kg/hm2,排放系數(shù)為0.04%,顯著低于其他處理,周年產(chǎn)量為11 493.65 kg/hm2,與FP100%處理相比沒有顯著差異(表6)。

    表6 不同灌溉施肥方式下夏玉米-冬小麥輪作產(chǎn)量、農(nóng)田N2O排放總量、平均排放通量和排放系數(shù)Table 6 Grain yields, cumulative emissions, average N2O emission flux and emission factors coefficient during winter wheat-summer maize growing season

    3 討論

    3.1 管理措施對(duì)農(nóng)田N2O排放的影響

    N2O是硝化反應(yīng)的重要副產(chǎn)物和反硝化反應(yīng)的中間產(chǎn)物,不同的灌溉方式能改變土壤結(jié)構(gòu)和理化性質(zhì),造成土壤通氣性等發(fā)生改變,從而對(duì)土壤硝化、反硝化以及 N2O排放產(chǎn)生影響[9]。已有研究顯示,漫灌條件下農(nóng)田土壤N2O排放主要來自于有氧條件下的硝化反應(yīng)[16],而溝灌條件下農(nóng)田土壤N2O排放主要來自于厭氧條件下的反硝化作用。也有研究認(rèn)為漫灌條件下N2O主要來源于反硝化反應(yīng),硝化反應(yīng)是滴灌條件下 N2O的主要來源[17]。本試驗(yàn)結(jié)果顯示,與FP100%相比,DN100%處理玉米田和小麥田平均排放通量分別增加了 237.5%和116.7%,這與郭樹芳等[18]的研究結(jié)果規(guī)律一致,這可能是因?yàn)榈喂喾绞较峦寥浪趾扛哂诼?,促進(jìn)了土壤微生物作用下的硝化和反硝化反應(yīng),從而增加了土壤N2O排放量;也可能由于撒施氮肥后進(jìn)行漫灌造成施入的氮肥淋入土壤深層,減少土壤表層多余氮素殘留,從而減少FP100%處理農(nóng)田土壤N2O排放。而Wu等[19]和Wang等[10]分別對(duì)新疆棉田和華北麥田的研究認(rèn)為,與漫灌相比,滴灌減少了28.8%和14.6%的農(nóng)田N2O排放量,并認(rèn)為滴灌方式下土壤含水量相對(duì)低,且土壤局部濕潤,改善了土壤通氣性,提高了O2含量,促進(jìn)硝化作用,在一定范圍內(nèi)抑制了反硝化反應(yīng)。

    另一方面,氮肥施用方式也影響農(nóng)田土壤N2O排放。與DN100%相比,F(xiàn)N100%處理的玉米田和小麥田N2O排放總量分別減少69.0%和30.8%,這說明在滴灌條件下,根據(jù)作物的生長發(fā)育需求隨水施肥比傳統(tǒng)的撒施氮肥有明顯的N2O減排效果,這可能是因?yàn)榉峙⒎制谑┓逝c單次施肥相比,可以保證施肥時(shí)間與養(yǎng)分需求的時(shí)間同步,提高作物對(duì)氮肥的吸收利用效率,減少了土壤氮?dú)埩?,降低了土壤硝化和反硝化作用需要的氮源,從而降低土?N2O排放量[20-21]。同時(shí),本試驗(yàn)結(jié)果顯示,與FP100%相比,F(xiàn)N100%處理玉米田和輪作農(nóng)田N2O排放總量沒有顯著性差異,麥田土壤N2O顯著增加了114.2%,而Tian等[11]對(duì)華北平原玉米田的研究證明,與漫灌撒施氮肥相比,滴灌施肥減少了7.7%的農(nóng)田土壤N2O排放量。這可能是由于本研究中翻耕及秸稈還田帶來土壤底物的差異,從而增加了灌溉施肥方式間土壤溫室氣體排放的差異。另外,氮肥施用量也顯著影響農(nóng)田土壤N2O排放,過量施氮會(huì)造成N2O排放量增加。本試驗(yàn)結(jié)果顯示,與FN100%相比,F(xiàn)N40%處理的玉米田和小麥田 N2O排放總量分別減少了56.7%和84.3%,周年排放總量顯著減少了 77.1%(P<0.05),而周年產(chǎn)量并沒有顯著降低,并且在滴灌水肥一體化下,增加施氮量會(huì)促進(jìn)農(nóng)田N2O排放。

    以上分析可以看出,不同灌溉施肥方式土壤含水量存在差異,并且通過影響土壤水分分布狀況影響無機(jī)氮空間運(yùn)移,影響作物吸收效率,造成對(duì)農(nóng)田土壤N2O的排放產(chǎn)生區(qū)別。綜上所述,在該地區(qū)夏玉米-冬小麥輪作制度下,若采用滴灌方式進(jìn)行灌溉,則根據(jù)作物需肥規(guī)律同時(shí)采用水肥一體化方式進(jìn)行施肥才有增產(chǎn),減少農(nóng)田N2O排放的效果,并且在該技術(shù)下,減少60%施氮量在不減產(chǎn)的同時(shí),有顯著的N2O減排效果。

    3.2 環(huán)境因子對(duì)農(nóng)田N2O排放的影響

    降水所引起的土壤快速再濕潤能激發(fā)N2O形成與排放。本試驗(yàn)結(jié)果顯示,夏玉米季常規(guī)施肥各處理均在 7月19日觀測(cè)到排放峰,這是由于7月16日-18日出現(xiàn)連續(xù)降水(26.4、28.3和81 mm),土壤含水量顯著增加,WFPS達(dá)到35.8%~42.2%,而此時(shí)的土壤溫度(25.6 ℃~27.3 ℃)相對(duì)較高,微生物活性較強(qiáng),土壤硝化與反硝化的共同作用促使土壤N2O的大量產(chǎn)生和排放。

    土壤溫度影響微生物的代謝活動(dòng),并通過影響其硝化和反硝化反應(yīng)速率來影響N2O的排放[22]。本試驗(yàn)得出,不同灌溉施肥處理土壤溫度與N2O 排放通量的相關(guān)性并不相同,玉米季DN100%處理土壤溫度與N2O排放通量相關(guān)性最好,呈極顯著正相關(guān)(P<0.01),小麥季FP100%處理相關(guān)性最好,呈極顯著負(fù)相關(guān)(P<0.01)。這可能由于不同灌溉方式下土壤含水率不同,影響了硝化和反硝化反應(yīng)的發(fā)生,從而影響N2O形成和排放。謝軍飛等[23]研究認(rèn)為,麥田N2O排放通量隨著土壤溫度的升高在不同程度上有一定的增加,這與本研究 CK、FN100%和FN40%處理的結(jié)果相符。

    土壤孔隙含水率對(duì)N2O的產(chǎn)生和擴(kuò)散也有很重要的影響。Granli等研究認(rèn)為,WFPS在70%~90%,N2O排放主要是由反硝化過程產(chǎn)生的,WFPS在30%~70%,則主要由硝化過程產(chǎn)生,最適宜N2O排放的土壤孔隙含水量范圍為60%~80%WFPS。當(dāng)土壤含水量低于飽和含水量時(shí),硝化作用產(chǎn)生的N2O占總產(chǎn)生量的61%~98%,N2O的排放量隨土壤水分的增加而增加,硝化作用是最基本來源[25]。本試驗(yàn)結(jié)果顯示,玉米田和小麥田N2O排放最高時(shí)期,各處理 WFPS分別為 32.2%~41.7%和31.8%~58.1%;在傳統(tǒng)灌溉施肥措施下,高于整個(gè)觀測(cè)期平均值(0.05 mg/m2·h)的N2O排放通量中有57.1%出現(xiàn)在WFPS在29%~35.5%的條件下,其余42.9%出現(xiàn)在WFPS<29%和WFPS>23.5%的水分條件下,這說明該地區(qū)29%~35.5%的濕度條件有利于N2O的產(chǎn)生和排放,而滴灌處理土壤孔隙含水率顯著高于常規(guī)灌溉,更有利于保持土壤水分。

    4 結(jié) 論

    1)整個(gè)觀測(cè)期內(nèi)農(nóng)田N2O排放峰出現(xiàn)在施肥、灌溉以及集中降雨后,一般持續(xù)5 d左右。夏玉米季,滴灌+傳統(tǒng)施肥和傳統(tǒng)灌溉施肥農(nóng)田N2O排放通量變化規(guī)律基本一致;冬小麥季,滴灌+傳統(tǒng)施肥、等氮量滴灌水肥一體化和減氮60%滴灌水肥一體化農(nóng)田N2O排放通量變化規(guī)律基本一致。小麥季農(nóng)田土壤N2O排放通量高于玉米季,夏玉米季土壤N2O階段排放峰值出現(xiàn)在拔節(jié)期和抽雄期;而冬小麥季土壤N2O階段排放峰值出現(xiàn)在冬前苗期和拔節(jié)期。

    2)土壤溫度和土壤濕度均影響N2O排放,但不同處理在夏玉米和冬小麥生長季與土壤溫度和土壤濕度的相關(guān)性并不相同。滴灌+不施氮肥和減氮60%滴灌水肥一體化處理輪作農(nóng)田N2O排放與土壤溫度表現(xiàn)出極顯著的正相關(guān)關(guān)系(P<0.01);傳統(tǒng)灌溉施肥處理玉米田、小麥田和輪作農(nóng)田N2O排放與土壤孔隙含水率(WFPS)均表現(xiàn)出極顯著的正相關(guān)關(guān)系(P<0.01)。

    3)從排放總量來看,不同處理夏玉米-冬小麥輪作農(nóng)田土壤N2O排放總量為0.98~8.84 kg/hm2,排放次序?yàn)椋篋N100%>FN100%>FP100%>FN40%>CK,排放系數(shù)為0.04%~1.72%。在該地區(qū)夏玉米-冬小麥輪作制度下,若采用滴灌方式進(jìn)行灌溉,則根據(jù)作物需肥規(guī)律同時(shí)采用水肥一體化方式進(jìn)行施肥才有增產(chǎn),減少農(nóng)田N2O排放的效果,并且在滴灌水肥一體化技術(shù)下,減少 60%施氮量在不減少周年產(chǎn)量的同時(shí),顯著降低農(nóng)田N2O排放。

    [1]Ravishankara A R, Daniel J S, Portmann R W. Nitrous oxide(N2O): The dominant ozone-depleting substance emitted in the 21st century[J]. Science, 2009, 326(5949): 123-125.

    [2]IPCC. Changes in Atmospheric Constituents and in Radioactive forcing[M]//Climate change 2007: The physical science basis, contribution of working group I to the Fourth assessment report of panel on climate change Cambridge University Press. Cambridge, United Kingdom and New York,NY, USA. 2007.

    [3]國家發(fā)展改革委員會(huì). 中華人民共和國氣候變化第二次國家信息通報(bào)[M]. 北京:中國計(jì)劃出版社,2011.

    [4]Du T S, Kang S Z, Sun J S, et al. An improved water use efficiency of cereals under temporal and spatial deficit irrigation in north China[J]. Agricultural Water Management,2010, 97: 66-74.

    [5]Zhao R F, Chen X P, Zhang F S, et al. Fertilization and nitrogen balance in a wheat-maize rotation system in North China[J]. Agronomy Journal, 2006, 98: 938-945.

    [6]Liu X J, Ju X T, Zhang F S, et al. Nitrogen dynamics and budgets in a winter wheat-maize cropping system in the North China Plain[J]. Field Crops Research, 2003, 83: 111-124.

    [7]李虎,王立剛,邱建軍. 基于DNDC模型的華北典型農(nóng)田氮素?fù)p失分析及綜合調(diào)控途徑及綜合調(diào)控途徑[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2012,20(4):414-421.Li Hu, Wang Ligang, Qiu Jianjun. Nitrate loss simulated with DNDC model and control technologies in typical cropland of North China[J]. Chinese Journal of Eco-Agriculture, 2012,20(4): 414-421. (in Chinese with English abstract)

    [8]Burt C, Connor K O, Ruehr T. Fertigation. Irrigation Training and Research Center[M]. California Polytechnic State University, San Luis Obispo, 1998.

    [9]Galbally I E. Factors controlling NOXemission from soils. In:Exchange of Trace Gases Between Terrestrial Ecosystems and the Atmosphere[M]. Andreae, M.O. and Schimel. D.S.,Eds. DahlemKonferenzen. Wiley, Chichester, 1989, 23-27.

    [10]Wang G S, Liang Y P, Zhang Q, et al. Mitigated CH4and N2O emissions and improved irrigation water use efficiency in winter wheat field with surface drip irrigation in the North China Plain[J]. Agricultural Water Management. 2016, 163:403-407.

    [11]Tian D, Zhang Y Y, Mu Y Y, et al. The effect of drip irrigation and drip fertigation on N2O and NO emissions,water saving and grain yields in a maize field in the North China Plain[J]. Science of Total Environment. 2017, 575:1034-1040.

    [12]山楠,趙同科,畢曉慶,等. 適宜施氮量降低京郊小麥-玉米農(nóng)田 N2O 排放系數(shù)增加產(chǎn)量[J]. 農(nóng)業(yè)工程學(xué)報(bào).2016,32(22):163-169.Shan Nan, Zhao Tongke, Bi Xiaoqing, et al. Suitable nitrogen application reducing N2O emission and improving grain yield in wheat-maize crop rotation system in Beijing suburb[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2016, 32(22): 163-170. (in Chinese with English abstract)

    [13]張婧,李虎,王立剛,等. 京郊典型設(shè)施蔬菜地土壤N2O排放特征[J]. 生態(tài)學(xué)報(bào),2014,34(14):4088-4098.Zhang Jing, Li Hu, Wang Ligang, et al. Characteristics of nitrous oxide emissions from typical greenhouse vegetable fields in Beijing suburbs[J]. Acta Ecologica Sinica, 2014,34(14): 4088-4098. (in Chinese with English abstract)

    [14]Vallejo A, Skiba U M, García-Torres L, et al. Nitrogen oxides emission from soils bearing a potato crop as influenced by fertilization with treated pig slurries and composts[J]. Soil Biology & Biochemistry, 2006, 38(9): 2782-2793.

    [15]Bouwman A F. Factors regulating nitrous oxide and nitric oxide emission[C]//Bouwman A F ,Boumans L J M , BatjesN H, eds. Global Estimates of Gaseous Emissions of NH3, NO and N2O from Agricultural Land. Rome, Italy: FAO and IFA.2001: 11-16.

    [16]Ju X T, Lu X, Gao Z, et al. Processes and factors controlling N2O production in an intensively managed low carbon calcareous soil under sub-humid monsoon conditions[J].Environmental Pollution. 2011, 159: 1007-1016.

    [17]Taryn L K, Emma C S. Reduced nitrous oxide emissions and increased yields in California tomato cropping systems under drip irrigation and fertigation[J]. Agriculture, Ecosystem &Environment. 2013. 170: 16-27.

    [18]郭樹芳,齊玉春,尹飛虎,等. 不同灌溉方式對(duì)華北平原冬小麥田土壤CO2和N2O排放通量的影響[J]. 環(huán)境科學(xué).2016,37(5):1880-1889.Guo Shufang, Qi Yuchun, Yin Feihu, et al. Effect of irrigation patterns on soil CO2and N2O emissions from winter wheat field in North China Plain[J]. Environmental Science, 2016, 37(5): 1880-1889. (in Chinese with English abstract)

    [19]Wu J, Guo W, Feng J F, et al. Greenhouse gas emissions from cotton field under different irrigation methods and fertilization regimes in arid northwestern China[J]. The Scientific World Journal, 2014(2014): 1-10.

    [20]鐘文輝,蔡祖聰,尹力初,等. 種植水稻和長期施用無機(jī)肥對(duì)紅壤氨氧化細(xì)菌多樣性和硝化作用的影響[J]. 土壤學(xué)報(bào),2008,45(1):105-111.Zhong Wenhui, Cai Zucong, Yin Lichu, et al. Effects of rice cultivation and long-term application of inorganic fertilizers on ammonium oxidizers diversity and nitrification of red soils[J]. Acta Pedologica Sinica, 2008, 45(1): 105-111. (in Chinese with English abstract)

    [21]Bouwman A F, Boumans L J M, Batjes N H. Emissions of N2O and NO from fertilized fi elds: Summary of available measurement data[J]. Global Biogeochemical Cycles. 2002,16(4): 1058.

    [22]鄭循華,王明星,王躍思,等. 溫度對(duì)農(nóng)田N2O產(chǎn)生與排放的影響[J]. 環(huán)境科學(xué),1997,18(5):1-5.Zheng Xunhua, Wang Mingxing, Wang Yuesi, et al. Impacts of temperature on N2O production and emission[J]. Chinese Journal of Environmental Science, 1997, 18(5): 1-5. (in Chinese with English abstract)

    [23]謝軍飛,李玉娥. 土壤溫度對(duì)北京旱地農(nóng)田N2O排放的影響[J]. 中國農(nóng)業(yè)氣象,2005,26(1):7-10.Xie Junfei, Li Yu¢e. Effect of soil temperature on N2O emission in upland farm of Beijing[J]. Chinese Journal of Agrometeorology, 2005, 26(1): 7-10. (in Chinese with English abstract)

    [24]Granli T, Bockman O C. Nitrous oxide production in aerobic soils under varying pH, temperature and water content[J].Soil Biology & Biochemistry. 1984. 16: 39-43.

    [25]齊玉春,董云社. 土壤氧化亞氮產(chǎn)生、排放及其影響因素[J]. 地理學(xué)報(bào),1999,54(6):534-542.Qi Yuchun, Dong Yunshe. Nitrous oxide emissions from soil and some influence factors[J]. Acta Geographicasinica, 1999,54(6): 534-542. (in Chinese with English abstract)

    亚洲真实伦在线观看| 搡女人真爽免费视频火全软件| 亚洲,欧美,日韩| 九色成人免费人妻av| 国产黄a三级三级三级人| 精品人妻视频免费看| 成人亚洲欧美一区二区av| 国产伦在线观看视频一区| 国产午夜精品论理片| 中文资源天堂在线| 国产精品乱码一区二三区的特点| 亚洲色图av天堂| 国产熟女欧美一区二区| ponron亚洲| 97人妻精品一区二区三区麻豆| av免费观看日本| 国产精品久久久久久精品电影| 日本猛色少妇xxxxx猛交久久| 国产高清国产精品国产三级 | 乱码一卡2卡4卡精品| 久久鲁丝午夜福利片| 婷婷色麻豆天堂久久 | 中国美白少妇内射xxxbb| 成人毛片a级毛片在线播放| 免费不卡的大黄色大毛片视频在线观看 | 午夜日本视频在线| 亚洲av免费在线观看| 午夜老司机福利剧场| 久久久久久久久久久丰满| 天天一区二区日本电影三级| 只有这里有精品99| 美女xxoo啪啪120秒动态图| 亚洲成人久久爱视频| 精品国内亚洲2022精品成人| 一级毛片我不卡| 欧美丝袜亚洲另类| 日韩在线高清观看一区二区三区| 蜜臀久久99精品久久宅男| a级毛片免费高清观看在线播放| 99在线视频只有这里精品首页| 国产一区有黄有色的免费视频 | 成人毛片60女人毛片免费| 最近最新中文字幕免费大全7| 99久久无色码亚洲精品果冻| 久久精品国产亚洲av天美| 亚洲国产精品合色在线| 亚洲精品成人久久久久久| 听说在线观看完整版免费高清| 最近的中文字幕免费完整| 男女那种视频在线观看| 三级国产精品欧美在线观看| 日韩大片免费观看网站 | 国产免费男女视频| 汤姆久久久久久久影院中文字幕 | 日日摸夜夜添夜夜爱| 干丝袜人妻中文字幕| 精品久久久久久久人妻蜜臀av| 久久热精品热| 成人漫画全彩无遮挡| 久久精品国产亚洲网站| 亚洲在久久综合| 日韩欧美精品v在线| 综合色丁香网| 欧美另类亚洲清纯唯美| 成人亚洲欧美一区二区av| 99热网站在线观看| 国产精品久久久久久精品电影小说 | 美女大奶头视频| 天美传媒精品一区二区| 国产在线一区二区三区精 | 免费大片18禁| 内地一区二区视频在线| 精品不卡国产一区二区三区| 长腿黑丝高跟| 久久99热6这里只有精品| av在线亚洲专区| 国产真实伦视频高清在线观看| 一级毛片久久久久久久久女| av线在线观看网站| 国产极品精品免费视频能看的| 爱豆传媒免费全集在线观看| 国产中年淑女户外野战色| 人人妻人人看人人澡| av在线观看视频网站免费| 深爱激情五月婷婷| 日产精品乱码卡一卡2卡三| 国产成人aa在线观看| 国产乱来视频区| 精品人妻一区二区三区麻豆| 看免费成人av毛片| 神马国产精品三级电影在线观看| 久久99蜜桃精品久久| 少妇人妻一区二区三区视频| 菩萨蛮人人尽说江南好唐韦庄 | 日韩高清综合在线| 国产精品伦人一区二区| 久久这里有精品视频免费| 在线观看66精品国产| 97人妻精品一区二区三区麻豆| 九色成人免费人妻av| 亚洲五月天丁香| 麻豆乱淫一区二区| 国产亚洲av片在线观看秒播厂 | 丝袜喷水一区| 内射极品少妇av片p| 国产免费一级a男人的天堂| 久久精品夜夜夜夜夜久久蜜豆| 激情 狠狠 欧美| 亚洲国产最新在线播放| 99久国产av精品| 99久久九九国产精品国产免费| 欧美bdsm另类| 大香蕉97超碰在线| 欧美日本视频| 日本色播在线视频| 亚洲四区av| av专区在线播放| 日韩精品青青久久久久久| 最近视频中文字幕2019在线8| 国产69精品久久久久777片| 国产成人freesex在线| 亚洲欧美成人综合另类久久久 | 精品人妻偷拍中文字幕| 久久精品久久精品一区二区三区| 成人美女网站在线观看视频| av播播在线观看一区| 两个人视频免费观看高清| 亚洲美女搞黄在线观看| 午夜福利高清视频| 亚洲av熟女| 亚洲国产精品国产精品| 久久久午夜欧美精品| 综合色丁香网| 亚洲人成网站高清观看| 国产精品福利在线免费观看| 禁无遮挡网站| 亚洲国产日韩欧美精品在线观看| 免费av毛片视频| 1024手机看黄色片| 岛国毛片在线播放| 色哟哟·www| 青春草国产在线视频| 国产探花极品一区二区| 国产精品一及| 人妻系列 视频| 亚洲,欧美,日韩| 日本色播在线视频| 免费观看的影片在线观看| 啦啦啦韩国在线观看视频| 亚洲精品久久久久久婷婷小说 | 天堂网av新在线| 亚洲精品亚洲一区二区| 久久精品久久精品一区二区三区| 中文在线观看免费www的网站| 我的老师免费观看完整版| 青春草亚洲视频在线观看| 国产熟女欧美一区二区| 麻豆国产97在线/欧美| 免费看av在线观看网站| 亚洲第一区二区三区不卡| 丰满乱子伦码专区| 国产在线一区二区三区精 | 嫩草影院入口| 人人妻人人看人人澡| 26uuu在线亚洲综合色| 99久久精品国产国产毛片| 国产成人午夜福利电影在线观看| 亚洲综合色惰| 久久久国产成人精品二区| 亚洲av二区三区四区| 一边亲一边摸免费视频| 欧美另类亚洲清纯唯美| 国产乱人视频| 能在线免费看毛片的网站| 麻豆一二三区av精品| 亚洲真实伦在线观看| 卡戴珊不雅视频在线播放| 日韩欧美国产在线观看| 免费看a级黄色片| 18禁在线无遮挡免费观看视频| 啦啦啦韩国在线观看视频| 国产一区二区三区av在线| 一区二区三区免费毛片| 尾随美女入室| av又黄又爽大尺度在线免费看 | 国内精品美女久久久久久| 亚洲av电影在线观看一区二区三区 | 午夜亚洲福利在线播放| 久久久精品大字幕| 久久精品夜夜夜夜夜久久蜜豆| 精品一区二区三区视频在线| 内射极品少妇av片p| 国产精品一及| 亚洲精品456在线播放app| 美女脱内裤让男人舔精品视频| 成人特级av手机在线观看| 亚洲国产成人一精品久久久| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| 国产av一区在线观看免费| 久久久亚洲精品成人影院| 99久久无色码亚洲精品果冻| av专区在线播放| 最新中文字幕久久久久| 国产淫片久久久久久久久| 一个人看视频在线观看www免费| 99久久中文字幕三级久久日本| 国产中年淑女户外野战色| 欧美最新免费一区二区三区| 深爱激情五月婷婷| 成人一区二区视频在线观看| 亚洲综合色惰| 国产视频内射| 一区二区三区高清视频在线| 97热精品久久久久久| 18+在线观看网站| 亚洲美女视频黄频| 人妻少妇偷人精品九色| 一级黄色大片毛片| 我要搜黄色片| 看片在线看免费视频| 26uuu在线亚洲综合色| 能在线免费观看的黄片| av.在线天堂| 少妇裸体淫交视频免费看高清| 日本一本二区三区精品| av国产免费在线观看| 老司机影院成人| 欧美激情久久久久久爽电影| 欧美性感艳星| 99久久精品一区二区三区| 男人舔女人下体高潮全视频| 99热这里只有是精品在线观看| 你懂的网址亚洲精品在线观看 | 亚洲经典国产精华液单| 听说在线观看完整版免费高清| 老师上课跳d突然被开到最大视频| 国产av在哪里看| 男人的好看免费观看在线视频| 特级一级黄色大片| 狂野欧美激情性xxxx在线观看| 国产综合懂色| 九九久久精品国产亚洲av麻豆| 秋霞在线观看毛片| 午夜福利在线观看免费完整高清在| 一夜夜www| 99热这里只有是精品50| 免费观看精品视频网站| 看黄色毛片网站| 老司机影院毛片| 日韩中字成人| 成人午夜高清在线视频| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 久久久成人免费电影| 美女cb高潮喷水在线观看| 啦啦啦观看免费观看视频高清| 精品人妻偷拍中文字幕| 久久久精品欧美日韩精品| 蜜桃亚洲精品一区二区三区| 成年版毛片免费区| 高清日韩中文字幕在线| 亚洲欧美日韩高清专用| 国产精品电影一区二区三区| 在线观看av片永久免费下载| 成人漫画全彩无遮挡| 亚洲图色成人| 少妇被粗大猛烈的视频| 女人被狂操c到高潮| 成人毛片60女人毛片免费| 欧美xxxx黑人xx丫x性爽| 国产伦理片在线播放av一区| 亚洲精品国产成人久久av| 十八禁国产超污无遮挡网站| 91av网一区二区| 亚洲欧美一区二区三区国产| 麻豆久久精品国产亚洲av| 欧美成人a在线观看| 亚洲性久久影院| 麻豆国产97在线/欧美| 国产 一区 欧美 日韩| 免费观看a级毛片全部| 午夜激情欧美在线| 国产乱人视频| 熟女人妻精品中文字幕| 大话2 男鬼变身卡| 亚洲欧美清纯卡通| 美女脱内裤让男人舔精品视频| 成人午夜高清在线视频| 国产探花极品一区二区| 观看美女的网站| 亚洲人成网站在线观看播放| 国产国拍精品亚洲av在线观看| 亚洲图色成人| 婷婷六月久久综合丁香| 18禁动态无遮挡网站| 久久99精品国语久久久| 免费黄色在线免费观看| 1000部很黄的大片| 禁无遮挡网站| 成年女人看的毛片在线观看| 亚洲va在线va天堂va国产| 人妻制服诱惑在线中文字幕| 国产黄色小视频在线观看| 青春草视频在线免费观看| 少妇被粗大猛烈的视频| 99久国产av精品国产电影| 中文字幕制服av| 久久久国产成人免费| 久久这里有精品视频免费| 久久久成人免费电影| 最近的中文字幕免费完整| 亚洲在久久综合| 国产亚洲91精品色在线| 2021少妇久久久久久久久久久| 精品人妻一区二区三区麻豆| 在线a可以看的网站| 丰满人妻一区二区三区视频av| 大香蕉97超碰在线| av专区在线播放| 久久99蜜桃精品久久| 色播亚洲综合网| 精品人妻一区二区三区麻豆| 久久久久久久久久久免费av| 久久这里只有精品中国| 亚洲精品国产av成人精品| 免费观看精品视频网站| 男人舔奶头视频| 日韩强制内射视频| 精品99又大又爽又粗少妇毛片| 欧美精品一区二区大全| 国产精品熟女久久久久浪| 色网站视频免费| 91在线精品国自产拍蜜月| 欧美日本视频| 狂野欧美激情性xxxx在线观看| 国产欧美日韩精品一区二区| 国产男人的电影天堂91| 国产免费又黄又爽又色| 免费不卡的大黄色大毛片视频在线观看 | 人妻少妇偷人精品九色| 国产高潮美女av| 国产亚洲av嫩草精品影院| 亚洲成色77777| 搡老妇女老女人老熟妇| 久久鲁丝午夜福利片| 国产精品精品国产色婷婷| 国产成人午夜福利电影在线观看| 舔av片在线| 午夜免费男女啪啪视频观看| 久久久精品欧美日韩精品| 亚洲国产最新在线播放| 搡女人真爽免费视频火全软件| av卡一久久| 久久久久久久亚洲中文字幕| 亚洲欧美日韩无卡精品| 大香蕉97超碰在线| 伦精品一区二区三区| 美女大奶头视频| 久久韩国三级中文字幕| 99热精品在线国产| 一区二区三区高清视频在线| 丰满乱子伦码专区| 久久草成人影院| 夜夜看夜夜爽夜夜摸| 91久久精品国产一区二区三区| 老司机福利观看| 在线天堂最新版资源| 久久精品综合一区二区三区| 国产精品蜜桃在线观看| 午夜福利在线观看免费完整高清在| 免费看美女性在线毛片视频| 亚洲美女搞黄在线观看| 身体一侧抽搐| 91久久精品国产一区二区三区| 国产乱人偷精品视频| 中文天堂在线官网| 日日摸夜夜添夜夜爱| 精品久久久噜噜| 精华霜和精华液先用哪个| 亚洲av免费高清在线观看| 九九在线视频观看精品| 99久久精品国产国产毛片| 成年女人看的毛片在线观看| 好男人在线观看高清免费视频| 麻豆久久精品国产亚洲av| 国产精品嫩草影院av在线观看| av在线蜜桃| 久久久久网色| 22中文网久久字幕| 国产精品一二三区在线看| av在线观看视频网站免费| 啦啦啦啦在线视频资源| 精华霜和精华液先用哪个| 午夜精品国产一区二区电影 | 国产免费又黄又爽又色| 国产黄a三级三级三级人| 免费人成在线观看视频色| 99在线人妻在线中文字幕| 啦啦啦韩国在线观看视频| 欧美一区二区国产精品久久精品| 在线播放无遮挡| 偷拍熟女少妇极品色| 麻豆成人av视频| 国产精品精品国产色婷婷| 搞女人的毛片| 一区二区三区四区激情视频| 欧美日韩综合久久久久久| 一区二区三区乱码不卡18| 精品熟女少妇av免费看| 精品一区二区三区视频在线| 汤姆久久久久久久影院中文字幕 | 精品人妻一区二区三区麻豆| 黑人高潮一二区| 成人一区二区视频在线观看| 色综合站精品国产| 亚洲,欧美,日韩| 国产黄片视频在线免费观看| 国语自产精品视频在线第100页| av卡一久久| 你懂的网址亚洲精品在线观看 | 人人妻人人澡欧美一区二区| 又爽又黄无遮挡网站| 亚洲成人中文字幕在线播放| 别揉我奶头 嗯啊视频| 美女黄网站色视频| 成人综合一区亚洲| 中国国产av一级| 91狼人影院| 狠狠狠狠99中文字幕| 国产黄片美女视频| 插阴视频在线观看视频| 一区二区三区高清视频在线| 免费在线观看成人毛片| 国产一区二区在线观看日韩| 国产精品永久免费网站| 亚洲18禁久久av| 欧美97在线视频| 亚洲最大成人av| 亚洲,欧美,日韩| 男女视频在线观看网站免费| 亚洲电影在线观看av| 精品人妻熟女av久视频| 久久99蜜桃精品久久| 中文字幕免费在线视频6| www.色视频.com| 丝袜喷水一区| 精品久久国产蜜桃| 久久久久久久久久成人| 日日摸夜夜添夜夜爱| 一个人看的www免费观看视频| 亚洲va在线va天堂va国产| 少妇的逼水好多| 亚洲av福利一区| 美女被艹到高潮喷水动态| 99九九线精品视频在线观看视频| 国产精品99久久久久久久久| 一边摸一边抽搐一进一小说| 国产女主播在线喷水免费视频网站 | 一级黄片播放器| 中文欧美无线码| 亚洲国产精品久久男人天堂| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合| 最近中文字幕高清免费大全6| 亚洲av二区三区四区| 美女cb高潮喷水在线观看| 91午夜精品亚洲一区二区三区| 精品无人区乱码1区二区| 亚洲国产精品sss在线观看| 内地一区二区视频在线| 日韩大片免费观看网站 | 日韩三级伦理在线观看| 国语对白做爰xxxⅹ性视频网站| 天堂中文最新版在线下载 | 狂野欧美激情性xxxx在线观看| 搞女人的毛片| 亚洲高清免费不卡视频| 在线a可以看的网站| 国产精品熟女久久久久浪| 欧美一区二区精品小视频在线| 女人被狂操c到高潮| 久久精品国产自在天天线| 成年av动漫网址| 亚洲av福利一区| 国产在线男女| 观看美女的网站| 美女脱内裤让男人舔精品视频| 国产精品嫩草影院av在线观看| 日产精品乱码卡一卡2卡三| 老司机福利观看| 精品久久久久久电影网 | 99九九线精品视频在线观看视频| 三级经典国产精品| 亚洲四区av| 日韩人妻高清精品专区| 久久久久久久国产电影| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 少妇裸体淫交视频免费看高清| 亚洲色图av天堂| 男插女下体视频免费在线播放| 1000部很黄的大片| 国产v大片淫在线免费观看| 久久99热这里只有精品18| 久久人妻av系列| 日韩强制内射视频| 久久久久久久久久久免费av| 一本一本综合久久| 午夜福利在线观看免费完整高清在| 丰满少妇做爰视频| .国产精品久久| 亚洲欧美清纯卡通| 最近中文字幕2019免费版| 91久久精品国产一区二区成人| 看黄色毛片网站| 一边亲一边摸免费视频| 欧美97在线视频| 久久久久久国产a免费观看| 久久久精品大字幕| 精品不卡国产一区二区三区| 美女黄网站色视频| 久久精品国产亚洲av天美| 精品人妻视频免费看| 麻豆精品久久久久久蜜桃| 国内少妇人妻偷人精品xxx网站| 免费观看的影片在线观看| 国产一级毛片七仙女欲春2| 一本久久精品| 免费av不卡在线播放| 亚洲aⅴ乱码一区二区在线播放| 日日摸夜夜添夜夜添av毛片| 欧美性猛交╳xxx乱大交人| 欧美日本视频| 国产视频首页在线观看| 日韩一区二区视频免费看| 97超视频在线观看视频| 天堂影院成人在线观看| 国产成人免费观看mmmm| 久久精品国产鲁丝片午夜精品| 26uuu在线亚洲综合色| 蜜桃亚洲精品一区二区三区| 91av网一区二区| 晚上一个人看的免费电影| 深夜a级毛片| 成年免费大片在线观看| 日韩制服骚丝袜av| 国产伦精品一区二区三区四那| 久久6这里有精品| 国产极品精品免费视频能看的| 国产精品伦人一区二区| 丰满人妻一区二区三区视频av| 免费av观看视频| 亚洲精品影视一区二区三区av| 边亲边吃奶的免费视频| 蜜桃亚洲精品一区二区三区| 国产精品一区二区三区四区久久| 国产午夜精品论理片| 亚洲最大成人手机在线| 国产精品国产三级专区第一集| 别揉我奶头 嗯啊视频| 97在线视频观看| 日韩,欧美,国产一区二区三区 | 夜夜看夜夜爽夜夜摸| 精品无人区乱码1区二区| 我的女老师完整版在线观看| 少妇的逼水好多| 久久韩国三级中文字幕| 欧美又色又爽又黄视频| 直男gayav资源| 欧美性猛交╳xxx乱大交人| 国产高清不卡午夜福利| 中文字幕精品亚洲无线码一区| 久久99热这里只有精品18| 午夜免费男女啪啪视频观看| 国产在视频线精品| 哪个播放器可以免费观看大片| 国产精品电影一区二区三区| 亚洲婷婷狠狠爱综合网| 国产精品福利在线免费观看| 国产成人免费观看mmmm| 国产精品久久久久久精品电影| 亚洲人成网站高清观看| 日韩强制内射视频| 网址你懂的国产日韩在线| 欧美成人a在线观看| av.在线天堂| 亚洲国产精品成人久久小说| 久久久久久九九精品二区国产| 麻豆av噜噜一区二区三区| 男女下面进入的视频免费午夜| 日本一本二区三区精品| 美女黄网站色视频| 成年av动漫网址| 久久人人爽人人爽人人片va| 亚洲无线观看免费| 亚洲av熟女| 亚洲欧美日韩东京热| 一本久久精品| 22中文网久久字幕| 国产一级毛片在线| 美女黄网站色视频| 99久久精品一区二区三区| 亚洲内射少妇av| 国产亚洲精品av在线| 99九九线精品视频在线观看视频| 女人久久www免费人成看片 | 三级经典国产精品| 人妻系列 视频| 国产一区亚洲一区在线观看| 免费无遮挡裸体视频| 秋霞在线观看毛片| 国产精品一区二区三区四区久久| 国产黄片美女视频| 国产精品一二三区在线看| 久久人人爽人人爽人人片va|