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

    太湖上游丘陵區(qū)典型水塘反硝化潛力變化特征及其影響因素*

    2022-11-09 03:35:02王慧亮張汪壽李恒鵬夏天宇陳永娟
    湖泊科學(xué) 2022年6期
    關(guān)鍵詞:過(guò)飽和水氣水塘

    王慧亮,曹 恒,,張汪壽,李恒鵬,何 鵬,夏天宇,,陳永娟,3**

    (1:鄭州大學(xué)水利科學(xué)與工程學(xué)院,鄭州 450001) (2:中國(guó)科學(xué)院南京地理與湖泊研究所,中國(guó)科學(xué)院流域地理學(xué)重點(diǎn)實(shí)驗(yàn)室,南京 210008) (3:滁州學(xué)院土木與建筑工程學(xué)院,滁州 239099)

    隨著人口的增長(zhǎng)和糧食需求的增加,在追求高產(chǎn)量的同時(shí),更多的農(nóng)業(yè)化肥被投入使用. 由于其利用率較低,全球超過(guò)一半的施用化肥通過(guò)水塘等匯入河流、湖泊和海洋,并造成嚴(yán)重的生態(tài)環(huán)境問(wèn)題[1]. 丘陵區(qū)的農(nóng)業(yè)發(fā)展模式主要以“坡地-溝塘”為主[2],近年來(lái),丘陵區(qū)的土地開(kāi)發(fā)呈現(xiàn)快速增長(zhǎng)的趨勢(shì)[3],其中坡地種植以茶樹(shù)為主,施氮量可達(dá)959 kg/(hm2·a),約是水田的2倍[4],造成了大量的氮素流失,嚴(yán)重威脅水質(zhì)安全,并導(dǎo)致水體富營(yíng)養(yǎng)化、生物多樣性降低等一系列環(huán)境問(wèn)題[5]. 水塘作為氮匯的重要場(chǎng)所,同時(shí)也多分布在流域上游,因此研究水塘反硝化作用對(duì)于流域源頭的水質(zhì)保護(hù)具有重要意義.

    丘陵區(qū)水塘是一種常見(jiàn)的連接農(nóng)業(yè)區(qū)與河流的重要水文通道[6],其能夠截留、過(guò)濾和凈化大量的污染物,是氮磷等污染物發(fā)生遷移轉(zhuǎn)化的重要場(chǎng)所,不僅保障了農(nóng)業(yè)用水,也促進(jìn)了養(yǎng)分的消納和循環(huán)利用[7]. 水塘中氮循環(huán)過(guò)程主要包括:氨化作用、礦化作用、植物吸收、硝化作用和反硝化作用等,其中氮素去除的機(jī)制主要有吸附沉降、植物吸收、反硝化作用,前兩者都難以從根本上將氮素去除,而后者使活性氮最終以惰性氮(N2)形式回歸大氣,將氮素永久去除[8]. 反硝化是氮循環(huán)的重要過(guò)程,也是水塘氮素去除的主要機(jī)制[9],對(duì)于減少氮素輸出、改善丘陵區(qū)水質(zhì)具有重要作用.

    反硝化作用受諸多因素的影響,如溫度、pH、硝態(tài)氮濃度、有機(jī)碳和溶解氧(DO)等[1]. 由于水塘水體受污染程度和污染類(lèi)型不同,導(dǎo)致反硝化作用的主要影響因素也存在差異. 反硝化作用主要發(fā)生在淹水環(huán)境中(如稻田、水塘等),是氮素有效去除的最佳途徑[10]. 然而,反硝化具有強(qiáng)烈的時(shí)空異質(zhì)性,精確、簡(jiǎn)便測(cè)定反硝化速率一直是一個(gè)世界性難題[11]. 這是由于反硝化主要產(chǎn)物N2在大氣中的背景值高達(dá)79%,在如此高的背景環(huán)境中一般難以直接測(cè)定[12]. 當(dāng)前對(duì)原位狀態(tài)下水塘反硝化過(guò)程的研究報(bào)道較少. 已經(jīng)建立的反硝化測(cè)定方法多以間接測(cè)定為主,如乙炔抑制法、15N同位素示蹤法、化學(xué)計(jì)量法等[9],但是這些方法普遍存在人為擾動(dòng)大、操作繁瑣、誤差大等不足,無(wú)法精確測(cè)定水土界面反硝化速率. 近年來(lái),基于N2∶Ar的膜進(jìn)樣質(zhì)譜法(MIMS)被廣泛用于直接測(cè)定淹水環(huán)境反硝化速率[13],其測(cè)定速度快(<3 min)、精度高(<0.03%)、進(jìn)樣少(<7 mL)、操作簡(jiǎn)便等,被廣泛用于不同水體環(huán)境反硝化作用研究[9]. 目前,MIMS不僅已大量應(yīng)用于河道[14]、河口以及稻田土壤反硝化速率直接測(cè)定[14-15],還可以與反硝化培養(yǎng)裝置串聯(lián),模擬室內(nèi)環(huán)境條件控制下反硝化速率的變化[16-17].

    天目湖位于太湖流域上游,由沙河和大溪兩大水庫(kù)組成,是溧陽(yáng)市重要的飲用水源地. 天目湖流域內(nèi)分布著1000多個(gè)大小不一的水塘,每1 km2約有7個(gè)[18]. 這些水塘主要由降雨沖刷和坡面產(chǎn)流匯集而成,一般多具有小流域結(jié)構(gòu). 由于其數(shù)量眾多,這些水塘構(gòu)成了丘陵坡地污染攔截和消納的第一道防線,其反硝化脫氮潛力大小將直接影響下游天目湖湖體輸入的氮量及水生態(tài)健康. 本文根據(jù)天目湖水塘集水區(qū)土地利用及水塘類(lèi)型的特點(diǎn),選擇天目湖流域內(nèi)4類(lèi)(茶園塘、村塘、養(yǎng)殖塘、林塘)共14個(gè)典型水塘,采用膜進(jìn)樣質(zhì)譜法(MIMS)測(cè)定水體中溶存的氮?dú)?N2)濃度,并結(jié)合水氣交換通量模型估算水體反硝化潛力,分析水塘反硝化潛力變化特征及其主要影響因素,揭示其時(shí)空變化規(guī)律,以期為太湖上游丘陵區(qū)氮污染源頭控制和水質(zhì)改善提供科學(xué)支撐.

    1 材料與方法

    1.1 研究區(qū)域概況

    天目湖流域位于江蘇溧陽(yáng)境內(nèi),為太湖上游低山丘陵區(qū),流域總面積246 km2,由沙河和大溪兩大水庫(kù)流域組成,面積分別為154和92 km2. 根據(jù)2017年天目湖流域土地利用遙感數(shù)據(jù)統(tǒng)計(jì)分析,天目湖流域土地類(lèi)型主要為:林地、耕地、茶園、水域等,其中林地是該流域最主要的土地類(lèi)型,面積約110 km2,占流域總面積44.6%;耕地面積僅次于林地,約為47 km2,占比18.9%;茶園面積29 km2,位居第三,占比11.9%. 天目湖流域?qū)賮啛釒Ъ撅L(fēng)氣候,四季分明,夏季和秋季溫度高,降水量大. 流域內(nèi)農(nóng)業(yè)種植包括茶樹(shù)、水稻、油菜、小麥、稻-麥和稻-油,氮肥施用量分別約為959、445、315、247、692和760 kg/(hm2·a)[4],農(nóng)業(yè)施肥也主要集中在夏季和秋季,受氣候條件和施肥的影響,反硝化作用在夏季和秋季較為劇烈.

    1.2 樣品采集與沉積物培養(yǎng)實(shí)驗(yàn)

    本文選擇天目湖流域內(nèi)14個(gè)水塘,根據(jù)其集水區(qū)土地利用類(lèi)型及用途不同劃分為4種類(lèi)型水塘,其中茶園塘3個(gè)、村塘4個(gè)、林塘4個(gè)、養(yǎng)殖塘3個(gè)(圖1). 在反硝化作用劇烈的夏季和秋季,即2020年6-11月期間,每月月末采集1次各水塘表層水樣,測(cè)定各水塘水體中溶存的N2濃度. 同時(shí),于2020年10月采集1次各水塘的柱狀沉積物,用于室內(nèi)反硝化的測(cè)定. 其中,由于林塘沉積物厚度較薄,無(wú)法采集完整柱狀沉積物,因此林塘的反硝化速率未進(jìn)行測(cè)定.

    圖1 天目湖流域采樣點(diǎn)位分布及不同類(lèi)型水塘Fig.1 Distribution of sampling points and diagram of different types of ponds in Lake Tianmu Basin

    采集樣品用于N2測(cè)定時(shí),應(yīng)全程動(dòng)作緩慢,避免氣泡產(chǎn)生,并減少與空氣接觸. 首先,將采樣器緩慢沒(méi)入水塘水體表面約10 cm后,緩慢提起采樣器. 打開(kāi)采樣器底部止水閥,將采樣器中的水樣通過(guò)硅膠管緩慢引流至帶橡膠墊片的螺口瓶(7 mL)的底部,當(dāng)水樣裝滿(mǎn)整個(gè)樣品瓶后應(yīng)繼續(xù)引流,讓其溢出樣品瓶體積3倍以上. 隨后,在樣品瓶液面下注射0.1 mL飽和ZnCl2溶液,以抑制微生物的活動(dòng). 最后,迅速擰緊瓶蓋,并將樣品瓶倒置觀察其中是否有氣泡產(chǎn)生,如果有氣泡則需要重新采集,每次樣品均取3個(gè)重復(fù). 同時(shí)采集水塘表層水樣約550 mL,用于測(cè)定水體中各項(xiàng)水質(zhì)指標(biāo).

    在樣品采集時(shí)的采樣現(xiàn)場(chǎng),利用便攜式多參數(shù)水質(zhì)檢測(cè)儀(YSI 6000,USA)對(duì)采樣點(diǎn)位的各項(xiàng)環(huán)境要素進(jìn)行現(xiàn)場(chǎng)測(cè)定,包括水溫、pH、溶解氧等指標(biāo).

    此外,還在2020年10月采集水塘無(wú)擾動(dòng)沉積物柱狀(10 cm)樣品,在各水塘不同區(qū)域采集3個(gè)平行沉積物柱狀樣品,用于室內(nèi)反硝化培養(yǎng)實(shí)驗(yàn),采樣結(jié)束后并采集各水塘原位水樣作為反硝化培養(yǎng)的上覆水. 室內(nèi)培養(yǎng)時(shí),首先將上覆水倒入培養(yǎng)裝置中,并將柱狀沉積物垂直放入培養(yǎng)裝置中,并控制上覆水水位高出培養(yǎng)柱約3~5 cm,調(diào)節(jié)溫度與采樣時(shí)的水溫一致,同時(shí)利用通氣泵保持上覆水通氣6 h以上,使得沉積物達(dá)到穩(wěn)定條件,并使得柱樣淹水環(huán)境溶解氧濃度與原位水體基本一致. 培養(yǎng)開(kāi)始前,首先關(guān)閉通氣閥,將柱蓋在水下擰緊,并連接柱蓋上的進(jìn)出水管;然后打開(kāi)進(jìn)出水管的閥門(mén),當(dāng)進(jìn)水管的原位水樣可在重力作用下進(jìn)入培養(yǎng)柱,同時(shí)出水管有水樣流出時(shí),說(shuō)明整個(gè)培養(yǎng)系統(tǒng)是密閉的. 同時(shí),要觀察進(jìn)出水管和密閉柱樣是否有氣泡殘留,若有氣泡則應(yīng)排空氣泡. 培養(yǎng)開(kāi)始后,開(kāi)啟培養(yǎng)裝置頂部的可調(diào)速電機(jī),轉(zhuǎn)速設(shè)置為75轉(zhuǎn)/min,以保證沉積物培養(yǎng)柱中水體的充分混勻,采樣第一個(gè)樣品作為0時(shí)刻樣品,同時(shí)依次在2、4、6、8、10、12、22、24 h時(shí)刻分別取樣,每個(gè)時(shí)空均采集3個(gè)5 mL的平行樣品,以減少誤差,采集的樣品均于48h內(nèi)測(cè)定完畢. 具體的室內(nèi)培養(yǎng)實(shí)驗(yàn)方法詳見(jiàn)李曉波等[9].

    1.3 理化指標(biāo)的室內(nèi)測(cè)定

    對(duì)于N2濃度,利用膜進(jìn)樣質(zhì)譜儀(MIMS)測(cè)定水樣(水塘表層水和室內(nèi)沉積物培養(yǎng)上覆水)中N2/Ar的比值來(lái)計(jì)算其氮?dú)鉂舛?

    1.4 反硝化速率計(jì)算

    水塘反硝化作用的估算包括2種方法,一是對(duì)水體的直接測(cè)定,采用MIMS測(cè)定水體N2濃度,計(jì)算N2的過(guò)飽和濃度,并結(jié)合水氣交換通量模型計(jì)算水體反硝化速率;二是沉積物室內(nèi)培養(yǎng)測(cè)定的反硝化速率,根據(jù)不同培養(yǎng)時(shí)間反硝化產(chǎn)物N2濃度和時(shí)間變化進(jìn)行線性擬合,計(jì)算反硝化速率. 第一種方法操作簡(jiǎn)便,對(duì)原位條件無(wú)擾動(dòng),能夠快速獲取野外原位條件下反硝化速率. 但同時(shí)存在較大的不確定性,這主要是由于當(dāng)前有多個(gè)備選的水氣交換通量估算模型,當(dāng)這些模型應(yīng)用至丘陵區(qū)溝塘反硝化速率測(cè)算中可能存在誤差不一的問(wèn)題;第二種方法是目前最為常用的方法,但一般在實(shí)驗(yàn)室控制條件下進(jìn)行,操作繁瑣且不可避免地改變了原位條件(比如光照、風(fēng)速擾動(dòng)等),對(duì)估算精度會(huì)產(chǎn)生一定的影響. 盡管如此,眾多學(xué)者認(rèn)為該方法估算的反硝化速率相對(duì)可靠,與真實(shí)的反硝化速率較為接近[8,13-14,17]. 本研究根據(jù)沉積物測(cè)算的反硝化速率來(lái)約束模型估算的結(jié)果,并優(yōu)選適合于丘陵溝塘的水氣交換通量模型,來(lái)進(jìn)一步提升反硝化測(cè)算的可靠性和精度.

    上述2種方法N2濃度的測(cè)定均采用MIMS法. MIMS法是利用氣體過(guò)濾膜系統(tǒng)直接將水樣中溶存的氣體吸入質(zhì)譜儀真空分析室,經(jīng)過(guò)離子化,由質(zhì)譜儀測(cè)出氣體組成和含量,對(duì)N2、Ar的測(cè)定精度分別小于0.5%、0.05%,但對(duì)于N2/Ar的測(cè)定精度小于0.03%,水環(huán)境中Ar溶解度一般只受溫度和鹽度的影響,其濃度非常穩(wěn)定.

    通過(guò)水氣交換通量模型來(lái)計(jì)算水體反硝化速率時(shí),首先需計(jì)算出過(guò)飽和濃度(Δ[N2],μmol/L),具體可通過(guò)下述公式得出:

    Δ[N2]=[N2]/[Ar]×[Ar]*-[N2]*

    (1)

    式中,[N2]/[Ar]為經(jīng)過(guò)質(zhì)譜校準(zhǔn)過(guò)的水樣中N2和Ar的濃度比值,[N2]*、[Ar]*為特定溫鹽條件下N2、Ar的理論平衡濃度(μmol/L).

    依據(jù)水-氣界面的分子擴(kuò)散模型和亨利定律,估算不同單位面積水塘單位時(shí)間反硝化N2的增量(F,mmol/(m2·d)),計(jì)算公式如下:

    F=k·Δ[N2]

    (2)

    式中,k為氣體擴(kuò)散系數(shù),由現(xiàn)場(chǎng)測(cè)定的溫度、風(fēng)速等參數(shù)結(jié)合氣體Sc數(shù)(schmidt number)計(jì)算. 考慮到風(fēng)速變化對(duì)估算k值的影響,本研究采用Cole和Caraco提出的CC98模型[19],該模型開(kāi)發(fā)時(shí)的風(fēng)力條件與本研究的水塘最為相似,估算k值的計(jì)算公式如下:

    CC98:k=[2.07+(0.215U101.7)](Sc/600)-2/3

    (3)

    式中,U10為水面上方10 m處的風(fēng)速(m/s),由流域內(nèi)安裝的氣象站獲得. 而Sc數(shù)為水的動(dòng)力黏度與待測(cè)氣體分子擴(kuò)散速率之比,對(duì)于特定氣體,Sc數(shù)與水溫、鹽度等物理參數(shù)有關(guān),Wanninkhof[20]提出N2氣體Sc數(shù)與水溫的關(guān)系式如下:

    Sc=1970.7-131.45t+4.1390t2-0.052106t3

    (4)

    式中,t為水溫(℃).

    除了CC98模型中的k值計(jì)算方法外,還有許多其他的方法來(lái)計(jì)算k值. 本文列舉了另外2種計(jì)算k值的經(jīng)驗(yàn)?zāi)P?,將不同模型?jì)算的結(jié)果進(jìn)行比較,以探究不同k值計(jì)算方法的選擇所帶來(lái)的不確定性. 模型方程如下:

    1)LM86模型: Liss and Merlivat (1986)[21]:

    k=0.17U10(Sc/600)-2/3(0k=(2.85U10-9.65)(Sc/600)-1/2(3.6k=(5.9U10-49.3)(Sc/600)-1/2(U10>13)

    (5)

    2)CW03模型: Crusius and Wanninkhof (2003)[22]:

    (6)

    在室內(nèi)沉積物培養(yǎng)實(shí)驗(yàn)中,根據(jù)反硝化產(chǎn)生的N2濃度變化與培養(yǎng)時(shí)間的關(guān)系來(lái)估算反硝化速率. 具體計(jì)算公式如下:

    Y=mX+b

    (7)

    式中,Y是不同培養(yǎng)時(shí)間對(duì)應(yīng)的反硝化產(chǎn)物N2濃度(μmol/L);斜率m是沉積物上覆水中N2濃度變化速率(μmol/(L·h));X是培養(yǎng)時(shí)間(h);b是初始上覆水中N2濃度(μmol/L);結(jié)合上覆水體積V(L)和沉積物柱樣橫截面積S(m2),即得反硝化速率(F,mmol/(m2·d)),計(jì)算公式如下:

    F=0.024mV/S

    (8)

    式中,0.024為單位換算系數(shù).

    2 結(jié)果與分析

    2.1 上覆水與沉積物理化特性

    表1 夏季和秋季各個(gè)水塘水體水質(zhì)參數(shù)(平均值±標(biāo)準(zhǔn)偏差)*

    表2 各水塘沉積物理化參數(shù)(平均值±標(biāo)準(zhǔn)偏差)*

    2.2 水塘水體N2過(guò)飽和濃度變化

    用膜進(jìn)樣質(zhì)譜法測(cè)定水塘水體反硝化作用,結(jié)果表明,不同類(lèi)型水塘反硝化作用存在一定的時(shí)空差異. 圖2a為不同類(lèi)型水塘夏季和秋季N2過(guò)飽和濃度(ΔN2)的變化. 14個(gè)水塘上覆水溶解N2過(guò)飽和濃度介于1.36~28.35 μmol/L之間,平均值為(8.23±6.04) μmol/L. 夏季4種類(lèi)型水塘N2過(guò)飽和濃度平均為(8.81±4.08) μmol/L,其中茶園塘N2過(guò)飽和濃度在所有類(lèi)型水塘中最高,達(dá)到(13.91±4.79) μmol/L;其次是村塘,為(8.28±1.36) μmol/L;林塘和養(yǎng)殖塘較為接近,分別為(6.54±1.34)、(6.51±1.79) μmol/L. 夏季水塘上覆水溶解N2過(guò)飽和濃度大小依次為:茶園塘>村塘>林塘>養(yǎng)殖塘.

    秋季4種類(lèi)型水塘N2過(guò)飽和濃度均值為(7.64±7.46) μmol/L,其中茶園塘ΔN2最高,為(19.39±5.48)μmol/L;養(yǎng)殖塘次之,為(4.61±1.57) μmol/L;村塘為(4.57±1.29) μmol/L;林塘ΔN2最低,為(2.00±0.42)μmol/L,秋季水塘上覆水ΔN2大小依次為:茶園塘>養(yǎng)殖塘>村塘>林塘. 總體上,水塘水體夏季ΔN2高于秋季,村塘、林塘和養(yǎng)殖塘均呈現(xiàn)相同的規(guī)律,而茶園塘夏季ΔN2低于秋季,可能與秋季茶園高強(qiáng)度施肥有關(guān).

    2.3 水塘水體反硝化潛力

    采用CC98模型中的水氣交換通量模型來(lái)估算不同類(lèi)型水塘水體反硝化潛力,由圖2b可知水塘反硝化潛力在時(shí)間和空間上均存在差異. 14個(gè)水塘的反硝化速率均值為(4.75±3.27) mmol/(m2·d). 其中茶園塘反硝化速率最高,為(9.36±1.44) mmol/(m2·d),村塘和養(yǎng)殖塘分別為(5.47±1.52)和(3.09±0.52)mmol/(m2·d),林塘最低為(1.09±0.48)mmol/(m2·d),茶園塘反硝化速率大約是村塘的1.7倍、林塘的8.6倍、養(yǎng)殖塘的3倍,其反硝化速率均值大小順序?yàn)椋翰鑸@塘>村塘>養(yǎng)殖塘>林塘. 反硝化速率與氮?dú)膺^(guò)飽和濃度的季節(jié)變化規(guī)律一致.

    室內(nèi)沉積物培養(yǎng)實(shí)驗(yàn)結(jié)果表明,不同類(lèi)型水塘沉積物的反硝化速率表現(xiàn)出一定的差異性(圖2b). 9個(gè)水塘的沉積物反硝化速率在2.21~10.34 mmol/(m2·d)之間,平均值為(5.73±2.65)mmol/(m2·d). 其中村塘和養(yǎng)殖塘反硝化速率分別為5.24和2.75 mmol/(m2·d),茶園塘的反硝化脫氮能力最強(qiáng),反硝化速率為9.19 mmol/(m2·d). 總體上,室內(nèi)培養(yǎng)試驗(yàn)測(cè)得的反硝化速率與CC98模型估算結(jié)果較為接近.

    圖2 夏季和秋季不同類(lèi)型水塘N2過(guò)飽和濃度(a)和水體反硝化速率(b)(b圖中“-”表示水塘沉積物反硝化速率)Fig.2 The excess concentrations of N2(a) and denitrification rate in water(b)in different types of ponds in summer and autumn(“-” in the figure b represents the measured denitrification rates from pond sediments)

    2.4 水塘反硝化作用與環(huán)境因子的關(guān)系

    圖3 水塘水體N2過(guò)飽和濃度與水體(a)和沉積物(b)理化指標(biāo)的相關(guān)性分析Fig.3 Correlation analysis of N2 supersaturated concentration in pond water with physicochemical indexes of water body (a) and sediment (b)

    沉積物反硝化速率與硝態(tài)氮和總氮濃度呈現(xiàn)顯著的正相關(guān)性(P<0.01,圖3b),而與容重存在顯著的負(fù)相關(guān)性(P<0.05),pH、含水率和氨氮與沉積物反硝化速率的相關(guān)性不顯著. 水體和沉積物反硝化作用均與硝態(tài)氮濃度呈顯著正相關(guān)關(guān)系.

    3 討論

    3.1 水塘沉積物反硝化估算及其不確定性

    通過(guò)膜進(jìn)樣質(zhì)譜法測(cè)定上覆水N2過(guò)飽和濃度,結(jié)合水氣交換通量模型估算上覆水反硝化速率在0.53~12.65 mmol/(m2·d)之間,平均值為(4.75±3.27)mmol/(m2·d). 與其他研究相比,本研究報(bào)道的水塘反硝化速率范圍明顯高于其他不同類(lèi)型水體的范圍. 比如,張波等對(duì)太湖流域典型河流沉積物反硝化速率研究發(fā)現(xiàn),利用15N同位素配對(duì)法測(cè)得25℃時(shí)反硝化速率為0.79 mmol/(m2·d),10℃時(shí)為0.26 mmol/(m2·d),明顯低于水塘沉積物反硝化速率[23]. Bonnett等利用乙炔抑制法測(cè)定英國(guó)西南部Tamar河流沉積物反硝化速率為0.64 mmol/(m2·d),低于本研究水塘反硝化速率[24]. Lazar等利用15N同位素配對(duì)技術(shù)測(cè)得美國(guó)東北部農(nóng)村海貍池塘沉積物反硝化速率為1.73~4.21 mmol/(m2·d),也總體低于本研究的結(jié)果[25]. 然而,一些河流的反硝化速率要高于本研究的結(jié)果. 比如,陳能汪等利用膜進(jìn)樣質(zhì)譜儀對(duì)九龍江河口區(qū)溶解N2濃度測(cè)定,并采用水氣交換模型估算了其反硝化速率,發(fā)現(xiàn)其均值為21.7 mmol/(m2·d),明顯高于本研究水塘反硝化速率[15]. 龍虹竹等利用改進(jìn)的乙炔抑制法對(duì)川中丘陵區(qū)農(nóng)業(yè)源頭溝渠沉積物反硝化速率測(cè)定,其反硝化速率均值為8.96 mmol/(m2·d),高于本研究水塘沉積物反硝化速率[26]. 李佳霖等同樣利用15N同位素配對(duì)法測(cè)得彌河沉積物反硝化速率為16.34 mmol/(m2·d),明顯高于天目湖流域水塘沉積物反硝化速率[27]. 通過(guò)這些對(duì)比發(fā)現(xiàn),本研究報(bào)道的水塘反硝化速率波動(dòng)范圍在全球不同類(lèi)型的水體中均有報(bào)道. 由于丘陵區(qū)水塘與坡地相連并直接受納坡地產(chǎn)生的污染物,本研究的結(jié)果證實(shí)了其相對(duì)較強(qiáng)的反硝化速率在削減氮負(fù)荷及減少下游污染方面發(fā)揮了重要的作用.

    本研究通過(guò)對(duì)水樣中氮?dú)獾闹苯訙y(cè)定,并通過(guò)水氣交換模型測(cè)算了反硝化速率. 與其他傳統(tǒng)的測(cè)定方法相比,膜進(jìn)樣質(zhì)譜法直接測(cè)定水體反硝化作用具有精度高、操作簡(jiǎn)便、結(jié)果可靠等特點(diǎn)[28-30],并且適用于無(wú)沉積物或少沉積物的水體. 本研究對(duì)不同類(lèi)型水塘水體反硝化速率直接測(cè)定,同時(shí)采集各水塘沉積物柱樣進(jìn)行室內(nèi)培養(yǎng)測(cè)定反硝化速率,兩種測(cè)定方法的結(jié)果具有可比性(圖4),并且前者實(shí)驗(yàn)周期更短.

    然而,本研究的估算存在一定的不確定性. 首先,反硝化估算受采樣頻率的影響. 例如每月采集一次水塘上覆水存在很大的不確定性,因?yàn)镹2濃度時(shí)刻在變化,低頻的觀測(cè)難以準(zhǔn)確獲得N2的動(dòng)態(tài)變化. 此外,本研究?jī)H估算了N2水氣通量,未考慮不完全反硝化產(chǎn)生的N2O排放. 盡管對(duì)于大多數(shù)水體而言,N2O水氣通量一般只占總通量的0.1%~1.0%,比如九龍江河口區(qū)反硝化作用產(chǎn)生的N2O水氣交換通量只占0.25%[15]. 但在某些情況下,最高可達(dá)到6%[31],這說(shuō)明我們?cè)谝欢ǔ潭壬系凸懒怂练聪趸瘽摿?,需要?duì)N2和N2O水氣通量同時(shí)評(píng)估. 此外,反硝化估算是否準(zhǔn)確與水氣交換模型的選擇有很大的關(guān)系,這是因?yàn)楫?dāng)前有很多不同的估算模型(見(jiàn)本文2.4節(jié)部分). 其中,Cole和Caraco[19]提出的水氣交換模型被應(yīng)用最廣泛、精度最高,且適用于丘陵區(qū)水文氣象條件. 該模型采用風(fēng)速和溫度作為模型的輸入,使其更接近原位環(huán)境條件,能更精準(zhǔn)地刻畫(huà)水氣交換過(guò)程. 本研究還進(jìn)一步采用沉積物柱狀培養(yǎng)測(cè)定反硝化速率來(lái)驗(yàn)證水氣交換通量的結(jié)果,研究進(jìn)一步證實(shí)CC98模型估算的反硝化速率與傳統(tǒng)室內(nèi)反硝化模擬試驗(yàn)測(cè)得的速率非常接近. 然而,對(duì)于其他2個(gè)常用的水氣交換模型,估算的結(jié)果遠(yuǎn)小于室內(nèi)沉積物培養(yǎng)測(cè)得的反硝化速率. 比如,當(dāng)采用LM86和CW03模型時(shí),估算的反硝化速率僅為實(shí)驗(yàn)室測(cè)得的速率的1/10~1/3,這說(shuō)明這些模型不適用于丘陵區(qū)水塘反硝化潛力的估算. 然而,水氣交換過(guò)程不僅受到大氣擾動(dòng),也受到水文和微氣象等因素的影響,在未來(lái)的工作中,應(yīng)原位測(cè)定k值,以期更準(zhǔn)確地估算反硝化速率.

    圖4 不同k值計(jì)算模型估算的水塘反硝化速率差異(LM86、CW03和CC98指采用相應(yīng)的公式計(jì)算得出的反硝化速率,具體公式見(jiàn)本文1.4節(jié)部分;圖中“-”代表水塘沉積物反硝化速率)Fig.4 The difference of denitrification rates of ponds estimated by different k calculation models

    3.2 影響水塘反硝化作用的主要因素

    溶解氧水平是影響反硝化脫氮的重要因素. 這是因?yàn)榉聪趸话銥閰捬踹^(guò)程,高濃度溶解氧的存在會(huì)顯著抑制反硝化功能微生物的活性,進(jìn)而抑制反硝化脫氮過(guò)程[37]. 本研究發(fā)現(xiàn),水塘反硝化速率與溶解氧濃度之間呈顯著負(fù)相關(guān)關(guān)系(r=-0.47,P<0.05). 在不同O2濃度水平下檢測(cè)土壤反硝化速率,其最高值出現(xiàn)在無(wú)氧條件下[38];丘陵區(qū)水塘水體溶解氧濃度夏季低于秋季,其反硝化速率夏季大于秋季;在不同類(lèi)型水塘上,茶園塘和村塘水體溶解氧濃度低于林塘和養(yǎng)殖塘,低溶解氧濃度有利于反硝化細(xì)菌的生長(zhǎng),從而促進(jìn)反硝化作用[39];在沉積物中較高的含水率會(huì)阻止氧氣的傳輸,從而形成厭氧條件,茶園塘沉積物含水率較高,其反硝化速率也較高,二者有一定的相關(guān)關(guān)系,但相關(guān)性不顯著.

    4 結(jié)論

    1)研究區(qū)夏、秋兩季4種類(lèi)型水塘水體N2過(guò)飽和濃度介于1.36~28.35 μmol/L之間,反硝化速率在0.53~12.65 mmol/(m2·d)之間,總體上夏季反硝化速率高于秋季.

    2)不同類(lèi)型水塘反硝化速率存在顯著空間差異,茶園塘反硝化速率最高,村塘和養(yǎng)殖塘次之,林塘最低.

    4)與水塘沉積物室內(nèi)培養(yǎng)測(cè)定反硝化速率相比,直接測(cè)定水體原位反硝化速率具有可行性,且操作簡(jiǎn)便,實(shí)驗(yàn)周期短.

    5)對(duì)不同水氣交換模型的計(jì)算結(jié)果進(jìn)行比較,發(fā)現(xiàn)CC98模型更適用于太湖上游丘陵區(qū)水塘反硝化速率的估算.

    猜你喜歡
    過(guò)飽和水氣水塘
    一個(gè)小水塘
    紊動(dòng)促進(jìn)過(guò)飽和總?cè)芙鈿怏w釋放研究
    遼中區(qū)患病草魚(yú)體內(nèi)嗜水氣單胞菌分離、鑒定與致病力測(cè)定
    海上邊水氣藏利用試井資料確定水侵狀況研究
    海洋石油(2021年3期)2021-11-05 07:42:26
    硼碳過(guò)飽和馬氏體激光熔覆層耐磨性研究
    醉在水塘
    基于MFD的城市區(qū)域過(guò)飽和交通信號(hào)優(yōu)化控制
    荒漠水塘
    風(fēng)速對(duì)過(guò)飽和總?cè)芙鈿怏w釋放速率的影響
    醫(yī)院感染嗜水氣單胞菌的臨床治療分析
    久久久久久九九精品二区国产 | 午夜日韩欧美国产| 俺也久久电影网| 亚洲欧美日韩无卡精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲色图 男人天堂 中文字幕| 国产精品九九99| 又粗又爽又猛毛片免费看| 国产成人精品久久二区二区免费| 桃色一区二区三区在线观看| 成人午夜高清在线视频| 国产精品久久久av美女十八| 国语自产精品视频在线第100页| www.精华液| 久久天躁狠狠躁夜夜2o2o| 男女床上黄色一级片免费看| 亚洲欧美日韩无卡精品| 国产探花在线观看一区二区| 午夜精品在线福利| 亚洲精品色激情综合| 一级黄色大片毛片| 亚洲国产中文字幕在线视频| 露出奶头的视频| 国产免费av片在线观看野外av| 麻豆久久精品国产亚洲av| 国产免费男女视频| 欧洲精品卡2卡3卡4卡5卡区| 久99久视频精品免费| 亚洲狠狠婷婷综合久久图片| 久久精品aⅴ一区二区三区四区| 久久久久久免费高清国产稀缺| 久久精品亚洲精品国产色婷小说| 国产99久久九九免费精品| 成人三级做爰电影| 国产成人一区二区三区免费视频网站| 天天添夜夜摸| 黄色女人牲交| 又粗又爽又猛毛片免费看| 激情在线观看视频在线高清| 亚洲国产精品合色在线| 桃色一区二区三区在线观看| 91av网站免费观看| 草草在线视频免费看| 国产精品久久久久久久电影 | 国产精品1区2区在线观看.| 好看av亚洲va欧美ⅴa在| 看片在线看免费视频| 亚洲成a人片在线一区二区| 免费人成视频x8x8入口观看| 国产一区二区在线av高清观看| 国产又色又爽无遮挡免费看| 又粗又爽又猛毛片免费看| 国产精品国产高清国产av| 成人精品一区二区免费| 久久精品国产99精品国产亚洲性色| 高清毛片免费观看视频网站| 色播亚洲综合网| 男人舔女人下体高潮全视频| 国产精品影院久久| 人妻夜夜爽99麻豆av| 我的老师免费观看完整版| 欧美 亚洲 国产 日韩一| 级片在线观看| 久久久久久亚洲精品国产蜜桃av| 妹子高潮喷水视频| 男人舔女人的私密视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av美国av| 听说在线观看完整版免费高清| 18禁国产床啪视频网站| 日韩精品青青久久久久久| 国产乱人伦免费视频| 国产精品亚洲美女久久久| 人人妻人人看人人澡| 非洲黑人性xxxx精品又粗又长| 一边摸一边抽搐一进一小说| 亚洲九九香蕉| 麻豆国产97在线/欧美 | 欧美成人性av电影在线观看| 我要搜黄色片| 亚洲人成77777在线视频| 精品久久久久久久毛片微露脸| 无遮挡黄片免费观看| 三级男女做爰猛烈吃奶摸视频| 色老头精品视频在线观看| 韩国av一区二区三区四区| 看片在线看免费视频| 高潮久久久久久久久久久不卡| 欧美中文综合在线视频| 精品第一国产精品| 蜜桃久久精品国产亚洲av| e午夜精品久久久久久久| 国产av又大| 久久精品91无色码中文字幕| 日本免费a在线| 国产一区二区在线av高清观看| 亚洲精品久久成人aⅴ小说| 亚洲精品一区av在线观看| 巨乳人妻的诱惑在线观看| 激情在线观看视频在线高清| 脱女人内裤的视频| 午夜福利成人在线免费观看| 国产视频内射| 亚洲成av人片免费观看| 两个人视频免费观看高清| 亚洲无线在线观看| 国产亚洲av嫩草精品影院| 午夜精品一区二区三区免费看| 色综合欧美亚洲国产小说| 在线看三级毛片| 1024手机看黄色片| 少妇的丰满在线观看| 桃红色精品国产亚洲av| 最近在线观看免费完整版| 又大又爽又粗| 欧美色视频一区免费| 非洲黑人性xxxx精品又粗又长| 亚洲欧美激情综合另类| 午夜精品在线福利| 禁无遮挡网站| 在线观看日韩欧美| 很黄的视频免费| 99久久国产精品久久久| 女人被狂操c到高潮| 欧美丝袜亚洲另类 | 国产成人av激情在线播放| 国产单亲对白刺激| 首页视频小说图片口味搜索| 精品一区二区三区四区五区乱码| 两性夫妻黄色片| 免费看十八禁软件| 久久久久九九精品影院| 青草久久国产| 色精品久久人妻99蜜桃| 成人av在线播放网站| bbb黄色大片| 桃红色精品国产亚洲av| 欧美一级毛片孕妇| 国产真实乱freesex| 给我免费播放毛片高清在线观看| 高清毛片免费观看视频网站| 黄片大片在线免费观看| www.www免费av| 两性夫妻黄色片| 久久午夜亚洲精品久久| a级毛片在线看网站| 国内毛片毛片毛片毛片毛片| 18禁美女被吸乳视频| 国内毛片毛片毛片毛片毛片| 99久久无色码亚洲精品果冻| 一个人免费在线观看电影 | 久久久久久九九精品二区国产 | a级毛片在线看网站| 制服丝袜大香蕉在线| 床上黄色一级片| 女生性感内裤真人,穿戴方法视频| 亚洲精品美女久久久久99蜜臀| 老汉色av国产亚洲站长工具| 一二三四社区在线视频社区8| 91av网站免费观看| 久久人妻福利社区极品人妻图片| 久久中文字幕一级| 日本a在线网址| 精品欧美国产一区二区三| 两人在一起打扑克的视频| 成年免费大片在线观看| 99久久精品国产亚洲精品| 很黄的视频免费| 国产区一区二久久| 女人被狂操c到高潮| av片东京热男人的天堂| 狂野欧美白嫩少妇大欣赏| 少妇的丰满在线观看| 亚洲aⅴ乱码一区二区在线播放 | 99riav亚洲国产免费| 毛片女人毛片| 亚洲成人久久爱视频| 国产亚洲欧美98| 国产成人欧美在线观看| 亚洲欧美精品综合久久99| 蜜桃久久精品国产亚洲av| 国产片内射在线| 50天的宝宝边吃奶边哭怎么回事| 日日夜夜操网爽| 中文字幕av在线有码专区| 亚洲一区二区三区不卡视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲av成人一区二区三| 婷婷亚洲欧美| 国产精品免费视频内射| 国产亚洲精品久久久久久毛片| 国产视频一区二区在线看| 日本成人三级电影网站| 成人国产一区最新在线观看| 午夜a级毛片| 久久精品国产清高在天天线| or卡值多少钱| 欧美成人性av电影在线观看| 成人手机av| 亚洲五月天丁香| 好男人在线观看高清免费视频| 法律面前人人平等表现在哪些方面| 国产在线精品亚洲第一网站| 成人永久免费在线观看视频| 欧美成人午夜精品| 久久久久久免费高清国产稀缺| 欧美在线黄色| 长腿黑丝高跟| 视频区欧美日本亚洲| 免费人成视频x8x8入口观看| 久久九九热精品免费| 桃红色精品国产亚洲av| 亚洲人成电影免费在线| 狂野欧美激情性xxxx| 脱女人内裤的视频| 亚洲va日本ⅴa欧美va伊人久久| 国产高清视频在线播放一区| 欧美日本亚洲视频在线播放| 欧美又色又爽又黄视频| 亚洲美女黄片视频| 国产免费av片在线观看野外av| 中文字幕人妻丝袜一区二区| 九九热线精品视视频播放| 国产亚洲精品久久久久5区| 中国美女看黄片| 成人av在线播放网站| 校园春色视频在线观看| 午夜福利成人在线免费观看| 成熟少妇高潮喷水视频| 精品午夜福利视频在线观看一区| 亚洲一区高清亚洲精品| 又爽又黄无遮挡网站| 成人三级做爰电影| 国产人伦9x9x在线观看| 最近最新中文字幕大全电影3| 久久国产精品人妻蜜桃| 午夜两性在线视频| 亚洲欧美日韩东京热| 91老司机精品| www.精华液| 午夜久久久久精精品| 国产在线精品亚洲第一网站| 久久久久久久久中文| 国产精品亚洲美女久久久| 亚洲成人免费电影在线观看| 欧美黑人精品巨大| 最近在线观看免费完整版| 好男人电影高清在线观看| 国产69精品久久久久777片 | 中文字幕av在线有码专区| 国产高清视频在线观看网站| 亚洲成人国产一区在线观看| 国产精品久久久人人做人人爽| 好男人电影高清在线观看| 男人舔女人的私密视频| 精品福利观看| 黄色毛片三级朝国网站| bbb黄色大片| 国产成人一区二区三区免费视频网站| 桃色一区二区三区在线观看| 国产区一区二久久| aaaaa片日本免费| 欧美一级a爱片免费观看看 | 久久人妻av系列| 亚洲成人久久性| 制服丝袜大香蕉在线| 特大巨黑吊av在线直播| 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| svipshipincom国产片| 日本成人三级电影网站| 国产一区二区在线av高清观看| 久久草成人影院| 欧美日韩黄片免| 在线观看美女被高潮喷水网站 | 国产男靠女视频免费网站| 1024视频免费在线观看| 国产免费男女视频| 黄色视频不卡| 国产精品永久免费网站| 国产三级在线视频| tocl精华| 欧美日韩瑟瑟在线播放| 久久中文字幕人妻熟女| 国产精品电影一区二区三区| 国产私拍福利视频在线观看| 亚洲国产中文字幕在线视频| 国产野战对白在线观看| 免费在线观看黄色视频的| 亚洲色图av天堂| 欧美av亚洲av综合av国产av| 18禁美女被吸乳视频| 久久久精品大字幕| 国产精品影院久久| 国产成人精品无人区| 日日夜夜操网爽| 亚洲片人在线观看| 可以免费在线观看a视频的电影网站| 日本黄大片高清| 欧美zozozo另类| 看片在线看免费视频| 男女做爰动态图高潮gif福利片| 桃色一区二区三区在线观看| www.自偷自拍.com| 成人av在线播放网站| 婷婷精品国产亚洲av| 俄罗斯特黄特色一大片| 亚洲 欧美一区二区三区| 午夜福利欧美成人| 老司机午夜十八禁免费视频| 日本黄色视频三级网站网址| 69av精品久久久久久| 在线国产一区二区在线| 国产精品九九99| 国产视频内射| 搞女人的毛片| 亚洲乱码一区二区免费版| 窝窝影院91人妻| 成在线人永久免费视频| avwww免费| 中文在线观看免费www的网站 | 久久婷婷成人综合色麻豆| 丁香欧美五月| 床上黄色一级片| 亚洲人成77777在线视频| www日本在线高清视频| 亚洲人成网站高清观看| 久久香蕉精品热| 久久亚洲精品不卡| 国产欧美日韩一区二区三| 久久精品aⅴ一区二区三区四区| 人成视频在线观看免费观看| 最近视频中文字幕2019在线8| 麻豆国产av国片精品| 午夜免费激情av| 免费高清视频大片| 在线a可以看的网站| 国产av在哪里看| 亚洲精品色激情综合| 亚洲av第一区精品v没综合| 国产成人av教育| 国产激情偷乱视频一区二区| 亚洲中文字幕日韩| 亚洲国产看品久久| 变态另类丝袜制服| 免费在线观看黄色视频的| 正在播放国产对白刺激| 国产亚洲精品久久久久久毛片| 中文字幕人成人乱码亚洲影| 国产精品九九99| 国产成人一区二区三区免费视频网站| 国产熟女xx| 日韩精品免费视频一区二区三区| 日韩av在线大香蕉| 狂野欧美白嫩少妇大欣赏| 欧美乱色亚洲激情| 五月伊人婷婷丁香| 精品久久久久久久末码| 国产av一区在线观看免费| 亚洲国产精品合色在线| 国产高清视频在线播放一区| 不卡av一区二区三区| 伦理电影免费视频| 18禁黄网站禁片午夜丰满| 亚洲一区二区三区不卡视频| 亚洲一码二码三码区别大吗| 国产成人精品无人区| 亚洲人成伊人成综合网2020| 久久精品综合一区二区三区| 人人妻人人澡欧美一区二区| 免费观看人在逋| 亚洲欧美日韩高清专用| 成年女人毛片免费观看观看9| 国产成人aa在线观看| 老鸭窝网址在线观看| 日日夜夜操网爽| 婷婷精品国产亚洲av| 国产亚洲欧美在线一区二区| 亚洲人成电影免费在线| 人妻丰满熟妇av一区二区三区| 国产视频一区二区在线看| 一二三四社区在线视频社区8| 99精品久久久久人妻精品| 91麻豆av在线| 成熟少妇高潮喷水视频| 国产精华一区二区三区| 日日夜夜操网爽| 欧美日本亚洲视频在线播放| 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 成人国产综合亚洲| 亚洲国产看品久久| 麻豆av在线久日| 国产精品 欧美亚洲| 亚洲色图av天堂| 老汉色av国产亚洲站长工具| 男插女下体视频免费在线播放| 国内精品久久久久精免费| 亚洲精品在线观看二区| 亚洲一码二码三码区别大吗| www.精华液| 久久精品91蜜桃| 99国产综合亚洲精品| 色av中文字幕| 在线观看www视频免费| 欧美一区二区国产精品久久精品 | 首页视频小说图片口味搜索| 午夜日韩欧美国产| 国产亚洲精品久久久久久毛片| 欧美日韩国产亚洲二区| 亚洲精品粉嫩美女一区| 大型av网站在线播放| 好男人在线观看高清免费视频| 国产精品久久电影中文字幕| 97碰自拍视频| 性色av乱码一区二区三区2| 精品国内亚洲2022精品成人| av国产免费在线观看| 嫁个100分男人电影在线观看| 国产精品亚洲av一区麻豆| 国语自产精品视频在线第100页| 91大片在线观看| 国产精品国产高清国产av| 欧美成狂野欧美在线观看| 欧美zozozo另类| 国产精品久久电影中文字幕| 久久久久久九九精品二区国产 | 精品久久蜜臀av无| 久久精品亚洲精品国产色婷小说| 天天躁狠狠躁夜夜躁狠狠躁| 18美女黄网站色大片免费观看| 精品国产乱子伦一区二区三区| 日韩成人在线观看一区二区三区| 久久人妻福利社区极品人妻图片| 99精品在免费线老司机午夜| 精华霜和精华液先用哪个| 免费看十八禁软件| 成人一区二区视频在线观看| 国产男靠女视频免费网站| 亚洲精华国产精华精| 成人精品一区二区免费| 精品国产美女av久久久久小说| 国产区一区二久久| 中文亚洲av片在线观看爽| 国产麻豆成人av免费视频| 99热6这里只有精品| 亚洲一区二区三区不卡视频| 国产精品一及| 亚洲第一欧美日韩一区二区三区| 国产高清视频在线播放一区| 一边摸一边抽搐一进一小说| 欧美乱色亚洲激情| 久久久久久久午夜电影| 最近最新中文字幕大全电影3| 天天躁夜夜躁狠狠躁躁| 俺也久久电影网| 激情在线观看视频在线高清| 一个人免费在线观看的高清视频| 黑人欧美特级aaaaaa片| av国产免费在线观看| 久久久久久亚洲精品国产蜜桃av| 啦啦啦观看免费观看视频高清| 亚洲片人在线观看| 别揉我奶头~嗯~啊~动态视频| 久久久久国产精品人妻aⅴ院| 成在线人永久免费视频| 亚洲男人的天堂狠狠| 正在播放国产对白刺激| 欧美av亚洲av综合av国产av| 午夜激情福利司机影院| 国产午夜福利久久久久久| 欧美成人性av电影在线观看| 亚洲一区二区三区不卡视频| 一个人免费在线观看的高清视频| 色综合婷婷激情| 久久香蕉国产精品| 日本在线视频免费播放| 国产三级中文精品| 麻豆国产av国片精品| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av熟女| 超碰成人久久| 一个人免费在线观看的高清视频| 日日夜夜操网爽| 国产精品一区二区三区四区久久| 久久午夜综合久久蜜桃| 一区福利在线观看| 久久久久国产一级毛片高清牌| 国产精品久久久久久亚洲av鲁大| 亚洲欧美日韩高清专用| 999久久久精品免费观看国产| 18禁裸乳无遮挡免费网站照片| 国产乱人伦免费视频| 国产成人精品久久二区二区91| 久久久久免费精品人妻一区二区| 九色成人免费人妻av| 午夜精品一区二区三区免费看| 亚洲男人天堂网一区| 亚洲人成网站高清观看| 丰满人妻一区二区三区视频av | 女人高潮潮喷娇喘18禁视频| 日本三级黄在线观看| 亚洲,欧美精品.| 亚洲av片天天在线观看| 99久久无色码亚洲精品果冻| 亚洲成人久久性| 国产单亲对白刺激| 久久久久久久久久黄片| 91九色精品人成在线观看| 亚洲欧美一区二区三区黑人| 欧美成人午夜精品| 国产日本99.免费观看| 美女大奶头视频| 久久久久精品国产欧美久久久| 亚洲精品国产精品久久久不卡| 特级一级黄色大片| 毛片女人毛片| 国产男靠女视频免费网站| 丰满人妻熟妇乱又伦精品不卡| АⅤ资源中文在线天堂| 亚洲最大成人中文| 欧美色视频一区免费| 欧美日本亚洲视频在线播放| 91大片在线观看| 亚洲人成77777在线视频| 亚洲国产精品合色在线| 麻豆成人午夜福利视频| 一级a爱片免费观看的视频| 亚洲第一电影网av| 国产午夜精品久久久久久| 757午夜福利合集在线观看| 亚洲天堂国产精品一区在线| 长腿黑丝高跟| 岛国在线观看网站| 麻豆国产av国片精品| 97人妻精品一区二区三区麻豆| 久久人妻av系列| 亚洲九九香蕉| 亚洲人成网站在线播放欧美日韩| 国产成人aa在线观看| avwww免费| 久久 成人 亚洲| cao死你这个sao货| 亚洲全国av大片| 精品久久久久久成人av| 久久久久久国产a免费观看| 亚洲av五月六月丁香网| 一进一出抽搐动态| avwww免费| 天堂av国产一区二区熟女人妻 | 免费一级毛片在线播放高清视频| 国产成人精品久久二区二区91| 12—13女人毛片做爰片一| 色综合欧美亚洲国产小说| 哪里可以看免费的av片| 婷婷亚洲欧美| 黄色片一级片一级黄色片| 久久精品夜夜夜夜夜久久蜜豆 | 91麻豆精品激情在线观看国产| 黑人操中国人逼视频| 亚洲精品久久国产高清桃花| 三级国产精品欧美在线观看 | 成人国产综合亚洲| 91国产中文字幕| 国产在线观看jvid| 亚洲 欧美一区二区三区| 人成视频在线观看免费观看| 又粗又爽又猛毛片免费看| 欧美一级毛片孕妇| 国产黄a三级三级三级人| 成人国产综合亚洲| 黄色片一级片一级黄色片| 色综合亚洲欧美另类图片| 在线观看舔阴道视频| 久久精品成人免费网站| av中文乱码字幕在线| 两个人看的免费小视频| 听说在线观看完整版免费高清| 精品少妇一区二区三区视频日本电影| 成年免费大片在线观看| 亚洲,欧美精品.| 亚洲va日本ⅴa欧美va伊人久久| 制服诱惑二区| 国内揄拍国产精品人妻在线| 精品高清国产在线一区| 麻豆av在线久日| 亚洲美女黄片视频| 亚洲人成网站高清观看| 岛国在线免费视频观看| 亚洲专区中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 极品教师在线免费播放| 色哟哟哟哟哟哟| 国产成人啪精品午夜网站| 天天添夜夜摸| 全区人妻精品视频| 可以免费在线观看a视频的电影网站| 变态另类丝袜制服| 变态另类成人亚洲欧美熟女| 国产视频内射| 99精品欧美一区二区三区四区| 成人国产综合亚洲| 精品无人区乱码1区二区| 男女做爰动态图高潮gif福利片| 中文资源天堂在线| 国产欧美日韩一区二区精品| 草草在线视频免费看| 在线国产一区二区在线| 给我免费播放毛片高清在线观看| 色噜噜av男人的天堂激情| 看黄色毛片网站| 99riav亚洲国产免费| 亚洲自拍偷在线| 中文在线观看免费www的网站 | 日韩中文字幕欧美一区二区|