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

    湖泊和水庫氧化亞氮通量分析

    2018-09-22 08:09:10謝燕紅李旭輝
    中國環(huán)境科學(xué) 2018年9期
    關(guān)鍵詞:緯度湖泊通量

    謝燕紅,張 彌*,肖 薇,王 偉,李旭輝,卞 航

    ?

    湖泊和水庫氧化亞氮通量分析

    謝燕紅1,2,張 彌1,2*,肖 薇1,2,王 偉1,2,李旭輝1,2,卞 航1,2

    (1.南京信息工程大學(xué)大氣環(huán)境中心,江蘇 南京 210044;2.南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,江蘇 南京 210044)

    內(nèi)陸水體中,湖泊與水庫作為重要的N2O潛在排放源,在氮循環(huán)中扮演著重要的角色.通過調(diào)研文獻(xiàn),在獲取全球479個(gè)湖泊和83個(gè)水庫N2O通量數(shù)據(jù)的基礎(chǔ)上,主要分析了湖泊和水庫所處的緯度、面積、深度以及營養(yǎng)狀況4個(gè)因素對(duì)N2O通量的影響,并估算全球湖泊和水庫N2O年總通量.結(jié)果表明,湖泊和水庫表現(xiàn)為N2O的源,湖泊N2O平均通量[(3.21±5.71)μmol/(m2?d)]低于水庫N2O平均通量[(20.82±113.94)μmol/(m2?d)].在全年尺度上,湖泊和水庫N2O通量隨緯度增加呈下降的趨勢(shì);湖泊N2O平均通量隨著面積的增加而增加,水庫N2O平均通量隨著面積的增加而下降;深水湖泊的N2O通量大于淺水湖泊;富營養(yǎng)湖泊及水庫的N2O通量高于貧營養(yǎng)湖泊及水庫.湖泊和水庫N2O年平均總通量分別為0.12和0.06Tg N/a,占到內(nèi)陸水體N2O排放的12%和7%.

    湖泊;水庫;N2O通量;影響因子;N2O年通量

    N2O作為第3大溫室氣體,其百年增溫潛勢(shì)是當(dāng)量CO2的298倍[1],同時(shí)還是消耗平流層臭氧的重要物質(zhì)[2].由于化學(xué)氮肥的大量使用,大量的氮(N)進(jìn)入水生生態(tài)系統(tǒng),使得水生生態(tài)系統(tǒng)儲(chǔ)存的N增多[3-4].全球湖泊和水庫的面積為3.0×106~4.8×106km2[5-7],占到內(nèi)陸水體的83%~89%.并且,這兩類水體作為內(nèi)陸水體最終的氮匯聚地,被認(rèn)為是N2O潛在的排放源.因此,湖泊和水庫向大氣中釋放的N2O受到了重點(diǎn)的關(guān)注.

    目前,對(duì)湖泊和水庫N2O通量進(jìn)行了大量的研究.Mccrackin等[4]在2011年首次利用模型和IPCC方法估算得出全球湖泊和水庫N2O平均排放量達(dá)到0.6Tg N/a,但該研究的估算方法只考慮了大氣N輸入對(duì)湖泊和水庫N2O排放的影響,并沒有包括農(nóng)田氮肥徑流和污水輸入的影響,因而估計(jì)值可能偏低.也有學(xué)者通過文獻(xiàn)調(diào)研的方式對(duì)全球湖泊或者水庫的N2O排放進(jìn)行了估計(jì),但由于并沒有考慮湖泊及水庫N2O的時(shí)間代表性問題,可能導(dǎo)致其估計(jì)值偏差較大[8-9].因此,在對(duì)全球湖泊及水庫N2O排放進(jìn)行估算時(shí)還存在不確定性,需要考慮湖泊、水庫的空間及時(shí)間的代表性.

    水體N2O主要是由微生物的硝化和反硝化過程產(chǎn)生,并通過擴(kuò)散途徑釋放到大氣中[10],在其產(chǎn)生和排放過程中會(huì)受到諸多因子的影響.對(duì)于單個(gè)湖泊和水庫,N2O主要受到水體物理化學(xué)參數(shù)的影響.研究表明,N2O濃度與水溫、含氮無機(jī)鹽(DIN、NO3-、NH4+)呈正比,與溶解氧(DO)、pH呈反比[11-12].在全球尺度上,不同湖泊和水庫的地理位置、形態(tài)特征以及營養(yǎng)狀況的不同差異可能決定著N2O的排放量.因此,在全球尺度上,得到影響湖泊與水庫N2O排放的因素成為準(zhǔn)確估算全球湖泊及水庫N2O排放的關(guān)鍵問題.

    本文基于文獻(xiàn)調(diào)研,在收集全球不同區(qū)域湖泊與水庫N2O排放通量數(shù)據(jù)及相關(guān)數(shù)據(jù)的基礎(chǔ)上分析湖泊和水庫所處緯度、形態(tài)特征(面積和深度)、營養(yǎng)狀況對(duì)N2O通量的影響;并對(duì)全球湖泊和水庫N2O排放量進(jìn)行估算,通過與河流、濕地等水生生態(tài)系統(tǒng)的N2O排放量進(jìn)行對(duì)比,進(jìn)而評(píng)價(jià)湖泊和水庫在內(nèi)陸水體N2O排放中的水平.從而為全球不同水體的N2O排放的準(zhǔn)確估算提供基礎(chǔ).

    1 材料與方法

    1.1 數(shù)據(jù)來源及特征描述

    本文通過調(diào)研已發(fā)表的文獻(xiàn),獲取湖泊和水庫N2O通量的直接測(cè)量數(shù)據(jù).共篩選了43篇文獻(xiàn),包括479個(gè)湖泊和83個(gè)水庫,其分布見圖1(詳細(xì)信息見表1和表2),可以看出選取的湖泊和水庫主要集中在北美、歐洲和中國這3個(gè)地區(qū),緯度跨度為70°S~70°N.調(diào)研的北半球的湖泊和水庫為549個(gè),南半球的湖泊和水庫為12個(gè).參考Soued等[9]對(duì)湖泊的緯度劃分,本文調(diào)研的湖泊和水庫位于熱帶地區(qū)(緯度<24°)的有23個(gè),占到調(diào)研湖泊和水庫總數(shù)的4%;位于溫帶地區(qū)(緯度在24°~55°之間)有451個(gè),占到調(diào)研湖泊和水庫總數(shù)的80%;位于高緯度地區(qū)(緯度>55o)有88個(gè),占到調(diào)研湖泊和水庫總數(shù)的16%.

    本文選取的文獻(xiàn)測(cè)定湖泊和水庫N2O通量的方法有2種:箱式法和頂空平衡法,分別占到40%和60%.篩選的湖泊和水庫中,觀測(cè)時(shí)間最早從1977年,最新至2017年;最短為1個(gè)月,最長(zhǎng)為7a,具有全年觀測(cè)的占到29%.空間采樣點(diǎn)最少有1個(gè),最多有29個(gè),其中,觀測(cè)點(diǎn)數(shù)>3的占到30%左右.

    圖1 本文選取湖泊及水庫的地理位置

    不同符號(hào)代表位置比較集中的湖泊

    1.2 數(shù)據(jù)提取與處理

    1.2.1 湖泊和水庫N2O通量數(shù)據(jù)提取及處理 本文篩選的湖泊和水庫的N2O通量既有全年的平均值,又有不同季節(jié)的平均值(表1、表2).具有全年觀測(cè)的湖泊有116個(gè),占到調(diào)研湖泊的24%.其中102個(gè)湖泊包含在12個(gè)區(qū)域內(nèi),文獻(xiàn)僅給出該12個(gè)區(qū)域?qū)?yīng)的N2O平均通量,故在之后的結(jié)果分析中用該12個(gè)平均觀測(cè)值代表對(duì)應(yīng)區(qū)域的湖泊[28];另外14個(gè)湖泊具有其各自的年N2O通量,其中太湖和Dendre stone pit lake的年N2O通量通過4個(gè)季節(jié)平均獲得[19,29],Muskegon Lake的年N2O通量通過月平均獲得[23],Lake Neusiedl的年N2O通量通過月份和空間平均獲得[27],其余均是文獻(xiàn)中給出的年平均值.對(duì)于進(jìn)行季節(jié)觀測(cè)的湖泊,Lake Huahu的N2O通量通過空間平均獲得[21],其余也均是文獻(xiàn)中給出的平均值.

    表1 湖泊的地理位置及觀測(cè)時(shí)間

    續(xù)表1

    續(xù)表1

    表2 水庫的地理位置及觀測(cè)時(shí)間

    續(xù)表2

    續(xù)表2

    全年觀測(cè)的水庫有49個(gè),占到調(diào)研水庫的60%.其中Lake Wivenhoe、Lake Baroon和Little Nerang Dam的年N2O通量通過4個(gè)季節(jié)平均獲得[35], Harsha Lake的年N2O通量通過月份和空間平均獲得[45],玉渡山水庫的年N2O通量通過季節(jié)和空間平均獲得[47],其余均是文獻(xiàn)中給出的年平均值.對(duì)于季節(jié)觀測(cè)的水庫,Lake Lungern、High C reservoir site、Medium C reservoir site和Low C reservoir site的N2O通量通過日平均獲得[49-50],Gold Creek Reservoir的N2O通量通過空間平均獲得[10],2nd Baoshan、Liyutan和Tsengwen的N2O通量通過月和空間平均獲得[40],其余是文獻(xiàn)中給出的平均值.

    在調(diào)研的湖泊中,Great Lakes和Lake Baldegg給出N2O通量范圍,本文用所給范圍值中最大值和最小值的平均值代替[24,26];Muskegon Lake、Lake Neusiedl和玉渡山水庫的通量數(shù)據(jù)以圖的方式給出[23,27,47],本文用Get Data Graph Digitizer軟件提取所需數(shù)據(jù);另有在加拿大Québec地區(qū)、中國天津市、北美五大湖地區(qū)、愛爾蘭地區(qū)以及英國的Upper Thurne catchment分別進(jìn)行了多個(gè)湖泊的觀測(cè),但只給出了區(qū)域的平均值,本文用平均值來代替這5個(gè)區(qū)域湖泊的N2O通量[9,22,24,28,30-31];有3個(gè)湖泊(太湖、鄱陽湖和Lake Daming)、加拿大Québec地區(qū)湖泊以及1個(gè)水庫(三峽水庫)出現(xiàn)了重復(fù)觀測(cè),本文選取采樣點(diǎn)更多和觀測(cè)時(shí)間更長(zhǎng)的來代表該湖泊(區(qū)域湖泊)或水庫的N2O通量,但是對(duì)于有不同季節(jié)觀測(cè)的,會(huì)保留以便于體現(xiàn)季節(jié)變化對(duì)N2O通量的影響[9,12-13,15-20,28,42];對(duì)于在湖濱帶測(cè)量得到的N2O通量,只選取了長(zhǎng)期淹水區(qū)域的數(shù)值[21,48].

    盡管具有季節(jié)觀測(cè)的湖泊和水庫占到71%,但由于N2O通量存在季節(jié)差異[17,29,43],因此,主要對(duì)具有全年觀測(cè)的湖泊和水庫N2O通量進(jìn)行探討,并用全年觀測(cè)得到的數(shù)值進(jìn)行全球估算.

    1.2.2 湖泊和水庫特征數(shù)據(jù)提取及處理 除了統(tǒng)計(jì)湖泊和水庫的N2O通量外,本文還收集了湖泊和水庫的地理位置信息、面積、容量、滯留時(shí)間、水庫庫齡、深度(包括平均深度和最大深度)、水溫、營養(yǎng)狀態(tài).對(duì)于文獻(xiàn)沒有提及的,通過其研究相關(guān)的文獻(xiàn)或相關(guān)網(wǎng)站獲取.對(duì)于仍無法獲取經(jīng)緯度的湖泊或水庫,選擇用其所在地區(qū)的經(jīng)緯度信息代替.

    基于Kolmogorov-Smirnov 檢驗(yàn)發(fā)現(xiàn)湖泊和水庫的N2O通量均不符合正態(tài)分布,因此在結(jié)果部分不同面積、深度以及營養(yǎng)狀態(tài)對(duì)N2O通量的影響均采用非參數(shù)檢驗(yàn).

    2 結(jié)果與分析

    2.1 湖泊和水庫N2O通量

    基于篩選,本文得到34個(gè)湖泊站點(diǎn)年N2O通量數(shù)據(jù)(表3),得出N2O通量的變化范圍為-5.32~ 21.28μmol/(m2?d),平均值為(3.21±5.71)μmol/(m2?d).其中,71%的數(shù)值為正,說明湖泊主要表現(xiàn)為N2O的源;29%的數(shù)值為負(fù),說明有小部分湖泊表現(xiàn)為N2O的匯,其原因可能是這部分湖泊的N2O處在極其厭氧的環(huán)境下被進(jìn)一步還原成N2,或是由于NO3-濃度水平受到限制,抑制反硝化速率,進(jìn)而影響N2O產(chǎn)生,N2O濃度在水中的濃度減小,從而導(dǎo)致N2O從大氣進(jìn)入水體.不同于自然湖泊,水庫作為人工水體,受到人類的干擾更大,基于50個(gè)水庫站點(diǎn)年數(shù)據(jù)(表4),N2O通量數(shù)值波動(dòng)較大,變化范圍為-2.50~ 815.50μmol/(m2?d),平均值(20.82±113.74)μmol/ (m2?d),98%是正通量,說明水庫也主要是N2O的源.綜上所述,水庫N2O通量變化范圍較湖泊的大,并且水庫的N2O平均通量也高于湖泊的N2O平均通量.

    表3 開展全年觀測(cè)的湖泊N2O平均通量

    表4 開展全年觀測(cè)的水庫N2O平均通量

    2.2 影響 N2O通量的因素

    2.2.1 緯度對(duì)N2O通量的影響 從圖2(a)可以看出:將文獻(xiàn)調(diào)研的所有數(shù)據(jù),包括全年、不同季節(jié)的數(shù)據(jù)與緯度建立關(guān)系,發(fā)現(xiàn)隨著緯度的增加,湖泊N2O通量的變化趨勢(shì)不明顯,N2O通量在-5.46~ 38.70μmol/(m2?d)范圍內(nèi).由于對(duì)湖泊N2O通量的觀測(cè)主要集中在夏季和全年,分析緯度對(duì)湖泊夏季和全年N2O通量的影響發(fā)現(xiàn):北半球湖泊夏季N2O通量與緯度呈線性減小(2=0.34,<0.05,=13)[圖2(b)].圖2(b)中圈出來的點(diǎn)是由74個(gè)湖泊的N2O通量數(shù)值平均獲得,這74個(gè)湖泊的N2O通量隨緯度增加呈指數(shù)遞減(2=0.13,<0.001,=74).基于湖泊全年的數(shù)值,在北半球N2O通量與緯度也存在著較好的負(fù)相關(guān)關(guān)系(2=0.21,<0.01,=33)[圖2(c)].

    (a)~(c)分別代表湖泊所有時(shí)間段、夏季、全年;(d)~(f)分別代表水庫所有時(shí)間段、夏季、全年

    對(duì)水庫的分析與湖泊相似,將所有調(diào)研數(shù)據(jù)與緯度建立關(guān)系時(shí),水庫N2O通量隨著緯度的增加也并無明顯的變化趨勢(shì),N2O通量的范圍為-3.18~ 815.50 μmol/(m2?d)[圖2(d)].相比湖泊,夏季水庫N2O通量并未隨緯度增加而有明顯的變化趨勢(shì)[圖2(e)];全年水庫N2O通量盡管隨緯度增加呈下降的趨勢(shì)[圖2(f)],但是在北緯40o附近,出現(xiàn)兩個(gè)較高的值.

    2.2.2 面積對(duì)N2O通量的影響 基于湖泊和水庫站點(diǎn)年數(shù)據(jù),湖泊和水庫面積大小范圍分別為0.03~3210km2和0.5~4318km2.參照Downing等[6]的研究對(duì)湖泊和水庫面積大小的分級(jí),分成6個(gè)等級(jí)(<0.1km2、0.1~1km2、1~10km2、10~100km2、100~ 1000km2和>1000km2).圖3(a)、3(b)所示,湖泊N2O通量最低出現(xiàn)在0.1~1km2范圍內(nèi),為(2.15± 3.06)μmol/(m2?d),最高值出現(xiàn)在面積>1000km2范圍內(nèi),為(7.77+8.65)μmol/(m2?d);湖泊面積大于0.1km2后,N2O通量隨湖泊面積增加而增加,但不同面積大小等級(jí)之間的N2O通量并無顯著差異.

    (a)湖泊,(b)水庫;不同字母表示差異顯著(<0.05),相同字母表示差異不顯著

    水庫N2O通量最低出現(xiàn)在面積>1000km2范圍內(nèi),為(1.25±2.10)μmol/(m2?d),最高值出現(xiàn)在1~10km2范圍內(nèi),為(173.72±321.33)μmol/(m2?d);水庫面積大于1km2后,隨著水庫面積的增加,N2O通量呈下降的趨勢(shì),但不同面積大小等級(jí)之間的N2O通量也并無顯著差異.

    2.2.3 深度對(duì)N2O通量的影響 根據(jù)明尼蘇達(dá)州污染控制機(jī)構(gòu)給出的標(biāo)準(zhǔn),認(rèn)為當(dāng)湖泊最大深度<4.6m時(shí),為淺水湖泊.基于獲取的湖泊和水庫站點(diǎn)年數(shù)據(jù),本文將最大深度£4.6m的湖泊及水庫歸為淺水湖泊及水庫,將最大深度>4.6m的湖泊及水庫歸為深水湖泊及水庫.深水湖泊的平均N2O通量[(5.22±6.41)μmol/(m2?d)]要大于淺水湖泊[(2.26± 2.64)μmol/(m2?d)], 但二者差異并不顯著(圖4);由于本文統(tǒng)計(jì)到的淺水水庫最大深度的樣本量小于3,因此并未對(duì)水庫進(jìn)行分析.

    圖4 湖泊N2O通量隨深度的變化

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

    2.2.4 營養(yǎng)狀態(tài)對(duì)N2O通量的影響 根據(jù)文獻(xiàn)調(diào)研的湖泊和水庫站點(diǎn)年數(shù)據(jù),湖泊和水庫營養(yǎng)狀態(tài)有4種:貧營養(yǎng)、中營養(yǎng)、中富營養(yǎng)和富營養(yǎng).由于在調(diào)研的湖泊與水庫中只有1個(gè)水庫是中富營養(yǎng),因此在后面的結(jié)果與分析中沒有考慮.如圖5(a)所示,雖然湖泊N2O平均通量是中營養(yǎng)湖泊[(11.62± 10.87)μmol/(m2?d)]>富營養(yǎng)湖泊[(3.89±5.24)μmol/ (m2?d)]>貧營養(yǎng)湖泊[(0.75±1.21)μmol/(m2?d)],但不同富營養(yǎng)程度之間無明顯差異.對(duì)于水庫 [圖5(b)], N2O通量平均值是富營養(yǎng)水庫[(109.63±285.62) μmol/(m2?d)]>中營養(yǎng)水庫[(3.28±0.64)μmol/(m2?d)]>貧營養(yǎng)水庫[(1.73±0.98)μmol/(m2?d)],但其相互之間也沒有顯著性差異.

    (a)湖泊, (b)水庫; 不同字母表示差異顯著(<0.05),相同字母表示差異不顯著

    2.3 估算全球湖泊和水庫N2O年通量

    基于湖泊和水庫站點(diǎn)年數(shù)據(jù)得到的N2O平均通量值,即(3.21±5.71)μmol/(m2?d)和(20.82±113.94)μmol/ (m2?d),再根據(jù)文獻(xiàn)得出的全球湖泊和水庫的面積,分別為2.7×106~4.5×106km2[5-7]和0.3×106km2[52],按平年(=365d)計(jì)算,最后估算得出全球湖泊和水庫N2O年總通量分別為0.09~0.15Tg N/a和0.06Tg N/a.

    3 討論

    3.1 湖泊和水庫N2O通量的對(duì)比

    水生生態(tài)系統(tǒng)儲(chǔ)存的氮主要有以下3個(gè)過程:反硝化過程、沉淀到底泥以及被水生植物吸收.其中,反硝化作用是最主要的機(jī)制[53],即大部分的N會(huì)以N2O或N2的形式釋放到大氣中.相對(duì)于湖泊,由于水庫具有更高的徑流系數(shù)和N表觀沉降速度,以及更大的N平均負(fù)荷率,使得水庫單位面積上每年儲(chǔ)存的N要大于湖泊的N[54].由于水庫有更多的可利用N,更有利于發(fā)生反硝化作用[55],進(jìn)而導(dǎo)致水庫N2O平均通量要大于湖泊的N2O平均通量.本研究結(jié)果與此一致,即水庫的N2O通量高于湖泊.

    3.2 湖泊和水庫N2O通量的影響因子

    對(duì)于湖泊夏季、全年以及水庫全年,N2O通量均隨著緯度的增加而呈下降的趨勢(shì).即緯度是N2O變化的影響因子,這主要是因?yàn)榫暥戎苯佑绊憸囟?溫度的改變導(dǎo)致N2O通量的變化.但本文調(diào)研的湖泊(包括水庫)N2O通量與水溫卻不相關(guān)(=0.04,=0.78,=62),這說明N2O排放還受到其他因素的影響.由于低緯度地區(qū)湖泊和水庫的N平均表觀沉降速度要高于高緯度地區(qū),表明低緯度地區(qū)湖泊和水庫保留的N比高緯度地區(qū)的多[54],再加上低緯度地區(qū)的溫度高于高緯度地區(qū),更有利于硝化和反硝化過程,進(jìn)而使得低緯度地區(qū)湖泊和水庫N2O通量要高于高緯度地區(qū).由于本文提取到的南半球湖泊和水庫較少,南半球湖泊和水庫N2O通量與緯度的關(guān)系并不明顯[圖2(a)、2(d)].

    湖泊及水庫面積的大小,決定了湖泊及水庫與周圍生態(tài)系統(tǒng)物質(zhì)和能量的交換量.在全球尺度上,小型湖泊和水庫單位面積上每年儲(chǔ)存的N要大于大型湖泊和水庫[54],并且每年由反硝化過程釋放的N也要大于大型湖泊和水庫[55],因此小型湖泊和水庫單位面積上產(chǎn)生的N2O可能要比大型湖泊和水庫的多.本文得出小型水庫的N2O平均通量高于大型水庫的結(jié)論與上述解釋一致.對(duì)于湖泊而言,由于面積在1~10km2范圍的湖泊中,Lake Baldegg的N2O通量主要由數(shù)值范圍的最大值和最小值平均獲得[26],面積在1~10km2和100~1000km2范圍的湖泊中,Muskegon Lake和Lake Neusiedl的N2O通量是通過提取數(shù)據(jù)的平均獲得[23,27],會(huì)導(dǎo)致結(jié)果產(chǎn)生一定的不確定性;面積>1000km2的湖泊主要來自于中國大型的富營養(yǎng)湖泊,根據(jù)2.2.4結(jié)果可知富營養(yǎng)湖泊具有較高的N2O通量,從而導(dǎo)致湖泊N2O通量與面積的關(guān)系與之前的研究不相符,因此湖泊N2O通量與水域面積的關(guān)系還會(huì)受到湖泊富營養(yǎng)化程度的影響.

    本文分析發(fā)現(xiàn),深水湖泊的N2O通量大于淺水湖泊,即深水湖泊有利于N2O的產(chǎn)生和排放.對(duì)于深、淺水湖泊,一個(gè)很重要的區(qū)別在于水體溶解氧含量,即湖泊深度越接近底泥,溶解氧越低.對(duì)于淺水湖泊,溶解氧含量較高,有利于硝化作用的發(fā)生,而由深水湖泊營造的厭氧環(huán)境則更有利于反硝化作用[56].盡管極厭氧環(huán)境下,N2O反而會(huì)被消耗而生成N2[57],且N的去除與湖泊的平均深度呈負(fù)相關(guān)[58],但在深水湖泊,N2O除了主要由底泥產(chǎn)生外[26],也有研究發(fā)現(xiàn)變溫層是N2O的第二個(gè)源,并且這部分N2O是由硝化作用產(chǎn)生[59],因此可能使得深水湖泊釋放的N2O通量較高.

    總氮(TN)、總磷(TP)和葉綠素a是評(píng)價(jià)湖泊和水庫富營養(yǎng)化狀態(tài)的重要指標(biāo).其中,與N2O直接相關(guān)的是TN.研究表明N2O濃度與TN呈正相關(guān)[11,32].但也有研究發(fā)現(xiàn)在TP濃度較低的情況下,反硝化速率反而更快[56].Knowles等[60]發(fā)現(xiàn)一個(gè)N:P比較低即TP濃度高的富營養(yǎng)化湖泊可能會(huì)使得細(xì)菌的反硝化代謝所需的N和P比例失衡,進(jìn)而影響反硝化速率.本文得出中營養(yǎng)湖泊的N2O通量具有相對(duì)較高的水平,說明這些湖泊TN和TP濃度比例可能更有利于N2O的產(chǎn)生.貧營養(yǎng)湖泊和水庫由于受到N的限制使得到的N2O通量最低,這與其他研究一致[32].

    3.3 與相關(guān)研究的比較

    本文調(diào)研得到的湖泊N2O通量值為(3.21± 5.71)μmol/(m2?d),年通量為0.12(0.09~0.15)Tg N/a (表5).水庫的N2O通量為(20.82±113.94)μmol/(m2?d),年通量為0.06Tg N/a,由于本文是在Deemer等[8]的研究基礎(chǔ)上再增加了水庫的樣本量,并且采用的是水庫站點(diǎn)年數(shù)據(jù),因而得到的N2O通量的平均值相對(duì)較大,估算的N2O排放量是Deemer等[8]得出的排放量值的兩倍.整合湖泊和水庫得到N2O通量的平均值為12.02μmol/(m2?d),年通量為0.4~0.6Tg N/a,與McCrackin等[4]得到的均值(0.6Tg N/a)比較相近,但由于其采用的是IPCC方法計(jì)算得到的結(jié)果,在使用IPCC方法時(shí),沒有針對(duì)湖泊的排放系數(shù),且給出的排放系數(shù)為0.0075,不確定性范圍較大(0.0005~0.025之間),得到的估算值可能也會(huì)存在較大的不確定性.本文得到的湖泊和水庫N2O年排放量與Soued等[9]得到的結(jié)果接近.

    根據(jù)表5,相比于其他水生生態(tài)系統(tǒng),湖泊N2O通量值相對(duì)偏低,尤其對(duì)比河流,湖泊N2O通量低了20倍左右,而水庫N2O通量則處于一個(gè)相對(duì)較高的水平.對(duì)于年排放量,湖泊和水庫N2O年排放量均低于濕地、河口、大陸架以及海洋的排放量.對(duì)于河流N2O排放量的計(jì)算既有單純文獻(xiàn)調(diào)研的方式得到的數(shù)值[9],又有通過模型或IPCC方法計(jì)算得到[3,61-66].由表5可知,本文得到的湖泊N2O排放量與Soued等[9]得到的河流排放的值相當(dāng),但與用模型或IPCC方法得到的排放量相比又處于較低水平.

    表5 不同水生生態(tài)系統(tǒng)N2O通量的比較

    基于表5不同水生生態(tài)系統(tǒng)的N2O年排放量(湖泊和水庫均用的是本文得到的值,河流采用的是平均值),本文計(jì)算得出湖泊排放的N2O占到內(nèi)陸水體排放的12%,水庫占7%.即在內(nèi)陸水體,湖泊和水庫N2O排放占19%.在計(jì)算總的內(nèi)陸水體排放量時(shí),忽略了池塘這部分小型水體的排放量,因而該比例被高估.如果考慮小型水體的排放量,湖泊和水庫所占比例可能會(huì)降低.此外,根據(jù)Dowing等[6]和Lehner等[52]的估算,全球有3.0×108個(gè)湖泊,1.67×108個(gè)水庫,而本文通過文獻(xiàn)調(diào)研到479個(gè)湖泊和83個(gè)水庫,約占全球湖泊的1.58×10-4%,全球水庫的4.97×10-5%.盡管所占的比例較小,但本文調(diào)研的湖泊和水庫位于中高緯度地區(qū)的占到調(diào)研總數(shù)的90%以上,且包括了熱帶及亞熱帶的湖泊,因而仍具有一定的代表性.由于文獻(xiàn)調(diào)研的方式存在很大的不確定性,因此需要建立模型得出更精確的結(jié)果.

    在本文調(diào)研的文獻(xiàn)中采用箱式法和頂空平衡法觀測(cè)N2O通量的比例比較接近,并且有觀測(cè)得到N2O在水體里面是以擴(kuò)散方式為主[10],因而本文在對(duì)數(shù)據(jù)進(jìn)行分析時(shí),并沒有考慮觀測(cè)方法的問題,而假定這兩種方法之間是沒有顯著差異的.但也有研究表明水平衡法得到的通量要低估箱式法得到的通量[70],因而觀測(cè)方法的不同仍會(huì)增加N2O通量及排放值估算的誤差.因此,在之后的研究中需要考慮觀測(cè)方法帶來的不確定性.

    4 結(jié)論

    4.1 湖泊和水庫均主要表現(xiàn)為N2O的源,湖泊N2O通量的變化范圍較水庫的小,其N2O通量的平均值[(3.21±5.71)μmol/(m2?d)]小于水庫的N2O通量[(20.82±113.94)μmol/(m2?d)].湖泊N2O年總通量為0.12Tg N/a,水庫N2O年總通量為0.06Tg N/a,湖泊和水庫N2O年排放量均低于其他水生生態(tài)系統(tǒng)的排放量;通過計(jì)算,湖泊和水庫N2O排放分別占到內(nèi)陸水體排放的12%和7%.

    4.2 在年尺度上,無論是湖泊還是水庫的N2O通量,都隨著緯度的增加而呈下降的趨勢(shì).湖泊N2O通量隨著面積大小等級(jí)的增加而增加,而水庫N2O通量隨著面積大小等級(jí)的增加而降低;深水湖泊N2O通量大于淺水湖泊;不同營養(yǎng)狀態(tài)湖泊和水庫的N2O平均通量均無顯著差異.

    [1] Climate change 2013-the physical science basis: contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [M]. Cambridge: University Press, 2013.

    [2] Ravishankara A R, Daniel J S, Portmann R W. Nitrous oxide (N2O): the dominant ozone-depleting substance emitted in the 21st century [J]. Journal of Transport Geography, 2009,19(19):807-820.

    [3] Seitzinger S P, Kroeze C, Styles R V. Global distribution of N2O emissions from aquatic systems: natural emissions and anthropogenic effects [J]. Chemosphere-Global Change Science, 2000,2(3/4):267- 279.

    [4] Mccrackin M L, Elser J J. Greenhouse gas dynamics in lakes receiving atmospheric nitrogen deposition [J]. Global Biogeochemical Cycles, 2011,25(4):327-336.

    [5] Raymond P A, Hartmann J, Lauerwald R, et al. Global carbon dioxide emissions from inland waters [J]. Nature, 2013,503(7476):355-359.

    [6] Downing J A, Prairie Y T, Cole J J, et al. The global abundance and size distribution of lakes, ponds, and impoundments [J]. The American Society of Limnology and Oceanography, 2006,51:2388-2397.

    [7] Verpoorter C, Kutser T, Seekell D A, et al. A global inventory of lakes based on high‐resolution satellite imagery [J]. Geophysical Research Letters, 2014,41(18):6396-6402.

    [8] Deemer B R, Harrison J A, Li S, et al. Greenhouse gas emissions from reservoir water surfaces: a new global synthesis [J]. Bioscience, 2016,66(11):949-964.

    [9] Soued C, Giorgio P A D, Maranger R. Nitrous oxide sinks and emissions in boreal aquatic networks in Quebec [J]. Nature Geoscience, 2015,9(2):116-120.

    [10] Sturm K, Yuan Z, Gibbes B, et al. Methane and nitrous oxide sources and emissions in a subtropical freshwater reservoir, South East Queensland, Australia [J]. Biogeosciences, 2014,11(18):5245-5258.

    [11] Liu X L, Liu C Q, Li S L, et al. Spatiotemporal variations of nitrous oxide (N2O) emissions from two reservoirs in SW China [J]. Atmospheric Environment, 2011,45(31):5458-5468.

    [12] Wang H X, Zhang L, Yao X L, et al. Dissolved nitrous oxide and emission relating to denitrification across the Poyang Lake aquatic continuum [J]. Journal of Environmental Sciences, 2017,52:130-140.

    [13] Liu Y S, Zhu R B, Ma D W, et al. Temporal and spatial variations of nitrous oxide fluxes from the littoral zones of three alga-rich lakes in coastal Antarctica [J]. Atmospheric Environment, 2011,45(7):1464- 1475.

    [14] Roland F A E, Darchambeau F, Morana C, et al. Nitrous oxide and methane seasonal variability in the epilimnion of a large tropical meromictic lake (Lake Kivu, East-Africa) [J]. Aquatic Sciences, 2016, 79(2):1-10.

    [15] 李香華.太湖水-氣界面溫室氣體通量及時(shí)空變化特征研究[D]. 南京:河海大學(xué), 2005.

    [16] Liu L X, Ming X, Mao L, et al. Spatial variability of greenhouse gas effluxes and their controlling factors in the Poyang Lake in China [J]. Polish Journal of Environmental Studies, 2013,22(3):749-758

    [17] 徐會(huì)顯,姜星宇,姚曉龍,等.鄱陽湖水體氧化亞氮排放特征及影響因素[J]. 湖泊科學(xué), 2016,28(5):972-981.

    [18] Wang S L, Liu C Q, Yeager K, et al. The spatial distribution and emission of nitrous oxide (N2O) in a large eutrophic lake in eastern China: anthropogenic effects [J]. Science of the Total Environment, 2009,407:3330-3337.

    [19] 肖啟濤.太湖水-氣界面溫室氣體(CO2、CH4、N2O)通量研究[D]. 南京:南京信息工程大學(xué), 2014.

    [20] 鄭小蘭,劉 敏,文帥龍,等.冬季太湖草、藻型湖區(qū)N2O的生成與排放特征 [J]. 中國環(huán)境科學(xué), 2018,38(2):683-693.

    [21] Zhu D, Wu Y, Wu N, et al. Nitrous oxide emission from infralittoral zone and pelagic zone in a shallow lake: implications for whole lake flux estimation and lake restoration [J]. Ecological Engineering, 2015, 82:368-375.

    [22] Liu X L, Bai L, Wang Z L, et al. Nitrous oxide emissions from river network with variable nitrogen loading in Tianjin, China [J]. Journal of Geochemical Exploration, 2015,157:153-161.

    [23] Salk K R, Ostrom P H, Biddanda B A, et al. Ecosystem metabolism and greenhouse gas production in a mesotrophic northern temperate lake experiencing seasonal hypoxia [J]. Biogeochemistry, 2016,131(3): 1-17.

    [24] Lemon E, Lemon D. Nitrous oxide in freshwaters of the great lakes basin [J]. Limnology & Oceanography, 1981,26(5):867-879.

    [25] Tangen B A, Finocchiaro R G, Gleason R A, et al. Greenhouse gas fluxes of a shallow lake in South-Central North Dakota, USA [J]. Wetlands, 2016,36(4):779-787.

    [26] Mengis M, G?chter R, Wehrli B. Nitrous oxide emissions to the atmosphere from an artificially oxygenated lake [J]. Limnology & Oceanography, 1996,41(3):548-553.

    [27] Soja G, Kitzler B, Soja A M. Emissions of greenhouse gases from Lake Neusiedl, a shallow steppe lake in Eastern Austria [J]. Hydrobiologia, 2014,731(1):125-138.

    [28] Tremblay A, Varfalvy L, Roehm C,et al. Greenhouse gas emissions- fluxes and processes [M]. New York: Springer Berlin Heidelberg, 2005:209-232.

    [29] Roland F A, Darchambeau F, Morana C, et al. Emission and oxidation of methane in a meromictic, eutrophic and temperate lake (Dendre, Belgium) [J]. Chemosphere, 2017,168:756-764.

    [30] Whitfield C J, Aherne J, Baulch H M. Controls on greenhouse gas concentrations in polymictic headwater lakes in Ireland [J]. Science of the Total Environment, 2011,s410-411(411):217-225.

    [31] Outram F N, Hiscock K M. Indirect nitrous oxide emissions from surface water bodies in a lowland arable catchment: a significant contribution to agricultural greenhouse gas budgets? [J]. Environmental Science & Technology, 2012,46(15):8156-8163.

    [32] Yang H, Andersen T, D?rsch P, et al. Greenhouse gas metabolism in Nordic boreal lakes [J]. Biogeochemistry, 2015,126(1/2):211-225.

    [33] Huttunen J T, Alm J, Liikanen A, et al. Fluxes of methane, carbon dioxide and nitrous oxide in boreal lakes and potential anthropogenic effects on the aquatic greenhouse gas emissions [J]. Chemosphere, 2003,52(3):609-621.

    [34] Miettinen H, Pumpanen J, Heiskanen J J, et al. Towards a more comprehensive understanding of lacustrine greenhouse gas dynamics - two-year measurements of concentrations and fluxes of CO2, CH4and N2O in a typical boreal lake surrounded by managed forests [J]. Boreal Environment Research, 2015,20(1):75-89.

    [35] Musenze R S, Grinham A, Werner U,et al. Assessing the spatial and temporal variability of diffusive methane and nitrous oxide emissions from subtropical freshwater reservoirs [J]. Environmental Science & Technology, 2014,48(24):14499-14507.

    [36] Guérin F, Abril G, Tremblay A, et al. Nitrous oxide emissions from tropical hydroelectric reservoirs [J]. Geophysical Research Letters, 2008,35(6):L06404.

    [37] Lima I B T, Victoria R L, Novo E M L M, et al. Methane, carbon dioxide and nitrous oxide emissions from two Amazonian reservoirs during high water table [J].Verhandlungen des Internationalen Verein Limnologie, 2002,28:438-442.

    [38] Deshmukh C. Greenhouse gas emissions (CH4, CO2and N2O) from a newly flooded hydroelectric reservoir in subtropical South Asia: the case of Nam Theun 2Reservoir, Lao PDR [D]. Ocean, Atmosphere. Universit′e Paul Sabatier - Toulouse III, 2013.

    [39] 張曉萌.廣東省17座大中型水庫及其入庫河流N2O分布與水-氣交換通量的研究[D]. 廣州:暨南大學(xué), 2012.

    [40] Wang Y H, Huang H H, Chu C P, et al. A preliminary survey of greenhouse gas emission from three reservoirs in Taiwan [J]. Sustainable Environment Research, 2013,23(3):215-225.

    [41] Chen N, Chen Z, Wu Y, et al. Understanding gaseous nitrogen removal through direct measurement of dissolved N2, and N2O in a subtropical river-reservoir system [J]. Ecological Engineering, 2014,70(5):56-67.

    [42] Zhao Y, Wu B F, Zeng Y. Spatial and temporal patterns of greenhouse gas emissions from Three Gorges Reservoir of China [J]. Biogeosciences, 2013,10(2):1219-1230.

    [43] Xia Y Q, Li Y F, Ti C P, et al. Is indirect N2O emission a significant contributor to the agricultural greenhouse gas budget? A case study of a rice paddy-dominated agricultural watershed in eastern China [J]. Atmospheric Environment, 2013,77(3):943-950.

    [44] 韓 洋,鄭有飛,吳榮軍,等.南京典型水體春季溫室氣體排放特征研究 [J]. 中國環(huán)境科學(xué), 2013,33(8):1360-1371.

    [45] Beaulieu J J, Smolenski R L, Nietch C T,et al. Denitrification alternates between a source and sink of nitrous oxide in the hypolimnion of a thermally stratified reservoir [J]. Limnology & Oceanography, 2015,59(2):495-506.

    [46] Jacinthe P A, Filippelli G M, Tedesco L P, et al. Carbon storage and greenhouse gases emission from a fluvial reservoir in an agricultural landscape [J]. Catena, 2012,94(9):53-63.

    [47] 李紅麗,楊 萌,張明祥,等.玉渡山水庫生長(zhǎng)季溫室氣體排放特征及其影響因素[J]. 生態(tài)學(xué)雜志, 2012,31(2):406-412.

    [48] Yang M, Geng X M, Grace J, et al. N2O fluxes from the littoral zone of a Chinese reservoir [J]. Biogeosciences, 2015,12(15):4711-4723.

    [49] Diem T, Koch S, Schwarzenbach S, et al. Greenhouse gas emissions (CO2, CH4, and N2O) from several perialpine and alpine hydropower reservoirs by diffusion and loss in turbines [J]. Aquatic Sciences, 2012,74:619-635.

    [50] Hendzel L L, Matthews C J D, Venkiteswaran J J, et al. Nitrous oxide fluxes in three experimental boreal forest reservoirs [J]. Environmental Science & Technology, 2005,39(12):4353-4360.

    [51] Huttunen J T, V?is?nen T S, Hellsten S K, et al. Fluxes of CH4, CO2, and N2O in hydroelectric reservoirs Lokka and Porttipahta in the northern boreal zone in Finland [J]. Global Biogeochemical Cycles, 2002,16(1):3-1-3–17.

    [52] Lehner B, Liermann C R, Revenga C, et al. High-resolution mapping of the world's reservoirs and dams for sustainable river-flow management [J]. Frontiers in Ecology & the Environment, 2011,9(9): 494-502.

    [53] Saunders D L, Kalff J. Nitrogen retention in wetlands, lakes and rivers [J]. Hydrobiologia, 2001,443(1-3):205-212.

    [54] Harrison J A, Maranger R J, Alexander R B, et al. The regional and global significance of nitrogen removal in lakes and reservoirs [J]. Biogeochemistry, 2009,93(1-2):143-157.

    [55] Seitzinger S, Harrison J A, B?hlke J K, et al. Denitrification across landscapes and waterscapes: a synthesis [J]. Ecological Applications a Publication of the Ecological Society of America, 2006,16(6):2064- 2090.

    [56] Pi?aochoa E, álvarezcobelas M. Denitrification in aquatic environments: a cross-system analysis [J]. Biogeochemistry, 2006, 81(1):111-130.

    [57] Mengis M, G?chter R, Wehrli B. Sources and sinks of nitrous oxide (N2O) in deep lakes [J]. Biogeochemistry, 1997,38(3):281-301.

    [58] Windolf J, Jeppesen E, Jensen J P, et al. Modelling of seasonal variation in nitrogen retention and in-lake concentration: a four-year mass balance study in 16shallow Danish lakes [J]. Biogeochemistry, 1996,33(1):25-44.

    [59] Knowles R, Chan Y K. Nitrous oxide concentrations in lakes: variations with depth and time [J]. Limnology & Oceanography, 1981, 26(5):855-866.

    [60] Knowles R. Denitrification [J]. Microbiological Reviews, 1982, 46(1):43-70.

    [61] Beaulieu J J, Shuster W D, Rebholz J A. Nitrous oxide emissions from a large, impounded river: the Ohio River [J]. Environmental Science & Technology, 2010,44(19):7527-7533.

    [62] Hu M, Chen D, Dahlgren R A. Modeling nitrous oxide emission from rivers: a global assessment [J]. Global Change Biology, 2016,22(11): 3566-3582.

    [63] Kroeze C, Dumont E, Seitzinger S P. New estimates of global emissions of N2O from rivers and estuaries [J]. Environmental Sciences, 2005,2(2/3):159-165.

    [64] Kroeze C, Dumont E, Seitzinger S P. Future trends in emissions of N2O from rivers and estuaries [J]. Journal of Integrative Environmental Sciences, 2010,7(sup1):71-78.

    [65] Mosier A, Kroeze C, Nevison C, et alClosing the global atmospheric N2O budget: nitrous oxide emissions through the agricultural nitrogen cycle [J]. Nutrient Cycling in Agroecosystems, 1998,52(2):225-248.

    [66] Seitzinger S P, Kroeze C. Global distribution of nitrous oxide production and N inputs in freshwater and coastal marine ecosystems [J]. Global Biogeochemical Cycles, 1998,12(1):93-113.

    [67] Melton J R, Wania R, Hodson E L, et al. Present state of global wetland extent and wetland methane modelling: conclusions from a model inter-comparison project (WETCHIMP) [J]. Geoscientific Model Development, 2013,6(3):617-641.

    [68] Tian H, Chen G, Lu C, et al. Global methane and nitrous oxide emissions from terrestrial ecosystems due to multiple environmental changes [J]. Ecosystem Health & Sustainability, 2015,1(1):1-20.

    [69] Nevison C D, Weiss R F, Erickson D J. Global oceanic emissions of nitrous oxide [J]. Journal of Geophysical Research Oceans, 1995,100(C8):15809-15820.

    [70] Duchemin E, Lucotte M, Canuel R. Comparison of static chamber and thin boundary layer equation methods for measuring greenhouse gas emissions from large water bodies§ [J]. Environmental Science & Technology, 1999,33(2):350-357.

    致謝:感謝University of Oslo的Dr. Yang Hong提供了北歐74個(gè)湖泊詳細(xì)的經(jīng)緯度、面積、深度、TN、TP和N2O通量數(shù)據(jù).

    Analysis of nitrous oxide flux from lakes and reservoirs.

    XIE Yan-hong1,2, ZHANG Mi1,2*, XIAO Wei1,2, WANG Wei1,2, LI Xu-hui1,2, BIAN Hang1,2

    (1.Yale-NUIST Center on Atmospheric Environment, Nanjing University of Information Science & Technology, Nanjing 210044, China;2.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China)., 2018,38(9):3481~3493

    Among inland water bodies, lakes and reservoirs are potential N2O sources and play important roles in nitrogen cycle. Based on the N2O flux data of 479 lakes and 83 reservoirs collected from published literatures, the influences of the latitude, area, depth and trophic status of the lakes and reservoirs were analyzed, and the global mean total N2O flux from lakes and reservoirs was estimated. The results indicated that most lakes and reservoirs acted as N2O sources. The mean N2O flux from lakes [(3.21±5.71)μmol/(m2?d)] was lower than that from reservoirs [(20.82±113.94)μmol/(m2?d)]. At the annual temporal scale, N2O flux from lakes and reservoirs decreased with increasing latitude. The N2O flux from lakes increased with increasing area, but opposite for reservoirs. The N2O flux over deep lakes were larger than those over shallow lakes. The N2O flux over eutrophic lakes and reservoirs were higher with comparison to oligotrophic lakes and reservoirs. The annual mean total N2O flux from lakes and reservoirs were 0.12 and 0.06Tg N/a, respectively, and accounted for 12% and 7% of the total N2O emission from inland waters.

    lakes;reservoirs;N2O flux;impact factors;annual N2O flux

    X511

    A

    1000-6923(2018)09-3481-13

    謝燕紅(1992-),女,廣西桂林人,南京信息工程大學(xué)碩士研究生,主要研究方向是地–氣溫室氣體交換.發(fā)表論文2篇.

    2018-01-10

    國家自然科學(xué)基金資助項(xiàng)目(41575141,41475147);國家自然科學(xué)基金青年基金資助項(xiàng)目(41505005)

    * 責(zé)任作者, 講師, zhangm.80@nuist.edu.cn

    猜你喜歡
    緯度湖泊通量
    冬小麥田N2O通量研究
    你相信嗎?湖泊也可以“生死輪回”
    “害羞”的湖泊
    緯度
    齊魯周刊(2017年29期)2017-08-08 06:28:15
    奇異的湖泊
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    基于時(shí)空緯度的國內(nèi)農(nóng)民工創(chuàng)業(yè)研究
    春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量
    常用緯度差異極值符號(hào)表達(dá)式
    江西省碳通量時(shí)空演變與模型構(gòu)建
    欧美xxxx黑人xx丫x性爽| 18禁黄网站禁片免费观看直播| 99热这里只有是精品50| 老司机福利观看| 亚洲一区高清亚洲精品| 国产69精品久久久久777片| 精品不卡国产一区二区三区| 欧美色欧美亚洲另类二区| av卡一久久| 亚洲高清免费不卡视频| 毛片女人毛片| 国产欧美日韩一区二区精品| 不卡一级毛片| 春色校园在线视频观看| 性插视频无遮挡在线免费观看| 免费大片18禁| 色综合站精品国产| 色噜噜av男人的天堂激情| 村上凉子中文字幕在线| 国产亚洲精品久久久久久毛片| 久久久久久久久久成人| 国产亚洲精品综合一区在线观看| 国产精品永久免费网站| 免费看av在线观看网站| 乱系列少妇在线播放| 美女cb高潮喷水在线观看| 一级毛片久久久久久久久女| 国产高清视频在线观看网站| 欧美bdsm另类| 亚洲国产精品sss在线观看| 亚洲欧美精品自产自拍| 久久久久国产网址| 91狼人影院| 丝袜喷水一区| 久久久国产成人精品二区| 国产真实乱freesex| 精品久久久久久久末码| 黄色视频,在线免费观看| 少妇被粗大猛烈的视频| 欧美一区二区精品小视频在线| 亚洲人成网站在线播放欧美日韩| 不卡视频在线观看欧美| 国产亚洲精品综合一区在线观看| 99热这里只有是精品50| 成年av动漫网址| 亚洲国产精品成人久久小说 | 少妇高潮的动态图| 成人精品一区二区免费| 国产精品99久久久久久久久| 久久久欧美国产精品| 亚洲第一电影网av| 中文字幕人妻熟人妻熟丝袜美| 老师上课跳d突然被开到最大视频| 美女内射精品一级片tv| 国产av不卡久久| 精品久久久久久久久久久久久| 中文字幕精品亚洲无线码一区| 美女 人体艺术 gogo| 女人被狂操c到高潮| 亚洲精品乱码久久久v下载方式| 中文字幕人妻熟人妻熟丝袜美| 久久久久精品国产欧美久久久| 老司机影院成人| 日本黄大片高清| 亚洲中文字幕日韩| 精品人妻一区二区三区麻豆 | 啦啦啦啦在线视频资源| 欧美成人免费av一区二区三区| 色尼玛亚洲综合影院| 黄色欧美视频在线观看| 97超碰精品成人国产| 黄色一级大片看看| 成年女人看的毛片在线观看| 人妻夜夜爽99麻豆av| 免费看a级黄色片| 精品国产三级普通话版| 波多野结衣高清作品| 亚洲性夜色夜夜综合| 亚洲aⅴ乱码一区二区在线播放| 日本欧美国产在线视频| 一区二区三区四区激情视频 | 色5月婷婷丁香| 成年版毛片免费区| 白带黄色成豆腐渣| 亚洲图色成人| 中文字幕熟女人妻在线| 91在线精品国自产拍蜜月| 99riav亚洲国产免费| 天天躁日日操中文字幕| 国产亚洲av嫩草精品影院| 亚洲国产精品成人久久小说 | 亚洲av一区综合| 国产成人a∨麻豆精品| 国产亚洲av嫩草精品影院| 97超级碰碰碰精品色视频在线观看| 日本在线视频免费播放| 99久久九九国产精品国产免费| 国产一区二区三区av在线 | 国产视频内射| 久久鲁丝午夜福利片| 色在线成人网| 91在线精品国自产拍蜜月| 国产精品一区www在线观看| 又爽又黄无遮挡网站| 欧美中文日本在线观看视频| 亚洲图色成人| 免费一级毛片在线播放高清视频| 熟妇人妻久久中文字幕3abv| 日本成人三级电影网站| 一区二区三区高清视频在线| 亚洲一级一片aⅴ在线观看| 日韩高清综合在线| 人妻丰满熟妇av一区二区三区| 精品一区二区三区视频在线观看免费| 国产黄a三级三级三级人| 精品欧美国产一区二区三| 亚洲欧美日韩高清在线视频| 12—13女人毛片做爰片一| АⅤ资源中文在线天堂| 久久精品夜色国产| 国产伦一二天堂av在线观看| 最新中文字幕久久久久| 俺也久久电影网| 亚洲美女搞黄在线观看 | 久久午夜亚洲精品久久| 午夜视频国产福利| av在线观看视频网站免费| 国产一区二区三区在线臀色熟女| 老司机午夜福利在线观看视频| 日本熟妇午夜| 久久久成人免费电影| 久久天躁狠狠躁夜夜2o2o| 国内久久婷婷六月综合欲色啪| 亚洲av二区三区四区| 久久精品国产亚洲av香蕉五月| 夜夜爽天天搞| 淫妇啪啪啪对白视频| 国产成人精品久久久久久| 亚洲熟妇中文字幕五十中出| 乱系列少妇在线播放| 久久久欧美国产精品| 蜜臀久久99精品久久宅男| 欧美性猛交╳xxx乱大交人| 99热6这里只有精品| 在线国产一区二区在线| eeuss影院久久| 六月丁香七月| 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 国产精品爽爽va在线观看网站| 夜夜夜夜夜久久久久| 久久精品国产亚洲av涩爱 | 亚洲欧美日韩卡通动漫| 国产色爽女视频免费观看| 国内精品美女久久久久久| 天堂√8在线中文| 日韩 亚洲 欧美在线| 亚洲精华国产精华液的使用体验 | 一级毛片久久久久久久久女| 国产视频内射| 禁无遮挡网站| 1000部很黄的大片| 亚洲精品成人久久久久久| 在线免费观看的www视频| 嫩草影院新地址| 精品人妻熟女av久视频| 国产久久久一区二区三区| 97超碰精品成人国产| 色哟哟哟哟哟哟| 精品国产三级普通话版| 亚洲欧美日韩高清专用| 欧美日韩国产亚洲二区| 免费av毛片视频| 99久国产av精品国产电影| 人妻久久中文字幕网| 18+在线观看网站| 欧美最新免费一区二区三区| 99久久久亚洲精品蜜臀av| 国产一区二区在线观看日韩| 长腿黑丝高跟| 嫩草影院入口| 国产亚洲欧美98| 中文在线观看免费www的网站| 欧美日本视频| 亚洲国产高清在线一区二区三| 亚洲国产精品sss在线观看| 国产av不卡久久| 岛国在线免费视频观看| 午夜激情欧美在线| 男女那种视频在线观看| 免费看光身美女| 日韩精品青青久久久久久| 亚洲自偷自拍三级| 亚洲欧美日韩卡通动漫| 日本黄色片子视频| 欧美极品一区二区三区四区| 亚洲久久久久久中文字幕| 日本在线视频免费播放| 欧美+日韩+精品| 91午夜精品亚洲一区二区三区| 男女做爰动态图高潮gif福利片| 国产视频内射| 成人美女网站在线观看视频| 亚洲经典国产精华液单| 午夜福利18| 男女边吃奶边做爰视频| 高清毛片免费观看视频网站| 如何舔出高潮| 日韩欧美在线乱码| 两个人的视频大全免费| 国产精品久久久久久精品电影| 美女xxoo啪啪120秒动态图| 国产淫片久久久久久久久| 亚洲国产精品国产精品| 日本熟妇午夜| av中文乱码字幕在线| 精品久久久久久久久av| 日本一二三区视频观看| 国产在视频线在精品| 精品少妇黑人巨大在线播放 | 国产精品一二三区在线看| 一级毛片aaaaaa免费看小| 国产亚洲精品久久久久久毛片| 最新中文字幕久久久久| 日韩成人av中文字幕在线观看 | 久久久午夜欧美精品| 哪里可以看免费的av片| 日日摸夜夜添夜夜爱| 欧美丝袜亚洲另类| 看片在线看免费视频| 国产极品精品免费视频能看的| 久久精品国产亚洲av香蕉五月| 两个人视频免费观看高清| 国产亚洲欧美98| 人人妻人人看人人澡| 久久热精品热| 国产成人福利小说| 国产片特级美女逼逼视频| av在线老鸭窝| 熟女电影av网| 别揉我奶头~嗯~啊~动态视频| 精品无人区乱码1区二区| 看黄色毛片网站| 亚洲图色成人| 亚洲三级黄色毛片| 亚洲一区二区三区色噜噜| 国产精品久久久久久久电影| 久久久久久伊人网av| 99国产精品一区二区蜜桃av| 亚洲一区高清亚洲精品| 欧美日韩国产亚洲二区| 中文字幕av在线有码专区| 国产一级毛片七仙女欲春2| 99热网站在线观看| 毛片女人毛片| 黄色欧美视频在线观看| 免费在线观看影片大全网站| 亚洲av.av天堂| 91在线精品国自产拍蜜月| 成人鲁丝片一二三区免费| 色5月婷婷丁香| 91av网一区二区| 99九九线精品视频在线观看视频| 99热全是精品| 免费av观看视频| 麻豆久久精品国产亚洲av| 久久精品国产亚洲av涩爱 | 熟女电影av网| 日产精品乱码卡一卡2卡三| 国产成人福利小说| 午夜爱爱视频在线播放| 国内精品久久久久精免费| 亚洲国产高清在线一区二区三| 国产高清视频在线播放一区| 国产成年人精品一区二区| 99久久精品国产国产毛片| 真实男女啪啪啪动态图| 欧美三级亚洲精品| 亚洲自偷自拍三级| 欧美色视频一区免费| 一个人看的www免费观看视频| 日日啪夜夜撸| 最近视频中文字幕2019在线8| 久久人人爽人人片av| 午夜视频国产福利| av天堂在线播放| av在线天堂中文字幕| 成年免费大片在线观看| 国产成人一区二区在线| 深夜a级毛片| 两个人视频免费观看高清| 欧洲精品卡2卡3卡4卡5卡区| 国产毛片a区久久久久| 精华霜和精华液先用哪个| 男女之事视频高清在线观看| 精品午夜福利视频在线观看一区| 不卡视频在线观看欧美| 久久精品国产自在天天线| 99久久无色码亚洲精品果冻| 日本三级黄在线观看| 在线a可以看的网站| 草草在线视频免费看| 91久久精品国产一区二区成人| 免费看av在线观看网站| 国产亚洲精品综合一区在线观看| 日韩精品青青久久久久久| 麻豆国产97在线/欧美| 久久6这里有精品| a级毛色黄片| 国产成人freesex在线 | 国产欧美日韩一区二区精品| 亚洲人成网站高清观看| 在线a可以看的网站| 亚洲精华国产精华液的使用体验 | 一区福利在线观看| 91久久精品国产一区二区成人| 高清毛片免费看| 伦理电影大哥的女人| 男女做爰动态图高潮gif福利片| 俺也久久电影网| 啦啦啦观看免费观看视频高清| 成人av一区二区三区在线看| 97在线视频观看| 国产日本99.免费观看| 欧美+日韩+精品| 色噜噜av男人的天堂激情| 美女内射精品一级片tv| 日本一本二区三区精品| 亚洲第一区二区三区不卡| 亚洲精品粉嫩美女一区| 久久久久性生活片| 国产午夜精品论理片| 日本与韩国留学比较| 免费看美女性在线毛片视频| 国产亚洲91精品色在线| 免费人成视频x8x8入口观看| 最后的刺客免费高清国语| 免费人成在线观看视频色| 午夜福利在线观看免费完整高清在 | 国产精品野战在线观看| 成人亚洲欧美一区二区av| 99久久精品热视频| 亚洲四区av| 国产在线男女| 午夜a级毛片| 一本一本综合久久| 黄色欧美视频在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲内射少妇av| 给我免费播放毛片高清在线观看| 99热精品在线国产| 国产精品久久久久久亚洲av鲁大| 麻豆国产97在线/欧美| 国产视频内射| 大香蕉久久网| 亚洲一级一片aⅴ在线观看| 国产精品一区www在线观看| 国产伦在线观看视频一区| 午夜a级毛片| 高清毛片免费观看视频网站| 啦啦啦韩国在线观看视频| 日本五十路高清| 一区二区三区四区激情视频 | 桃色一区二区三区在线观看| 一个人看视频在线观看www免费| 看非洲黑人一级黄片| 插阴视频在线观看视频| 国产激情偷乱视频一区二区| 亚洲国产精品sss在线观看| 插阴视频在线观看视频| 久久久久久久久久久丰满| 精品一区二区三区视频在线观看免费| 免费高清视频大片| 国产成人aa在线观看| 国产欧美日韩精品一区二区| 国产精品福利在线免费观看| 午夜亚洲福利在线播放| 特大巨黑吊av在线直播| 欧美最新免费一区二区三区| 亚洲熟妇熟女久久| 中文字幕精品亚洲无线码一区| 亚洲,欧美,日韩| 日日撸夜夜添| 久久精品夜色国产| 国产精品国产高清国产av| 97超视频在线观看视频| 狂野欧美白嫩少妇大欣赏| 少妇熟女aⅴ在线视频| 国产精品一及| 中文字幕精品亚洲无线码一区| 国产免费男女视频| 97超视频在线观看视频| 别揉我奶头~嗯~啊~动态视频| 91在线观看av| 男女做爰动态图高潮gif福利片| 亚洲国产精品成人久久小说 | 精品福利观看| 变态另类丝袜制服| 欧美色欧美亚洲另类二区| 变态另类成人亚洲欧美熟女| 亚洲精品粉嫩美女一区| 亚洲欧美精品综合久久99| 亚洲国产高清在线一区二区三| 激情 狠狠 欧美| 久久精品国产鲁丝片午夜精品| 老熟妇仑乱视频hdxx| 狠狠狠狠99中文字幕| 男女那种视频在线观看| 国产精品福利在线免费观看| 欧美成人一区二区免费高清观看| 小蜜桃在线观看免费完整版高清| 色噜噜av男人的天堂激情| 久久久国产成人免费| 黄色配什么色好看| АⅤ资源中文在线天堂| 精品午夜福利视频在线观看一区| 一个人看的www免费观看视频| 99在线人妻在线中文字幕| 欧美丝袜亚洲另类| 免费观看在线日韩| 观看美女的网站| 国产高清三级在线| 婷婷色综合大香蕉| av女优亚洲男人天堂| 亚洲国产高清在线一区二区三| 高清毛片免费观看视频网站| 亚洲第一区二区三区不卡| 美女cb高潮喷水在线观看| 你懂的网址亚洲精品在线观看 | 九九热线精品视视频播放| 日日摸夜夜添夜夜添小说| 国产成人福利小说| 人人妻人人看人人澡| 亚洲精品国产成人久久av| 久久久欧美国产精品| 噜噜噜噜噜久久久久久91| 18+在线观看网站| 狂野欧美激情性xxxx在线观看| 午夜免费男女啪啪视频观看 | 看片在线看免费视频| 国产精品久久久久久精品电影| 97超级碰碰碰精品色视频在线观看| 亚洲av二区三区四区| 国国产精品蜜臀av免费| 九九爱精品视频在线观看| 日韩人妻高清精品专区| 精品久久久久久久末码| 国产精品不卡视频一区二区| 色吧在线观看| 女人被狂操c到高潮| 亚洲精品在线观看二区| 欧美中文日本在线观看视频| 欧美日韩精品成人综合77777| 色av中文字幕| 色综合亚洲欧美另类图片| 亚洲欧美成人精品一区二区| 五月伊人婷婷丁香| 91精品国产九色| www日本黄色视频网| 午夜福利成人在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 两个人视频免费观看高清| 国产精品99久久久久久久久| 欧美在线一区亚洲| 99视频精品全部免费 在线| 性色avwww在线观看| 日本-黄色视频高清免费观看| 亚洲欧美日韩高清专用| 国产高清视频在线播放一区| 久久人人爽人人片av| 午夜影院日韩av| 中出人妻视频一区二区| 久久婷婷人人爽人人干人人爱| 中文在线观看免费www的网站| 黑人高潮一二区| 99精品在免费线老司机午夜| 丰满乱子伦码专区| 欧美激情国产日韩精品一区| 免费搜索国产男女视频| 99视频精品全部免费 在线| 日本欧美国产在线视频| 97碰自拍视频| 久久99热6这里只有精品| 精华霜和精华液先用哪个| 蜜桃久久精品国产亚洲av| 国产精品一及| 国产综合懂色| 久久久久久久久大av| 中文字幕免费在线视频6| 亚洲av中文av极速乱| 久久久久九九精品影院| 亚洲最大成人中文| 中国国产av一级| 欧美最黄视频在线播放免费| 久久久久久伊人网av| 午夜福利成人在线免费观看| 在线观看av片永久免费下载| 一个人免费在线观看电影| 国产色婷婷99| 国产精品久久久久久精品电影| 欧美中文日本在线观看视频| 色哟哟哟哟哟哟| 亚洲婷婷狠狠爱综合网| 久久精品国产亚洲av涩爱 | 亚洲av一区综合| 亚洲人成网站在线观看播放| 久久中文看片网| 免费黄网站久久成人精品| 亚洲真实伦在线观看| 国产精品嫩草影院av在线观看| 不卡视频在线观看欧美| 日韩强制内射视频| 精品一区二区三区av网在线观看| 国内揄拍国产精品人妻在线| 美女高潮的动态| 97碰自拍视频| 老司机影院成人| 国产熟女欧美一区二区| 69av精品久久久久久| 欧美日韩国产亚洲二区| 久久精品综合一区二区三区| 亚洲欧美日韩高清在线视频| 一进一出好大好爽视频| 校园春色视频在线观看| 色在线成人网| 韩国av在线不卡| 悠悠久久av| 亚洲四区av| 国产老妇女一区| 欧美区成人在线视频| 搡老妇女老女人老熟妇| 欧美高清成人免费视频www| 日本精品一区二区三区蜜桃| 99久久成人亚洲精品观看| 国产不卡一卡二| 男人狂女人下面高潮的视频| 久久久午夜欧美精品| 可以在线观看的亚洲视频| 国产成人精品久久久久久| 一本久久中文字幕| .国产精品久久| 99久久无色码亚洲精品果冻| 日韩,欧美,国产一区二区三区 | av.在线天堂| 中国美女看黄片| 99久久九九国产精品国产免费| 校园人妻丝袜中文字幕| 亚洲性久久影院| 露出奶头的视频| 看片在线看免费视频| 少妇的逼好多水| 免费观看精品视频网站| 国产精品人妻久久久影院| 久久人妻av系列| 日韩制服骚丝袜av| 亚洲人成网站在线播| www日本黄色视频网| 亚洲精品456在线播放app| 亚洲精品一区av在线观看| 少妇熟女欧美另类| 日日啪夜夜撸| 欧美色欧美亚洲另类二区| 日韩精品青青久久久久久| 最近2019中文字幕mv第一页| 91久久精品电影网| 99久久精品国产国产毛片| 亚洲中文字幕日韩| 高清毛片免费观看视频网站| 国产亚洲精品久久久久久毛片| 18禁裸乳无遮挡免费网站照片| 黄色配什么色好看| 久久久久久久久久成人| 日日撸夜夜添| 精品久久久久久成人av| 免费在线观看影片大全网站| 麻豆成人午夜福利视频| 高清午夜精品一区二区三区 | 高清毛片免费观看视频网站| av福利片在线观看| 日韩高清综合在线| 欧美日韩一区二区视频在线观看视频在线 | 午夜影院日韩av| 国产亚洲av嫩草精品影院| 国产精品1区2区在线观看.| 免费人成视频x8x8入口观看| 免费看光身美女| 男女之事视频高清在线观看| 国产精品99久久久久久久久| 少妇裸体淫交视频免费看高清| 内射极品少妇av片p| 69人妻影院| 白带黄色成豆腐渣| av卡一久久| 亚洲国产精品合色在线| 一区二区三区免费毛片| 国产三级中文精品| 亚洲精品乱码久久久v下载方式| 国产一区二区亚洲精品在线观看| 欧美不卡视频在线免费观看| 国产人妻一区二区三区在| 91av网一区二区| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久久久大av| 久久人妻av系列| 免费不卡的大黄色大毛片视频在线观看 | 熟女人妻精品中文字幕| 国产爱豆传媒在线观看| 99热6这里只有精品| 村上凉子中文字幕在线| 如何舔出高潮| 久久亚洲精品不卡| 级片在线观看| 五月玫瑰六月丁香| 免费一级毛片在线播放高清视频| 国产精品日韩av在线免费观看|