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

    沱江浮游植物群落特征及水質(zhì)評(píng)價(jià)*

    2016-01-15 08:21:25謝碧文齊澤民楊品紅
    海洋與湖沼 2016年4期
    關(guān)鍵詞:沱江豐水期浮游

    陶 敏 謝碧文 齊澤民 楊品紅 李 斌①

    (1. 內(nèi)江師范學(xué)院生命科學(xué)學(xué)院 長(zhǎng)江上游魚(yú)類資源保護(hù)與利用四川省重點(diǎn)實(shí)驗(yàn)室 內(nèi)江 641112; 2. 湖南文理學(xué)院 水產(chǎn)高效健康生產(chǎn)湖南省協(xié)同創(chuàng)新中心 環(huán)洞庭湖水產(chǎn)健康養(yǎng)殖及加工湖南省重點(diǎn)實(shí)驗(yàn)室 動(dòng)物學(xué)湖南省高校重點(diǎn)實(shí)驗(yàn)室 常德415000)

    水質(zhì)是水生態(tài)系統(tǒng)健康的重要指標(biāo), 進(jìn)行水質(zhì)評(píng)價(jià)是水體開(kāi)發(fā)利用、保護(hù)與生態(tài)恢復(fù)的前提與基礎(chǔ)(焦珂?zhèn)サ? 2015)。早期的水質(zhì)評(píng)價(jià)多針對(duì)理化指標(biāo)(Bozzetti et al, 2004)。然而由于環(huán)境與生物之間的復(fù)雜相互作用, 僅僅用理化指標(biāo)無(wú)法完整體現(xiàn)水體的健康狀況。作為水域的初級(jí)生產(chǎn)者, 浮游植物是整個(gè)水生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動(dòng)的基礎(chǔ)(王愛(ài)愛(ài)等,2014), 在維持水生態(tài)系統(tǒng)平衡中起著重要作用。浮游植物群落結(jié)構(gòu)對(duì)水環(huán)境條件變化敏感, 其種類組成和分布格局對(duì)環(huán)境變化具有重要指示作用, 已被廣泛應(yīng)用于水質(zhì)評(píng)價(jià)中(況琪軍等, 2005; 王愛(ài)愛(ài)等,2014; 茍婷等, 2015; 楊敏等, 2015)。

    沱江為長(zhǎng)江上游重要支流, 流域內(nèi)人口密度高,工農(nóng)業(yè)發(fā)達(dá), 污染嚴(yán)重, 富營(yíng)養(yǎng)化程度在四川省主要河流中最高(劉佳等, 2012)。2003年的特大水污染事件(孫莉等, 2005)引起了學(xué)者們對(duì)沱江水質(zhì)的密切關(guān)注, 從各種理化指標(biāo)的角度對(duì)沱江水質(zhì)進(jìn)行評(píng)價(jià)(謝賢健等, 2009; 肖玖金等, 2011; 周實(shí)際等, 2011; 劉建平, 2013; 陳雨艷等, 2015; 周開(kāi)錫等, 2015), 而生物學(xué)監(jiān)測(cè)方面的研究仍停留在20世紀(jì)80、90年代(李進(jìn), 1986; 楊昌述, 1987)。本研究對(duì)沱江中下游10個(gè)斷面開(kāi)展采樣調(diào)查, 從理化指標(biāo)和浮游植物兩個(gè)角度出發(fā)進(jìn)行水質(zhì)評(píng)價(jià), 為沱江水體健康狀況提供重要的基礎(chǔ)數(shù)據(jù), 也為水體的利用及河流生態(tài)修復(fù)提供參考資料。

    1 材料與方法

    1.1 樣品采集與分析

    1.1.1 采樣點(diǎn)布設(shè) 根據(jù)研究江段周邊環(huán)境特點(diǎn),在城市區(qū)(蓮花山、資陽(yáng)一橋、西林渡口、沱橋)、工業(yè)區(qū)或建筑工地附近(麻柳壩工業(yè)園、資州大橋)、城郊上游或鄉(xiāng)鎮(zhèn)(順河場(chǎng)鎮(zhèn)、五里店水電站、銀山鎮(zhèn)、二水廠)布設(shè)采樣斷面 10個(gè)(圖 1)。各采樣點(diǎn)環(huán)境描述參考陶敏等(2016)。

    圖1 采樣點(diǎn)分布圖Fig.1 Map of sampling sites

    1.1.2 樣品采集 于2013年8月(豐水期)與11月(枯水期)在研究江段進(jìn)行采樣和調(diào)查。表層水樣以辛德勒采水器(Schindler-Patalas trap)采集。水化學(xué)及葉綠素樣品直接裝入干凈的樣品瓶密封, 冷藏帶回實(shí)驗(yàn)室。同時(shí)將1L浮游植物定量樣品裝入樣品瓶并立即加入魯哥氏試劑(Lugol’s reagent) 現(xiàn)場(chǎng)固定。浮游植物定性樣品用 25號(hào)浮游生物網(wǎng)在水面下作“∞”字型撈取并過(guò)濾收集, 加入4%甲醛現(xiàn)場(chǎng)固定。

    1.1.3 浮游植物樣品的處理與分析 將浮游植物定量樣品在浮游生物沉淀器中避光靜置沉降 48h, 以虹吸法去除上層清液后, 濃縮定容至 50mL, 加入數(shù)滴甲醛溶液保存。計(jì)數(shù)方法參照章宗涉等(1991), 種類鑒定參考胡鴻鈞等(2006)。

    1.1.4 理化指標(biāo)的分析測(cè)定 使用塞氏盤(pán)(Secchi disc)、便攜式溶解氧測(cè)量?jī)x(YSI-550A)、電導(dǎo)率儀(Mettler Toledo SevenGo SG3)等設(shè)備對(duì)研究區(qū)域透明度(Secchi Depth, SD)、水溫、pH、溶解氧(DO)和電導(dǎo)率(Conductivity)等進(jìn)行現(xiàn)場(chǎng)測(cè)定。總氮(TN)、氨氮(NH4-N)、硝酸鹽氮(NO3-N)、總磷(TP)、正磷酸鹽(PO4-P)、高錳酸鹽指數(shù)(CODMn)、葉綠素a(Chl a)含量在實(shí)驗(yàn)室中按照《水和廢水監(jiān)測(cè)分析方法》(國(guó)家環(huán)境保護(hù)總局, 2002)進(jìn)行分析。

    1.2 評(píng)價(jià)方法

    1.2.1 優(yōu)勢(shì)種與優(yōu)勢(shì)度 浮游植物優(yōu)勢(shì)種依據(jù)各種的優(yōu)勢(shì)度指數(shù)(Y)值進(jìn)行鑒別(劉成士等, 2015):

    式中, Pi為第i種浮游植物豐度與樣品中浮游植物總豐度的比值(ni/N), fi為第i種浮游植物在各采樣點(diǎn)出現(xiàn)的頻度。Y≥0.02的物種為優(yōu)勢(shì)種。

    1.2.2 綜合營(yíng)養(yǎng)狀態(tài)指數(shù) 以總氮(TN)、總磷(TP)、高錳酸鹽指數(shù)(CODMn)、葉綠素a(Chl a)和透明度(SD)為評(píng)價(jià)指標(biāo), 采用湖泊(水庫(kù))富營(yíng)養(yǎng)化評(píng)價(jià)方法及分級(jí)技術(shù)規(guī)定(中國(guó)環(huán)境監(jiān)測(cè)總站, 2001)中的綜合營(yíng)養(yǎng)狀態(tài)指數(shù) TLI(∑)評(píng)價(jià)方法與標(biāo)準(zhǔn)對(duì)沱江營(yíng)養(yǎng)狀態(tài)進(jìn)行分析。

    式中, Wj為第 j種參數(shù)的營(yíng)養(yǎng)狀態(tài)指數(shù)的相關(guān)權(quán)重;TLI(j)為第 j種參數(shù)的營(yíng)養(yǎng)狀態(tài)指數(shù)。TLI(∑)≤30時(shí)為貧營(yíng)養(yǎng)狀態(tài); 30<TLI(∑)≤50為中營(yíng)養(yǎng)狀態(tài);TLI(∑)>50 則為富營(yíng)養(yǎng)化狀態(tài), 其中 50<TLI(∑)≤60為輕度富營(yíng)養(yǎng)化, 60<TLI(∑)≤70為中度富營(yíng)養(yǎng)化,TLI(∑)>70 為重度富營(yíng)養(yǎng)化。

    1.2.3 多樣性指數(shù) 同時(shí)運(yùn)用 Shannon-Wiener多樣性指數(shù)(H′)、Margalef豐富度指數(shù)(d)、Pielou均勻度指數(shù)(J)進(jìn)行評(píng)價(jià)。計(jì)算公式分別為

    式中, S為樣品中物種種類總數(shù), n為樣品中第i種浮游植物的豐度, N為樣品中浮游植物總豐度。H’>3時(shí)為輕或無(wú)污染, 1—3為中污染, 其中2—3為β-中污,1—2為 α-中污, 在 0—1之間為重污染。d>5時(shí)為清潔, 在4—5之間為寡污型, 3—4之間為β-中污, 1—3為α-中污, 0—1為重污染。J值在0.8—1.0之間時(shí)為清潔, 0.5—0.8為寡污型, 0.4—0.5為β-中污, 0.3—0.4為α-中污, 0—0.3為重污染(況琪軍等, 2005; 周廣杰等, 2006; 郭春燕等, 2010; Spatharis et al, 2011; 劉成士等, 2015; 楊敏等, 2015; 茍婷等, 2015)。

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

    本文統(tǒng)計(jì)分析采用 SPSS 19.0與 Salford Predictive Modeler (SPM) 7.0軟件實(shí)現(xiàn); 基本圖形繪制在Origin 9.0軟件中完成。

    1.3.1 分類回歸樹(shù) 為分析營(yíng)養(yǎng)狀態(tài)與生物多樣性與浮游植物密度之間的關(guān)系, 以H′、J、d、TLI(∑)以及采樣點(diǎn)類型(城市、農(nóng)村、工業(yè)區(qū))為預(yù)測(cè)指標(biāo), 采用分類回歸樹(shù)(classification and regression tree, CART)評(píng)價(jià)這5個(gè)預(yù)測(cè)指標(biāo)的對(duì)浮游植物密度的重要性。對(duì)浮游植物密度進(jìn)行常用對(duì)數(shù)(lg)轉(zhuǎn)換, 采用R2評(píng)價(jià)CART預(yù)測(cè)效果。CART分析采用Salford Predictive Modeler (SPM) 7.0。

    1.3.2 聚類分析 為探尋浮游植物的空間分布差異, 根據(jù)浮游植物種類和密度, 利用SPSS 19.0軟件對(duì)不同采樣斷面的浮游植物進(jìn)行聚類分析。浮游植物密度經(jīng)過(guò)平方根轉(zhuǎn)換。

    2 結(jié)果與討論

    2.1 沱江環(huán)境參數(shù)與營(yíng)養(yǎng)狀態(tài)

    2.1.1 主要環(huán)境參數(shù) 從季節(jié)來(lái)看, 沱江水體中的浮游植物密度、葉綠素a、TN與TP數(shù)值總體呈現(xiàn)出豐水期高于枯水期的趨勢(shì); 與之相反, 透明度則表現(xiàn)為枯水期較高(圖2、圖3)。這一現(xiàn)象可能與豐水期較大的地表徑流量帶來(lái)的豐富氮磷(夏品華等, 2011)及該時(shí)期適宜的水溫、光照有關(guān)。其余參數(shù)無(wú)明顯的季節(jié)差異。

    圖2 葉綠素a、透明度、高錳酸鹽指數(shù)、總氮和總磷的空間分布Fig.2 Spatial distribution of chlorophyll a, Secchi depth (SD), CODMn, TN and TP

    圖3 浮游植物種類數(shù)與密度分布Fig.3 Spatial distribution of phytoplankton species and density

    從空間分布來(lái)看, 豐水期研究區(qū)域上游河段(資陽(yáng)市區(qū)) 的葉綠素 a、浮游植物密度與 TP等參數(shù)值均高于中、下游河段(資中、內(nèi)江市區(qū)), 見(jiàn)圖2、圖3。這可能與沱江上游磷礦開(kāi)采有關(guān)(施澤明等, 2012),在豐水期高氮含量的情況下, 浮游植物的生長(zhǎng)主要受到磷的限制(田時(shí)彌等, 2015), 故而上游水域因地表徑流帶來(lái)較多的磷, 浮游植物與葉綠素含量均較高??菟竟?jié)葉綠素a、浮游植物密度與TN在中、下河段較高(圖2、圖3), 說(shuō)明枯水期由于地表徑流沖刷力度降低, 帶入沱江中的氮素逐減少并成為藻類生長(zhǎng)限制因子。在兩個(gè)季節(jié)中, 透明度變化與浮游植物密度變化規(guī)律不相符(圖2、圖3), 可能與沱江水體中泥沙含量有關(guān)。

    2.1.2 綜合營(yíng)養(yǎng)狀態(tài)指數(shù) 根據(jù)綜合營(yíng)養(yǎng)狀態(tài)指數(shù) TLI(∑)的大小對(duì)研究江段水體進(jìn)行富營(yíng)養(yǎng)化狀態(tài)評(píng)價(jià)(圖4a), 沱江水質(zhì)總體處于中-富營(yíng)養(yǎng)狀態(tài)。與四川省其它主要河流如岷江、嘉陵江等相比, 沱江的TLI(∑)較高, 富營(yíng)養(yǎng)化較為嚴(yán)重(劉佳等, 2012)??菟竟?jié)水體整體為中營(yíng)養(yǎng)狀態(tài)(<50), 水質(zhì)良好; 豐水期多存在各種程度的污染, 尤其是位于資陽(yáng)市區(qū)的 1—3號(hào)采樣點(diǎn), 富營(yíng)養(yǎng)化程度較高(>60)。豐水期較大的降水沖刷使得更多的氮磷營(yíng)養(yǎng)進(jìn)入水體, 尤其是較接近上游磷礦的江段, 富營(yíng)養(yǎng)化水平較高, 而中、下游水體的自凈作用在一定程度上緩解了富營(yíng)養(yǎng)化進(jìn)程。

    圖4 各種參數(shù)的空間分布情況Fig.4 Spatial distribution of various indicators

    2.2 浮游植物群落結(jié)構(gòu)分析與水質(zhì)評(píng)價(jià)

    2.2.1 浮游植物群落組成與優(yōu)勢(shì)種 調(diào)查期間共檢出浮游植物7門(mén)88種(屬), 其中綠藻種類數(shù)最多為37種(屬), 硅藻23種(屬), 藍(lán)藻15種(屬), 其它門(mén)類13種(屬)。浮游植物種類數(shù)在空間分布上變化不大,但在時(shí)間分布上有較大差異(圖 3), 該差異主要由藍(lán)藻和綠藻引起。豐水期各斷面種類數(shù)變化范圍為30—52種, 以6#最低, 4#最高; 枯水期為12—28種,3#最低, 9#最高。根據(jù)圖3, 豐水期浮游植物密度主要由藍(lán)藻(35%)、硅藻(35%)和綠藻(24%)組成; 枯水期硅藻成為優(yōu)勢(shì)類群(44%), 隱藻數(shù)量增多(18%), 藍(lán)藻和綠藻水平相當(dāng)(15%)。浮游植物群落結(jié)構(gòu)是評(píng)價(jià)水體富營(yíng)養(yǎng)化的重要參數(shù)。一般認(rèn)為甲藻、隱藻和硅藻易在中營(yíng)養(yǎng)型水體中占優(yōu)勢(shì), 綠藻和藍(lán)藻則易在富營(yíng)養(yǎng)型水體中占優(yōu)勢(shì)(金相燦等, 1990)。沱江豐、枯水季節(jié)浮游植物群落結(jié)構(gòu)分別為綠藻+硅藻+藍(lán)藻和硅藻+綠藻+隱藻, 由此推斷該江段目前應(yīng)處于中到富營(yíng)養(yǎng)狀態(tài), 尤其是豐水期, 富營(yíng)養(yǎng)化程度更高, 這與前文營(yíng)養(yǎng)狀態(tài)分析結(jié)果一致。

    研究江段浮游植物優(yōu)勢(shì)種較多且優(yōu)勢(shì)度均不高(表 1), 表明其群落結(jié)構(gòu)較復(fù)雜也較為完整。硅藻在兩個(gè)季節(jié)均有優(yōu)勢(shì), 其中梅尼小環(huán)藻(Cyclotellameneghiniana)和花環(huán)小環(huán)藻(C. operculata)存在豐水期高峰, 脆桿藻(Fragilariasp.)、線形舟形藻(Navicula graciloides)和針桿藻(Synedrasp.)則在枯水期占優(yōu)勢(shì)。藍(lán)藻和隱藻有明顯的季節(jié)差異: 藍(lán)藻為豐水期優(yōu)勢(shì)類群, 優(yōu)勢(shì)種為湖泊假魚(yú)腥藻(Pseudanabaena limnetica)、拉氏擬柱孢藻(Cylindrospermopsis raciborskii)、微小平裂藻(Merismopedia tenuissima)和細(xì)小平裂藻(M. minima); 隱藻則主要在枯水季節(jié)占優(yōu)勢(shì), 包括尖尾藍(lán)隱藻(Chroomonas acuta)、卵形隱藻(Cryptomonas ovata)和嚙蝕隱藻(C. erosa)。其它類群如綠藻和金藻優(yōu)勢(shì)種類少且優(yōu)勢(shì)度均較低。豐水期優(yōu)勢(shì)度最高的梅尼小環(huán)藻和微小平裂藻均為 α-中污指示種(沈蘊(yùn)芬等, 1990), 且富營(yíng)養(yǎng)化水體中常見(jiàn)的絲狀藍(lán)藻湖泊假魚(yú)腥藻(黃成等, 2011)等出現(xiàn)頻度亦較高, 說(shuō)明該季節(jié)沱江已處于富營(yíng)養(yǎng)化狀態(tài)??菟竟?jié)優(yōu)勢(shì)種則反映出寡污-中污狀態(tài)(沈蘊(yùn)芬等, 1990)。

    表1 沱江浮游植物優(yōu)勢(shì)種及優(yōu)勢(shì)度Tab.1 Dominant species and dominance degree of phytoplankton in Tuojiang River

    2.2.2 浮游植物群落聚類分析 聚類分析(圖 5)結(jié)果顯示, 研究江段中游3個(gè)位于五里店水電站壩下的斷面(5#、6#、7#)聚集為集合 1; 下游天宮堂水電站壩上 3個(gè)斷面(8#、9#、10#)與中上游的 4#斷面共同組成集合2; 上游南津驛水電站壩上3個(gè)斷面(1#、2#、3#)則聚為集合3。3個(gè)集合的藍(lán)藻、綠藻、硅藻、隱藻比例分別為15.6%、19.4%、50%、8.7%和33.7%、20.3%、30.5%、8.7%以及21.5%、17.7%、41.4、16.9%,硅藻在壩下集合 1中所占比例最大, 藍(lán)藻+綠藻在壩上集合2中最多, 隱藻則在壩上集合3中最多。由于水電壩的攔截作用, 大壩上下游之間生態(tài)條件會(huì)發(fā)生較大變化, 壩上湖泊區(qū)由于流速和透明度的變化其浮游植物群落結(jié)構(gòu)也顯著區(qū)別與壩下區(qū)域(韓博平,2010; 楊敏等, 2015), 因此與壩下集合 1相比, 壩上集合2與集合3中湖泊型藍(lán)藻、綠藻類及隱藻增多,河流型硅藻類有所下降(邱光勝等, 2011), 這一變化趨勢(shì)與嘉陵江等鄰近水體接近(楊敏等, 2015)。

    單因素方差分析(One-way ANOVA)顯示, 浮游植物密度與NH+4數(shù)值均表現(xiàn)出從集合1到集合3明顯升高的趨勢(shì)(P<0.05, 表 2), 其余指標(biāo)差異不顯著(P>0.05, 表 2), 說(shuō)明水電壩的建設(shè)與空間分布格局對(duì)水體多種營(yíng)養(yǎng)鹽濃度作用較小, 但對(duì)與浮游植物群落結(jié)構(gòu)已有影響。由于水電壩的攔截作用, 大壩上游附近區(qū)域流速減緩, NH+4等營(yíng)養(yǎng)積累(趙晨等, 2014), 適宜浮游植物增殖, 故壩上集合2與3中浮游植物密度與NH+4均高于壩下集合1; 同為壩上采樣點(diǎn), 較上游的集合3中NH+4濃度最高可能與沱江上游的企業(yè)污水排放有關(guān)(孫莉等, 2005),喜好NH+4的浮游植物(de Aldaet al, 1996)密度也隨之更高。

    圖5 采樣點(diǎn)聚類分析樹(shù)狀圖Fig.5 Cluster analysis for the sampling sites

    表2 沱江不同采樣斷面集群水體特征比較Tab.2 Comparison of biological and physiochemical parameters among sample sites in 3 groups

    2.2.3 多樣性指數(shù) 沱江浮游植物多樣性指數(shù)變化趨勢(shì)如圖 4b—d。豐水期中各指數(shù)變化較為一致,物種多樣性指數(shù)H′平均3.66, 大致處于清潔狀態(tài), 最高值出現(xiàn)在 7#斷面(4.14), 最低為 6#斷面(1.98); 物種豐富度指數(shù)d平均6.21, 整體為清潔, 1#斷面最高(7.47), 6#斷面最低(4.39); 均勻度指數(shù)J平均0.69, 多為清潔-寡污狀態(tài), 與多樣性指數(shù)H′一樣, 最高與最低值分別出現(xiàn)在6#(0.84)和7#(0.40)。6#與10#等斷面物種多樣性與均勻度指數(shù)較低呈中度污染狀態(tài), 是由于細(xì)胞很小的群體藍(lán)藻, 如平裂藻屬(Merismopediaspp.)的優(yōu)勢(shì)造成。盡管該類藻細(xì)胞數(shù)量較大, 但對(duì)浮游植物生物量和水質(zhì)的影響卻較小, 形成水華的可能性也較低, 根據(jù)此結(jié)果豐水期沱江處于清潔-寡污狀態(tài)。枯水期各指數(shù)一致性較差: 多樣性指數(shù)H′(3.13—4.21)與均勻度指數(shù)J(0.80—0.93)各點(diǎn)差異不明顯; 物種豐富度指數(shù)d最低值出現(xiàn)在 3#(3.12), 最高值為 8#斷面(5.40), 多數(shù)點(diǎn)處于輕度到中度污染狀態(tài)。三個(gè)指數(shù)中, 物種豐富度指數(shù)d(圖4)與浮游植物種類數(shù)(圖 3)變化趨勢(shì)一致性較好, 說(shuō)明該指數(shù)對(duì)沱江中浮游植物群落物種數(shù)變化較為敏感, 因此按物種豐富度指數(shù)d評(píng)價(jià)標(biāo)準(zhǔn), 枯水季節(jié)沱江各斷面多處于輕度到中度污染狀態(tài)。但是應(yīng)當(dāng)看到, 盡管從多樣性來(lái)說(shuō)沱江污染程度較低, 其實(shí)營(yíng)養(yǎng)水平和浮游植物群落結(jié)構(gòu)卻呈現(xiàn)出明顯的富營(yíng)養(yǎng)化狀態(tài)。由于環(huán)境的復(fù)雜性, 從不同角度評(píng)價(jià)水體狀態(tài)均有其局限性(易斌等, 2007), 綜合各種指標(biāo)進(jìn)行判定, 沱江污染程度仍不容樂(lè)觀。

    2.3 浮游植物密度與營(yíng)養(yǎng)和物種多樣性的關(guān)系

    浮游植物群落結(jié)構(gòu)受外界環(huán)境(營(yíng)養(yǎng)水平)與內(nèi)部因素(種內(nèi)與種間關(guān)系, 即群落結(jié)構(gòu))的共同作用,從數(shù)據(jù)上來(lái)說(shuō), 浮游植物密度與營(yíng)養(yǎng)狀態(tài)指數(shù)以及多樣性指數(shù)有某種關(guān)聯(lián)。以浮游植物密度為因變量,綜合營(yíng)養(yǎng)狀態(tài)指數(shù)和多樣性指數(shù)為自變量所作的分類回歸樹(shù)模型R2=0.93, 說(shuō)明這兩類指數(shù)能較好地解釋浮游植物密度的變化。該模型將浮游植物密度數(shù)據(jù)作了3次細(xì)分(圖6)。第一個(gè)分支將4個(gè)物種豐富度較高(d>6.56)并且浮游植物密度較高(未進(jìn)行 lg轉(zhuǎn)換的平均值為 4.6×106cell/L)的樣品與其它樣品分離開(kāi)來(lái); 第二個(gè)分支又將 5個(gè)具有較低營(yíng)養(yǎng)水平(TLI≤43.9)且浮游植物密度較低(平均值為0.4 ×106cell/L)的樣品從數(shù)據(jù)集中分離出來(lái); 第三個(gè)分支將剩余的 11個(gè)樣品以J=0.85為標(biāo)準(zhǔn)分為兩部分, 均勻度低(≤0.85)的 8個(gè)樣品浮游植物密度相對(duì)較高(平均值為2.6×106cell/L), 均勻度較高(>0.85)的 3 個(gè)樣品浮游植物密度較低(平均 1.0 ×106cell/L)。Shannon-Wiener多樣性指數(shù)(H′)與采樣點(diǎn)位置類型未被納入模型中。該模型顯示, 營(yíng)養(yǎng)水平與物種豐富度較高, 但均勻度偏低的水體中, 浮游植物密度易出現(xiàn)高值。在營(yíng)養(yǎng)物質(zhì)富足的情況下, 浮游植物種類較多, 但有一部分為偶見(jiàn)種(俞秋佳, 2014), 少數(shù)優(yōu)勢(shì)種的大量繁殖降低了群落均勻度, 這些優(yōu)勢(shì)種若繼續(xù)發(fā)展, 水體就有發(fā)生水華的潛在風(fēng)險(xiǎn)(馮波等, 2008)。對(duì)漢江硅藻水華的研究表明, 較豐富的物種多樣性有可能是水華初期的特點(diǎn)(潘曉潔等, 2014)。因此在營(yíng)養(yǎng)水平較高的水體中,為預(yù)警水華的發(fā)生, 豐富的浮游植物種類不僅不能成為放松警惕的原因, 更應(yīng)成為密切監(jiān)測(cè)的理由。

    圖6 浮游植物密度與綜合營(yíng)養(yǎng)狀態(tài)指數(shù)、多樣性指數(shù)間的關(guān)系Fig.6 Classification and regression tree of phytoplankton abundance versus trophic level index and diversity index x為未進(jìn)行l(wèi)g轉(zhuǎn)換的浮游植物密度平均值; n為樣本數(shù)量

    3 結(jié)論

    根據(jù)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)判斷沱江水體處于中-富營(yíng)養(yǎng)狀態(tài), 富營(yíng)養(yǎng)化程度較高。由于雨水淋溶沖刷作用和上游磷礦的影響, 豐水期污染程度高于枯水期,上游斷面高于較下游斷面。

    從浮游植物群落組成來(lái)看, 沱江為中-富營(yíng)養(yǎng)型水體, 其中豐水期為綠藻+硅藻+藍(lán)藻的形式, 枯水期為硅藻+綠藻+隱藻的形式, 仍表現(xiàn)出豐水期富營(yíng)養(yǎng)化程度高于枯水期的狀態(tài)。優(yōu)勢(shì)種類的組成顯示, 沱江中浮游植物群落結(jié)構(gòu)較為復(fù)雜和完整, 豐水期為α-中污染, 枯水期則為寡污-中污狀態(tài)。由于沿江梯級(jí)大壩的影響, 沱江水文條件與 NH4+分布具有明顯空間異質(zhì)性, 導(dǎo)致了大壩上下游之間、上游斷面與下游斷面之間浮游植物群落結(jié)構(gòu)與密度的差異, 整體表現(xiàn)出大壩上游湖泊型藻類數(shù)量增加、上游斷面浮游植物密度高于下游斷面的趨勢(shì)。盡管多樣性指數(shù)顯示沱江污染狀態(tài)較輕, 綜合考慮各類指標(biāo)的局限性, 沱江污染程度仍不容樂(lè)觀。

    沱江中浮游植物密度受到營(yíng)養(yǎng)和群落結(jié)構(gòu)的顯著影響: 在營(yíng)養(yǎng)水平高, 物種豐富度較大但均勻度較低的情況下, 容易暴發(fā)水華。因此對(duì)于營(yíng)養(yǎng)物質(zhì)含量豐富并且浮游植物種類數(shù)較多的水體, 不能放松警惕, 應(yīng)加強(qiáng)監(jiān)控。

    王愛(ài)愛(ài), 馮 佳, 謝樹(shù)蓮, 2014. 汾河中下游浮游藻類群落特征及水質(zhì)分析. 環(huán)境科學(xué), 5(3): 915—923

    中國(guó)環(huán)境監(jiān)測(cè)總站, 2001. 湖泊(水庫(kù))富營(yíng)養(yǎng)化評(píng)價(jià)方法及分級(jí)技術(shù)規(guī)定. 北京: 中國(guó)環(huán)境科學(xué)出版社

    田時(shí)彌, 楊 揚(yáng), 喬永民等, 2015. 珠江流域東江干流浮游植物葉綠素 a時(shí)空分布及與環(huán)境因子的關(guān)系. 湖泊科學(xué),27(1): 31—37

    馮 波, 陳俊輝, 盧伙勝, 2008. 茂名市大放雞島海域浮游植物群落的季節(jié)變化. 海洋漁業(yè), 30(1): 67—73

    劉 佳, 徐 亮, 張秋勁等, 2012. 四川省水體富營(yíng)養(yǎng)狀況時(shí)空分布調(diào)查. 中國(guó)環(huán)境監(jiān)測(cè), 28(5): 6—8

    劉成士, 黃 燕, 劉 瑜等, 2015. 青山水庫(kù)輪蟲(chóng)群落結(jié)構(gòu)特征的時(shí)空差異對(duì)水質(zhì)變化的影響. 海洋與湖沼, 46(3):694—702

    劉建平, 2013. 沱江富順段水質(zhì)分析及污染狀況研究. 四川環(huán)境, 32(2): 23—26

    孫 莉, 朱鴻斌, 張成云等, 2005. 一起工業(yè)廢水污染沱江水源水事故的調(diào)查. 環(huán)境與健康雜志, 22(3): 192—193

    李 進(jìn), 1986. 沱江水質(zhì)污染的微生物監(jiān)測(cè)評(píng)價(jià). 四川環(huán)境,5(3): 26—32

    楊 敏, 張 晟, 劉朔孺, 2015. 草街水庫(kù)蓄水后嘉陵江浮游植物群落特征及水質(zhì)評(píng)價(jià). 環(huán)境科學(xué), 36(7): 2480—2486

    楊昌述, 1987. 沱江干流底棲動(dòng)物調(diào)查及水質(zhì)生物學(xué)評(píng)價(jià). 環(huán)境科學(xué), 8(5): 50—56

    肖玖金, 譚周亮, 李旭東等, 2011. 基于物元分析法的沱江地表水水質(zhì)評(píng)價(jià). 安全與環(huán)境學(xué)報(bào), 11(1): 104—108

    邱光勝, 胡 圣, 葉 丹等, 2011. 三峽庫(kù)區(qū)支流富營(yíng)養(yǎng)化及水華現(xiàn)狀研究. 長(zhǎng)江流域資源與環(huán)境, 20(3): 311—316

    況琪軍, 馬沛明, 胡征宇等, 2005. 湖泊富營(yíng)養(yǎng)化的藻類生物學(xué)評(píng)價(jià)與治理研究進(jìn)展. 安全與環(huán)境學(xué)報(bào), 5(2): 87—91

    沈韞芬, 章宗涉, 龔循矩等, 1990. 微型生物監(jiān)測(cè)新技術(shù). 北京: 中國(guó)建筑工業(yè)出版社, 1—524

    陳雨艷, 余 恒, 向秋實(shí)等, 2015. 沱江流域水環(huán)境質(zhì)量分析.四川環(huán)境, 34(2): 85—89

    茍 婷, 馬千里, 許振成等, 2015. 貝江浮游藻類群落特征及富營(yíng)養(yǎng)化風(fēng)險(xiǎn)分析. 環(huán)境科學(xué), 36(3): 946—954

    國(guó)家環(huán)境保護(hù)總局, 2002. 水和廢水監(jiān)測(cè)分析方法(第四版).北京: 中國(guó)環(huán)境科學(xué)出版社

    易 斌, 吳玲玲, 林 端等, 2007. 營(yíng)養(yǎng)指數(shù)法和浮游植物多樣性指數(shù)法在評(píng)價(jià)富營(yíng)養(yǎng)化上的差異. 見(jiàn): 中國(guó)海洋學(xué)會(huì)赤潮研究與防治專業(yè)委員會(huì)第二屆學(xué)術(shù)研討會(huì). 上海: 中國(guó)海洋學(xué)會(huì)

    金相燦, 屠清瑛, 1990. 湖泊富營(yíng)養(yǎng)化調(diào)查規(guī)范(第二版). 北京: 中國(guó)環(huán)境科學(xué)出版社

    周廣杰, 況琪軍, 胡征宇等, 2006. 三峽庫(kù)區(qū)四條支流藻類多樣性評(píng)價(jià)及“水華”防治. 中國(guó)環(huán)境科學(xué), 26(3): 337—341

    周開(kāi)錫, 李 敏, 吳 建, 2015. 沱江河內(nèi)江段主要監(jiān)控?cái)嗝嫠|(zhì)變化趨勢(shì)分析. 見(jiàn): 2015年中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第一卷). 深圳: 中國(guó)環(huán)境科學(xué)學(xué)會(huì)

    周實(shí)際, 朱 俊, 張 科等, 2011. 沱江內(nèi)江城區(qū)段水質(zhì)分析.環(huán)境科學(xué)與管理, 36(10): 114—116

    趙 晨, 董世魁, 劉世梁等, 2014. 漫灣大壩上下游沉積物重金屬與營(yíng)養(yǎng)元素分布特征及環(huán)境風(fēng)險(xiǎn)評(píng)價(jià). 環(huán)境科學(xué)學(xué)報(bào), 34(9): 2417—2425

    胡鴻鈞, 魏印心, 2006. 中國(guó)淡水藻類—系統(tǒng)、分類及生態(tài). 北京: 科學(xué)出版社, 23—915

    俞秋佳, 2014. 蘇州河浮游植物群落變化特征及其影響因子研究. 上海: 華東師范大學(xué)碩士學(xué)位論文

    施澤明, 倪師軍, 張成江等, 2012. 沱江流域磷礦開(kāi)采和加工過(guò)程中放射性環(huán)境問(wèn)題探討. 地球科學(xué)進(jìn)展, 27(10):1134—1139

    夏品華, 李秋華, 林 陶等, 2011. 貴州高原百花湖水庫(kù)湖沼學(xué)變量特征及環(huán)境效應(yīng). 環(huán)境科學(xué)學(xué)報(bào), 31(8):1660—1669

    郭春燕, 馮 佳, 謝樹(shù)蓮, 2010. 山西晉陽(yáng)湖浮游藻類分布的時(shí)空格局及水質(zhì)分析. 湖泊科學(xué), 22(2): 251—255

    陶 敏, 王永明, 謝碧文等, 2016. 沱江浮游生物群落時(shí)空分布及相關(guān)環(huán)境因子分析. 水生生物學(xué)報(bào), 40(2): 301—312

    黃 成, 侯 偉, 顧繼光等, 2011. 珠江三角洲城市周邊典型中小型水庫(kù)富營(yíng)養(yǎng)化與藍(lán)藻種群動(dòng)態(tài). 應(yīng)用與環(huán)境生物學(xué)報(bào), 17(3): 295—302

    章宗涉, 黃祥飛, 1991. 淡水浮游生物研究方法. 北京: 科學(xué)出版社

    韓博平, 2010. 中國(guó)水庫(kù)生態(tài)學(xué)研究的回顧與展望. 湖泊科學(xué),22(2): 151—160

    焦珂?zhèn)? 周啟星, 2015. 基于水質(zhì)與生物指標(biāo)的松花江流域水生態(tài)健康評(píng)價(jià). 生態(tài)學(xué)雜志, 34(6): 1731—1737

    謝賢健, 蘭代萍, 2009. 基于因子分析法的沱江流域地表水水質(zhì)的綜合評(píng)價(jià). 安徽農(nóng)業(yè)科學(xué), 37(3): 1304—1306

    潘曉潔, 朱愛(ài)民, 鄭志偉等, 2014. 漢江中下游春季浮游植物群落結(jié)構(gòu)特征及其影響因素. 生態(tài)學(xué)雜志, 33(1): 33—40

    Bozzetti M, Schulz U H, 2004. An index of biotic integrity based on fish assemblages for subtropical streams in southern Brazil. Hydrobiologia, 529(1): 133—144

    de Alda J A G O, Tapia M I, Franck F et al, 1996. Changes in nitrogen source modify distribution of excitation energy in the cyanobacterium Phormidium laminosum. Physiologia Plantarum, 97(1): 69—78

    Spatharis S, Roelke D L, Dimitrakopoulos P G et al, 2011.Analyzing the (mis) behavior of Shannon index in eutrophication studies using field and simulated phytoplankton assemblages. Ecological Indicators, 11(2): 697—703

    猜你喜歡
    沱江豐水期浮游
    豐水期巢湖水體與沉積物重金屬污染及分配特征研究
    牛欄江流域浮游動(dòng)物群落結(jié)構(gòu)特征及驅(qū)動(dòng)因子
    膠東沿海某地區(qū)生活飲用水總α、總β放射性水平調(diào)查分析
    龍龍草的詩(shī)
    浞河浮游藻類的調(diào)查研究與水質(zhì)評(píng)價(jià)
    《浮游》
    流行色(2019年10期)2019-12-06 08:13:26
    恐龍科普秀
    ——多棘沱江龍
    星河映沱江 煙雨落鳳凰
    又唱沱江謠
    黃河之聲(2018年17期)2018-11-20 10:43:40
    浮游藻類在水體PAHs富集與降解中的研究
    国产成人精品婷婷| 成年美女黄网站色视频大全免费 | 成人无遮挡网站| 一本一本综合久久| av在线老鸭窝| h日本视频在线播放| 中文欧美无线码| 国产视频首页在线观看| 黄色视频在线播放观看不卡| 秋霞伦理黄片| 91久久精品电影网| 丰满乱子伦码专区| 国产视频内射| 日日摸夜夜添夜夜添av毛片| 国产 一区精品| 国产高清国产精品国产三级 | 国产精品成人在线| 一级二级三级毛片免费看| 日韩 亚洲 欧美在线| 国产亚洲av片在线观看秒播厂| 三级经典国产精品| 亚洲av男天堂| 日本一二三区视频观看| 日韩视频在线欧美| 国产精品一区二区性色av| 18禁裸乳无遮挡动漫免费视频| 18禁在线无遮挡免费观看视频| 欧美极品一区二区三区四区| 成人亚洲精品一区在线观看 | a级毛色黄片| 人妻 亚洲 视频| 寂寞人妻少妇视频99o| 99久久中文字幕三级久久日本| 卡戴珊不雅视频在线播放| 如何舔出高潮| 国产在线视频一区二区| 婷婷色综合大香蕉| 18禁在线无遮挡免费观看视频| 成年av动漫网址| 国产成人精品福利久久| 国产成人91sexporn| 毛片一级片免费看久久久久| 老女人水多毛片| 九九在线视频观看精品| 久久精品国产亚洲av天美| 一个人看视频在线观看www免费| 午夜福利在线在线| 久久ye,这里只有精品| 精品人妻熟女av久视频| 18+在线观看网站| 亚洲av中文av极速乱| 日韩强制内射视频| 国产亚洲精品久久久com| 观看免费一级毛片| av网站免费在线观看视频| 国产男人的电影天堂91| 少妇人妻精品综合一区二区| 欧美人与善性xxx| 国产精品偷伦视频观看了| 国产精品三级大全| 国产精品一二三区在线看| 久久亚洲国产成人精品v| 久久99蜜桃精品久久| 久久久精品94久久精品| 青春草国产在线视频| 国产精品不卡视频一区二区| 亚洲熟女精品中文字幕| 亚洲欧美日韩卡通动漫| 中文在线观看免费www的网站| 99精国产麻豆久久婷婷| 日韩大片免费观看网站| 久久精品夜色国产| 免费看光身美女| 日韩精品有码人妻一区| videossex国产| 97在线人人人人妻| 91久久精品电影网| 国产精品精品国产色婷婷| 欧美丝袜亚洲另类| 看免费成人av毛片| 少妇人妻久久综合中文| 国产精品久久久久成人av| 成人综合一区亚洲| 自拍偷自拍亚洲精品老妇| 天堂8中文在线网| 国产一级毛片在线| 国产av精品麻豆| 99久久中文字幕三级久久日本| 国内精品宾馆在线| 国产色婷婷99| 一区二区三区精品91| 深爱激情五月婷婷| 国产深夜福利视频在线观看| 在线看a的网站| 一级a做视频免费观看| 亚洲人与动物交配视频| www.av在线官网国产| 亚洲av成人精品一二三区| 好男人视频免费观看在线| 内地一区二区视频在线| 欧美日韩综合久久久久久| 各种免费的搞黄视频| 久久久色成人| 久久av网站| 精品99又大又爽又粗少妇毛片| 黄色欧美视频在线观看| 国产黄频视频在线观看| 久久午夜福利片| 国产亚洲av片在线观看秒播厂| 制服丝袜香蕉在线| 欧美一区二区亚洲| 欧美高清性xxxxhd video| 欧美区成人在线视频| 各种免费的搞黄视频| 男人和女人高潮做爰伦理| 免费人成在线观看视频色| 亚洲精品日韩在线中文字幕| 日韩成人av中文字幕在线观看| 熟妇人妻不卡中文字幕| 99视频精品全部免费 在线| 国产一区亚洲一区在线观看| 建设人人有责人人尽责人人享有的 | 国产又色又爽无遮挡免| av国产免费在线观看| 欧美成人a在线观看| av在线老鸭窝| 免费高清在线观看视频在线观看| 久久午夜福利片| 中文字幕久久专区| 五月伊人婷婷丁香| 在现免费观看毛片| 国产毛片在线视频| 亚洲一级一片aⅴ在线观看| 久久久久精品性色| 国产欧美亚洲国产| 欧美精品国产亚洲| 干丝袜人妻中文字幕| 女性生殖器流出的白浆| 麻豆国产97在线/欧美| 韩国高清视频一区二区三区| 高清午夜精品一区二区三区| 噜噜噜噜噜久久久久久91| 免费看日本二区| 久久久久久久久久人人人人人人| 成人影院久久| 97在线人人人人妻| 久久热精品热| av国产精品久久久久影院| 高清在线视频一区二区三区| 久久99精品国语久久久| 亚洲va在线va天堂va国产| 秋霞伦理黄片| 黄片wwwwww| 春色校园在线视频观看| 久久毛片免费看一区二区三区| 3wmmmm亚洲av在线观看| 国产美女午夜福利| 夫妻性生交免费视频一级片| 噜噜噜噜噜久久久久久91| 纯流量卡能插随身wifi吗| 国产熟女欧美一区二区| 免费久久久久久久精品成人欧美视频 | 成人毛片a级毛片在线播放| av.在线天堂| 日本免费在线观看一区| 男女免费视频国产| 久久6这里有精品| 国产大屁股一区二区在线视频| 久久热精品热| 精品久久久精品久久久| 亚洲性久久影院| 亚洲伊人久久精品综合| 妹子高潮喷水视频| 亚洲电影在线观看av| 激情五月婷婷亚洲| 国产大屁股一区二区在线视频| 亚洲精品亚洲一区二区| 18禁裸乳无遮挡免费网站照片| 成人午夜精彩视频在线观看| 日本一二三区视频观看| 亚洲性久久影院| 亚洲av男天堂| 欧美成人精品欧美一级黄| 国产精品久久久久久久电影| 国国产精品蜜臀av免费| 五月天丁香电影| 亚洲av综合色区一区| 如何舔出高潮| 新久久久久国产一级毛片| 久久精品国产亚洲av涩爱| 天堂中文最新版在线下载| 国产成人91sexporn| 一本久久精品| 亚洲av成人精品一区久久| 久久久久久九九精品二区国产| 国产亚洲最大av| 成年免费大片在线观看| 国产伦在线观看视频一区| 最黄视频免费看| 91精品国产国语对白视频| 直男gayav资源| 夜夜看夜夜爽夜夜摸| 免费不卡的大黄色大毛片视频在线观看| 成年人午夜在线观看视频| 免费人成在线观看视频色| 亚洲精品视频女| 多毛熟女@视频| 熟女人妻精品中文字幕| av视频免费观看在线观看| 国产永久视频网站| 夫妻午夜视频| 国产乱来视频区| 伊人久久精品亚洲午夜| 99热国产这里只有精品6| 一区二区三区免费毛片| av在线app专区| 一级毛片黄色毛片免费观看视频| 一级毛片电影观看| 亚洲国产成人一精品久久久| 久久久久久人妻| 亚洲色图av天堂| 国产欧美另类精品又又久久亚洲欧美| 少妇人妻精品综合一区二区| 国产精品一区二区三区四区免费观看| 久久久久精品久久久久真实原创| 中文乱码字字幕精品一区二区三区| 亚洲,欧美,日韩| 亚洲伊人久久精品综合| 最近的中文字幕免费完整| 国产成人精品福利久久| 又粗又硬又长又爽又黄的视频| 国产成人aa在线观看| 久久久欧美国产精品| 国产精品无大码| 国产av国产精品国产| 亚洲图色成人| 熟女电影av网| 啦啦啦在线观看免费高清www| 国产精品99久久99久久久不卡 | 午夜福利在线观看免费完整高清在| 夫妻午夜视频| 亚洲精品成人av观看孕妇| 18+在线观看网站| 日本午夜av视频| 成人影院久久| 欧美极品一区二区三区四区| 日本-黄色视频高清免费观看| 国产男女内射视频| a级毛色黄片| 免费大片黄手机在线观看| 精品国产乱码久久久久久小说| 国产亚洲最大av| 2021少妇久久久久久久久久久| 精品一区二区三卡| 久久热精品热| 精品一区二区免费观看| 日韩强制内射视频| 日韩亚洲欧美综合| 精品亚洲成a人片在线观看 | 高清av免费在线| 99热全是精品| 中文资源天堂在线| av网站免费在线观看视频| 国产亚洲欧美精品永久| 国产91av在线免费观看| 免费观看a级毛片全部| 国产成人aa在线观看| 精品国产露脸久久av麻豆| 嫩草影院新地址| 亚洲人成网站在线观看播放| 夫妻午夜视频| 精品视频人人做人人爽| 国产一区亚洲一区在线观看| 国产大屁股一区二区在线视频| 久久人人爽av亚洲精品天堂 | 国产精品一区二区在线不卡| 成人一区二区视频在线观看| 国产精品人妻久久久影院| 久久久久久久国产电影| 五月伊人婷婷丁香| 免费看日本二区| 亚洲精品乱码久久久久久按摩| 婷婷色综合www| 久久鲁丝午夜福利片| 精品一区在线观看国产| 秋霞伦理黄片| 麻豆国产97在线/欧美| 精品亚洲成国产av| 国产成人91sexporn| 欧美 日韩 精品 国产| 亚洲精品久久久久久婷婷小说| 日本欧美国产在线视频| 国产有黄有色有爽视频| 国语对白做爰xxxⅹ性视频网站| 国产成人freesex在线| 激情五月婷婷亚洲| 赤兔流量卡办理| 亚洲av福利一区| 直男gayav资源| av国产免费在线观看| 麻豆精品久久久久久蜜桃| 91久久精品电影网| a 毛片基地| 久久人妻熟女aⅴ| 美女福利国产在线 | 国产亚洲午夜精品一区二区久久| 免费人成在线观看视频色| 一级毛片 在线播放| 中文字幕亚洲精品专区| 日本色播在线视频| 欧美 日韩 精品 国产| 只有这里有精品99| 色视频www国产| 国产爱豆传媒在线观看| 亚洲四区av| av福利片在线观看| 亚洲国产最新在线播放| 人体艺术视频欧美日本| 午夜福利影视在线免费观看| 你懂的网址亚洲精品在线观看| 六月丁香七月| 免费观看无遮挡的男女| 亚洲久久久国产精品| 久久久精品94久久精品| 国产又色又爽无遮挡免| 国产精品爽爽va在线观看网站| 精品视频人人做人人爽| 亚洲欧洲日产国产| 夜夜爽夜夜爽视频| 亚洲真实伦在线观看| 大片免费播放器 马上看| 久久精品久久久久久久性| 亚洲精品乱久久久久久| 午夜激情福利司机影院| 久久婷婷青草| 99热国产这里只有精品6| 久久久久国产精品人妻一区二区| 欧美最新免费一区二区三区| 欧美精品国产亚洲| 丝瓜视频免费看黄片| 久久久亚洲精品成人影院| 免费观看无遮挡的男女| 人妻 亚洲 视频| 久久久久久伊人网av| 狂野欧美白嫩少妇大欣赏| av天堂中文字幕网| 精品一区在线观看国产| av在线app专区| 亚洲欧美精品专区久久| 免费观看的影片在线观看| 久久久久国产网址| 国产亚洲午夜精品一区二区久久| 亚洲精品乱码久久久v下载方式| 亚洲av国产av综合av卡| 夜夜骑夜夜射夜夜干| 国产成人freesex在线| 一级a做视频免费观看| 97在线人人人人妻| 免费观看a级毛片全部| 国产人妻一区二区三区在| 亚洲色图av天堂| 青春草亚洲视频在线观看| 老女人水多毛片| 在线观看免费日韩欧美大片 | 日韩人妻高清精品专区| 亚洲高清免费不卡视频| 成人综合一区亚洲| 久久久午夜欧美精品| 亚洲欧美日韩卡通动漫| 精品久久久噜噜| 街头女战士在线观看网站| 日本一二三区视频观看| 人妻一区二区av| 各种免费的搞黄视频| 久久久久性生活片| 成人亚洲欧美一区二区av| 噜噜噜噜噜久久久久久91| 亚洲精品一区蜜桃| 亚洲精品视频女| 男男h啪啪无遮挡| 成人一区二区视频在线观看| 婷婷色综合大香蕉| 成年人午夜在线观看视频| 日本黄色日本黄色录像| 成年免费大片在线观看| 国产乱人偷精品视频| 日韩人妻高清精品专区| 18禁裸乳无遮挡免费网站照片| 另类亚洲欧美激情| 国产精品一及| 如何舔出高潮| 久久99蜜桃精品久久| 国产精品嫩草影院av在线观看| 少妇高潮的动态图| 校园人妻丝袜中文字幕| 久久婷婷青草| www.色视频.com| 如何舔出高潮| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 久久久久久久久久成人| 国产精品免费大片| 久久久久久久国产电影| 久久久午夜欧美精品| 六月丁香七月| 少妇的逼水好多| 国产精品一及| 天堂俺去俺来也www色官网| 一级毛片黄色毛片免费观看视频| 黑人猛操日本美女一级片| 美女高潮的动态| 伊人久久精品亚洲午夜| 久久综合国产亚洲精品| 网址你懂的国产日韩在线| 少妇猛男粗大的猛烈进出视频| 国产男女超爽视频在线观看| 我要看黄色一级片免费的| 亚洲欧美精品自产自拍| 亚洲欧美一区二区三区国产| 91aial.com中文字幕在线观看| 国产精品蜜桃在线观看| 精品一区二区免费观看| 亚洲三级黄色毛片| 我要看日韩黄色一级片| 国产成人精品福利久久| 夜夜骑夜夜射夜夜干| 永久免费av网站大全| 在线亚洲精品国产二区图片欧美 | 国产一区二区三区综合在线观看 | 纵有疾风起免费观看全集完整版| 国产黄色视频一区二区在线观看| 中文欧美无线码| 欧美极品一区二区三区四区| 18禁裸乳无遮挡动漫免费视频| 精品视频人人做人人爽| 深爱激情五月婷婷| 一个人看的www免费观看视频| 一本久久精品| 日韩av不卡免费在线播放| 国产69精品久久久久777片| 亚洲成人中文字幕在线播放| 尤物成人国产欧美一区二区三区| 国产亚洲欧美精品永久| 色5月婷婷丁香| 亚洲,一卡二卡三卡| .国产精品久久| 我要看黄色一级片免费的| 国产爱豆传媒在线观看| 国产一级毛片在线| 高清日韩中文字幕在线| 80岁老熟妇乱子伦牲交| 美女xxoo啪啪120秒动态图| 少妇高潮的动态图| 在现免费观看毛片| 久久久午夜欧美精品| 国产精品蜜桃在线观看| 97超视频在线观看视频| 人人妻人人看人人澡| 最近中文字幕高清免费大全6| 国产黄片美女视频| 亚洲国产精品成人久久小说| 国产一区二区三区av在线| 丝瓜视频免费看黄片| 国产精品秋霞免费鲁丝片| 极品教师在线视频| 欧美bdsm另类| 日本欧美视频一区| 欧美国产精品一级二级三级 | 黑丝袜美女国产一区| 久久久久网色| 国产精品无大码| 亚洲图色成人| 97在线视频观看| 亚洲最大成人中文| 亚洲av免费高清在线观看| 亚洲精品亚洲一区二区| 免费黄网站久久成人精品| 大码成人一级视频| 欧美老熟妇乱子伦牲交| 亚洲av福利一区| 成年女人在线观看亚洲视频| 日本av手机在线免费观看| 最近中文字幕2019免费版| 乱系列少妇在线播放| 日本色播在线视频| 国产v大片淫在线免费观看| 国产极品天堂在线| 日韩欧美精品免费久久| 国产人妻一区二区三区在| 黄色日韩在线| 少妇裸体淫交视频免费看高清| 人人妻人人看人人澡| 777米奇影视久久| 新久久久久国产一级毛片| 亚洲怡红院男人天堂| av黄色大香蕉| 97在线视频观看| 黄色一级大片看看| 日韩亚洲欧美综合| 哪个播放器可以免费观看大片| 久久久久久久久久久丰满| 亚洲国产欧美在线一区| 国产精品99久久久久久久久| 高清欧美精品videossex| 亚洲一区二区三区欧美精品| 九九爱精品视频在线观看| 九色成人免费人妻av| 伦精品一区二区三区| 一区二区三区乱码不卡18| 亚洲一区二区三区欧美精品| 老熟女久久久| av一本久久久久| 成人免费观看视频高清| 熟女人妻精品中文字幕| 日韩精品有码人妻一区| 国产伦精品一区二区三区四那| 久久97久久精品| 三级国产精品欧美在线观看| 亚洲精品中文字幕在线视频 | 亚洲伊人久久精品综合| 午夜免费男女啪啪视频观看| 一区二区三区免费毛片| 汤姆久久久久久久影院中文字幕| 尤物成人国产欧美一区二区三区| 我要看日韩黄色一级片| 国产av码专区亚洲av| 亚洲精品成人av观看孕妇| 九九久久精品国产亚洲av麻豆| 国产欧美日韩精品一区二区| 午夜福利视频精品| 亚洲激情五月婷婷啪啪| 亚洲精品一区蜜桃| 噜噜噜噜噜久久久久久91| 久久国产精品大桥未久av | 日韩成人伦理影院| 一区二区三区乱码不卡18| 亚洲真实伦在线观看| 国产亚洲av片在线观看秒播厂| 成人漫画全彩无遮挡| 夜夜看夜夜爽夜夜摸| 欧美3d第一页| 啦啦啦在线观看免费高清www| 久久久久国产网址| 国产精品99久久99久久久不卡 | 精品人妻视频免费看| 一个人免费看片子| 超碰av人人做人人爽久久| 欧美xxⅹ黑人| 久久久久久久大尺度免费视频| 日韩av免费高清视频| 久久久久国产网址| 国产美女午夜福利| 亚洲国产av新网站| 国产伦理片在线播放av一区| 国内精品宾馆在线| 少妇被粗大猛烈的视频| 成人午夜精彩视频在线观看| 国模一区二区三区四区视频| 亚洲精品,欧美精品| 国产高清三级在线| 观看美女的网站| 国产精品一区二区性色av| 狠狠精品人妻久久久久久综合| 免费av中文字幕在线| av天堂中文字幕网| 欧美3d第一页| 丰满人妻一区二区三区视频av| 成人午夜精彩视频在线观看| 97热精品久久久久久| 街头女战士在线观看网站| 夜夜爽夜夜爽视频| 国产成人a∨麻豆精品| 这个男人来自地球电影免费观看 | 亚洲不卡免费看| 搡老乐熟女国产| 免费观看av网站的网址| 蜜桃亚洲精品一区二区三区| 午夜激情久久久久久久| 人人妻人人澡人人爽人人夜夜| 人妻系列 视频| 国产精品人妻久久久影院| 欧美激情极品国产一区二区三区 | 99久久精品热视频| 人妻一区二区av| 五月天丁香电影| 日韩不卡一区二区三区视频在线| 91精品一卡2卡3卡4卡| 美女cb高潮喷水在线观看| 久久国产乱子免费精品| 国产伦精品一区二区三区四那| 欧美日韩国产mv在线观看视频 | tube8黄色片| 久久久久久久精品精品| 天天躁夜夜躁狠狠久久av| 老师上课跳d突然被开到最大视频| 精品少妇久久久久久888优播| 国内少妇人妻偷人精品xxx网站| av黄色大香蕉| 久久精品久久精品一区二区三区| 中文字幕av成人在线电影| 观看免费一级毛片| 久久 成人 亚洲| 只有这里有精品99| 精品一区二区三卡| 好男人视频免费观看在线| 国产 一区精品| av视频免费观看在线观看| 春色校园在线视频观看| 在线观看一区二区三区| 熟女av电影| 男人爽女人下面视频在线观看| 久久99热这里只有精品18| 国产探花极品一区二区| 欧美xxxx性猛交bbbb| 黑人高潮一二区|