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

    華北平原農(nóng)田N2O排放通量的高頻動(dòng)態(tài)觀測(cè)*

    2018-02-08 02:05:02李發(fā)東張秋英古叢珂冷佩芳喬云峰黃勇彬施生錦
    關(guān)鍵詞:農(nóng)田降雨大氣

    李發(fā)東, 杜 錕, 張秋英, 古叢珂, 冷佩芳, 喬云峰, 朱 農(nóng), 郝 帥, 黃勇彬, 施生錦

    ?

    華北平原農(nóng)田N2O排放通量的高頻動(dòng)態(tài)觀測(cè)*

    李發(fā)東1,2,3, 杜 錕1,2,3, 張秋英4, 古叢珂1,2,3, 冷佩芳1,2,3, 喬云峰1,2,3, 朱 農(nóng)1,2,3, 郝 帥1,2,3, 黃勇彬4, 施生錦4

    (1. 中國(guó)科學(xué)院地理科學(xué)與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測(cè)與模擬重點(diǎn)實(shí)驗(yàn)室 北京 100101; 2. 中國(guó)科學(xué)院禹城綜合試驗(yàn)站 禹城 251200; 3. 中國(guó)科學(xué)院大學(xué)資源與環(huán)境學(xué)院 北京 100190; 4. 北京雨根科技有限公司 北京 100094)

    N2O是主要源自農(nóng)田的重要溫室氣體之一, 可破壞臭氧層而導(dǎo)致全球增溫。目前對(duì)N2O的原位高頻觀測(cè)尚不多。為完善N2O的觀測(cè)方法, 為華北地區(qū)N2O變化研究提供參考, 本研究以華北平原典型農(nóng)田為研究對(duì)象, 利用新型的N2O測(cè)定儀器TGA200A, 進(jìn)行實(shí)時(shí)、自動(dòng)、晝夜連續(xù)地觀測(cè)中國(guó)科學(xué)院禹城綜合試驗(yàn)站農(nóng)田大氣N2O的日動(dòng)態(tài)變化。本次觀測(cè)自2017年6月中旬玉米播種后開(kāi)始, 持續(xù)至2017年9月(8月份儀器調(diào)試)。結(jié)果顯示: 1)晴朗天氣下, 農(nóng)田大氣N2O呈現(xiàn)出夜晚(0:00—6:00、18:00—24:00)高(0.618~1.171 mg·m-3)、白天(6:00—18:00)低(0.526~1.145 mg·m-3)的趨勢(shì), 而白天高溫又促進(jìn)農(nóng)田N2O排放, 在午后15:00—17:00大多出現(xiàn)1次峰值, 表明溫度的促進(jìn)作用存在滯后性。2)降雨天氣下, 農(nóng)田N2O在適當(dāng)?shù)挠炅肯轮饾u增加(3 h內(nèi)增加0.033 mg·m-3), 且存在累積效應(yīng), 但過(guò)度淹水后N2O表現(xiàn)出逐漸降低的趨勢(shì)。3)大風(fēng)天氣下, N2O的濃度產(chǎn)生變化, 但規(guī)律并不明顯。研究結(jié)果表明, 利用TGA200A可以實(shí)現(xiàn)對(duì)溫室氣體N2O的實(shí)時(shí)、連續(xù)、動(dòng)態(tài)的自動(dòng)觀測(cè), 觀測(cè)結(jié)果具有較高可信度, 可以反映出當(dāng)前華北地區(qū)農(nóng)田N2O在不同環(huán)境要素(溫度、降水及大風(fēng))下的動(dòng)態(tài)變化趨勢(shì)。

    華北平原; 農(nóng)田生態(tài)系統(tǒng); N2O排放; 通量監(jiān)測(cè); 高頻監(jiān)測(cè)

    現(xiàn)今全球氣候變暖已成為主要的氣候變化趨勢(shì), 由此引發(fā)的溫室效應(yīng)等問(wèn)題已日益成為當(dāng)前研究熱點(diǎn)之一。公認(rèn)的產(chǎn)生溫室效應(yīng)氣體主要包括CO2、CH4、N2O, 其對(duì)溫室效應(yīng)的貢獻(xiàn)率分別為60%、15%、5%, 綜合貢獻(xiàn)率高達(dá)80%[1], 又以N2O的溫室效應(yīng)最大, 其增溫效應(yīng)最高可達(dá)CO2的270倍[2]??諝庵械腘2O與其他物質(zhì)反應(yīng), 對(duì)臭氧層產(chǎn)生危害, 可增強(qiáng)地球表面的紫外線輻射, 進(jìn)一步對(duì)地球表面增溫產(chǎn)生影響。因此, 對(duì)于N2O的觀測(cè)和研究對(duì)于摸清當(dāng)前植被-土壤-大氣系統(tǒng)中的N2O濃度, 解釋當(dāng)前大氣中N2O的變化機(jī)制并以此促進(jìn)N2O的減排具有重要意義。目前, 前人關(guān)于N2O排放研究表明, 土壤對(duì)大氣N2O的貢獻(xiàn)為80%~90%, 是重要的N2O排放源[3]。對(duì)于不同生態(tài)系統(tǒng)的研究表明, 植被究竟是源還是匯與植被類(lèi)型和地區(qū)有關(guān)[4], 且不同生態(tài)系統(tǒng)排放的N2O量存在差異[5-7]。目前, 生態(tài)系統(tǒng)中N2O的觀測(cè)主要以靜態(tài)箱-氣相色譜法較為常見(jiàn)[5,8-9], 但該方法的測(cè)定時(shí)間間隔相對(duì)較長(zhǎng), 實(shí)現(xiàn)實(shí)時(shí)連續(xù)測(cè)定的難度較大, 且難以捕捉到降雨、風(fēng)等異常天氣下的排放情況。

    農(nóng)田生態(tài)系統(tǒng)受人類(lèi)影響, 土壤中氮元素循環(huán)較快, N2O排放較大, 是公認(rèn)的重要的N2O排放源[10]。華北平原地區(qū)是我國(guó)重要的糧食產(chǎn)地, 耕地面積廣, 耕作強(qiáng)度大, 屬于典型的具有代表性的農(nóng)田生態(tài)系統(tǒng)。閆翠萍等[11]研究表明, 華北平原小麥玉米輪作農(nóng)田是N2O的排放源, 小麥季土壤排放N2O總量為0.83~1.80 kg(N2O-N)·hm-2, 玉米季土壤排放N2O總量1.21~2.07 kg(N2O-N)·hm-2, 土壤溫度升高和降水量增加以及干濕交替加劇均會(huì)造成N2O排放速率增加。同時(shí)在持續(xù)干燥和低溫的冬季不施氮處理觀測(cè)到土壤對(duì)N2O 的吸收現(xiàn)象[12]。然而, 這些研究多是利用靜態(tài)箱-氣相色譜法測(cè)定, 測(cè)定頻度多為每周1次, 每次測(cè)定2~3 h, 總排放量通常采用插值法計(jì)算。這可能導(dǎo)致總排放量的計(jì)算誤差, 且由于需人工測(cè)定, 異常天氣下測(cè)定難度較大, 進(jìn)一步影響了測(cè)定精度。

    TGA200A(Trace Gas Analyzer, Campbell Scientific Ltd., USA)儀器是一種新型的TGA系列的溫室氣體觀測(cè)儀器, 具有響應(yīng)時(shí)間短、測(cè)定迅速的特點(diǎn), 可以每秒10 Hz的速度對(duì)大氣中的N2O進(jìn)行實(shí)時(shí)分析測(cè)定, 該類(lèi)儀器目前在國(guó)內(nèi)應(yīng)用尚未見(jiàn)報(bào)道。因此, 本研究采用TGA200A儀器, 配以渦度相關(guān)儀器和技術(shù), 以位于華北平原的中國(guó)科學(xué)院禹城綜合試驗(yàn)站農(nóng)田為研究對(duì)象, 觀測(cè)農(nóng)田N2O排放的日變化特征, 并利用相關(guān)環(huán)境要素如溫度、風(fēng)速、濕度等, 對(duì)華北農(nóng)田玉米季N2O變化及其影響要素進(jìn)行研究, 以期完善N2O的觀測(cè)方法, 并為華北地區(qū)N2O變化研究提供參考。

    1 材料和方法

    1.1 試驗(yàn)地概況

    試驗(yàn)地位于中國(guó)科學(xué)院禹城綜合試驗(yàn)站(36°40′~37°12′N(xiāo), 116°22′~116°45′E), 地處山東省西北部、黃河下游。該地土壤為黃河沖積物母質(zhì)發(fā)育, 土壤類(lèi)型主要為潮土和鹽化潮土, 是華北平原地區(qū)的典型農(nóng)田生態(tài)系統(tǒng)。試驗(yàn)地屬溫帶半干旱季風(fēng)氣候區(qū), 年均氣溫13.4 ℃, 年均降水量567 mm, 夏季為主要降雨期; 多年無(wú)霜期220 d, 大于0 ℃積溫為4 951 ℃, 日照時(shí)間超過(guò)2 600 h[13], 水熱條件優(yōu)良, 利于作物種植和生長(zhǎng)。

    該地農(nóng)田采用冬小麥-夏玉米輪作的種植方式, 小麥?zhǔn)斋@后于6月10—15日使用機(jī)器播種玉米, 不作翻耕處理, 生長(zhǎng)期使用地下水適度灌溉1~2次, 僅7月底施肥1次。

    1.2 TGA200A測(cè)定系統(tǒng)

    于2016年5月初安裝TGA200A系統(tǒng), 主要組成部分包括支架、TGA分析測(cè)定主機(jī)箱(Campbell Scientific Ltd., USA, 體積約0.5 m3)、CR6 datalogger數(shù)采器(Campbell Scientific Ltd., USA)、抽真空高壓泵、電腦、N2O標(biāo)氣瓶。此外, 配以渦度相關(guān)技術(shù)(eddy covariance technique)的CAST-3(Campbell Scientific Ltd., USA)、Li-7500(Li-COR Ltd., USA), 為防止相互干擾, 兩者中軸線不在同一水平面上, 存在一定夾角。在CAST-3超聲裝置外側(cè)架設(shè)N2O的進(jìn)樣口, 高2.5 m(圖1)。

    TGA主機(jī)箱內(nèi)配置以電流和溫度控制的激光發(fā)射管及其相關(guān)測(cè)定裝置, 通過(guò)掃描和分析痕量氣體(trace gas)對(duì)激光能量的線性吸收波長(zhǎng), 實(shí)現(xiàn)對(duì)痕量氣體的監(jiān)測(cè)[14]。當(dāng)真空泵抽低壓后, 待測(cè)氣體從進(jìn)樣口進(jìn)入管道, 經(jīng)過(guò)過(guò)濾和干燥后, 進(jìn)入主機(jī)箱待測(cè)。與此同時(shí), 標(biāo)氣瓶的N2O標(biāo)準(zhǔn)氣體進(jìn)入主機(jī)箱(濃度已知), 利用激光技術(shù)同時(shí)測(cè)定樣品和標(biāo)準(zhǔn)氣體的N2O的吸光值, 以此測(cè)定樣品的N2O濃度值。設(shè)定測(cè)定頻次為10 Hz, 每秒內(nèi)各次N2O均值即為大氣實(shí)時(shí)N2O濃度。

    渦度相關(guān)技術(shù)通過(guò)測(cè)定和分析溫度、H2O、CO2等的脈動(dòng)與垂直風(fēng)速脈動(dòng)協(xié)方差, 對(duì)生物圈(森林、農(nóng)田等)與大氣間湍流通量進(jìn)行非破壞性的標(biāo)準(zhǔn)測(cè)定方法, 具有適用范圍廣、可信度高等特點(diǎn)[15-16]。利用渦度相關(guān)技術(shù)中的CAST-3、Li-7500設(shè)備測(cè)定大氣風(fēng)速、溫度、濕度, 相關(guān)測(cè)定和分析方法見(jiàn)文獻(xiàn)[17-18]。

    圖1 測(cè)定農(nóng)田生態(tài)系統(tǒng)N2O濃度的TGA200A測(cè)定系統(tǒng)示意圖

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

    對(duì)半小時(shí)內(nèi)N2O濃度取算術(shù)平均值, 得每半小時(shí)的大氣N2O濃度, 利用Microsoft Excel軟件對(duì)數(shù)據(jù)進(jìn)行計(jì)算和作圖。

    2 結(jié)果與分析

    2.1 晴朗天氣下農(nóng)田N2O日動(dòng)態(tài)變化

    2.1.1 玉米出苗期N2O日動(dòng)態(tài)變化

    6月10—15日玉米播種后1周內(nèi)屬于出苗期, 選擇6月18日(晴朗天氣)的數(shù)據(jù)進(jìn)行分析(圖2)發(fā)現(xiàn), N2O濃度在0.611~0.669 mg·m-3范圍內(nèi), 整體上呈現(xiàn)出白天(6:00—18:00, 下同)低(多數(shù)低于0.630 mg·m-3), 夜晚(0:00—6:00、18:00—24:00, 下同)高(多數(shù)高于0.630 mg·m-3)的特征。自凌晨開(kāi)始至日出階段, 農(nóng)田N2O隨大氣溫度降低而逐漸降低, 僅在日出前后短暫升高。N2O在8:30—18:00整體上表現(xiàn)出先升后降的趨勢(shì), 即高溫時(shí)段15:30出現(xiàn)一次峰值0.627 mg·m-3; 在18:00前后逐漸升高, 至23:00保持平穩(wěn)。結(jié)果表明, 溫度對(duì)農(nóng)田N2O排放產(chǎn)生重要影響, 白天高溫可提高N2O濃度, 但這種作用存在一定的滯后性。大氣濕度和風(fēng)速對(duì)大氣溫度的變化表現(xiàn)出相應(yīng)的變化, 即溫度逐漸升高, 濕度隨后逐漸降低, 風(fēng)速在白天表現(xiàn)出升高的趨勢(shì), 進(jìn)一步佐證了環(huán)境要素對(duì)N2O排放的協(xié)同作用及其相互影響。

    圖2 玉米出苗期農(nóng)田N2O濃度及環(huán)境要素的日動(dòng)態(tài)變化

    NC: N2O濃度; AMC: 大氣含水量; T: 溫度; WS: 風(fēng)速。NC: nitrous oxide concentration; AMC:atmospheric water content; T: temperature; WS: wind speed.

    2.1.2 玉米三葉期N2O日動(dòng)態(tài)變化

    選擇玉米三葉期晴朗天氣的7月8日數(shù)據(jù)進(jìn)行分析。從圖3可知, 氣溫在日出前后最低, 下午14:00—15:00高溫, 甚至在17:00之前大氣溫度仍高于36 ℃; 濕度則隨溫度變化表現(xiàn)出相反的變化, 風(fēng)速在白天溫度較高時(shí)段達(dá)到較大水平。在該種溫度主導(dǎo)的環(huán)境要素下, N2O濃度呈現(xiàn)出夜晚高白天低的特征, 即在凌晨至8:00表現(xiàn)出降低趨勢(shì), 之后隨溫度升高N2O濃度逐漸升高, 但其降低的時(shí)間晚于苗期的6月18日, 這是因?yàn)樵?月份白天積溫較高, 相對(duì)高溫時(shí)段較長(zhǎng), 延遲了N2O降低的時(shí)間節(jié)點(diǎn)。

    東半壁店小流域內(nèi)垃圾主要為村莊、設(shè)施大棚區(qū)及工業(yè)企業(yè)區(qū)的生活、生產(chǎn)垃圾。隨著近年新農(nóng)村建設(shè)的大力投入,流域內(nèi)實(shí)施垃圾密閉化、無(wú)害化管理,對(duì)垃圾進(jìn)行分類(lèi)處理,生產(chǎn)垃圾統(tǒng)一運(yùn)到鎮(zhèn)垃圾填埋場(chǎng)回填;生活垃圾統(tǒng)一運(yùn)到鎮(zhèn)垃圾壓縮轉(zhuǎn)運(yùn)站進(jìn)行壓縮后運(yùn)往安定垃圾消納場(chǎng)消納。

    綜合分析以上6月18日和7月8日的數(shù)據(jù)可以發(fā)現(xiàn), 兩個(gè)時(shí)間的N2O濃度變化特點(diǎn)進(jìn)一步驗(yàn)證了溫度對(duì)植被-土壤系統(tǒng)N2O濃度存在重要影響, 且N2O濃度和溫度變化的節(jié)點(diǎn)不同, 溫度對(duì)N2O的作用存在滯后作用, 該種作用隨季節(jié)(高溫時(shí)長(zhǎng))而存在差異。

    圖3 玉米三葉期農(nóng)田N2O濃度及環(huán)境要素日動(dòng)態(tài)變化

    NC: N2O濃度; AMC: 大氣含水量; T: 溫度; WS: 風(fēng)速。NC: nitrous oxide concentration; AMC:atmospheric water content; T: temperature“WS: wind speed”.

    2.1.3 玉米花粒期N2O日動(dòng)態(tài)變化

    選擇玉米花粒期9月9日的數(shù)據(jù)進(jìn)行分析發(fā)現(xiàn)(圖4), 溫度、濕度、風(fēng)速均表現(xiàn)出白天高、夜晚低的特征。濕度在溫度最高的階段呈降低的趨勢(shì), 而風(fēng)速與溫度的變化存在相似趨勢(shì)。N2O在日出前后達(dá)到1.171 mg·m-3, 遠(yuǎn)高于玉米生長(zhǎng)前期平均水平, 而白天的N2O濃度較低, 其均值低于0.600 mg·m-3, 遠(yuǎn)低于相對(duì)應(yīng)的其他月份的濃度。

    圖4 玉米花粒期農(nóng)田N2O濃度及環(huán)境要素日動(dòng)態(tài)變化

    NC: N2O濃度; AMC: 大氣含水量; T: 溫度; WS: 風(fēng)速。NC: nitrous oxide concentration; AMC:atmospheric water content; T: temperature; WS: wind speed.

    2.2 降雨天氣下農(nóng)田N2O日動(dòng)態(tài)變化

    選擇降雨天氣(6月22日)的數(shù)據(jù)進(jìn)行分析發(fā)現(xiàn)(圖5), 降雨天氣下, 大氣溫度的全天變化較小, 為20~25 ℃, 顯著低于晴天(6月18日)的溫度(圖2); 濕度隨降雨時(shí)長(zhǎng)累積逐漸增大, 風(fēng)速在白天時(shí)段基本保持平穩(wěn), 夜間風(fēng)速逐漸降低。該種環(huán)境下, 農(nóng)田N2O表現(xiàn)出隨降雨累積逐漸增大的趨勢(shì)(9:00—12:30), 但在12:30過(guò)后表現(xiàn)出短時(shí)間降低的趨勢(shì), 之后急劇上升至夜間最高值0.762 mg·m-3。降雨(6月22日)與晴天(6月18日)的N2O濃度比較, 凌晨至日出階段以晴天較高, 而白天至深夜均以降雨天氣較高, 且受降雨影響, 夜晚階段的最高值遠(yuǎn)高于晴朗天氣, 表明降雨可促進(jìn)植被-土壤系統(tǒng)釋放N2O, 增加大氣N2O濃度, 但過(guò)度降雨抑制N2O排放。

    圖5 雨天玉米農(nóng)田N2O濃度及環(huán)境要素日動(dòng)態(tài)變化

    NC: N2O濃度; AMC: 大氣含水量; T: 溫度; WS: 風(fēng)速。NC: nitrous oxide concentration; AMC:atmospheric water content; T: temperature; WS: wind speed.

    2.3 大風(fēng)天氣下農(nóng)田N2O日動(dòng)態(tài)變化

    晴朗強(qiáng)風(fēng)天氣7月8日(圖3)和晴朗弱風(fēng)7月3日(圖6)數(shù)據(jù)進(jìn)行分析比較發(fā)現(xiàn), 兩者的溫度范圍較為相近, 濕度差距較小, 但風(fēng)速前者是后者3~20倍。對(duì)比N2O數(shù)據(jù)發(fā)現(xiàn), 風(fēng)速較低的7月3日的N2O濃度顯著高于風(fēng)速較高的7月8日。7月3日的N2O濃度隨時(shí)間呈先升后降趨勢(shì), 最高點(diǎn)出現(xiàn)在6:00前后; 而7月8日則表現(xiàn)出隨時(shí)間逐漸降低趨勢(shì)。凌晨至日出前后的溫度和風(fēng)速均以7月3日較低, 但其N(xiāo)2O濃度反而相對(duì)較高。8:00—12:00, 兩者溫度極度相近, 但濕度以7月3日相對(duì)較低, N2O反而較為相近。推測(cè)風(fēng)速對(duì)N2O的影響可能不是主導(dǎo)因素, 是否有其他因素的交互作用影響, 還需要進(jìn)一步研究。

    圖6 晴朗弱風(fēng)天氣玉米農(nóng)田N2O濃度及其環(huán)境要素日動(dòng)態(tài)變化

    NC: N2O濃度; AMC: 大氣含水量; T: 溫度; WS: 風(fēng)速。NC: nitrous oxide concentration; AMC:atmospheric water content; T: temperature; WS: wind speed.

    3 討論

    3.1 溫度對(duì)農(nóng)田N2O排放的影響

    6月中旬玉米播種后, TGA對(duì)農(nóng)田大氣6—9月份(8月中下旬進(jìn)行了相關(guān)調(diào)試)N2O濃度進(jìn)行了實(shí)時(shí)連續(xù)的觀測(cè), TGA監(jiān)測(cè)功能基本保持穩(wěn)定。分析晴朗天氣下玉米生長(zhǎng)期的N2O濃度及其環(huán)境要素發(fā)現(xiàn), 白天農(nóng)田N2O變化具有一定的規(guī)律性, 與葉欣等[19]和鄭循華等[20]研究結(jié)論相近。日出之后, 日照增強(qiáng)導(dǎo)致大氣溫度和植被-土壤系統(tǒng)溫度升高, 植物產(chǎn)生更多O2進(jìn)入到土壤中, 土壤參與N循環(huán)的相關(guān)微生物活性提高, 氮元素有機(jī)物的礦化作用增強(qiáng), 土壤排放到大氣中的N2O增多。中午過(guò)后, 溫度到達(dá)峰值后逐漸下降, 植物產(chǎn)生O2降低, 土壤含氧量減少, 土壤氮循環(huán)相關(guān)微生物活性降低, 土壤向大氣中排放N2O也隨之降低[20]。杜睿等[21]利用原狀土開(kāi)展相關(guān)增溫試驗(yàn), 得出了溫度顯著影響土壤N2O釋放的結(jié)論; 自然環(huán)境下進(jìn)行的農(nóng)田增溫試驗(yàn)同樣得出溫度升高可促進(jìn)土壤N2O的結(jié)論[13,22]。同時(shí), 對(duì)比不同月份發(fā)現(xiàn), 下午高溫時(shí)長(zhǎng)增加(高于30 ℃)進(jìn)一步促進(jìn)農(nóng)田N2O的釋放, 這有可能是高溫累積效應(yīng)導(dǎo)致的。

    3.2 降雨天氣對(duì)農(nóng)田大氣N2O排放的影響

    降雨引起農(nóng)田土壤含水量加大, 提高了土壤氮循環(huán)相關(guān)微生物活性, 同時(shí)土壤濕度的增大導(dǎo)致進(jìn)入土壤中的O2降低, 氮循環(huán)過(guò)程產(chǎn)生的N2O作為替代O2進(jìn)行還原作用, 因此硝化過(guò)程產(chǎn)N2O增多, 排放到大氣的N2O進(jìn)而增大[24]。隨著降雨的進(jìn)行, 土壤含水量逐漸增大, 形成一定的厭氧環(huán)境, 促進(jìn)反硝化作用[25], 表現(xiàn)出N2O在降雨發(fā)生后一定時(shí)間持續(xù)增長(zhǎng), 因此適當(dāng)?shù)慕涤昕纱龠M(jìn)土壤N2O排放[26], 這也與灌溉促進(jìn)N2O釋放具有相似的機(jī)制[27-28]。但白天過(guò)量的降水和灌溉導(dǎo)致土壤過(guò)度淹水, 影響土壤O2進(jìn)入和N2O排放, 導(dǎo)致N2O隨降水的過(guò)量累積表現(xiàn)出降低趨勢(shì)[5]。因此, 在降水和灌溉影響下土壤含水量增大的快慢和變化梯度對(duì)土壤N2O排放存在影響。

    3.3 大風(fēng)天氣對(duì)農(nóng)田N2O排放的影響

    大氣的水熱條件與農(nóng)田風(fēng)力存在較為密切的相互影響, 難以在野外單獨(dú)控制溫度和含水量(降水)條件以此研究風(fēng)力對(duì)農(nóng)田產(chǎn)生N2O的影響。本文僅選擇不同風(fēng)力下一定時(shí)段的N2O進(jìn)行對(duì)比分析, 發(fā)現(xiàn)大風(fēng)天氣可能降低N2O排放。相同溫度下, 一方面, 農(nóng)田風(fēng)速較大可促進(jìn)地表水分蒸發(fā), 降低表土含水量, 改變土壤與大氣的氣體(如O2、N2O)交換速率; 另一方面, 強(qiáng)風(fēng)導(dǎo)致植株葉片翻動(dòng), 葉片正面接受光照面積降低, 降低光合作用[29], 對(duì)溫室氣體產(chǎn)量造成影響, 而強(qiáng)風(fēng)甚至可引起植株受損, 導(dǎo)致農(nóng)田排放的溫室氣體產(chǎn)生變化。

    3.4 TGA監(jiān)測(cè)系統(tǒng)的可利用性

    目前, 國(guó)外少數(shù)學(xué)者利用TGA對(duì)生物圈與大氣間溫室氣體的CO2中13C同位素[30]、CH4[31]、NO2[32]進(jìn)行了一定的研究, 得出了TGA的應(yīng)用具有可行性的結(jié)論。本研究所采用的TGA為最新一代的系列產(chǎn)品, 不再使用液氮進(jìn)行主機(jī)箱降溫, 操作更加方便簡(jiǎn)單。本研究中晴朗天氣下, 白天時(shí)間TGA測(cè)定的晴朗天氣下農(nóng)田N2O排放量多數(shù)在0.60mg·m-3左右, 與相關(guān)研究結(jié)果相近[1-2]; 且本研究中晴朗天氣下, 大氣N2O變化趨勢(shì)與溫度變化趨勢(shì)相似, 在一定程度上可以相互印證, 表明應(yīng)用TGA系統(tǒng)監(jiān)測(cè)華北地區(qū)N2O濃度變化具有可行性。當(dāng)前, 農(nóng)田生態(tài)系統(tǒng)的溫室氣體排放監(jiān)測(cè)多是單一土壤或植被系統(tǒng)的研究, 直接在植被-土壤系統(tǒng)冠層高度觀測(cè)大氣N2O的研究極少, 本研究提供了一種在華北農(nóng)田觀測(cè)溫室氣體的新方法。由于人類(lèi)活動(dòng)影響, 農(nóng)田的下墊面相對(duì)單一, TGA系統(tǒng)可以和渦度相關(guān)技術(shù)集合利用, 但TGA系統(tǒng)技術(shù)在國(guó)內(nèi)森林、草原和山地等復(fù)雜狀況下的研究是否可行, 目前尚不明了。

    TGA系統(tǒng)主機(jī)箱內(nèi)部更換相關(guān)配件后, 可用于測(cè)定農(nóng)田CO2、CH4。同時(shí), 對(duì)于進(jìn)樣口的改造, 可以對(duì)土壤或植被釋放的溫室氣體進(jìn)行自動(dòng)采集和測(cè)定, 具有較高的科研利用價(jià)值。但TGA造價(jià)高昂, 在全國(guó)乃至全世界的推廣具有一定難度。本文利用TGA觀測(cè)農(nóng)田大氣N2O僅僅是初步研究, 測(cè)定周期較短, 對(duì)于農(nóng)田N2O在長(zhǎng)期氣候變化下的響應(yīng)還不明了, 尚需開(kāi)展長(zhǎng)期觀測(cè)。

    4 結(jié)論

    1)晴朗天氣下, 溫度是影響玉米農(nóng)田N2O排放的主要環(huán)境因素, N2O排放呈現(xiàn)出夜晚高、白天低的特征, 溫度對(duì)N2O的影響存在滯后性。

    2)適量的降雨可促進(jìn)農(nóng)田N2O排放, 但過(guò)度降水抑制農(nóng)田N2O排放。

    3)大風(fēng)天氣影響下, 農(nóng)田N2O排放發(fā)生變化, 但該種變化不具有明顯規(guī)律性, 大風(fēng)天氣的作用機(jī)制尚需進(jìn)一步探究。

    4)本研究為國(guó)內(nèi)首次利用TGA系統(tǒng)對(duì)溫室氣體通量進(jìn)行原位高頻監(jiān)測(cè), 證實(shí)TGA系統(tǒng)測(cè)定華北農(nóng)田的N2O具有可行性, 可實(shí)現(xiàn)實(shí)時(shí)、動(dòng)態(tài)、晝夜連續(xù)監(jiān)測(cè)。

    [1] Climate Change 2014: Synthesis Report. Contribution of working groups Ⅰ, Ⅱ and Ⅲ to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[R]. Geneva, Switzerland: IPCC, 2014

    [2] IPCC. Climate Change 2007: The Physical Science Basis[M]. Cambridge: Cambridge University Press, 2007

    [3] MELILLO J M, STEUDLER P A, ABER J D, et al. Soil warming and carbon-cycle feedbacks to the climate system[J]. Science, 2002, 298(5601): 2173-2176

    [4] GIL J, PéREZ T, BOERING K, et al. Mechanisms responsible for high N2O emissions from subarctic permafrost peatlands studied via stable isotope techniques[J]. Global Biogeochemical Cycles, 2017, 31(1): 172-189

    [5] 趙光影, 劉景雙, 王洋. 陸地生態(tài)系統(tǒng)N2O排放源研究進(jìn)展[J]. 土壤通報(bào), 2008, 39(5): 1192-1197 ZHAO G Y, LIU J S, WANG Y. Research progresses on nitrous oxide emission in terrestrial ecosystem[J]. Journal of Soil Science, 2008, 39(5): 1192-1197

    [6] 韓琳, 王鴿, 王偉, 等. 全球森林土壤N2O排放通量的影響因子[J]. 生態(tài)學(xué)雜志, 2012, 31(2): 446-452 HAN L, WANG G, WANG W, et al. Factors affecting global forest soil N2O emission flux[J]. Chinese Journal of Ecology, 2012, 31(2): 446-452

    [7] LENKA N K, LAL R. Soil aggregation and greenhouse gas flux after 15 years of wheat straw and fertilizer management in a no-till system[J]. Soil and Tillage Research, 2013, 126: 78-89

    [8] 閻佩云, 劉建亮, 沈玉芳, 等. 黃土旱塬旱作覆膜春玉米農(nóng)田N2O排放通量及影響因素研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2013, 32(11): 2278-2285 YAN P Y, LIU J L, SHEN Y F, et al. Nitrous oxide emissions and its influencing factors in spring corn fields mulched with plastic film on the Loess Plateau[J]. Journal of Agriculture Environment Science, 2013, 32(11): 2278-2285

    [9] 花莉, 唐志剛, 解井坤, 等. 生物質(zhì)炭對(duì)農(nóng)田溫室氣體排放的作用效應(yīng)及其影響因素探討[J]. 生態(tài)環(huán)境學(xué)報(bào), 2013, 22(6): 1068-1073 HUA L, TANG Z G, XIE J K, et al. Effect and its influencing factors of biochar on agricultural greenhouse gases emissions[J]. Ecology and Environmental Sciences, 2013, 22(6): 1068-1073

    [10] KROEZE C, MOSIER A, BOUWMAN L. Closing the global N2O budget: a retrospective analysis 1500-1994[J]. Global Biogeochemical Cycles, 1999, 13(1): 1-8

    [11] 閆翠萍, 張玉銘, 胡春勝, 等. 不同耕作措施下小麥–玉米輪作農(nóng)田溫室氣體交換及其綜合增溫潛勢(shì)[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2016, 24(6): 704-715 YAN C P, ZHANG Y M, HU C S, et al. Greenhouse gas exchange and comprehensive global warming potential under different wheat-maize rotation patterns[J]. Chinese Journal of Eco-Agriculture, 2016, 24(6): 704-715

    [12] 王玉英, 胡春勝. 施氮水平對(duì)太行山前平原冬小麥-夏玉米輪作體系土壤溫室氣體通量的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2011, 19(5): 1122-1128 WANG Y Y, HU C S. Soil greenhouse gas emission in winter wheat/summer maize rotation ecosystem as affected by nitrogen fertilization in the Piedmont Plain of Mount Taihang, China[J].Chinese Journal of Eco-Agriculture, 2011, 19(5): 1122-1128

    [13] TU C, LI F D. Responses of greenhouse gas fluxes to experimental warming in wheat season under conventional tillage and no-tillage fields[J]. Journal of Environmental Sciences, 2017, 54(4): 314-327

    [14] BOWLING D R, SARGENT S D, TANNER B D, et al. Tunable diode laser absorption spectroscopy for stable isotope studies of ecosystem-atmosphere CO2exchange[J]. Agricultural & Forest Meteorology, 2003, 118(1/2): 1-19

    [15] 于貴瑞, 孫曉敏. 陸地生態(tài)系統(tǒng)通量觀測(cè)的原理與方法[M]. 北京: 高等教育出版社, 2006: 205-206 YU G R, SUN X M. Principle of Flux Measurement in Terrestrial Ecosystem[M]. Beijing: Higher Education Press, 2014, 119: 205-206

    [16] CAMPBELL Scientific, Inc. Open Path Eddy Covariance System Operator’s Manual Csat, Li7500, and KH20[G]. Campbell Scientific, Inc, 2009

    [17] BURBA G G, MCDERMITT D K, GRELLE A, et al. Addressing the influence of instrument surface heat exchange on the measurements of CO2flux from open-path gas analyzers[J]. Global Change Biology, 2010, 14(8): 1854-1876

    [18] OECHEL W C, LLAKOWSKI C A, BURBA G, et al. Annual patterns and budget of CO2flux in an arctic tussock tundra ecosystem[J]. Journal of Geophysical Research: Biogeosciences, 2014, 119: 323-339

    [19] 葉欣, 李俊, 王迎紅, 等. 華北平原典型農(nóng)田土壤氧化亞氮的排放特征[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2005, 24(6): 1186-1191 YE X, LI J, WANG Y H. Characterization of emissions of nitrous oxide from soils of typical crop fields in North China Plain[J]. Journal of Agriculture Environment Science, 2005, 24(6): 1186-1191

    [20] 鄭循華, 王明星, 王躍思, 等. 溫度對(duì)農(nóng)田N2O產(chǎn)生與排放的影響[J]. 環(huán)境科學(xué), 1997, 18(5): 1-6 ZHENG X H, WANG M X, WANG Y S, et al. Impacts of temperature on N2O production and emission[J]. Environmental Science, 1997, 18(5): 1-6

    [21] 杜睿, 呂達(dá)仁, 王庚辰, 等. 內(nèi)蒙古典型草原土壤N2O排放的影響[J]. 自然科學(xué)進(jìn)展, 2003, 13(1): 64-68 DU R, LYU D R, WANG G C. Study on the mechanism of N2O production in typical grassland soils of Inner Mongolia[J]. Progress in Natural Science, 2003, 13(1): 64-68

    [22] 劉艷. 模擬增溫對(duì)農(nóng)田土壤呼吸、硝化及反硝化作用的影響[D]. 南京: 南京信息工程大學(xué), 2013 LIU Y. Effects of simulated warming on soil nitrification and denitrification in a farmland[D]. Nanjing: Nanjing University of Information Science and Technology, 2013

    [23] 胡立峰. 不同耕法對(duì)麥玉兩熟及雙季稻農(nóng)田溫室氣體排放的影響[D]. 北京: 中國(guó)農(nóng)業(yè)大學(xué), 2006HU L F. Effects of different tillage on greenhouse gas emissions of wheat-maize double cropping and double rice cropping system[D]. Beijing: China Agricultural University, 2006

    [24] AZAM F, MüLLER C, WEISKE A, et al. Nitrification and denitrification as sources of atmospheric nitrous oxide —Role of oxidizable carbon and applied nitrogen[J]. Biology & Fertility of Soils, 2002, 35(1): 54-61

    [25] 孫志強(qiáng), 郝慶菊, 江長(zhǎng)勝, 等. 農(nóng)田土壤N2O的產(chǎn)生機(jī)制及其影響因素研究進(jìn)展[J]. 土壤通報(bào), 2010, 41(6): 1524-1530 SUN Z Q, HAO Q J, JIANG C S. Advances in the study of nitrous oxide production mechanism and its influencing factors in agricultural soils[J]. Journal of Soil Science, 2010, 41(6): 1524-1530

    [26] LAURA S M, ANA M, LOURDES G T, et al. Combination of drip irrigation and organic fertilizer for mitigating emissions of nitrogen oxides in semiarid climate[J]. Agriculture Ecosystems & Environment, 2010, 137(1): 99-107

    [27] 裴淑瑋, 張圓圓, 劉俊鋒, 等. 華北平原玉米-小麥輪作農(nóng)田N2O交換通量的研究[J]. 環(huán)境科學(xué), 2012, 33(10): 3641-3646 PEI S W, ZHANG Y Y, LIU J F. N2O exchange fluxes from wheat-maize crop rotation system in the North China Plain[J]. Environmental Science, 2012, 33(10): 3641-3646

    [28] FIERER N, SCHIMEL J P. Effects of drying-rewetting frequency on soil carbon and nitrogen transformations[J]. Soil Biology & Biochemistry, 2002, 34(6): 777-787

    [29] 郝再彬. 植物生理生化[M]. 哈爾濱: 哈爾濱出版社, 2002: 207-208 HAO Z B. Plant Physiology and Biochemistry[M]. Harbin: Harbin Press, 2002: 207-208

    [30] BARBOUR M M, EVANS J R, SIMONIN K A, et al. Online CO2and H2O oxygen isotope fractionation allows estimation of mesophyll conductance in C4plants, and reveals that mesophyll conductance decreases as leaves age in both C4and C3plants[J]. New Phytologist, 2016, 210(3): 875-889

    [31] LAURILA T, TUOVINEN J P, LOHILA A, et al. Measuring methane emissions from a landfill using a cost-effective micrometeorological method[J]. Geophysical Research Letters, 2005, 32(19): 312-321

    [32] RINNE J, PIHLATIE M, LOHILA A, et al. Nitrous oxide emissions from a municipal landfill[J]. Environmental Science & Technology, 2005, 39(20): 7790-7793

    李發(fā)東, 杜錕, 張秋英, 古叢珂, 冷佩芳, 喬云峰, 朱農(nóng), 郝帥, 黃勇彬, 施生錦. 華北平原農(nóng)田N2O排放通量的高頻動(dòng)態(tài)觀測(cè)[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2018, 26(2): 195-202

    LI F D, DU K, ZHANG Q Y, GU C K, LENG P F, QIAO Y F, ZHU N, HAO S, HUANG Y B, SHI S J. High-frequency dynamic observation of N2O emission flux from cropland in the North China Plain[J]. Chinese Journal of Eco-Agriculture, 2018, 26(2): 195-202

    High-frequency dynamic observation of N2O emission flux from cropland in the North China Plain*

    LI Fadong1,2,3, DU Kun1,2,3, ZHANG Qiuying4, GU Congke1,2,3, LENG Peifang1,2,3, QIAO Yunfeng1,2,3, ZHU Nong1,2,3, HAO Shuai1,2,3, HUANG Yongbin4, SHI Shengjin4

    (1. Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China; 2. Yucheng Comprehensive Experiment Station, Chinese Academy of Sciences, Yucheng 251200, China; 3. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100190, China; 4. Rainroot Scientific Limited, Beijing 100094, China)

    Nitrous oxide (N2O) is one of the primary greenhouse gases, which depletes ozone sphere and results in strong greenhouse effects. Nowadays, a series of studies on greenhouse gases emissions of vegetation-soil-atmosphere system has been carried out domestic and overseas with the measurement method of static chamber/gas chromatograph, in which gas sampling time is concentrated in a fixed period during daytime rather than around the clock. Moreover, its’ sampling frequency is usually longer than one minute inducing difficulty to implement a real-time, high-frequency and continuous determination of greenhouse gas emissionsfor a long time. As an important source of N2O emission, farmland ecosystem is a human-disturbed system with fast N recycle and large N2O emission. To improve N2O emission observation methods and provide a reference for N2O change study in the North China Plain, we chose a typical farmland ecosystem, a maize field at the Yucheng Comprehensive Experiment Station, Chinese Academy of Sciences, and monitored N2O concentration change with a new monitor instrument TGA200A (Trace Gas Analyzer 200A) to achieve a automatically real-time monitoring of N2O emission for day and night. The TGA200A was equipped with a laser launch tube simultaneously controlled by current and temperature and some related measuring devices. When sample and reference gases (with a known concentration) entered the analyzer synchronously, the target gas concentration was determined by through scanning and comparing the linear absorption wavelengths of laser energy between sample and reference gases. The monitoring was from the middle of June 2017 to September 2017. The results showed, firstly, the N2O emission was higher in night from 0:00 to 6:00 and 18:00 to 24:00 (0.618-1.171 mg·m-3) than in daytime from 6:00 to 18:00 (0.526-1.145 mg·m-3) in fine weather. N2O emission was facilitated under higher atmospheric temperature in daytime, but the emission peak appeared in 15:00 to 17:00 indicating a significant hysteresis of temperature effect. Secondly, in rainy day, maize field N2O emissions were increased by 0.033 mg·m-3in 3 hours with a proper rainfall and presented an accumulative effect. But if the rain was too heavy for a long time, the N2O emission would be reduced. Thirdly, N2O emission was affected by strong wind, but this conclusion needed more verification because such result was not regular. This study demonstrated the data determined by the TGA200A was useful and reliable for study on dynamic change of N2O emission in different weathers. TGA200A was available in an automatically and real-time monitor of N2O emission for day and night with different environment elements (temperature, rain, wind) and made it possible to reduce human costs and errors in greenhouse gas flux observation.

    North China Plain; Farmland ecosystem; Nitrous oxide emission; Flux observation; High-frequency observation

    , LI Fadong, E-mail: lifadong@igsnrr.ac.cn

    Dec. 19, 2017;

    Dec. 25, 2017

    10.13930/j.cnki.cjea.171209

    X51; S-3

    A

    1671-3990(2018)02-0195-08

    2017-12-19

    2017-12-25

    * This study was supported by the National Natural Science Foundation of China (41771292, 31170414), and the National Key Research and Development Project of China (2016YFD0800301).

    * 國(guó)家自然科學(xué)基金項(xiàng)目(41771292, 31170414)和國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFD0800301)資助

    李發(fā)東, 主要從事生態(tài)水文與水環(huán)境研究。E-mail: lifadong@igsnrr.ac.cn

    猜你喜歡
    農(nóng)田降雨大氣
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    達(dá)爾頓老伯的農(nóng)田
    滄州市2016年“7.19~7.22”與“8.24~8.25”降雨對(duì)比研究
    農(nóng)田創(chuàng)意秀
    大氣古樸揮灑自如
    大氣、水之后,土十條來(lái)了
    紅黏土降雨入滲的定量分析
    農(nóng)田搞養(yǎng)殖需辦哪些證
    農(nóng)田制作所
    南方降雨不斷主因厄爾尼諾
    免费看日本二区| 97超视频在线观看视频| 国产亚洲精品av在线| 高清午夜精品一区二区三区| 美女cb高潮喷水在线观看| 精品一区二区三区视频在线| 人人妻人人看人人澡| 免费观看无遮挡的男女| 久久久精品94久久精品| 亚洲第一区二区三区不卡| 欧美激情在线99| 午夜福利高清视频| 91久久精品国产一区二区三区| 国产黄片美女视频| av女优亚洲男人天堂| 国产精品麻豆人妻色哟哟久久 | 亚洲最大成人手机在线| 三级经典国产精品| 有码 亚洲区| 天天躁夜夜躁狠狠久久av| 免费观看在线日韩| 18禁在线播放成人免费| 十八禁国产超污无遮挡网站| 国产精品日韩av在线免费观看| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 国产精品人妻久久久久久| 免费在线观看成人毛片| 国产伦精品一区二区三区视频9| 国产男女超爽视频在线观看| 精品一区二区三区视频在线| 久久久久久国产a免费观看| 少妇人妻一区二区三区视频| 欧美激情久久久久久爽电影| 麻豆精品久久久久久蜜桃| 亚洲欧美精品专区久久| 国产精品国产三级国产专区5o| 啦啦啦啦在线视频资源| 晚上一个人看的免费电影| 免费看美女性在线毛片视频| 久久99蜜桃精品久久| 永久免费av网站大全| 久久鲁丝午夜福利片| 亚洲欧美日韩东京热| 国产乱人偷精品视频| 午夜福利视频1000在线观看| 69人妻影院| 国产精品爽爽va在线观看网站| 国产真实伦视频高清在线观看| 91久久精品电影网| 国产精品一区二区在线观看99 | 99热这里只有是精品在线观看| a级毛色黄片| 免费观看的影片在线观看| 性插视频无遮挡在线免费观看| 卡戴珊不雅视频在线播放| 卡戴珊不雅视频在线播放| 男女边吃奶边做爰视频| 久久精品夜夜夜夜夜久久蜜豆| 黄色配什么色好看| 国产成人精品一,二区| 中文字幕制服av| 国产精品久久久久久av不卡| 日韩三级伦理在线观看| 黄色一级大片看看| 午夜精品国产一区二区电影 | 午夜福利成人在线免费观看| 18禁在线无遮挡免费观看视频| 亚洲精品视频女| 国产久久久一区二区三区| 亚洲国产精品专区欧美| 国产成人91sexporn| 国内揄拍国产精品人妻在线| 国产伦精品一区二区三区视频9| 国产av不卡久久| 国产 一区 欧美 日韩| 毛片女人毛片| 丰满乱子伦码专区| 我的女老师完整版在线观看| 在线观看一区二区三区| 亚洲真实伦在线观看| 尤物成人国产欧美一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 国产精品美女特级片免费视频播放器| 一区二区三区四区激情视频| 亚洲最大成人av| 成人漫画全彩无遮挡| 日日摸夜夜添夜夜爱| 午夜激情久久久久久久| 国产精品无大码| 免费高清在线观看视频在线观看| 亚洲欧美精品自产自拍| 欧美日韩综合久久久久久| 国产毛片a区久久久久| 精品国产露脸久久av麻豆 | 国内精品一区二区在线观看| 中文字幕av在线有码专区| 成年av动漫网址| 99久国产av精品| 嫩草影院新地址| 一级片'在线观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲图色成人| 99热全是精品| 国产中年淑女户外野战色| 国语对白做爰xxxⅹ性视频网站| 久久这里只有精品中国| 婷婷色麻豆天堂久久| 日日摸夜夜添夜夜添av毛片| 麻豆乱淫一区二区| 高清av免费在线| 精品久久久久久电影网| 色网站视频免费| 观看免费一级毛片| 国产在视频线精品| 久久久久精品久久久久真实原创| 天堂av国产一区二区熟女人妻| 又爽又黄无遮挡网站| 欧美成人a在线观看| 69人妻影院| 精品久久久久久久久av| 国产亚洲91精品色在线| 九色成人免费人妻av| 日本wwww免费看| 搡老乐熟女国产| 搡老乐熟女国产| 九色成人免费人妻av| 亚洲av成人av| 国国产精品蜜臀av免费| 成人漫画全彩无遮挡| 91在线精品国自产拍蜜月| av卡一久久| 婷婷色综合www| 建设人人有责人人尽责人人享有的 | 男人舔女人下体高潮全视频| 亚洲欧美成人综合另类久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产永久视频网站| 菩萨蛮人人尽说江南好唐韦庄| 深爱激情五月婷婷| 最后的刺客免费高清国语| 狂野欧美白嫩少妇大欣赏| 国产亚洲av片在线观看秒播厂 | 亚洲av国产av综合av卡| 深夜a级毛片| 午夜激情久久久久久久| 中文精品一卡2卡3卡4更新| 一夜夜www| 99热全是精品| 十八禁国产超污无遮挡网站| 国内精品一区二区在线观看| 91久久精品电影网| 亚洲欧美日韩无卡精品| 久久久久久久久久黄片| 校园人妻丝袜中文字幕| 欧美成人一区二区免费高清观看| 亚洲欧美清纯卡通| 欧美成人精品欧美一级黄| 大又大粗又爽又黄少妇毛片口| 久久精品国产亚洲av涩爱| 最近中文字幕2019免费版| 少妇丰满av| 成人漫画全彩无遮挡| 在线免费十八禁| 国内精品宾馆在线| 男女视频在线观看网站免费| 亚洲欧美一区二区三区国产| av在线老鸭窝| 一个人看的www免费观看视频| 久久国产乱子免费精品| 一个人免费在线观看电影| 在线观看免费高清a一片| 日韩一区二区三区影片| 成人欧美大片| 亚洲精品亚洲一区二区| 精品久久久久久久久av| 熟女人妻精品中文字幕| 欧美潮喷喷水| 国产伦精品一区二区三区视频9| 精品一区二区三区视频在线| 亚洲av二区三区四区| 欧美精品国产亚洲| 狠狠精品人妻久久久久久综合| 寂寞人妻少妇视频99o| 舔av片在线| 免费黄网站久久成人精品| 两个人的视频大全免费| 91精品伊人久久大香线蕉| 欧美潮喷喷水| 国产精品一及| 天堂俺去俺来也www色官网 | 美女内射精品一级片tv| 国产伦理片在线播放av一区| 亚洲av福利一区| 青春草国产在线视频| 波野结衣二区三区在线| av在线天堂中文字幕| 99久久人妻综合| 男的添女的下面高潮视频| 国产av国产精品国产| 日韩电影二区| 我的女老师完整版在线观看| 女人被狂操c到高潮| 免费观看在线日韩| 成人午夜高清在线视频| 国产又色又爽无遮挡免| 久久久久久伊人网av| 99热全是精品| 纵有疾风起免费观看全集完整版 | 午夜视频国产福利| 国产精品人妻久久久影院| 麻豆av噜噜一区二区三区| 免费观看av网站的网址| 日韩av在线大香蕉| 在线播放无遮挡| 99久久中文字幕三级久久日本| 国产亚洲一区二区精品| 国产女主播在线喷水免费视频网站 | 国产色婷婷99| 精品一区二区免费观看| 国产单亲对白刺激| 久久久亚洲精品成人影院| 91久久精品国产一区二区成人| 美女被艹到高潮喷水动态| 少妇被粗大猛烈的视频| 久久99热6这里只有精品| 欧美日韩综合久久久久久| 我要看日韩黄色一级片| 中文字幕久久专区| 国产成人91sexporn| 色视频www国产| 肉色欧美久久久久久久蜜桃 | 干丝袜人妻中文字幕| 夜夜看夜夜爽夜夜摸| 观看美女的网站| 观看免费一级毛片| 国产伦精品一区二区三区视频9| 韩国av在线不卡| 国产精品久久久久久久电影| 91精品伊人久久大香线蕉| 国产午夜精品一二区理论片| 午夜激情久久久久久久| 丝袜喷水一区| 日韩制服骚丝袜av| 欧美 日韩 精品 国产| 国产精品一区二区性色av| 蜜桃亚洲精品一区二区三区| 一级毛片黄色毛片免费观看视频| 校园人妻丝袜中文字幕| 久久久欧美国产精品| 久久久久久久久大av| 99热这里只有是精品在线观看| 日本免费在线观看一区| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 精品国产三级普通话版| 国产精品国产三级国产专区5o| 国产熟女欧美一区二区| 日日干狠狠操夜夜爽| 国产伦精品一区二区三区四那| 亚洲欧美成人综合另类久久久| 久久人人爽人人片av| av.在线天堂| 日日摸夜夜添夜夜爱| 免费看不卡的av| 国产一区亚洲一区在线观看| 久久精品国产鲁丝片午夜精品| 汤姆久久久久久久影院中文字幕 | 国产 一区 欧美 日韩| 婷婷色综合大香蕉| 国产单亲对白刺激| 黄片无遮挡物在线观看| 国产视频内射| 久久久亚洲精品成人影院| 午夜免费激情av| 国产成人freesex在线| 亚洲欧洲日产国产| 成人美女网站在线观看视频| 色综合亚洲欧美另类图片| 嘟嘟电影网在线观看| av网站免费在线观看视频 | 少妇熟女aⅴ在线视频| 综合色av麻豆| 可以在线观看毛片的网站| 亚洲伊人久久精品综合| 精品99又大又爽又粗少妇毛片| 国产精品国产三级专区第一集| 中文乱码字字幕精品一区二区三区 | 精品久久国产蜜桃| 啦啦啦韩国在线观看视频| 亚洲av.av天堂| 一级毛片 在线播放| av在线播放精品| 中文字幕免费在线视频6| 国产午夜福利久久久久久| 欧美激情在线99| 91狼人影院| 丰满乱子伦码专区| 国产在线男女| 久久久精品94久久精品| 国产在线男女| 精品一区在线观看国产| 赤兔流量卡办理| 又粗又硬又长又爽又黄的视频| 国产久久久一区二区三区| 日本午夜av视频| 在线免费观看不下载黄p国产| 大话2 男鬼变身卡| 国产老妇伦熟女老妇高清| 在线a可以看的网站| 69人妻影院| 亚洲性久久影院| 国国产精品蜜臀av免费| 小蜜桃在线观看免费完整版高清| 男女国产视频网站| 欧美变态另类bdsm刘玥| 综合色av麻豆| kizo精华| 美女xxoo啪啪120秒动态图| 日韩一区二区视频免费看| 国产精品一区二区三区四区久久| 午夜精品一区二区三区免费看| 亚洲av.av天堂| 99久久人妻综合| 久久精品夜夜夜夜夜久久蜜豆| 嘟嘟电影网在线观看| 亚洲乱码一区二区免费版| 精品人妻熟女av久视频| 中文字幕av成人在线电影| 亚洲国产精品成人久久小说| 国内少妇人妻偷人精品xxx网站| 久久久色成人| 91久久精品电影网| 国产黄色视频一区二区在线观看| 久久精品国产自在天天线| 成人鲁丝片一二三区免费| 亚洲精品视频女| 久99久视频精品免费| 午夜福利网站1000一区二区三区| 亚洲国产色片| 久久热精品热| 秋霞在线观看毛片| 中文欧美无线码| 国产亚洲一区二区精品| 精品人妻视频免费看| 色视频www国产| 2021天堂中文幕一二区在线观| 一级黄片播放器| 22中文网久久字幕| 日韩 亚洲 欧美在线| ponron亚洲| 男人爽女人下面视频在线观看| 九九在线视频观看精品| 99久国产av精品国产电影| 美女脱内裤让男人舔精品视频| 美女内射精品一级片tv| 插逼视频在线观看| 久久久色成人| 在线观看av片永久免费下载| 99热全是精品| 少妇人妻一区二区三区视频| 欧美极品一区二区三区四区| 久久精品国产自在天天线| 晚上一个人看的免费电影| 亚洲成人中文字幕在线播放| 熟妇人妻不卡中文字幕| www.色视频.com| 伊人久久国产一区二区| 黑人高潮一二区| 真实男女啪啪啪动态图| 久久久国产一区二区| 啦啦啦韩国在线观看视频| 一边亲一边摸免费视频| 午夜免费男女啪啪视频观看| 美女xxoo啪啪120秒动态图| 亚洲无线观看免费| 天美传媒精品一区二区| 白带黄色成豆腐渣| 久久午夜福利片| 精品不卡国产一区二区三区| 午夜精品在线福利| 寂寞人妻少妇视频99o| 亚洲精品日本国产第一区| 亚洲精品aⅴ在线观看| 五月天丁香电影| 日韩欧美 国产精品| 熟女人妻精品中文字幕| 久久99精品国语久久久| 亚洲精品中文字幕在线视频 | 日韩一区二区三区影片| 美女大奶头视频| 晚上一个人看的免费电影| 久久久久精品久久久久真实原创| 精品久久久久久久久久久久久| 七月丁香在线播放| 免费观看精品视频网站| 一级毛片电影观看| 国产精品99久久久久久久久| 91在线精品国自产拍蜜月| 色综合色国产| 最近手机中文字幕大全| 成人毛片a级毛片在线播放| 国产人妻一区二区三区在| ponron亚洲| 秋霞在线观看毛片| 婷婷色综合大香蕉| 国产在视频线在精品| 我的女老师完整版在线观看| 可以在线观看毛片的网站| 国产亚洲最大av| 亚洲精品一二三| 日本av手机在线免费观看| 亚洲国产日韩欧美精品在线观看| 天美传媒精品一区二区| 国产男女超爽视频在线观看| 亚洲国产精品成人久久小说| 亚洲精品,欧美精品| 在线观看www视频免费| 观看av在线不卡| 久久精品国产a三级三级三级| 亚洲综合色惰| 久久婷婷青草| 我的亚洲天堂| 亚洲人成77777在线视频| 欧美日韩国产mv在线观看视频| 亚洲第一区二区三区不卡| 欧美日韩精品网址| 制服诱惑二区| 亚洲图色成人| 欧美日韩亚洲国产一区二区在线观看 | 我的亚洲天堂| 老司机影院毛片| 亚洲精品国产av成人精品| 纵有疾风起免费观看全集完整版| 韩国高清视频一区二区三区| 搡女人真爽免费视频火全软件| 亚洲少妇的诱惑av| 国产亚洲午夜精品一区二区久久| 国产伦理片在线播放av一区| 91在线精品国自产拍蜜月| 夫妻性生交免费视频一级片| 亚洲国产色片| 制服诱惑二区| videos熟女内射| 亚洲av福利一区| 久久久欧美国产精品| 日本av手机在线免费观看| 久久女婷五月综合色啪小说| 国产精品熟女久久久久浪| 视频在线观看一区二区三区| 国产欧美日韩一区二区三区在线| 久久久久视频综合| 久久毛片免费看一区二区三区| 日韩成人av中文字幕在线观看| 亚洲经典国产精华液单| 久久久久网色| 如日韩欧美国产精品一区二区三区| 亚洲美女搞黄在线观看| 亚洲av免费高清在线观看| 久久久久精品人妻al黑| 99九九在线精品视频| 国产女主播在线喷水免费视频网站| 韩国av在线不卡| 少妇被粗大猛烈的视频| 国产日韩欧美在线精品| 久久ye,这里只有精品| 宅男免费午夜| av网站在线播放免费| 1024香蕉在线观看| 久久 成人 亚洲| 午夜福利在线观看免费完整高清在| 人妻少妇偷人精品九色| 日本av免费视频播放| 国产乱人偷精品视频| 国产男女内射视频| 不卡视频在线观看欧美| 人妻人人澡人人爽人人| 亚洲精品日韩在线中文字幕| 九九爱精品视频在线观看| 欧美精品国产亚洲| 亚洲在久久综合| 王馨瑶露胸无遮挡在线观看| 久久免费观看电影| 亚洲精品国产av蜜桃| 精品久久蜜臀av无| 高清视频免费观看一区二区| 午夜福利影视在线免费观看| 日韩av在线免费看完整版不卡| 久久久久久久大尺度免费视频| 爱豆传媒免费全集在线观看| 尾随美女入室| av电影中文网址| 精品人妻偷拍中文字幕| 久久久久久久国产电影| 日韩av免费高清视频| 精品少妇黑人巨大在线播放| 中文字幕精品免费在线观看视频| 免费人妻精品一区二区三区视频| 97在线人人人人妻| 亚洲欧洲精品一区二区精品久久久 | 蜜桃在线观看..| 啦啦啦在线免费观看视频4| 亚洲欧美精品自产自拍| 成人二区视频| av又黄又爽大尺度在线免费看| 18+在线观看网站| 国产免费视频播放在线视频| 伊人亚洲综合成人网| 在线亚洲精品国产二区图片欧美| 2018国产大陆天天弄谢| 久久久久精品人妻al黑| 老司机影院成人| 国产 一区精品| 精品亚洲成国产av| 久久久久人妻精品一区果冻| 成人亚洲精品一区在线观看| 一二三四在线观看免费中文在| 国产国语露脸激情在线看| 一区福利在线观看| 啦啦啦在线免费观看视频4| 满18在线观看网站| 另类精品久久| 精品午夜福利在线看| 在线精品无人区一区二区三| 午夜福利影视在线免费观看| 国产精品成人在线| 国产精品久久久av美女十八| 亚洲精品久久午夜乱码| 欧美xxⅹ黑人| 欧美 日韩 精品 国产| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产亚洲av涩爱| av不卡在线播放| 亚洲天堂av无毛| 免费少妇av软件| 久久青草综合色| 午夜福利网站1000一区二区三区| 午夜精品国产一区二区电影| 亚洲欧洲日产国产| 五月开心婷婷网| 亚洲欧美一区二区三区国产| 日韩制服骚丝袜av| 欧美老熟妇乱子伦牲交| 国产在视频线精品| 制服诱惑二区| 久久狼人影院| 日日摸夜夜添夜夜爱| 老鸭窝网址在线观看| 日韩制服丝袜自拍偷拍| 国产亚洲最大av| 精品国产露脸久久av麻豆| 天堂俺去俺来也www色官网| 免费不卡的大黄色大毛片视频在线观看| 亚洲,欧美,日韩| 亚洲av免费高清在线观看| av免费观看日本| 日韩精品免费视频一区二区三区| 日日摸夜夜添夜夜爱| 99精国产麻豆久久婷婷| 久久久精品区二区三区| 亚洲av福利一区| 国产成人精品福利久久| 久久久久国产精品人妻一区二区| 久久久久久伊人网av| 人人妻人人爽人人添夜夜欢视频| 色婷婷久久久亚洲欧美| 视频区图区小说| 亚洲精华国产精华液的使用体验| 婷婷色综合www| 一区在线观看完整版| 精品一品国产午夜福利视频| 精品卡一卡二卡四卡免费| 婷婷色麻豆天堂久久| 久久久a久久爽久久v久久| 日本爱情动作片www.在线观看| 精品福利永久在线观看| 亚洲av在线观看美女高潮| 国产一区二区 视频在线| 大话2 男鬼变身卡| 中文字幕亚洲精品专区| 亚洲欧洲精品一区二区精品久久久 | 欧美日本中文国产一区发布| 美女国产高潮福利片在线看| 搡女人真爽免费视频火全软件| 日韩在线高清观看一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 一二三四在线观看免费中文在| 久久久久久伊人网av| 精品少妇黑人巨大在线播放| 水蜜桃什么品种好| av免费在线看不卡| 欧美日本中文国产一区发布| 一级a爱视频在线免费观看| 少妇人妻久久综合中文| 精品少妇黑人巨大在线播放| 欧美精品亚洲一区二区| 99久久中文字幕三级久久日本| 久久久久精品人妻al黑| 最近2019中文字幕mv第一页| 国产精品亚洲av一区麻豆 | 日本wwww免费看| 亚洲国产精品999| 在现免费观看毛片| 久久久精品国产亚洲av高清涩受| 人妻一区二区av| 国产片特级美女逼逼视频| 极品人妻少妇av视频| 国产探花极品一区二区| 最近的中文字幕免费完整| 中文字幕色久视频| 成人亚洲精品一区在线观看| 另类精品久久| 久久久精品94久久精品| 制服丝袜香蕉在线| 久久精品国产a三级三级三级|