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

    官廳水庫水質(zhì)特征及水體微生物多樣性的響應(yīng)

    2015-11-19 07:43:32孫寓姣丁愛中趙曉輝張惠淳北京師范大學(xué)水科學(xué)研究院北京00875中國國際工程咨詢有限公司資源與環(huán)境業(yè)務(wù)部北京00048
    中國環(huán)境科學(xué) 2015年5期
    關(guān)鍵詞:官廳水庫水河水華

    孫寓姣,陳 程,2,丁愛中*,趙曉輝,張惠淳(.北京師范大學(xué)水科學(xué)研究院,北京 00875;2.中國國際工程咨詢有限公司資源與環(huán)境業(yè)務(wù)部,北京 00048)

    官廳水庫水質(zhì)特征及水體微生物多樣性的響應(yīng)

    孫寓姣1,陳 程1,2,丁愛中1*,趙曉輝1,張惠淳1(1.北京師范大學(xué)水科學(xué)研究院,北京 100875;2.中國國際工程咨詢有限公司資源與環(huán)境業(yè)務(wù)部,北京 100048)

    針對北京市原水源地官廳水庫,選擇了上下游共6個采樣點進行了不同季節(jié)的水質(zhì)及環(huán)境因子特征分析.總體看來,4個季節(jié)水體氮磷碳營養(yǎng)物質(zhì)濃度依次為:夏季>秋季>春季>冬季,夏秋季出現(xiàn)明顯水華.上游媯水河水質(zhì)要好于下游水庫庫區(qū)水質(zhì),尤其以夏秋季節(jié)最為明顯.夏季上游媯水河的總碳、總氮、氨氮、總磷平均濃度為26.5, 0.95, 0.55, 0.077mg/L,而水庫庫區(qū)響應(yīng)平均濃度分別為111.47, 4.27, 3.16,0.25mg/L.結(jié)合PCR-DGGE對不同時期各點水體微生物進行了群落分析,發(fā)現(xiàn)上下游水體中細菌群落結(jié)構(gòu)差異很大,水華嚴重區(qū)的細菌群落多樣性和豐度要低于比水質(zhì)較好區(qū)域.CANOCA軟件分析了夏秋季水華時水體微生物群落變化與環(huán)境因子的響應(yīng)關(guān)系,發(fā)現(xiàn)細菌群落結(jié)構(gòu)與環(huán)境因子(T、pH值、DO、NFR、TP、Fe)之間存在較強的關(guān)聯(lián).水質(zhì)較好的上游水體夏秋細菌群落結(jié)構(gòu)主要受溫度、DO影響較大;而水質(zhì)較差的庫區(qū)水體夏秋季細菌群落結(jié)構(gòu)變化與固氮速率、pH值、總磷的相關(guān)性較好.

    官廳水庫;水華;環(huán)境因子;水質(zhì);微生物多樣性

    官廳水庫位于北京市西北部,是新中國成立后興建的第一座大型水庫.從20世紀70年代起由于受上游工農(nóng)業(yè)污染,官廳水庫水質(zhì)惡化,富營養(yǎng)化嚴重,1997年退出北京市應(yīng)用水體系[1-3].而先前的監(jiān)測數(shù)據(jù)顯示,總磷(TP)、氨氮(NH4+-N)、總氮(TN)是官廳水庫夏季富營養(yǎng)化的主要貢獻因子[4].近年來,由于北京市飲用短缺,官廳水庫的水體修復(fù)即成為亟待解決的問題之一.水體中氮營養(yǎng)元素的來源主要有人為源及天然源,目前針對官廳水庫氮含量超標的人為來源研究很多,但是對其天然來源的研究則很少.

    作為水生生態(tài)系統(tǒng)的一個重要組成部分,細菌數(shù)量巨大,具有很高的能量轉(zhuǎn)換效率,是水體生態(tài)系統(tǒng)能量流動、營養(yǎng)物質(zhì)交換的重要環(huán)節(jié).在水體環(huán)境發(fā)生變化時,細菌群落結(jié)構(gòu)也會發(fā)生相應(yīng)變化[5-6].因此細菌群落對水環(huán)境質(zhì)量相互關(guān)系十分密切.目前大量研究工作都是針對湖泊中的浮游動植物進行的[7],而水華過程中水體細菌群落的變化應(yīng)是研究水體性質(zhì)變化的一個重要切入點.鄭小紅等[8]用PCR-DGGE技術(shù)對玄武湖的水體進行了細菌群落的變化分析,發(fā)現(xiàn)在水華暴發(fā)期細菌優(yōu)勢種有16種,生物多樣性高,而進入衰退期后優(yōu)勢菌減少,多樣性降低, 但細菌數(shù)量增大.

    本文對官廳水庫不同季節(jié)水體水質(zhì)理化指標進行跟蹤監(jiān)測,探索水庫水質(zhì)年度變化情況.監(jiān)測的水質(zhì)指標主要包括:溫度、pH、DO、、、、TN、TP、BOD5、COD和葉綠素a,以及自然界內(nèi)源的氮輸入方式固氮作用的速率. 分析水體氮磷營養(yǎng)元素含量變化的原因.利用PCR-DGGE技術(shù)擴增細菌16srDNA基因來對官廳水庫不同區(qū)域水體中細菌的群落變化進行了研究,以期發(fā)現(xiàn)水體富營養(yǎng)化過程中微生物群落結(jié)構(gòu)的響應(yīng)關(guān)系.進而為水華水質(zhì)變化及水體中微觀生態(tài)代謝過程深入討論提供理論依據(jù).

    1 材料與方法

    1.1 采樣點設(shè)置

    水樣的采集于7月20日水華發(fā)生時進行,本次采樣共設(shè)6個采樣點(S1~S6)如圖1所示.其中S1、S2、S3位于官廳水庫上游媯水河,S4、S5、S6位于水庫庫區(qū),采樣過程中使用GPS系統(tǒng)對每個采樣點進行精確定位,采樣點位置S1~S6經(jīng)、緯度依次為115°46′5″,40°21′58″; 115°42′24″,40°20′51″;115°42′9″,40°20′15;115°36′75″,40°14′1 3;15°35′57″,40°14′7;115°25′39″, 40°14′58.

    1.2 水樣采集

    有機玻璃采樣器采集水深0.5m左右處水樣.用于測定固氮速率的水樣注入經(jīng)鹽酸浸泡洗凈、水樣反復(fù)潤洗的玻璃瓶中,低溫避光帶回實驗室.

    圖1 官廳水庫采樣點示意Fig.1 Location of sampling sites in Guanting Reservoir

    1.3 水樣分析

    1.3.1 水體理化參數(shù)的測定 水溫、pH值、溶解氧(DO)濃度現(xiàn)場測定,理化指標當(dāng)天送樣品至北京譜尼測試公司測定,包括-N,-N,總氮(TN), 總磷(TP), CODcr, Fe等.葉綠素a(Chl-a)濃度測定在北京理化測試中心完成.

    1.3.2 水體固氮速率測定 采用乙炔還原法測定水體固氮速率[9-10].取1500mL水樣8500×g離心20min 后棄上清135mL即將水樣濃縮10倍, 取濃縮水樣50mL于65mL磨口三角瓶中,用注射器抽取頂端15%的空氣,并注入等量乙炔氣體,同時用無菌水設(shè)置對照.27℃溫育72h后加入0.1mL50%TCA終止反應(yīng).取100μL頂端氣體,氣相色譜檢測乙炔含量,換算成氮含量表示固氮速率.

    1.3.3 總DNA提取及固氮16s基因擴增 水樣于24h內(nèi)經(jīng)0.22μm醋酸纖維素濾膜過濾,濃縮生物樣品.使用Omega Water DNA Kit提取水體微生物的總DNA.

    使用16SrDNA V3區(qū)通用引物341F-GC和907R(合成于上海生工生物公司)對水樣基因組總DNA擴增.在上游引物341f的5'端前面加了一段“GC發(fā)卡”結(jié)構(gòu)使后續(xù)的DGGE能有更好的分辨率(Poly et al., 2001).引物序列(GC_341f:5'-CGC CCG CCG CGC GCG GCG GGC GGG GCG GGG GCA CGG GGG GCC TAC GGG AGG CAG CAG-3';907r:5'-CCG TCA ATT CMT TTG AGT TT-3').以1μL DNA樣品為模板,PCR擴增體系含有:5.0μL 10×PCR緩沖液(Takara),4.0μL 2.5mM dNTP混合液,20mM引物各1μL,0.25μL (0.5U/μL) Taq酶 (Takara),無菌水補齊至50μL.采用降落PCR技術(shù)擴增,擴增程序為:95℃/5min;94℃/30s,63℃/60s,72℃/90s,8個循環(huán),每個循環(huán)降低1℃至56℃;94℃/30s,56℃/60s,72℃/ 90s; 25個循環(huán),72℃/7min.用Omega PCR產(chǎn)物純化試劑盒對產(chǎn)物進行切膠回收[6].

    1.3.4 16s基因擴增片段變性梯度凝膠電泳(DGGE) 聚丙烯酰胺的變性梯度范圍為35%~65%,DNA擴增產(chǎn)物上樣量約為200ng,60℃恒溫,80V恒壓條件下電泳13~15h,電泳完畢后采用銀染技術(shù)染色,白光成像.隨后使用Quantity One軟件對DGGE圖譜進行分析,分析各泳道條帶數(shù)目及灰度. 16s基因的多樣性使用Shannon-Weaver指數(shù)H進行表征.計算公式為:

    式中: Pi為每條泳道中第i條條帶灰度占該樣品總灰度的比率;S為泳道中含有的條帶總數(shù).

    所有實驗數(shù)據(jù)均在SPSS for windows (15.0)統(tǒng)計軟件上進行處理.

    2 結(jié)果與分析

    2.1 基本理化指標全年變化

    官廳水庫水體溫度全年變化較大,春季水體平均溫度為4.9℃,夏季為25.9℃,秋季降低到21.9℃,冬季冰面下水體平均溫度為3.6℃.水體pH值全年呈堿性,平均值在8.7~9之間,其中秋季水體堿性最大.水體溶解氧春、冬季平均值較高,分別為7.9mg/L,8.7mg/L.,夏季為6.6mg/L,為全年最低,其中媯水河溶解氧平均濃度為7.16mg/L,水庫庫區(qū)僅為6.04mg/L.秋季水體溶解氧濃度相對夏季略有上升,媯水河溶解氧平均濃度上升到8.01mg/L,庫區(qū)平均濃度上升到6.89mg/L.

    2.2 水體葉綠素a全年變化

    在飲用水源地和景觀娛樂水體,發(fā)生水華的界限通常為葉綠素a(Chl-a)含量超過10mg/m3或藻類細胞密度大于20000個/mL[11].通過對Chl-a濃度的監(jiān)測,反映了水庫全年藻類的生長情況.由圖2可以看到,官廳水庫4個季節(jié)Chl-a的平均濃度分別為1.87,41.79,29.82,2.86mg/m3,夏季秋季水庫Chl-a含量超高,水華嚴重.

    圖2 官廳水庫葉綠素a全年變化Fig.2 Chl-a concentration of different water samples

    2.3 有機物濃度全年變化

    圖3 官廳水庫全年COD變化Fig.3 COD concentration of different water samples

    從時間上看,官廳水庫四季水質(zhì)變化較大(如圖3),尤其是夏季,各項水質(zhì)指標均嚴重超標.其中夏季6個采樣點 CODcr平均值為68.99mg/L,庫區(qū)S5采樣點處CODcr值最高達到245mg/L,明顯高于其他3個季節(jié).從空間上看,媯水河流域采樣點(S1、S2、S3)和庫區(qū)采樣點(S4、S5、S6)水質(zhì)差異明顯.庫區(qū)CODcr濃度遠高于媯水河流域.其中夏季媯水河流域采樣點平均CODcr濃度為26.5mg/L,而庫區(qū)采樣點濃度值高達111.47mg/L.

    2.4 氮、磷營養(yǎng)濃度全年變化

    圖4 官廳水庫總氮,氨氮及總磷全年變化Fig.4 Concentration of total nitrogen, ammonia, and total phosphorus of different water samples in Guanting resevior

    官廳水庫已經(jīng)長期處于富營養(yǎng)狀態(tài),氮、磷是水庫水體主要營養(yǎng)污染物(圖4).通過全年監(jiān)測發(fā)現(xiàn),春夏秋冬水體總氮的平均濃度分別為:1.14,2.61, 1.29, 0.62mg/L,冬季總氮含量最低,夏季最高,已超過地表V類水標準(圖4A).乙炔還原法對官廳水庫4個季節(jié)所采水樣進行了固氮速率的檢測,只有夏季6個采樣點監(jiān)測到固氮速率,水華爆發(fā)的采樣點S5處固氮速率最高,達到9.35nmolN/(m3·d),其他5個采樣點的固氮速率都小于1nmolN/(m3·d).相比之下,富營養(yǎng)化水平低的媯水河流域固氮速率較庫區(qū)低很多,最小值0.172nmolN/(m3·d)出現(xiàn)在媯水河S3處.

    官廳水庫無機氮主要以氨氮形式存在,春季水體氨氮平均濃度為0.33mg/L,夏季上升至1.86mg/L(圖4B),最大值為8.18mg/L在S5點處測得;秋季水體氨氮平均濃度又降至0.63mg/L,最大值2.77mg/L仍在S5測得;冬季氨氮濃度平均濃度降低為0.16mg/L.

    通過對水體總磷濃度的監(jiān)測可知(圖4C),官廳水庫4個季節(jié)各采樣點總磷平均濃度分別為0.06, 0.16, 0.14, 0.058mg/L.即使在冬季,水庫濃度都超過地表Ⅲ類水標準,夏季總磷濃度超過地表Ⅳ類水標準.4個季節(jié)水體氮磷營養(yǎng)物濃度依次為:夏季>秋季>春季>冬季.另外從空間分布看,庫區(qū)采樣點總氮、氨氮、總磷濃度明顯高于媯水河各點.例如,在三者平均濃度最高的夏季,上游媯水河3個采樣點的總氮、氨氮、總磷平均濃度為0.95,0.55,0.077mg/L,表明媯水河水質(zhì)介于地表湖庫Ⅱ類-Ⅲ類水之間.而水庫庫區(qū)3個采樣點的總氮、氨氮、總磷平均濃度則分別達到了4.27,3.16, 0.25mg/L,均超過地表湖庫Ⅴ類水標準.這表明,官廳水庫庫區(qū)水體水質(zhì)及營養(yǎng)狀態(tài)比媯水河嚴重,尤其以夏秋季最為明顯.

    2.5 水樣DNA的提取及16s基因PCR-DGGE結(jié)果擴增分析

    為了增加DNA的提取效率,在使用試劑盒提取水體DNA樣品的過程中增加了機械力細胞破碎的步驟使水體微生物DNA盡量多的釋放.提取出的DNA樣品經(jīng)過1.2%的瓊脂糖凝膠電泳檢驗,片段長度均在23kb左右, OD260/280在1.6~1.8之間.細菌16S rDNA PCR目標片段的條帶清晰,位置正確,擴增產(chǎn)物可用于后續(xù)DGGE分析.樣品的PCR產(chǎn)物經(jīng)DGGE分析后,官廳水庫細菌群落DGGE結(jié)果見圖5.

    圖5 官廳水庫四季細菌群落DGGE圖譜Fig.5 DGGE band profiles of 16s gene iamplified using four seasons genomic DNA

    從圖5可以看到,21個樣品的PCR擴增產(chǎn)物經(jīng)DGGE分析后,各泳道中均分離出了不同的條帶,表明DGGE可以很好的將樣品中不同菌種區(qū)分開來.泳道條帶數(shù)目代表各采樣點微生物豐度,21個采樣點形成了數(shù)目不同的條帶,且條帶的亮度和遷移速率不同,表示官廳水庫不同采樣點細菌豐度和群落結(jié)構(gòu)不同.

    Quantity One 軟件對DGGE所得圖譜經(jīng)過泳道識別、去除背景、識別條帶、高斯模型分析之后得到各泳道的條帶數(shù)量(即豐度S),同時得出各條帶的光密度(即條帶灰度),最后根據(jù)Shannon-Weaver公式可以計算出多樣性指數(shù)H.得到官廳水庫全年細菌豐度和多樣性的變化如圖6.各采樣點微生物多樣性的計算結(jié)果顯示,各點微生物多樣性差別較大.H的數(shù)值在1.47~2.79之間.

    圖6 官廳水庫細菌豐度與多樣性變化Fig.6 Richness and diversities of bacteria in different water samples of Guanting resevior

    夏季隨著水體氮磷營養(yǎng)的升高,水體細菌豐度和多樣性上升,總體看來,水質(zhì)較好的媯水河水體中細菌的多樣性和豐度都好于水華嚴重的庫區(qū)水體.9月份秋季水華期間大多數(shù)點的細菌豐度和水體多樣性均有所下降.

    2.6 夏、秋季細菌群落與環(huán)境因子相關(guān)性

    在獲得官廳水庫全年細菌群落變化的基礎(chǔ)上,重點分析7月和9月水華期間細菌群落結(jié)構(gòu)的變化,并利用CANOCO軟件分析其與環(huán)境因子的相關(guān)性.結(jié)果見圖7.

    由圖7可見,媯水河采樣點(S1,S2,S3)在夏、秋季細菌群落結(jié)構(gòu)非常相似,被歸為一支,6個樣品之間的相似度大于57.4%.其中夏季媯水河S2、S3采樣點細菌群落結(jié)構(gòu)更為相似(75.3%),秋季媯水河S1、S3采樣點細菌群落結(jié)構(gòu)更為相似.庫區(qū)采樣點夏、秋兩季細菌群落結(jié)構(gòu)較為相似,夏季S4、S5細菌群落結(jié)構(gòu)相似度為48.6%.秋季S5、S6群落相似度為48%.由此可見夏季官廳水庫S6點細菌群落結(jié)構(gòu)較為特殊,秋季S4點細菌群落結(jié)構(gòu)較為特殊.

    圖7 官廳水庫夏季、秋季細菌群落DGGE指紋圖譜UPGMA聚類分析Fig.7 UPGMA of DGGE fringerprint of summer and autumn water samples of Guanting resevior

    總體來看,夏秋季水華期間媯水河和庫區(qū)細菌群落結(jié)構(gòu)差異一直很大,而各自的細菌結(jié)構(gòu)在兩次水華期間也有所不同,尤其庫區(qū)采樣點7月和9月的細菌群落結(jié)構(gòu)差異很大.

    對7月、9月官廳水庫水華期間細菌群落進行聚類分析之后,利用CANOCO for Windows 4.5軟件對DGGE結(jié)果和環(huán)境因子進行多變量統(tǒng)計分析.首先利用DCA對DGGE圖譜代表的細菌群落信息進行分析,得到的最長的梯度長度值為2.639,表明應(yīng)選擇線性模型RDA對樣品進行排序分析.水樣理化指標剔除膨脹因子大于30的指標之后,進行了標準化處理[12].典型對應(yīng)分析結(jié)果概括于表1.

    根據(jù)膨脹因子選擇了溫度(T)、pH值、DO、固氮速率(NFR)、總磷(TP)、Fe作為環(huán)境因子表征水體環(huán)境變化情況,Monte Carlo檢驗顯示6個環(huán)境參數(shù)與AX1軸(P=0.026)及全部排序軸(P=0.022)均有顯著的相關(guān)性.具體對應(yīng)結(jié)果(表1)表明第1排序軸解釋了樣本中24.4%的變異,第2和第3排序軸分別解釋了樣本中41.2%和49.3%的變異,前4個排序軸共解釋了55.5%的樣本總變異,表明排序圖能很好的解釋官廳水庫夏季和秋季的細菌變異情況.從物種-環(huán)境的相關(guān)性看,第1和第2排序軸的種?環(huán)境相關(guān)系數(shù)分別為0.970和0.984,同時,前4個排序軸的物種?環(huán)境累積百分率高達85.9%,這說明水樣中細菌群落結(jié)構(gòu)與環(huán)境因子間存在較強的關(guān)聯(lián).水樣中微生物群落結(jié)構(gòu)主要與DO、T、TP、NFR相關(guān)聯(lián),其中DO(R2=0.5306)、T(R2=0.3708)與第1排序軸正相關(guān),TP(R2=-0.4699)與第1排序軸負相關(guān);T(R2=0.7490)、NFR(R2=0.6848)與第2排序軸正相關(guān),pH(R2=-0.6005)與第2排序軸負相關(guān),鐵與第1、2個排序軸均有一定的正相關(guān)性(R2= 0.1466, 0.3207),但相關(guān)程度不高.

    表1 官廳水庫夏季、秋季細菌群落與環(huán)境因子典型對應(yīng)分析結(jié)果Table 1 Data of cooresponding relations between bacterial community and environmental factors

    圖8 官廳水庫夏季、秋季細菌群落與環(huán)境因子典型對應(yīng)分析Fig.8 Typical cooresponding relations between bacterial community and environmental factors

    物種環(huán)境排序圖中2個樣品距離的遠近表明樣品物種群落的相似程度(圖8),環(huán)境因子的箭頭方向與采樣點之間的夾角表示其和該點物種群落之間的相關(guān)性,夾角越小表示相關(guān)性越高,兩者箭頭若同向,表明正相關(guān),反向表明負相關(guān).從圖8可以看到,夏季官廳水庫上游媯水河3個采樣點距離較近,而其和庫區(qū)3個采樣點的距離遠,表明7月水庫上游媯水河和下游庫區(qū)的細菌群落結(jié)構(gòu)差異很大,這和DGGE聚類分析的結(jié)果一致.由于夏季水溫較高,營養(yǎng)物質(zhì)較豐富的環(huán)境中微生物生長較旺盛,會導(dǎo)致不同環(huán)境出現(xiàn)較大的微生物群落差異.從采樣點和環(huán)境因子的夾角來看,夏季7月媯水河細菌群落受水溫的影響較大,庫區(qū)細菌群落和NFR相關(guān)性較好.秋季9月水庫6個采樣點較為分散,表明9月水庫細菌群落結(jié)構(gòu)差異大,其中媯水河細菌群落受DO影響較大,庫區(qū)細菌群落受pH值和TP影響較大.可見在水質(zhì)較好地水域,非人為因素(溫度等)是影響水體中的微生物群落的關(guān)鍵;而存在營養(yǎng)化污染區(qū)域的水體,TP、pH值等相關(guān)人為因素的指標,成為影響微生態(tài)群落變化的主導(dǎo)因子.

    3 結(jié)論

    3.1 官廳水庫四季水質(zhì)變化較大,庫區(qū)水質(zhì)劣于上游媯水河水質(zhì),夏季媯水河水質(zhì)介于地表湖庫Ⅱ類-Ⅲ類之間,而庫區(qū)水質(zhì)超過地表湖庫Ⅴ類水標準.夏季和秋季庫區(qū)均存在大量的水華.

    3.2 用綜合營養(yǎng)狀態(tài)指數(shù)法對官廳水庫的富營養(yǎng)狀況進行評價,結(jié)果表明官廳水庫全年水體均已經(jīng)富營養(yǎng)化.水質(zhì)較好的媯水河水體中細菌的多樣性和豐度都好于水華嚴重的庫區(qū)水體. 9月秋季水華期間大多數(shù)點的細菌豐度和水體多樣性均比夏季有所下降.

    3.3 對7月和9月水華期間細菌群落和環(huán)境因子進行典型對應(yīng)分析,結(jié)果表明細菌群落結(jié)構(gòu)和環(huán)境因子(T、pH值、DO、NFR、TP、Fe)之間存在較強的關(guān)聯(lián).樣品中微生物群落結(jié)構(gòu)主要與DO、T、TP、NFR相關(guān),pH值、鐵與細菌群落結(jié)構(gòu)關(guān)聯(lián)度最小.7月水華期間媯水河和庫區(qū)的細菌群落結(jié)構(gòu)差異很大,媯水河細菌群落受水溫的影響較大.9月水庫各點的細菌群落較為分散,群落差異大,其中媯水河細菌群落受DO影響較大,庫區(qū)細菌群落受pH值和TP影響較大.

    [1]杜桂森,王建廳,張為華,等.官廳水庫水體營養(yǎng)狀況分析 [J]. 湖泊科學(xué), 2004,(3):277-281.

    [2]孫 峰,郝芳華.基于GIS的官廳水庫流域非點源污染負荷計算研究 [J]. 北京水利, 2004,(1):16-18.

    [3]Dai R H, Liu H J, Qu J H, et al. Cyanobacteria and their toxins in Guanting Reservoir of Beijing, China [J]. Journal of Hazard Mater, 2008,153:470-477.

    [4]楊大杰.官廳水庫水體氮污染特征分析 [J]. 中國水利, 2008,(9):51-53.

    [5]Riemann L, Steward G F, Azam F. Dynamics of bacterial community composition and activity during a mesocosm diatom bloom [J]. Applied and Environmenatl Microbiology, 2000,66: 578-587.

    [6]Fandino L B, Riemann L, Steward G F, et al. Variations in bacterial community structure during a dinoflagellate bloom analyzed by DGGE and 16SrDNA sequencing [J]. Aquatic Microbial Ecology, 2001,23:119-130.

    [7]Eiler A, Bertilsson S. Composition of freshwater bacterial communities associated with cyanobacterial blooms in four Swedish lakes [J]. Environmental Microbiology, 2004,6:1228-1243.

    [8]鄭小紅,肖 琳,任 晶,等.玄武湖微囊藻水華暴發(fā)及衰退期細菌群落變化分析 [J]. 環(huán)境科學(xué), 2008,29(10):2956-2962.

    [9]Stewart W D, Fitzgerald G P, Burris R H. In situ studies on N2fixation using the acetylene reduction technique [J]. PNAS,1967,58(5):2071-2078.

    [10]Présing M, V-Balogh K., V?r?s L, et al. Relative nitrogen deficiency without occurrence of nitrogen fixing bluegreen algae in a hypertrophic reservoir [J]. Hydrobiologia, 1997,342:55-61.

    [11]Oliver R L, Ganf G G. Freshwater blooms [M]//Whitton B A,Potts M. (Eds). The ecology of cyanobacteria. NL. Kluwer Academic Publisher, 2000,149-194.

    [12]Yannarell A C, Triplett E W. Geographic and environmental sources of variation in lake bacterial community composition [J]. Applied and Environmental Microbiology, 2005,71(1):227-239.

    The corresponding of microbial diversity on water quality and environmental variables of Guanting Reservoir.

    SUN Yu-jiao1, CHEN Cheng1,2, DING Ai-zhong1*, ZHAO Xiao-hui1, ZHANG Hui-chun1(1.College of Water Science,Beijing Normal University, Beijing 100875, China;2.China International Engineering Consulting Corporation Resources and Environment Business Department, Beijing 100048, China). China Environmental Science, 2015,35(5):1547~1553

    The water quality and environmental factors of Beijing Guanting Reservoir in different seasons were analyzed based on the characterization of the data from six sampling sites. Generally, the concentrations of the nutrients including carbon, nitrogen, and phosphorus were in the orders of: summer > autumn > spring > winter, and the outbreaks of water-blooms occurred in summer and autumn. The water quality in the upstream (Gui River), was better than that of the downstream (reservoir area), especially in summer and autumn. In summer, the average concentrations of total carbon,total nitrogen, ammonia nitrogen, total phosphorus in upstream were 26.5, 0.95, 0.55, 0.077mg/L respectively, while that of downstream were 111.47, 4.27, 3.16, 0.25mg/L respectively. Microbial community of the water was analyzed employing PCR-DGGE technique. The result showed that bacterial community structures varied considerably between the upstream and the downstream. The upstream water behaved the higher level of bacterial community diversity and richer degree than that of downstream water serious bloom area. In addition, CANOCA software was used to explore the response relationship between microbial community and environmental factors in summer and autumn bloom periods. The result showed strong correlation in response of bacterial community structure with environmental factors (T, pH, DO,NFR, TP, Fe). For upstream water, the better water quality was in the bacterial community structures mainly affected by the temperature and DO factor; while for the downstream, the poor water quality was with the bacteria community structures effected mainly by environmental factors of pH, TP and NFR.

    Guanting Reservoir;water bloom;environmental factors;water quality;microbiological diversity

    X703.5

    A

    1000-6923(2015)05-1547-07

    孫寓姣(1975-),女,哈爾濱人,副教授,博士,主要從事環(huán)境微生態(tài)學(xué)及環(huán)境污染生物處理技術(shù)研究.發(fā)表論文30余篇.

    2014-09-31

    國家自然科學(xué)基金(51178048,51378064);北京師范大學(xué)自主基金(2014KJJCB22)

    * 責(zé)任作者, 教授, ading@bnu.edu.cn

    猜你喜歡
    官廳水庫水河水華
    我愛家鄉(xiāng)的宜水河
    藻類水華控制技術(shù)及應(yīng)用
    鑒水河之歌
    歌海(2021年2期)2021-06-22 02:25:59
    我的水河
    南美白對蝦養(yǎng)殖池塘藍藻水華處理舉措
    南美白對蝦養(yǎng)殖池塘藍藻水華處理舉措
    我的洣水河
    官廳水庫夏季葉綠素a濃度短時分布特征
    官廳水庫跨界流域管理體制研究
    綠色科技(2017年6期)2017-04-20 10:28:01
    油酸酰胺去除藍藻水華的野外圍隔原位試驗
    色综合亚洲欧美另类图片| 亚洲国产精品sss在线观看| 国产精品综合久久久久久久免费| av福利片在线观看| 真人做人爱边吃奶动态| 一卡2卡三卡四卡精品乱码亚洲| 精品一区二区三区四区五区乱码| 精品熟女少妇八av免费久了| 亚洲国产欧美一区二区综合| 一个人免费在线观看电影 | 亚洲成av人片免费观看| 又大又爽又粗| 亚洲精品久久国产高清桃花| 欧美性猛交╳xxx乱大交人| 男女午夜视频在线观看| 熟妇人妻久久中文字幕3abv| 啦啦啦韩国在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜激情福利司机影院| 十八禁人妻一区二区| 亚洲精品久久国产高清桃花| 中亚洲国语对白在线视频| 91在线精品国自产拍蜜月 | 99在线人妻在线中文字幕| 国产精品久久久人人做人人爽| 日日夜夜操网爽| 久久久国产成人精品二区| 最新在线观看一区二区三区| 国产成人aa在线观看| 国产美女午夜福利| 可以在线观看的亚洲视频| 欧美成狂野欧美在线观看| 国产亚洲精品久久久com| 午夜福利视频1000在线观看| 性色avwww在线观看| 国产精品久久久久久人妻精品电影| 床上黄色一级片| 成年女人永久免费观看视频| 亚洲av美国av| 国产一区二区激情短视频| 成人特级av手机在线观看| 久久久久久国产a免费观看| 国产av不卡久久| 国产一区二区激情短视频| 一二三四社区在线视频社区8| 午夜免费观看网址| 亚洲美女视频黄频| 亚洲成人中文字幕在线播放| 最好的美女福利视频网| 女人高潮潮喷娇喘18禁视频| 91字幕亚洲| 99国产极品粉嫩在线观看| 午夜福利在线在线| 午夜免费成人在线视频| av中文乱码字幕在线| 亚洲自拍偷在线| av国产免费在线观看| 黄色女人牲交| 床上黄色一级片| 在线观看美女被高潮喷水网站 | 精品99又大又爽又粗少妇毛片 | 亚洲精品乱码久久久v下载方式 | 久久久色成人| 亚洲成av人片免费观看| 午夜精品在线福利| cao死你这个sao货| 最近最新免费中文字幕在线| 免费看十八禁软件| 亚洲精华国产精华精| 长腿黑丝高跟| 国产精品自产拍在线观看55亚洲| 欧美日韩中文字幕国产精品一区二区三区| 九九热线精品视视频播放| 俄罗斯特黄特色一大片| 日本a在线网址| 中亚洲国语对白在线视频| 国产精品 国内视频| 亚洲五月天丁香| 最近最新免费中文字幕在线| 国产精华一区二区三区| 精品熟女少妇八av免费久了| 在线看三级毛片| 国产亚洲欧美98| 欧美另类亚洲清纯唯美| 法律面前人人平等表现在哪些方面| 欧美国产日韩亚洲一区| 国产又黄又爽又无遮挡在线| 两个人看的免费小视频| 99久久精品一区二区三区| 超碰成人久久| 精品日产1卡2卡| 精品午夜福利视频在线观看一区| 特大巨黑吊av在线直播| 国产av一区在线观看免费| 91老司机精品| 人人妻人人看人人澡| 成人永久免费在线观看视频| 精品福利观看| 精品电影一区二区在线| 人人妻,人人澡人人爽秒播| 中文字幕人成人乱码亚洲影| 国产蜜桃级精品一区二区三区| 丰满的人妻完整版| 哪里可以看免费的av片| 极品教师在线免费播放| 夜夜爽天天搞| 成人性生交大片免费视频hd| 久久这里只有精品19| 久久久国产成人精品二区| 熟女少妇亚洲综合色aaa.| 级片在线观看| 亚洲欧美日韩高清在线视频| 中文资源天堂在线| 国产精品国产高清国产av| 亚洲人成网站在线播放欧美日韩| 91字幕亚洲| 2021天堂中文幕一二区在线观| 无遮挡黄片免费观看| 国产精品一及| 偷拍熟女少妇极品色| 全区人妻精品视频| 亚洲中文字幕一区二区三区有码在线看 | 黑人欧美特级aaaaaa片| 无限看片的www在线观看| 国产亚洲精品久久久久久毛片| 亚洲人成网站在线播放欧美日韩| 此物有八面人人有两片| 色噜噜av男人的天堂激情| 亚洲av成人av| 亚洲av第一区精品v没综合| 一边摸一边抽搐一进一小说| 午夜精品久久久久久毛片777| 色av中文字幕| 国内精品久久久久精免费| 亚洲午夜精品一区,二区,三区| 日本成人三级电影网站| 12—13女人毛片做爰片一| 亚洲美女黄片视频| 国产91精品成人一区二区三区| 国产午夜精品久久久久久| 亚洲av第一区精品v没综合| 1024手机看黄色片| 性欧美人与动物交配| 久久久精品大字幕| 国产亚洲精品av在线| 亚洲成人免费电影在线观看| 男人舔女人的私密视频| 色噜噜av男人的天堂激情| 欧美黑人欧美精品刺激| av在线天堂中文字幕| 美女cb高潮喷水在线观看 | 日韩免费av在线播放| 一a级毛片在线观看| 久久精品影院6| 18禁美女被吸乳视频| 久久午夜亚洲精品久久| e午夜精品久久久久久久| 日本在线视频免费播放| 精品久久久久久久人妻蜜臀av| 国产精品一区二区三区四区免费观看 | 黄色女人牲交| 国产精品日韩av在线免费观看| 夜夜爽天天搞| 99在线视频只有这里精品首页| 色综合欧美亚洲国产小说| 亚洲人成网站高清观看| 国产精品亚洲美女久久久| 精品国产超薄肉色丝袜足j| 国产精品久久久av美女十八| 久久中文看片网| av在线天堂中文字幕| 丝袜人妻中文字幕| 午夜激情福利司机影院| 久久午夜综合久久蜜桃| 一个人看视频在线观看www免费 | 午夜两性在线视频| 人人妻,人人澡人人爽秒播| 18禁国产床啪视频网站| av在线天堂中文字幕| 国产精华一区二区三区| 夜夜夜夜夜久久久久| 两个人看的免费小视频| 又黄又粗又硬又大视频| 黑人巨大精品欧美一区二区mp4| 一级黄色大片毛片| 国产精品影院久久| 久久精品影院6| 免费搜索国产男女视频| 麻豆国产av国片精品| 中出人妻视频一区二区| 99热这里只有是精品50| 香蕉久久夜色| 日本黄大片高清| 成人三级黄色视频| 久久香蕉国产精品| 欧美日韩福利视频一区二区| 欧美最黄视频在线播放免费| 此物有八面人人有两片| 亚洲国产精品sss在线观看| 黄片小视频在线播放| 欧美日韩中文字幕国产精品一区二区三区| 最近视频中文字幕2019在线8| 国产亚洲精品久久久com| 国内久久婷婷六月综合欲色啪| 欧美绝顶高潮抽搐喷水| 一区二区三区激情视频| 亚洲成人中文字幕在线播放| 国产精品爽爽va在线观看网站| 成人性生交大片免费视频hd| 国产av一区在线观看免费| 麻豆成人av在线观看| 久久中文看片网| 日本一二三区视频观看| 亚洲人成网站高清观看| 老熟妇仑乱视频hdxx| 在线观看美女被高潮喷水网站 | 久久午夜综合久久蜜桃| 天堂av国产一区二区熟女人妻| 色在线成人网| 亚洲国产欧美一区二区综合| 男女那种视频在线观看| 一区二区三区激情视频| 床上黄色一级片| 成年免费大片在线观看| 老熟妇乱子伦视频在线观看| 九九久久精品国产亚洲av麻豆 | 色噜噜av男人的天堂激情| 久久久久久久久久黄片| 欧美日韩一级在线毛片| 亚洲av片天天在线观看| 1000部很黄的大片| 99热只有精品国产| 一夜夜www| x7x7x7水蜜桃| 啦啦啦韩国在线观看视频| 黄色成人免费大全| 亚洲成人中文字幕在线播放| 天堂动漫精品| 天天躁日日操中文字幕| 特大巨黑吊av在线直播| 国产伦精品一区二区三区四那| 亚洲成人免费电影在线观看| 18禁黄网站禁片免费观看直播| 18禁美女被吸乳视频| 亚洲人成伊人成综合网2020| 久久人妻av系列| 女同久久另类99精品国产91| 美女大奶头视频| 99国产综合亚洲精品| 一本一本综合久久| 久久久久国产一级毛片高清牌| 国产成+人综合+亚洲专区| 国产精品影院久久| 九色成人免费人妻av| 欧美日韩福利视频一区二区| 天堂av国产一区二区熟女人妻| 午夜免费成人在线视频| 免费观看人在逋| 最新在线观看一区二区三区| 成人欧美大片| 999精品在线视频| 欧美午夜高清在线| 国产伦人伦偷精品视频| bbb黄色大片| 亚洲欧美日韩高清在线视频| 亚洲人成网站在线播放欧美日韩| 亚洲片人在线观看| 最近最新中文字幕大全免费视频| 很黄的视频免费| 国产激情欧美一区二区| 色哟哟哟哟哟哟| 亚洲av美国av| 这个男人来自地球电影免费观看| 亚洲国产精品久久男人天堂| www日本在线高清视频| 色在线成人网| 一夜夜www| 久久草成人影院| 两性午夜刺激爽爽歪歪视频在线观看| 午夜影院日韩av| 一个人免费在线观看的高清视频| 欧美日韩一级在线毛片| 国产欧美日韩精品一区二区| 日本黄色视频三级网站网址| 美女黄网站色视频| 午夜精品在线福利| 免费观看精品视频网站| 老汉色∧v一级毛片| 久久久久久久久久黄片| 一边摸一边抽搐一进一小说| av黄色大香蕉| 又紧又爽又黄一区二区| 国产精品久久久人人做人人爽| 精品久久久久久久末码| 欧美国产日韩亚洲一区| 亚洲欧美精品综合一区二区三区| 91麻豆精品激情在线观看国产| 青草久久国产| 精品一区二区三区av网在线观看| 熟女电影av网| 久久久久久久久久黄片| 无人区码免费观看不卡| 99久久精品一区二区三区| 无遮挡黄片免费观看| 99国产综合亚洲精品| 老司机午夜福利在线观看视频| 欧美高清成人免费视频www| 午夜福利在线观看免费完整高清在 | 午夜两性在线视频| 久久精品国产亚洲av香蕉五月| 18禁美女被吸乳视频| 狂野欧美激情性xxxx| 欧美绝顶高潮抽搐喷水| 最新中文字幕久久久久 | 最新在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看 | 97人妻精品一区二区三区麻豆| 制服丝袜大香蕉在线| 色噜噜av男人的天堂激情| 国产精品日韩av在线免费观看| aaaaa片日本免费| 亚洲国产欧美网| 一a级毛片在线观看| 久久精品91无色码中文字幕| 国产aⅴ精品一区二区三区波| 午夜福利在线观看免费完整高清在 | www.自偷自拍.com| 一个人免费在线观看的高清视频| 亚洲精品456在线播放app | 亚洲 欧美 日韩 在线 免费| АⅤ资源中文在线天堂| 国产麻豆成人av免费视频| 国产伦在线观看视频一区| 欧美在线黄色| 午夜福利在线观看吧| 日韩欧美国产在线观看| 亚洲中文字幕日韩| 长腿黑丝高跟| 在线免费观看不下载黄p国产 | 日韩欧美国产一区二区入口| 国产av不卡久久| 午夜激情欧美在线| 久久久久免费精品人妻一区二区| 午夜福利18| 亚洲黑人精品在线| 真实男女啪啪啪动态图| 老熟妇仑乱视频hdxx| 亚洲av日韩精品久久久久久密| 又爽又黄无遮挡网站| 在线观看66精品国产| 日韩三级视频一区二区三区| 久久精品国产99精品国产亚洲性色| 欧美不卡视频在线免费观看| 欧美黑人欧美精品刺激| xxxwww97欧美| 激情在线观看视频在线高清| 美女黄网站色视频| 99riav亚洲国产免费| 亚洲av成人精品一区久久| av在线蜜桃| 亚洲av第一区精品v没综合| 国产av一区在线观看免费| 99久久久亚洲精品蜜臀av| 国产精品自产拍在线观看55亚洲| 精品久久久久久成人av| 色哟哟哟哟哟哟| 亚洲狠狠婷婷综合久久图片| 亚洲第一欧美日韩一区二区三区| 一区二区三区激情视频| 90打野战视频偷拍视频| 欧美成人免费av一区二区三区| 我的老师免费观看完整版| 国产精品女同一区二区软件 | 国产 一区 欧美 日韩| 国产高清视频在线播放一区| 91av网站免费观看| 日韩欧美在线二视频| 国产高清视频在线播放一区| 日日干狠狠操夜夜爽| 久久久久精品国产欧美久久久| 99久久久亚洲精品蜜臀av| 免费大片18禁| 成人三级做爰电影| 国语自产精品视频在线第100页| 丰满人妻一区二区三区视频av | 变态另类丝袜制服| 国产69精品久久久久777片 | 日韩av在线大香蕉| 制服人妻中文乱码| 亚洲 国产 在线| 最近最新中文字幕大全免费视频| 久久热在线av| 精品人妻1区二区| a在线观看视频网站| 久久精品国产综合久久久| 久久亚洲精品不卡| 别揉我奶头~嗯~啊~动态视频| 激情在线观看视频在线高清| 波多野结衣高清作品| 很黄的视频免费| 精品不卡国产一区二区三区| 琪琪午夜伦伦电影理论片6080| 欧美中文日本在线观看视频| 精品不卡国产一区二区三区| 国产精品爽爽va在线观看网站| 精品午夜福利视频在线观看一区| 999久久久精品免费观看国产| 一个人免费在线观看电影 | 久久亚洲精品不卡| 别揉我奶头~嗯~啊~动态视频| 欧美激情久久久久久爽电影| 一二三四社区在线视频社区8| 亚洲欧美精品综合久久99| 熟女电影av网| 午夜亚洲福利在线播放| 淫妇啪啪啪对白视频| 午夜免费成人在线视频| 一边摸一边抽搐一进一小说| 日韩欧美在线乱码| 国产精品av视频在线免费观看| 国产成人aa在线观看| 黄色成人免费大全| 成年女人永久免费观看视频| 成人亚洲精品av一区二区| 精品一区二区三区视频在线观看免费| 国产精品av久久久久免费| 国产99白浆流出| 久久人人精品亚洲av| 俺也久久电影网| 亚洲最大成人中文| 亚洲精品粉嫩美女一区| 99热这里只有精品一区 | 老汉色av国产亚洲站长工具| 亚洲欧美日韩高清专用| 日韩欧美国产在线观看| 不卡一级毛片| 一本久久中文字幕| 亚洲欧美精品综合一区二区三区| 制服丝袜大香蕉在线| 中文资源天堂在线| 中文字幕人成人乱码亚洲影| 黄频高清免费视频| 给我免费播放毛片高清在线观看| 欧美黄色淫秽网站| 国产激情久久老熟女| 精品久久久久久久久久免费视频| 国产高清videossex| 久久久久久九九精品二区国产| 狠狠狠狠99中文字幕| 日日摸夜夜添夜夜添小说| 国产伦精品一区二区三区视频9 | 叶爱在线成人免费视频播放| 精品久久久久久久久久免费视频| 男人舔女人下体高潮全视频| 国内精品久久久久久久电影| 12—13女人毛片做爰片一| 丰满的人妻完整版| 亚洲欧美精品综合久久99| 久久国产乱子伦精品免费另类| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| 欧美黑人巨大hd| 国产熟女xx| 国产一区二区激情短视频| 最近最新中文字幕大全电影3| 啦啦啦免费观看视频1| e午夜精品久久久久久久| 好男人电影高清在线观看| 最近在线观看免费完整版| 亚洲国产高清在线一区二区三| 精品久久久久久久末码| 欧美不卡视频在线免费观看| 91av网站免费观看| 国产亚洲精品一区二区www| 嫩草影院精品99| 午夜日韩欧美国产| 久久久国产精品麻豆| 国产1区2区3区精品| 日韩三级视频一区二区三区| 欧美一级a爱片免费观看看| 精品乱码久久久久久99久播| www.熟女人妻精品国产| 少妇丰满av| 特大巨黑吊av在线直播| 亚洲国产精品999在线| 久久久国产成人精品二区| 亚洲成人精品中文字幕电影| xxxwww97欧美| 日本五十路高清| 夜夜夜夜夜久久久久| 久久精品aⅴ一区二区三区四区| 国内精品久久久久精免费| 99国产综合亚洲精品| 黄色成人免费大全| 国产精品美女特级片免费视频播放器 | 麻豆av在线久日| 免费在线观看影片大全网站| 欧美中文日本在线观看视频| 黄片小视频在线播放| 欧美高清成人免费视频www| 亚洲aⅴ乱码一区二区在线播放| 国产一区二区三区视频了| 色视频www国产| 久久久成人免费电影| 国产精品野战在线观看| 国产成人精品无人区| 欧美日韩黄片免| 夜夜夜夜夜久久久久| 精品国产美女av久久久久小说| 真人做人爱边吃奶动态| 99久久无色码亚洲精品果冻| ponron亚洲| 亚洲欧美精品综合久久99| 中文字幕人妻丝袜一区二区| 国产免费男女视频| 久久性视频一级片| 亚洲美女黄片视频| 亚洲人成网站高清观看| 观看美女的网站| 国产精品1区2区在线观看.| 熟女人妻精品中文字幕| 国产不卡一卡二| 怎么达到女性高潮| а√天堂www在线а√下载| 999久久久精品免费观看国产| 成人鲁丝片一二三区免费| 黄片小视频在线播放| 欧美乱码精品一区二区三区| 在线观看午夜福利视频| 亚洲真实伦在线观看| 中文字幕最新亚洲高清| 极品教师在线免费播放| 午夜两性在线视频| 亚洲一区高清亚洲精品| 亚洲精品色激情综合| 特级一级黄色大片| 亚洲欧美日韩高清专用| 美女午夜性视频免费| 啦啦啦韩国在线观看视频| av视频在线观看入口| 国产高清视频在线播放一区| 精品久久久久久成人av| 视频区欧美日本亚洲| 久久精品亚洲精品国产色婷小说| 性色av乱码一区二区三区2| 中出人妻视频一区二区| 少妇熟女aⅴ在线视频| 日本 欧美在线| av在线蜜桃| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 成人鲁丝片一二三区免费| 国产免费av片在线观看野外av| 天堂网av新在线| 亚洲av电影在线进入| 在线十欧美十亚洲十日本专区| 久久精品91无色码中文字幕| 悠悠久久av| 欧美中文日本在线观看视频| 在线永久观看黄色视频| 亚洲欧美日韩高清在线视频| 亚洲精品粉嫩美女一区| 久久久久免费精品人妻一区二区| 最近最新中文字幕大全免费视频| 国产亚洲精品久久久久久毛片| 99久久国产精品久久久| 搡老岳熟女国产| 欧美国产日韩亚洲一区| www.999成人在线观看| 亚洲中文日韩欧美视频| 制服丝袜大香蕉在线| 一个人看视频在线观看www免费 | 亚洲九九香蕉| 少妇熟女aⅴ在线视频| 男女视频在线观看网站免费| 日本成人三级电影网站| 精品欧美国产一区二区三| 亚洲av成人不卡在线观看播放网| 欧美黄色淫秽网站| 久久精品夜夜夜夜夜久久蜜豆| 观看免费一级毛片| 久久香蕉国产精品| 一本综合久久免费| 国产成人精品无人区| 九九在线视频观看精品| 毛片女人毛片| 国产精品久久久久久亚洲av鲁大| 国产v大片淫在线免费观看| 丰满的人妻完整版| 精品国内亚洲2022精品成人| 天天添夜夜摸| 欧美日韩精品网址| 好男人在线观看高清免费视频| 人人妻,人人澡人人爽秒播| 欧美另类亚洲清纯唯美| 岛国视频午夜一区免费看| 少妇的逼水好多| av视频在线观看入口| 成年女人看的毛片在线观看| 在线看三级毛片| 亚洲片人在线观看| 可以在线观看毛片的网站| 国产精品av久久久久免费| 岛国在线观看网站| 久99久视频精品免费| 国产精品久久久久久精品电影| 欧美绝顶高潮抽搐喷水| 18禁黄网站禁片午夜丰满| 美女黄网站色视频| 手机成人av网站| xxxwww97欧美| 日韩三级视频一区二区三区| 国产视频一区二区在线看| 亚洲第一欧美日韩一区二区三区| 欧美日韩中文字幕国产精品一区二区三区|