• <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)建
    国产伦理片在线播放av一区| 亚洲成人一二三区av| 国产 一区 欧美 日韩| 99re6热这里在线精品视频| 国内精品美女久久久久久| 3wmmmm亚洲av在线观看| 久久久精品94久久精品| 精品一区二区三区视频在线| 国产在线一区二区三区精| 欧美激情久久久久久爽电影| 日韩精品青青久久久久久| 亚州av有码| 午夜老司机福利剧场| 蜜桃久久精品国产亚洲av| 国产探花极品一区二区| 国产三级在线视频| 激情五月婷婷亚洲| 亚洲欧美中文字幕日韩二区| 亚洲人成网站高清观看| 免费无遮挡裸体视频| 久久国产乱子免费精品| 久久综合国产亚洲精品| 久久久久网色| 国产成人精品一,二区| 美女脱内裤让男人舔精品视频| 亚洲精品视频女| 久久精品国产自在天天线| 一级毛片 在线播放| 九九久久精品国产亚洲av麻豆| 国产伦在线观看视频一区| 亚洲欧美精品自产自拍| 免费观看a级毛片全部| 亚洲18禁久久av| 亚洲av一区综合| 能在线免费看毛片的网站| 亚洲国产色片| 成年免费大片在线观看| 丝袜喷水一区| 欧美人与善性xxx| 99re6热这里在线精品视频| 国产白丝娇喘喷水9色精品| 黄色配什么色好看| 欧美最新免费一区二区三区| 免费av不卡在线播放| 免费av不卡在线播放| 成年版毛片免费区| 老司机影院毛片| 高清午夜精品一区二区三区| 欧美xxⅹ黑人| 天堂中文最新版在线下载 | 免费看av在线观看网站| 别揉我奶头 嗯啊视频| 色综合色国产| 国产一区二区在线观看日韩| 一区二区三区四区激情视频| 成人亚洲精品av一区二区| 久久鲁丝午夜福利片| 亚洲熟女精品中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 久久综合国产亚洲精品| 婷婷色综合www| 亚洲精品中文字幕在线视频 | 国产成人福利小说| 成人性生交大片免费视频hd| 国产欧美日韩精品一区二区| 日韩av不卡免费在线播放| 亚洲欧洲日产国产| 亚洲人与动物交配视频| 亚洲人成网站在线观看播放| 亚洲精品自拍成人| 搡老乐熟女国产| 免费黄色在线免费观看| 精品久久久噜噜| 国内精品美女久久久久久| 99热网站在线观看| 观看美女的网站| 欧美激情久久久久久爽电影| 国产男女超爽视频在线观看| 国产成人精品久久久久久| 免费观看a级毛片全部| 国产 亚洲一区二区三区 | 婷婷六月久久综合丁香| av卡一久久| 久久精品国产亚洲av天美| h日本视频在线播放| 亚洲av成人精品一二三区| av国产久精品久网站免费入址| 亚洲精品一区蜜桃| 插逼视频在线观看| 国产精品人妻久久久影院| 国产精品久久久久久精品电影小说 | 97超视频在线观看视频| 午夜精品在线福利| 欧美潮喷喷水| 亚洲精品乱码久久久久久按摩| 精品一区在线观看国产| 久久99热6这里只有精品| 男女那种视频在线观看| 全区人妻精品视频| 亚洲熟妇中文字幕五十中出| 水蜜桃什么品种好| 久久久久精品久久久久真实原创| 热99在线观看视频| 观看免费一级毛片| 波多野结衣巨乳人妻| 亚洲人成网站高清观看| 亚洲成色77777| 日本与韩国留学比较| 一边亲一边摸免费视频| 国模一区二区三区四区视频| 又爽又黄无遮挡网站| 成人午夜精彩视频在线观看| 人妻制服诱惑在线中文字幕| 欧美三级亚洲精品| 国产精品国产三级专区第一集| 国产不卡一卡二| 一级毛片 在线播放| 成人漫画全彩无遮挡| 69av精品久久久久久| 亚洲av免费高清在线观看| 亚洲精品中文字幕在线视频 | 18禁动态无遮挡网站| 成人午夜精彩视频在线观看| 国产午夜精品久久久久久一区二区三区| 69人妻影院| 精品久久久久久久久久久久久| 天堂中文最新版在线下载 | 一区二区三区四区激情视频| 天堂影院成人在线观看| 日韩欧美 国产精品| 中国国产av一级| 特大巨黑吊av在线直播| 高清午夜精品一区二区三区| 麻豆国产97在线/欧美| 成年免费大片在线观看| 伦理电影大哥的女人| 久久99蜜桃精品久久| 91在线精品国自产拍蜜月| 国产在线一区二区三区精| 寂寞人妻少妇视频99o| a级毛片免费高清观看在线播放| 成人二区视频| 国产男女超爽视频在线观看| 日韩av免费高清视频| 久久久久性生活片| 黄色配什么色好看| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕亚洲精品专区| av福利片在线观看| 午夜激情久久久久久久| 久久99热6这里只有精品| 免费无遮挡裸体视频| 五月玫瑰六月丁香| 99热全是精品| 精品酒店卫生间| 国产一区亚洲一区在线观看| 亚洲精品第二区| 99热全是精品| 色综合色国产| 又粗又硬又长又爽又黄的视频| 麻豆精品久久久久久蜜桃| 九草在线视频观看| 亚洲精品乱久久久久久| 亚洲精品第二区| 99久久九九国产精品国产免费| 色吧在线观看| 三级毛片av免费| 国产成人精品福利久久| xxx大片免费视频| 久久人人爽人人片av| 国产探花极品一区二区| 最近最新中文字幕免费大全7| 少妇人妻精品综合一区二区| www.av在线官网国产| 亚洲欧美成人综合另类久久久| 国产黄a三级三级三级人| 国产午夜福利久久久久久| 亚洲av免费在线观看| 精品久久久精品久久久| 国产亚洲av片在线观看秒播厂 | 欧美日韩在线观看h| 亚洲精品久久午夜乱码| 欧美极品一区二区三区四区| 99久久九九国产精品国产免费| 免费播放大片免费观看视频在线观看| 亚洲av一区综合| 国产v大片淫在线免费观看| 国产毛片a区久久久久| 国内少妇人妻偷人精品xxx网站| 日韩大片免费观看网站| 国语对白做爰xxxⅹ性视频网站| 久久精品国产亚洲av天美| av福利片在线观看| 少妇人妻一区二区三区视频| 国产高清国产精品国产三级 | 久久久亚洲精品成人影院| 看非洲黑人一级黄片| 国产精品日韩av在线免费观看| 一级毛片我不卡| 国产精品久久久久久久久免| 免费av观看视频| 日韩欧美 国产精品| 精品一区二区免费观看| 女人十人毛片免费观看3o分钟| 国产高清不卡午夜福利| 午夜福利在线在线| 亚洲精品国产成人久久av| 国产精品嫩草影院av在线观看| 国产亚洲av片在线观看秒播厂 | 老司机影院成人| 精品少妇黑人巨大在线播放| 中文精品一卡2卡3卡4更新| 干丝袜人妻中文字幕| av.在线天堂| 深爱激情五月婷婷| 又爽又黄a免费视频| 午夜日本视频在线| 国产精品久久视频播放| 午夜免费男女啪啪视频观看| eeuss影院久久| av女优亚洲男人天堂| 亚洲最大成人中文| 中文精品一卡2卡3卡4更新| 晚上一个人看的免费电影| 永久网站在线| 两个人视频免费观看高清| 联通29元200g的流量卡| 看非洲黑人一级黄片| 成人性生交大片免费视频hd| 久久久久久久久久久免费av| 综合色av麻豆| 免费黄频网站在线观看国产| 精品人妻视频免费看| 99久国产av精品国产电影| 久久精品国产亚洲av涩爱| 久久久国产一区二区| 亚洲av电影不卡..在线观看| 久久热精品热| 一区二区三区免费毛片| 尾随美女入室| 亚洲国产欧美在线一区| 久久久久久久大尺度免费视频| 一区二区三区高清视频在线| 欧美成人a在线观看| 日韩在线高清观看一区二区三区| 97超视频在线观看视频| 女的被弄到高潮叫床怎么办| 亚洲精品自拍成人| 午夜福利视频1000在线观看| 伊人久久精品亚洲午夜| 亚洲欧洲日产国产| 精品久久久久久久久久久久久| 中文欧美无线码| 久久久久网色| 国产高清不卡午夜福利| 男人舔奶头视频| 亚洲欧美日韩卡通动漫| 波多野结衣巨乳人妻| 久久久久久伊人网av| 成人鲁丝片一二三区免费| 边亲边吃奶的免费视频| 日韩欧美精品免费久久| 日本熟妇午夜| 91久久精品国产一区二区三区| 极品少妇高潮喷水抽搐| 日韩欧美精品v在线| 国产精品一区www在线观看| 成人亚洲欧美一区二区av| 久久综合国产亚洲精品| a级一级毛片免费在线观看| 99热网站在线观看| 欧美日韩亚洲高清精品| 日产精品乱码卡一卡2卡三| 免费观看的影片在线观看| 天天一区二区日本电影三级| 亚洲最大成人av| 国产精品麻豆人妻色哟哟久久 | 久久久久久久亚洲中文字幕| 26uuu在线亚洲综合色| 少妇熟女欧美另类| 亚洲精品aⅴ在线观看| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩东京热| 国产视频首页在线观看| 99热这里只有是精品50| 日本wwww免费看| 一级av片app| 91午夜精品亚洲一区二区三区| 春色校园在线视频观看| 大香蕉久久网| 一个人免费在线观看电影| 秋霞伦理黄片| 免费不卡的大黄色大毛片视频在线观看 | 永久免费av网站大全| 六月丁香七月| 亚洲最大成人中文| 国产精品爽爽va在线观看网站| 日本一本二区三区精品| 精品一区二区免费观看| 中文字幕制服av| 免费高清在线观看视频在线观看| 在线免费十八禁| 肉色欧美久久久久久久蜜桃 | 特大巨黑吊av在线直播| 日本与韩国留学比较| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区黑人 | 91午夜精品亚洲一区二区三区| 亚洲最大成人手机在线| 人人妻人人澡人人爽人人夜夜 | 免费看不卡的av| 天堂√8在线中文| 久久久a久久爽久久v久久| 草草在线视频免费看| 国产一级毛片在线| freevideosex欧美| 亚洲人与动物交配视频| 两个人的视频大全免费| 婷婷色综合www| 夫妻性生交免费视频一级片| 91精品伊人久久大香线蕉| 熟妇人妻久久中文字幕3abv| 日韩强制内射视频| 波多野结衣巨乳人妻| 高清av免费在线| 久久久久久久亚洲中文字幕| 日日撸夜夜添| 亚洲av在线观看美女高潮| 日韩伦理黄色片| 嘟嘟电影网在线观看| 午夜福利在线在线| 国产中年淑女户外野战色| 边亲边吃奶的免费视频| 成年人午夜在线观看视频 | 毛片一级片免费看久久久久| 岛国毛片在线播放| 最近中文字幕高清免费大全6| 色5月婷婷丁香| 国产高清国产精品国产三级 | 成人美女网站在线观看视频| 色吧在线观看| 日韩欧美国产在线观看| 免费av毛片视频| 亚洲美女搞黄在线观看| 黄片wwwwww| 国产精品久久久久久久久免| 精品人妻视频免费看| 色视频www国产| 亚洲精品自拍成人| 国产中年淑女户外野战色| 中文字幕av成人在线电影| 亚洲国产精品成人久久小说| 国产女主播在线喷水免费视频网站 | 人妻制服诱惑在线中文字幕| 日韩成人av中文字幕在线观看| kizo精华| 国产精品精品国产色婷婷| 波多野结衣巨乳人妻| 国产一区二区亚洲精品在线观看| 天堂av国产一区二区熟女人妻| 亚洲内射少妇av| 久久久成人免费电影| 久久久久久久大尺度免费视频| 亚洲成人中文字幕在线播放| 欧美成人精品欧美一级黄| av播播在线观看一区| 国产视频内射| av播播在线观看一区| 日本黄色片子视频| 国产午夜精品久久久久久一区二区三区| 麻豆精品久久久久久蜜桃| 亚洲精品久久久久久婷婷小说| 非洲黑人性xxxx精品又粗又长| 日韩精品青青久久久久久| 午夜激情福利司机影院| 免费高清在线观看视频在线观看| 日本黄色片子视频| 久久国产乱子免费精品| 中文字幕亚洲精品专区| 婷婷色综合大香蕉| 女的被弄到高潮叫床怎么办| 26uuu在线亚洲综合色| 超碰97精品在线观看| 国产成人精品久久久久久| 欧美性感艳星| 男女那种视频在线观看| 在线观看一区二区三区| 亚洲人与动物交配视频| 国产黄片美女视频| 97在线视频观看| 成人亚洲精品av一区二区| 成人毛片60女人毛片免费| 最近中文字幕高清免费大全6| 免费少妇av软件| 久久精品国产自在天天线| 乱系列少妇在线播放| 看十八女毛片水多多多| 成人特级av手机在线观看| 国产在线一区二区三区精| 国产av在哪里看| 18禁裸乳无遮挡免费网站照片| 一个人看视频在线观看www免费| 色综合亚洲欧美另类图片| 国产精品国产三级国产专区5o| 夜夜看夜夜爽夜夜摸| 视频中文字幕在线观看| av专区在线播放| 日韩精品有码人妻一区| 三级国产精品片| 91aial.com中文字幕在线观看| 99久国产av精品| eeuss影院久久| 亚洲精品乱码久久久v下载方式| 婷婷六月久久综合丁香| 丝袜喷水一区| 精品久久久久久久久av| 欧美高清成人免费视频www| 免费大片18禁| 成年版毛片免费区| 亚洲精品国产av成人精品| 欧美三级亚洲精品| 日韩 亚洲 欧美在线| 国产av在哪里看| 91精品国产九色| 亚洲av成人精品一二三区| 乱码一卡2卡4卡精品| 日韩三级伦理在线观看| 秋霞伦理黄片| 偷拍熟女少妇极品色| 看免费成人av毛片| 最后的刺客免费高清国语| 精品一区二区三区人妻视频| 亚洲精品成人av观看孕妇| 中文欧美无线码| 精品99又大又爽又粗少妇毛片| 精品一区二区免费观看| 精品一区二区三卡| 乱人视频在线观看| 天天一区二区日本电影三级| 成人无遮挡网站| 我的老师免费观看完整版| 精品人妻熟女av久视频| 亚洲成人精品中文字幕电影| 久久精品国产亚洲av天美| 爱豆传媒免费全集在线观看| 少妇的逼水好多| 天美传媒精品一区二区| 亚洲经典国产精华液单| 五月天丁香电影| 神马国产精品三级电影在线观看| 99久国产av精品| 免费av毛片视频| 国产淫语在线视频| 免费无遮挡裸体视频| 国产综合精华液| 免费无遮挡裸体视频| 美女cb高潮喷水在线观看| 成人毛片a级毛片在线播放| 麻豆精品久久久久久蜜桃| 久久精品国产自在天天线| 乱码一卡2卡4卡精品| 麻豆国产97在线/欧美| 亚洲性久久影院| 欧美激情在线99| 日本免费在线观看一区| 免费黄色在线免费观看| 日韩精品有码人妻一区| 免费观看精品视频网站| 观看美女的网站| 欧美日韩综合久久久久久| 国产人妻一区二区三区在| 亚洲自拍偷在线| 精品亚洲乱码少妇综合久久| 色哟哟·www| 国产午夜精品论理片| 欧美成人a在线观看| 国产免费一级a男人的天堂| 韩国高清视频一区二区三区| 色综合亚洲欧美另类图片| 成年女人在线观看亚洲视频 | 超碰97精品在线观看| 色综合亚洲欧美另类图片| 可以在线观看毛片的网站| 国产精品女同一区二区软件| 人人妻人人澡人人爽人人夜夜 | 国产av码专区亚洲av| 精品人妻熟女av久视频| 看非洲黑人一级黄片| 99热这里只有是精品50| 深爱激情五月婷婷| 在线天堂最新版资源| 午夜激情福利司机影院| 草草在线视频免费看| 欧美潮喷喷水| 久久6这里有精品| 亚洲欧美一区二区三区黑人 | 日韩大片免费观看网站| 亚洲精华国产精华液的使用体验| 欧美成人a在线观看| 边亲边吃奶的免费视频| 18禁动态无遮挡网站| 亚洲自偷自拍三级| 亚洲性久久影院| 精品久久久精品久久久| 亚洲精品久久久久久婷婷小说| 高清毛片免费看| 又粗又硬又长又爽又黄的视频| 国产伦一二天堂av在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲av中文av极速乱| 亚洲精品乱码久久久v下载方式| 亚洲内射少妇av| 日韩,欧美,国产一区二区三区| 午夜免费激情av| 婷婷色麻豆天堂久久| 别揉我奶头 嗯啊视频| 男女边吃奶边做爰视频| www.av在线官网国产| 成年人午夜在线观看视频 | 秋霞在线观看毛片| 国产伦精品一区二区三区视频9| 可以在线观看毛片的网站| 国产av国产精品国产| 人妻夜夜爽99麻豆av| 99久国产av精品| 亚洲av成人av| 国产成人精品福利久久| 成人性生交大片免费视频hd| 男女下面进入的视频免费午夜| 精品国内亚洲2022精品成人| 日日摸夜夜添夜夜爱| 一区二区三区高清视频在线| 国产91av在线免费观看| 国产一区二区亚洲精品在线观看| 亚洲av二区三区四区| 亚洲人与动物交配视频| 色综合亚洲欧美另类图片| 18+在线观看网站| videossex国产| 国产精品国产三级专区第一集| 精品人妻一区二区三区麻豆| 搡老乐熟女国产| 国产av码专区亚洲av| 国产亚洲一区二区精品| 亚洲av免费在线观看| 久久精品国产亚洲av涩爱| 在现免费观看毛片| 少妇被粗大猛烈的视频| 午夜免费观看性视频| 亚洲高清免费不卡视频| 国产精品国产三级专区第一集| 国产亚洲av片在线观看秒播厂 | 80岁老熟妇乱子伦牲交| 六月丁香七月| 亚洲真实伦在线观看| 十八禁国产超污无遮挡网站| 搡老妇女老女人老熟妇| 中文天堂在线官网| av线在线观看网站| 我要看日韩黄色一级片| 在线天堂最新版资源| 亚洲精品中文字幕在线视频 | 久久热精品热| 亚洲精品一区蜜桃| 激情 狠狠 欧美| 肉色欧美久久久久久久蜜桃 | 色综合色国产| 少妇人妻精品综合一区二区| 七月丁香在线播放| 亚洲在久久综合| 中文乱码字字幕精品一区二区三区 | 国产美女午夜福利| 热99在线观看视频| 午夜视频国产福利| 日韩欧美精品免费久久| 亚洲欧美日韩无卡精品| 亚洲丝袜综合中文字幕| 一区二区三区免费毛片| 亚洲自拍偷在线| 青春草视频在线免费观看| 国产免费视频播放在线视频 | 99久久人妻综合| 18禁在线播放成人免费| 我的老师免费观看完整版| 欧美成人午夜免费资源| av在线天堂中文字幕| 亚洲av日韩在线播放| 蜜桃久久精品国产亚洲av| 最近视频中文字幕2019在线8| 午夜福利在线观看免费完整高清在| 久久久久久久大尺度免费视频| 丰满乱子伦码专区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美中文字幕日韩二区| 九九久久精品国产亚洲av麻豆| 亚洲性久久影院| 成人性生交大片免费视频hd| 久久久久久久久大av| 亚洲在久久综合| 国产黄色视频一区二区在线观看| 最近最新中文字幕免费大全7| 久久久久国产网址| 免费黄色在线免费观看| 黑人高潮一二区| 成年版毛片免费区| 韩国高清视频一区二区三区| 高清av免费在线| 色网站视频免费| 99久国产av精品| 国产精品爽爽va在线观看网站| 亚洲欧美精品自产自拍| 亚洲人成网站在线播| 欧美性感艳星| ponron亚洲| 丝袜喷水一区|