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

    廣東韓江潮州段浮游植物群落結(jié)構(gòu)特征與水質(zhì)評價

    2023-04-29 00:00:00林小植胡苑玲王瑞旋李冬梅林鴻生查廣才溫茹淑吳曉瓊
    水生態(tài)學雜志 2023年4期

    摘要:為探究廣東韓江潮州市飲用水源保護江段浮游植物群落結(jié)構(gòu)特征及水質(zhì)富營養(yǎng)化狀況,于2019年2月至2020年1月在潮州段設(shè)立4個采樣點,進行浮游植物調(diào)查和水質(zhì)監(jiān)測。結(jié)果表明,韓江潮州共鑒定出浮游植物8門103屬290種,其中綠藻門42屬99種,硅藻門21屬87種,藍藻門24屬70種。全年有9個月的種類組成以硅藻門為絕對優(yōu)勢種,綠藻門僅在6、9、10月種類數(shù)多于硅藻門,表明該江段浮游植物組成為硅藻-綠藻型。浮游植物豐度月變化在0.83×106~1.99×106個/L,月均豐度為0.90×106 個/L,全年以硅藻門豐度最高,占47.01%。優(yōu)勢種11種,其中硅藻門5種,隱藻門4種,藍藻門和裸藻門各1種,顆粒溝鏈藻(Melosira granulata)是絕對優(yōu)勢種。Shannon-Wiener多樣性指數(shù)(H)為3.14~5.08,平均值4.40;Pielou均勻度指數(shù)(J)為0.69~0.99,平均值0.87;Margalef豐富度指數(shù)(D)為2.60~6.85,平均值5.08,表明水質(zhì)為輕度污染。Cluster聚類和NMDS排序分析顯示,物種構(gòu)成可分為2個類群,類群I包括2019年2-7月,類群II包括2019年8-12月和2020年1月。Pearson相關(guān)性和冗余分析表明,影響浮游植物優(yōu)勢種的水環(huán)境因子主要是透明度和亞硝態(tài)氮,影響顆粒溝鏈藻分布的是鹽度和水溫,影響隱藻門分布的主要是透明度、pH和水溫。參照地表水環(huán)境質(zhì)量標準,韓江潮州段為Ⅱ類水,多樣性指數(shù)顯示為輕度污染,浮游植物群落結(jié)構(gòu)和優(yōu)勢物種顯示水體為偏中營養(yǎng)型。

    關(guān)鍵詞:浮游生物;群落結(jié)構(gòu);水環(huán)境;季節(jié)動態(tài);韓江潮州段

    中圖分類號:Q145,X824 " " " "文獻標志碼:A " " " "文章編號:1674-3075(2023)04-0052-09

    由于人類經(jīng)濟社會活動的影響,水環(huán)境生態(tài)系統(tǒng)和水資源污染日益嚴重,許多河流和湖泊的水質(zhì)富營養(yǎng)化正在快速惡化中(Douterelo et al, 2004;況琪軍等,2005)。20世紀初,韓博平等(2003)對廣東省各流域19座水庫進行研究,發(fā)現(xiàn)大部分水庫水質(zhì)狀況為中度污染,3座污染嚴重,水質(zhì)已超富營養(yǎng)化。浮游植物是水生生態(tài)系統(tǒng)中最重要的初級生產(chǎn)者,其生長速度快,對水生態(tài)系統(tǒng)敏感,能夠?qū)Νh(huán)境干擾作出反應(yīng),被認為是水體營養(yǎng)狀態(tài)的良好指標,也是反映水體污染狀況最好的指示生物(Paerl et al, 2010;Jiang et al,2014);其群落特征及動態(tài)變化是多個環(huán)境因子的時空綜合作用結(jié)果,與水環(huán)境中各種理化因子密切相關(guān),是評價水質(zhì)狀況及其變化趨勢的重要指標,可用于評價河流、湖泊富營養(yǎng)化(Kamenir et al,2004;況琪軍等,2005;Suikkanen et al,2007)。已有研究表明,通過浮游生物群落結(jié)構(gòu)特征結(jié)合水環(huán)境因子動態(tài)情況進行多樣性、Pearson相關(guān)性和RDA分析,多角度對水質(zhì)進行評價,能夠更好地反映水質(zhì)綜合狀況(吳慧雙等,2020)。

    韓江(115°13′~117°09′E,23°17′~26°05′N)是廣東省第二大河流,也是粵東1 400萬居民的飲用水源地,曾被評為“全國十大最美家鄉(xiāng)河”之一,具有重要的飲用水、灌溉、航運等綜合功能。上游由梅江和汀江匯合而成,中游為梅州的三河壩至潮州水利樞紐段,下游流經(jīng)潮州分為東溪、西溪、北溪,最終注入南海。近年來,由于經(jīng)濟快速發(fā)展,排入江河的污染物增多,水質(zhì)污染事件時有發(fā)生(張李源清和王玲, 2010),有時水質(zhì)為Ⅲ級(廣東省水利廳,2021)。有關(guān)韓江富營養(yǎng)狀況和浮游植物的研究甚少,僅見零星調(diào)查報告(陳霓彤,2014)。2019年2月至2020年1月,對韓江潮州段進行浮游植物組成和周年變化調(diào)查,并檢測水質(zhì)指標,分析其種群結(jié)構(gòu)、優(yōu)勢物種、污染指示物種、多樣性指數(shù)等變化,旨在為水生態(tài)學評價和水環(huán)境監(jiān)測提供依據(jù),為建立韓江飲用水源健康安全環(huán)境管理機制提供參考。

    1 " 材料與方法

    1.1 " 采樣地點

    本次調(diào)查共設(shè)置4個采樣點,分別是韓江西岸的竹竿山水廠(S1)和青龍廟(S4);東岸的橋東水廠(S2)和韓江大橋(S3),該江段為潮州市一級飲用水源保護江段(圖1)。每月采樣1次,共12次。

    1.2 " 采樣方法

    浮游生物定性、定量采樣方法參照《湖泊富營養(yǎng)化調(diào)查規(guī)范》進行(金相燦和屠清瑛,1990)。因采樣點的分布和距離相對較小(lt;5 km),均位于廣東省潮州水利樞紐上游的類湖泊環(huán)境中,統(tǒng)計數(shù)據(jù)取該江段4個采樣點的平均值。

    定量采樣用2 L有機玻璃采水器取水6 L,各采樣點設(shè)3個重復,使用25號浮游生物網(wǎng)(64 μm)過濾,以0.1 mL浮游生物計數(shù)框計數(shù)(章宗涉和黃祥飛,1991)。定性采樣用25號浮游生物網(wǎng),參照《湖泊富營養(yǎng)化調(diào)查規(guī)范》,在Nikon E100光學顯微鏡進行拍照并分類鑒定(胡鴻鈞和魏印心,2006)。

    1.3 " 水環(huán)境指標測定

    用哈希HQ40d雙通道多參數(shù)分析儀在采樣點測定水體溶解氧、鹽度、水溫、電導率和pH;用塞氏盤測定透明度,用LGY-II型智能流速儀測定流速。

    取表層水用于大腸桿菌、NH3-N、NO2-N、高錳酸鹽指數(shù)的測定,采用哈納HI83203微電腦多參數(shù)離子濃度測定儀測定水中NH3-N和NO2-N,高錳酸鹽指數(shù)、大腸菌群檢測參照《水和廢水監(jiān)測分析方法》(國家環(huán)保總局, 2002)。

    1.4 " 數(shù)據(jù)分析

    在Excel統(tǒng)計物種數(shù)、豐度、出現(xiàn)頻度,DPS 9.50計算浮游植物多樣性指數(shù),SPSS 19.0分析浮游植物與水環(huán)境因子的Pearson相關(guān)關(guān)系,采用Canoco5.0對全年浮游植物優(yōu)勢種與水環(huán)境因子的相關(guān)關(guān)系進行冗余分析。相對豐富度(Q)、優(yōu)勢度(Y)指數(shù)、Shannon-Wiener多樣性指數(shù)(H)、Pielou均勻度指數(shù)(J)和Margalef豐富度指數(shù)(D)計算公式(吳朝和張慶國,2009;王麗卿等,2011;Siham et al, 2019)如下:

    Q = (ni/N)×100% " " " " " " " " " " " "①

    Y = (ni/N)×fi " " " " " " " " " " " " " "②

    H = -[i=1snin] log2[nin] " " " " " ③

    J = H/log2S " " " " " " " " " " " " " ④

    D = (S - 1)/log2N " " " " " " " " " " " ⑤

    式中:ni為每月第i種浮游植物的個體數(shù),N為全年浮游植物總個體數(shù),fi為第i種浮游植物在12次采樣中的出現(xiàn)的頻率;n為每月浮游植物總個體數(shù);S為每月浮游植物總種數(shù)。

    浮游生物多樣性指數(shù)的污染狀態(tài)水質(zhì)評價標準(鄭丙輝等,2007)見表1。

    1.5 " 季節(jié)變化聚類

    以12個月浮游植物的種類和數(shù)量組成原始數(shù)據(jù)矩陣,用Bray-Curtis相似性系數(shù)為基礎(chǔ)構(gòu)建月份之間的相似性矩陣。采用等級聚類(UPGMA)的分類方法和非度量多維標度排序法(NMDS)分析潮州段浮游植物群落的季節(jié)變化特征。本文四季劃分為春季(2-4月)、夏季(5-7月)、秋季(8-10月)、冬季(11-1月),建立聚類分析圖和NMDS平面圖。任意兩樣本間的距離代表其相似程度;聚類越近,相似程度越高。以協(xié)強系數(shù)(Stress)檢驗NMDS分析結(jié)果的優(yōu)劣(Clarke amp; Warwick,2001;林小植等,2016)。

    通過相似性分析(ANOSIM)檢驗不同矩陣間的差異顯著性,組內(nèi)所有位點的平均秩相似性與組間所有位點間的平均秩相似性差異用R表示,域值為[-1, 1] (林小植等, 2016) 。應(yīng)用SIMPER (Similarity Percentages)方法分析每種浮游植物群落組內(nèi)相似性和組間相異性的平均貢獻率,將貢獻率≥5%的浮游植物定義為主要特征種;對全年優(yōu)勢物種浮游植物進行聚類和排序分析。所有群落結(jié)構(gòu)分析用PRIMER 5.0軟件包(Plymouth Marine Laboratory, UK)完成。

    2 " 結(jié)果

    2.1 " 浮游植物種類組成

    2019年2月至2020年1月的采樣調(diào)查,韓江潮州段共鑒定出浮游植物8門103屬290種。其中綠藻門最多,達42屬99種,硅藻門21屬87種,藍藻門24屬70種,裸藻門16種,甲藻門7種,黃藻門6種,隱藻門4種,金藻門1種,種類組成及其占比見圖2。

    從時間分布上看(圖3),浮游植物物種數(shù)出現(xiàn)最多在2019年5月,有127種,最少的在2019年6月,僅44種。

    各月的物種構(gòu)成中,硅藻門、綠藻門、藍藻門的物種在所有物種中占比較大。硅藻門的物種數(shù)在2019年2、3、4、5、7、8、11、12月和2020年1月中占比均最大,綠藻門在6、9、10月物種數(shù)大于硅藻門(圖4)。因此認為韓江潮州段浮游植物的種類組成特征為硅藻-綠藻型。

    2019年6月韓江潮州段僅由硅藻門、藍藻門、綠藻門構(gòu)成,其余各月均由3~7門浮游植物構(gòu)成;金藻門只有分歧錐囊藻(Dinobryon drivergens)1種,僅在2020年1月采到(圖3)。

    2.2 " 優(yōu)勢種組成和豐度

    全年Y≥0.02的優(yōu)勢種共有11種,其中硅藻門5種,顆粒溝鏈藻為絕對優(yōu)勢種;隱藻門4種,綠藻門、藍藻門各1種(表2)。

    2.3 " 浮游植物豐度構(gòu)成和季節(jié)變化

    韓江潮州段浮游植物月豐度在0.83×106~1.99×106 個/L,月均豐度0.90×106個/L。全年硅藻門豐度為5.06×106個/L,占全部類群的47.01%,全年的綠藻門、隱藻門、藍藻門占比也較大,分別占20.97%、16.87%、11.38%(圖5)。

    浮游植物全年豐度最高值出現(xiàn)在2019年10月,為1.99×106個/L,月均豐度0.83×106 個/L;豐度最低值在2019年2月,為0.20×106 個/L(圖6)。硅藻門全年都有出現(xiàn),且數(shù)量占比大,隱藻門則在秋末冬初出現(xiàn),金藻門和黃藻門在定量采樣中未采集到。根據(jù)湖泊富營養(yǎng)類型評價的細胞密度標準(況琪軍等,2005),潮州市江段水質(zhì)處于貧營養(yǎng)與貧中營養(yǎng)型。

    2.4 " 浮游植物多樣性指數(shù)

    韓江潮州段浮游植物Shannon-Wiener多樣性指數(shù)(H)均值為4.40(3.14~5.08),2019年3月最低,為3.14;Pielou均勻度指數(shù)(J)均值為0.87(0.69~0.99);Margalef豐富度指數(shù)(D)均值為5.08(2.60~6.85)(圖7)。

    2.5 " 理化指標

    根據(jù)地表水環(huán)境質(zhì)量標準(GB3838-2002),韓江潮州段的高錳酸鹽指數(shù)、溶解氧、氨氮、糞大腸桿菌各月平均值均在Ⅱ類水范圍內(nèi)。水體環(huán)境理化指標見表3。

    2.6 " 浮游植物與水環(huán)境因子的相關(guān)性及冗余分析

    浮游植物豐度與水環(huán)境因子的Pearson相關(guān)性分析顯示,隱藻門豐度與流速(FR)呈顯著正相關(guān)(r=0.768,P=0.011),與透明度(SD)呈極顯著正相關(guān)(r=0.872,P=0.002)。

    以浮游植物優(yōu)勢種(Y≥0.02)的生物量與各月份的相關(guān)水環(huán)境因子數(shù)據(jù)進行冗余(RDA)分析(圖8)。結(jié)果表明,RDA第1、2軸的特征值分別為0.4632、0.2777,共解釋了韓江潮州段浮游植物總變異數(shù)的98.4%,第1軸和第2軸浮游植物優(yōu)勢種與水環(huán)境因子的典型相關(guān)都達到了0.999,表明具有很強的相關(guān)性,透明度(SD)和亞硝態(tài)氮(NO2-N)對浮游植物優(yōu)勢種影響最大,且呈顯著相關(guān)(Plt;0.05),大部分優(yōu)勢種與透明度呈正相關(guān)關(guān)系,與亞硝態(tài)氮呈現(xiàn)負相關(guān)關(guān)系。以隱藻門為主(S8、S9、S10、S11)的優(yōu)勢種與FR、pH、SD呈正相關(guān)關(guān)系,與水溫(WT)呈負相關(guān)關(guān)系。硅藻門中顆粒溝鏈藻(S2)與鹽度(SAL)呈正相關(guān)關(guān)系,與WT呈負相關(guān)關(guān)系。

    2.7 " 聚類分析及群落特征

    等級聚類分析表明,可將12個調(diào)查月分成2個類群,類群I包括2019年2-7月;類群II包括2019年8-12月和2020年1月(圖9-a)。類群I中對群落內(nèi)貢獻率≥5%的種類有克洛脆桿藻(Fragilaria crotonensis)、鈍脆桿藻、尖針桿藻(Synedra acus)、游絲藻(Planctonema lauterbornii)、顆粒溝鏈藻、變異直鏈藻;類群II中對群落內(nèi)貢獻率≥5%的種類有顆粒溝鏈藻、顆粒溝鏈藻最窄變種、變異直鏈藻、嚙蝕隱藻、吻狀隱藻、尖尾藍隱藻。類群I與II的差異系數(shù)為86.95%,對類群間差異貢獻率≥5%的種類有顆粒溝鏈藻、窄異極藻延長變種、吻狀隱藻。

    NMDS排序的協(xié)強系數(shù)為0.14,表明NMDS排序的二維點陣圖結(jié)果可以接受。聚類分析的2個組群在NMDS圖上能夠清晰地表現(xiàn)在2個不同的區(qū)域中,標識為類群I和類群II(圖9-b)。ANOSM差異顯著性檢驗結(jié)果顯示,r=0.887 (lt;1), P=0.002 (lt;0.01),表明2個類群之間差異極顯著。

    將2019年8-12月和2020年1月與2019年2-7月的多樣性指數(shù)(H、D)進行獨立樣本T檢驗,發(fā)現(xiàn)2個類群的H和D之間均存在極顯著差異(P=0.007和P=0.002),與聚類分析結(jié)果一致。

    對11種全年優(yōu)勢物種進行Cluster聚類(圖10-a)和NMDS排序分析(圖10-b),NMDS排序的協(xié)強系數(shù)為0.08,表明二維點陣圖可以很好地解釋浮游植物物種間相似關(guān)系。11種浮游植物可分為4個類群,ANOSM顯著性檢驗結(jié)果顯示,r=0.907 (lt;1),P=0.001 (lt;0.01),表明2個群聚存在極顯著差異。類群A包括卵形隱藻、吻狀隱藻、嚙蝕隱藻、尖尾藍隱藻,其群落特點是在秋冬季出現(xiàn)。類群B包括細微顫藻,在4-5月出現(xiàn)。類群C包括鈍脆桿藻、變異直鏈藻、顆粒溝鏈藻、顆粒溝鏈藻最窄變種、小球藻,其群落特點是四季常見種。類群D包括窄異極藻延長變種,在10月中出現(xiàn)。

    3 " 討論

    3.1 " 浮游植物的群落結(jié)構(gòu)特征與季節(jié)演替

    浮游植物是重要的生產(chǎn)者,其群落構(gòu)成和特征是水生態(tài)系統(tǒng)的綜合表現(xiàn)(Jiang et al,2014)。韓江潮州段12個月的浮游植物構(gòu)成可分為2個類群,客觀反映了該江段浮游植物組成的季節(jié)性波動特征。韓江地處廣東東部,表層水溫19.8~30.1℃,適合大多數(shù)耐高溫藻類生長,特別是硅藻門、綠藻門和藍藻門。在類群I和類群II中,對群落貢獻率≥5%的共有物種是顆粒溝鏈藻和變異直鏈藻,可見在浮游植物群落結(jié)構(gòu)中,同一水體不同時間的類群也存在群落組成的邊界模糊性;況琪軍等(2005)研究中也有類似闡述,即不同類型水體的指示物種均有營養(yǎng)過渡性藻類存在,有的藻類適應(yīng)性很強,能夠生存在不同營養(yǎng)水平的水體中(金相燦和屠清瑛,1990;況琪軍等, 2005) 。

    浮游植物群落結(jié)構(gòu)和豐度受到季節(jié)環(huán)境因子、水動力學和營養(yǎng)鹽的共同作用(韓博平等,2003;Jiang et al,2014 )。從潮州段浮游植物豐度看,豐度最小的是2月,正處于冬春季節(jié)的交替期;與陳霓彤(2014 )對韓江下游梅溪河段的研究結(jié)果一致。2月雨量少,水溫最低,水的流動性較小,其浮游植物構(gòu)成主要是能夠適應(yīng)低溫生活的硅藻、綠藻和裸藻,其中硅藻數(shù)量最多,占59.61%,裸藻和綠藻各占10.10%和30.29%。這種浮游植物豐度構(gòu)成大致與同處亞熱帶地區(qū)的廣東省新豐江水庫相類似,其浮游植物豐度最低也出現(xiàn)在2月(趙孟緒等,2005);而珠江廣州段的周年調(diào)研發(fā)現(xiàn),其豐度最低出現(xiàn)在6月,由于該江段為重度富營養(yǎng)化水體,其浮游植物豐度構(gòu)成中藍藻、綠藻的占比較大(張俊逸等, 2011) 。本研究中,浮游植物豐度變化在3-5月出現(xiàn)一個小高峰,9-12月出現(xiàn)最高峰,其原因主要與韓江潮州段的水文環(huán)境因子和季節(jié)有關(guān)。冬春交替季節(jié)(2月)的水溫逐漸上升,浮游植物初級生產(chǎn)力會加速消耗冬季積累下來的營養(yǎng)鹽,加之降水量少、農(nóng)業(yè)灌溉、生活用水和自然蒸騰作用,使得浮游植物豐度經(jīng)過冬季的峰值之后逐漸降低到較低水平;6月是夏季汛期降雨量最大的月份,水流量大,流速較快,大量顆粒物沉降到底部,營養(yǎng)水平相對較低,水溫較高,不利于浮游植物的漂浮生活,因此浮游植物的豐度和初級生產(chǎn)力在此階段也較低;秋冬季水流相對平穩(wěn),對流垂直混合使得表層營養(yǎng)鹽得到補充,陽光充足,水溫適宜,因而浮游植物的生長得到了較好的環(huán)境支持,初級生產(chǎn)力增長最大。

    3.2 " 浮游植物群落結(jié)構(gòu)與環(huán)境因子的關(guān)系

    透明度是評價水質(zhì)的重要指標,是表征水體中懸浮物多少的重要依據(jù),也可以反映水的渾濁程度。當透明度降低時,會影響光在水下的穿透和傳播,導致浮游植物的光合作用受阻,不利于藻類生長;在淮河流域西淝河、里下河大縱湖也發(fā)現(xiàn)浮游植物與透明度存在顯著相關(guān)關(guān)系(鄧建明等,2010;王順天等,2020;張珊等, 2020) 。

    亞硝態(tài)氮是氨氮氧化或硝態(tài)氮還原過程的中間產(chǎn)物,環(huán)境污染加劇導致亞硝態(tài)氮含量上升(王九龍等,2012) ;氮濃度變化會直接影響浮游植物生長,不適合耐低氮的硅藻生長,但有助于耐高氮的藻類成為優(yōu)勢種(Elser et al, 2009) ,與盧云黎等(2021) 報道的漢江集家嘴受到亞硝態(tài)氮的影響一致。硅藻門全年數(shù)量分布廣,在秋末冬初時豐度很高,21℃為其適宜生存溫度,硅藻門的顆粒溝鏈藻是絕對優(yōu)勢種,直鏈藻對水環(huán)境變化非常敏感,容易受到水溫影響,而且其有一定的適鹽性,相比其他藻類,更能適應(yīng)高鹽度,在富營養(yǎng)化水體及污染水體中極易形成優(yōu)勢種(Wang et al, 2009;陳威等, 2018) 。

    隱藻門的繁殖主要集中在秋冬季,由于其特殊的光敏感度和營養(yǎng)型,在全年溫度較低的時期為優(yōu)勢藻種(陳霓彤,2014) ;而高pH值的堿性環(huán)境有助于藻類進行光合作用并促進其增殖(徐暢,2020) 。本研究中,隱藻門的豐度與透明度呈正相關(guān),與水溫呈顯著負相關(guān),符合隱藻門的生長特點。

    2019年2月至2020年1月,潮州段浮游植物種類組成的變化與韓江地區(qū)的氣候有關(guān),受海洋性東南季風的影響,4-9月為每年的汛期,陰天多,降水量占全年降雨量的70%以上。浮游植物的群落結(jié)構(gòu)與光照有關(guān),陰天會影響光合作用(林少君等, 2005) ,物種構(gòu)成較單一,造成2019年6月的物種數(shù)較少。

    3.3 " 不同指標的韓江潮州段水質(zhì)評價

    理化指標單因子評價法是我國水質(zhì)監(jiān)測公報的綜合水質(zhì)評價方法,參照地表水環(huán)境質(zhì)量標準(GB3838-2002),選擇單項最差指標所屬類別來確定水域綜合水質(zhì)類別(戴潤泉等, 2004;尹海龍和徐祖信,2008;叢銘等, 2021)。結(jié)合表4,本次監(jiān)測表明韓江潮州段水質(zhì)屬于II類,水體為堿性。

    多樣性指數(shù)通常是研究浮游植物群落結(jié)構(gòu)特征和動態(tài)規(guī)律的測度工具,同時也是判斷水體富營養(yǎng)狀況的常用指標。H、J、D指數(shù)值越大,其群落結(jié)構(gòu)越復雜,群落所包含的物種多度和豐度信息量也越大,穩(wěn)定性越強,水質(zhì)越好;但多樣性指數(shù)減小時,群落結(jié)構(gòu)的物種多度和豐度會趨于簡單,穩(wěn)定性變?nèi)?,水質(zhì)變差(金相燦和屠清瑛,1990;況琪軍等, 2005) 。根據(jù)多樣性指數(shù)評價標準(表1),韓江潮州段水質(zhì)總體處于輕度污染狀態(tài),對比廣東省部分河湖的多樣性指數(shù)(表5),潮州段水質(zhì)污染程度比小東江下游、連江、孔江、陸豐海域、珠海海域、珠江廣州段低,總體水質(zhì)較好。

    湖泊的富營養(yǎng)化與浮游植物群落密切相關(guān),其結(jié)構(gòu)特征及優(yōu)勢種變化趨勢是分析和預警水體富營養(yǎng)化和藍藻水華的依據(jù)之一(金相燦和屠清瑛,1990)。研究表明,浮游藻類的生物量及種類組成,對于評價河流湖泊富營養(yǎng)化狀況是有效的(江源等,2013;Jiang et al, 2014)。韓江潮州段浮游植物種類組成為硅藻-綠藻型(物種33.90%、29.79%,豐度20.24%、46.69%),根據(jù)不同藻類組成類型的指示作用顯示,韓江潮州江段為寡-中營養(yǎng)型至中營養(yǎng)型,從藍藻的種類和數(shù)量判斷(物種23.97%、豐度11.52%),該江段為偏中營養(yǎng)型。

    藻類作為指示物種,其種類和數(shù)量可以反映水體營養(yǎng)狀況,特別是對于優(yōu)勢物種的存在,其指示等級可劃分可為超富營養(yǎng)型、富營養(yǎng)型、中營養(yǎng)型和寡營養(yǎng)型(況琪軍等,2005;張麗娜和于洪賢,2010);通過對浮游植物的優(yōu)勢種進行等級劃分(表2),大多數(shù)優(yōu)勢種的指示等級為β-ms型和α~β-ms型,表明水體營養(yǎng)型為中營養(yǎng)型至中-中富營養(yǎng)型。

    參考文獻

    陳霓彤,2014. 韓江下游梅溪河段浮游藻類動態(tài)監(jiān)測及分析[J]. 城鎮(zhèn)供水, (1):32-34.

    陳威,魏南,金小偉,等,2018. 松花江哈爾濱段藻類植物分布及其與環(huán)境因子的關(guān)系[J]. 中國環(huán)境監(jiān)測, 34(4):102-110.

    叢銘,陽輝,張曉靜,等,2021. 單因子法與可變模糊法在水質(zhì)評價中的應(yīng)用[J]. 南水北調(diào)與水利科技, 19(4):720-728.

    戴潤泉,臧小平,邱光勝,2004. 三峽水庫蓄水前庫區(qū)水質(zhì)狀況研究[J]. 長江流域資源與環(huán)境,(2):124-127.

    鄧建明,蔡永久,陳宇煒,等,2010. 洪湖浮游植物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系[J]. 湖泊科學, 22(1):70-78.

    廣東省水利廳,2021. 2008-2020水資源公報[N/OL]. http://slt.gd.gov.cn/tjsj8778/index.html.

    國家環(huán)??偩郑?002. 水和廢水監(jiān)測分析方法[M]. 4版. 北京:中國環(huán)境科學出版社.

    韓博平,李鐵,林旭鈿,2003. 廣東省大中型水庫富營養(yǎng)化現(xiàn)狀與防治對策研究[M]. 北京:科學出版社.

    何森華,丘泓欣,鄧婉璐,等,2021. 廣東孔江國家濕地公園水體中的浮游植物群落結(jié)構(gòu)特征[J]. 濕地科學, 19(2):218-225.

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

    黃報遠,李秋華,陳桐生,2009. 廣東連江梯級電站開發(fā)后春季河流浮游生物群落結(jié)構(gòu)特征[J]. 安徽農(nóng)業(yè)科學, 37(7):3137-3140.

    江源,彭秋志,廖劍宇,等,2013. 浮游藻類與河流生境關(guān)系研究進展與展望[J]. 資源科學, 35(3):461-472.

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

    況琪軍,馬沛明,胡征宇,等,2005. 湖泊富營養(yǎng)化的藻類生物學評價與治理研究進展[J]. 安全與環(huán)境學報, 5(2):87-91.

    林少君,顧繼光,魏鵬,等,2005. 廣東省公平水庫與星湖生態(tài)特征的對比分析[J]. 生態(tài)學雜志,(7):773-779.

    林小植,李冬梅,劉煥章,等,2016. 廣東韓江潮州江段魚類多樣性及季節(jié)變化[J]. 生物多樣性, 24(2):185-194.

    盧云黎,王琪,羅春艷,2021. 漢江集家嘴2017-2019年浮游植物變化研究[J]. 水利水電快報, 42(5):43-48.

    粟麗,黃梓榮,陳作志,2017. 廣東沿岸不同海域浮游植物群落結(jié)構(gòu)特征的比較分析[J]. 海洋環(huán)境科學, 36(1):61-65.

    王九龍,張明明,王金葉,2012. 青島農(nóng)業(yè)大學生活飲用水亞硝態(tài)氮檢測與情況分析[J]. 青島農(nóng)業(yè)大學學報(自然科學版), 29(2):111-114.

    王麗卿,施榮,季高華,等,2011. 淀山湖浮游植物群落特征及其演替規(guī)律[J]. 生物多樣性, 19(1):48-56.

    王順天,雷俊山,賈海燕,等,2020. 三峽水庫浮游植物群落特征及水體富營養(yǎng)化評價[J]. 三峽生態(tài)環(huán)境監(jiān)測, 5(1):32-41.

    吳朝,張慶國,2009. 淡水浮游生物多樣性及數(shù)量分析方法研究進展[J]. 安徽農(nóng)學通報, 15(12):41-42.

    吳慧雙,胡泓,李正炎,等,2020. 2018年秋季大沽河浮游生物群落特征及其與水環(huán)境因子的關(guān)系[J]. 海洋湖沼通報,(5):148-154.

    徐暢,2020. 微山湖浮游生物群落時空變化特征及其驅(qū)動因子分析[D]. 濟南:濟南大學.

    尹海龍,徐祖信,2008. 河流綜合水質(zhì)評價方法比較研究[J]. 長江流域資源與環(huán)境, (5):729-733.

    張才學,詹冬玲,陳春亮,等, 2010. 小東江茂名段浮游生物及污染狀況分析[J]. 廣東海洋大學學報, 30(4):22-28.

    張俊逸,蔣江巒,劉擎,等,2011. 珠江廣州段微生物和浮游植物群落與水質(zhì)特征研究[J]. 水生態(tài)學雜志, 32(2):38-46.

    張李源清,王玲,2010. 紫金礦業(yè)污水入粵致漁業(yè)受損[N]. 中國經(jīng)濟時報, 2010-07-23(1).

    張麗娜,于洪賢,2010. 牡丹江海浪河夏季浮游植物群落結(jié)構(gòu)和水質(zhì)評價[J]. 水產(chǎn)學雜志, 23(1):40-46.

    張珊,商乃萱,張金羽,等,2020. 淮河流域西淝河浮游植物群落結(jié)構(gòu)特征[J]. 水生態(tài)學雜志, 41(4):81-89.

    章宗涉,黃祥飛,1991. 淡水浮游生物研究方法[M]. 北京:科學出版社.

    趙孟緒,雷臘梅,韓博平,2005. 亞熱帶水庫浮游植物群落季節(jié)變化及其影響因素分析——以湯溪水庫為例[J]. 熱帶亞熱帶植物學報, (5):386-392.

    鄭丙輝,田自強,張雷,等, 2007. 太湖西岸湖濱帶水生生物分布特征及水質(zhì)營養(yǎng)狀況[J]. 生態(tài)學報,(10):4214-4223.

    Clarke K R, Warwick R M, 2001. Change in Marine Communities: an An Approach to Statistical Analysis and Interpretation[M]. PRIMER-E Ltd., Plymouth.

    Douterelo I, Perona E, Mateo P, 2004. Use of cyanobacteria to assess water quality in running waters[J]. Enviromental Poluution, 127:377-384.

    Elser J J, Andersen T, Baron J S, et al, 2009. Shifts in Lake N:P Stoichiometry and Nutrient Limitation Driven by Atmospheric Nitrogen Deposition[J]. Science, 326:835-837.

    Jiang Y J, He W, Liu W X, et al, 2014. The seasonal and spatial variations of phytoplankton community and their correlation with environmental factors in a large eutrophic Chinese lake (Lake Chaohu) [J]. Ecological Indicatiors, 40:58-67.

    Kamenir Y, Dubinsky Z, Zohary T, 2004. Phytoplankton size structure stability in a meso-eutrophic subtropical lake[J]. Hydrobiologia, 520(1):89-104.

    Paerl H W, Rossignol K L, Hall S N, et al, 2010. Phytoplankton Community Indicators of Short-and Long-term Ecological Change in the Anthropogenically and Climatically Impacted Neuse River Estuary, North Carolina, USA[J]. Estuaries and Coasts, 33(2):485-497.

    Siham A, Somia H, Abdessamad R M, et al, 2019. Seasonal variation of water quality and phytoplankton dynamics and diversity in the surface water of Boukourdane Lake, Algeria[J]. Arabian Journal of Geosciences, DOI: 10.1007/s12517-018-4164-4

    Suikkanen S, Laamanen M, Huttunen M, 2007. Long-term changes in summer phytoplankton communities of the open northern Baltic Sea Estuarine[J]. Coastal and Shelf Science, 71(3/4): 580-592.

    Wang C, Li X H, Lai Z N, et al, 2009. Seasonal variations of Aulacoseira granulata population abundance in the Pearl River Estuary[J]. Estuarine, Coastal and Shelf Science, 85(4):585-592.

    (責任編輯 " 萬月華)

    Phytoplankton Community Structure and Water Quality Assessment

    of the Chaozhou Section of Hanjiang River

    LIN Xiao‐zhi1, HU Yuan‐ling1, WANG Rui‐xuan1, LI Dong‐mei1, LIN Hong‐sheng1,

    ZHA Guang‐cai1, WEN Ru‐shu2, WU Xiao‐qiong1

    (1. Hanshan Normal University, Chaozhou " 521041, P. R. China;

    2. Jiaying University, Meizhou " 514015, P. R. China)

    Abstract:Hanjiang River is the second largest river in Guangdong Province and serves as the drinking water source for eastern Guangdong. In this study, we characterized the phytoplankton community structure and eutrophication status of the Chaozhou section of Hanjiang River, which is the drinking water source for Chaozhou City. From February 2019 to January 2020, phytoplankton and water quality monitoring was carried out monthly at 4 sampling sites, focusing on phytoplankton community structure, dominant species, diversity indices, pollution indicator species and 11 water environment parameters. Our aim was to provide data for aquatic ecosystem evaluation and water quality monitoring of Hanjiang River, and to provide a reference for developing a management strategy for maintaining a healthy and safe environment of drinking water source. A total of 290 phytoplankton species from 103 genera and 8 phyla were identified, including Chlorophyta (99 species, 42 genera), "Bacillariophyta (87 species, 21 genera) and Cyanophyta (70 species, 24 genera). While species richness of the diatom phylum was absolutely dominant for nine months of the year, Chlorophyta had higher species richness in June, September and October. Thus, the phytoplankton community of the river section was of the diatom-green algae type. The monthly abundance of phytoplankton varied from 0.83×106 to 1.99×106 cells/L, with an average abundance of 0.90×106 cells/L. The abundance of diatoms was the highest, accounting for 47.01% of the total. There were 11 dominant species identified during the investigation, consisting of 5 Bacillariophyta species, 4 Cryptophyta species, 1 Cyanophyta species and 1 Euglena species, among which Melosira granulata was the absolutely dominant species. The Shannon-Wiener diversity (H), Pielou evenness (J) and Margalef richness (D) indices were in the range of 3.14-5.08, 0.69-0.99 and 2.60-6.85, with the average values of 4.40, 0.87 and 5.08, respectively, indicating a slightly polluted condition. Cluster analysis and NMDS ordination results show that phytoplankton communities from February to July 2019 was clustered into one group, and communities from August to December 2019 and January 2020 was clustered into another group. Pearson correlation and redundancy analysis show that transparency and nitrite nitrogen were the primary water environmental factors affecting the dominant species of phytoplankton. The factors affecting the distribution of Melosira granulata were salinity and water temperature, and the factors affecting the distribution of Cryptophyta were transparency, pH and water temperature. According to the environmental quality standard of surface water, the Chaozhou section of Hanjiang River met the standard of Class II surface water. The diversity indices also indicated that the water was lightly polluted. Analysis of phytoplankton community structure and dominant species indicated that the water is mesotrophic.

    Key words: plankton; community structure; water environment; seasonal dynamics; Chaozhou section of Hanjiang River

    av黄色大香蕉| av中文乱码字幕在线| 在线观看一区二区三区| av视频在线观看入口| 国产成人av教育| 欧美高清性xxxxhd video| 啦啦啦韩国在线观看视频| 波野结衣二区三区在线| 99久久精品热视频| 极品教师在线免费播放| 欧美3d第一页| 啦啦啦韩国在线观看视频| 熟女电影av网| 九九久久精品国产亚洲av麻豆| 超碰av人人做人人爽久久| 久久九九热精品免费| 久久精品国产亚洲网站| 在线播放无遮挡| 一区二区三区激情视频| 九九久久精品国产亚洲av麻豆| 欧美zozozo另类| 男插女下体视频免费在线播放| 搡女人真爽免费视频火全软件 | 最近在线观看免费完整版| videossex国产| 亚洲aⅴ乱码一区二区在线播放| 国产精品自产拍在线观看55亚洲| 国产亚洲91精品色在线| 久久久久久国产a免费观看| 丰满的人妻完整版| 日日撸夜夜添| 国产精品一及| 亚洲av免费在线观看| 亚洲18禁久久av| 亚洲精品国产成人久久av| 噜噜噜噜噜久久久久久91| 五月伊人婷婷丁香| 一进一出好大好爽视频| 色综合亚洲欧美另类图片| 深夜精品福利| 欧美xxxx黑人xx丫x性爽| 亚洲久久久久久中文字幕| 内射极品少妇av片p| 精品人妻1区二区| 麻豆国产97在线/欧美| 桃红色精品国产亚洲av| 精品久久久久久久久久久久久| 国内精品一区二区在线观看| 草草在线视频免费看| 少妇裸体淫交视频免费看高清| 无人区码免费观看不卡| 亚洲性久久影院| 国产真实伦视频高清在线观看 | 色综合婷婷激情| 日本成人三级电影网站| 97碰自拍视频| 成人毛片a级毛片在线播放| av在线蜜桃| 久久久午夜欧美精品| 亚洲 国产 在线| 久久人人精品亚洲av| 97人妻精品一区二区三区麻豆| 日日撸夜夜添| 色哟哟·www| 国产黄色小视频在线观看| 啪啪无遮挡十八禁网站| 日本欧美国产在线视频| 99热6这里只有精品| 91精品国产九色| 国产av一区在线观看免费| 亚洲人成网站高清观看| 免费看光身美女| 午夜免费成人在线视频| 最近中文字幕高清免费大全6 | 久久久久久九九精品二区国产| 国产精品人妻久久久影院| 高清毛片免费观看视频网站| 成人综合一区亚洲| 久久精品人妻少妇| 国内毛片毛片毛片毛片毛片| 老司机福利观看| 黄色配什么色好看| 香蕉av资源在线| 最近最新免费中文字幕在线| 国产成人a区在线观看| 精品一区二区三区人妻视频| 欧美性猛交黑人性爽| 18禁黄网站禁片午夜丰满| 成人欧美大片| av中文乱码字幕在线| 18禁在线播放成人免费| 成人鲁丝片一二三区免费| 日日摸夜夜添夜夜添av毛片 | 日本在线视频免费播放| 国产午夜福利久久久久久| 在线免费观看的www视频| 最后的刺客免费高清国语| 观看免费一级毛片| 日韩一本色道免费dvd| 日韩强制内射视频| 男人狂女人下面高潮的视频| .国产精品久久| 久久久成人免费电影| 国产精品久久视频播放| 国产大屁股一区二区在线视频| 小说图片视频综合网站| 久久久久久伊人网av| 日韩欧美精品免费久久| 免费人成在线观看视频色| 国产精品综合久久久久久久免费| 成人特级av手机在线观看| 18禁黄网站禁片午夜丰满| 欧美高清成人免费视频www| 村上凉子中文字幕在线| 国产精品久久久久久av不卡| 午夜精品一区二区三区免费看| 12—13女人毛片做爰片一| 狂野欧美激情性xxxx在线观看| 在线看三级毛片| 欧美绝顶高潮抽搐喷水| 香蕉av资源在线| 日韩亚洲欧美综合| 尾随美女入室| 美女免费视频网站| 一本久久中文字幕| 亚洲无线观看免费| 久久99热这里只有精品18| 亚洲黑人精品在线| 午夜精品一区二区三区免费看| 久久久久久九九精品二区国产| 日韩强制内射视频| 国产精品久久电影中文字幕| 亚洲 国产 在线| av在线天堂中文字幕| 很黄的视频免费| 国产精品亚洲美女久久久| 国产不卡一卡二| 一边摸一边抽搐一进一小说| 真人做人爱边吃奶动态| 成人国产麻豆网| 午夜老司机福利剧场| 人妻丰满熟妇av一区二区三区| 三级国产精品欧美在线观看| 免费看美女性在线毛片视频| 嫁个100分男人电影在线观看| 18禁裸乳无遮挡免费网站照片| 久久人人精品亚洲av| 99热网站在线观看| 国产国拍精品亚洲av在线观看| 俄罗斯特黄特色一大片| 嫩草影院入口| 真人一进一出gif抽搐免费| 色av中文字幕| 午夜福利成人在线免费观看| 美女被艹到高潮喷水动态| 久久久久久久久久成人| 精品一区二区三区人妻视频| 搡老岳熟女国产| 国产一区二区在线av高清观看| 欧美日韩乱码在线| 日韩精品有码人妻一区| 噜噜噜噜噜久久久久久91| 午夜老司机福利剧场| 老司机福利观看| 国产精品不卡视频一区二区| 不卡一级毛片| 国产伦人伦偷精品视频| 蜜桃久久精品国产亚洲av| 国产午夜福利久久久久久| 亚洲精品国产成人久久av| 如何舔出高潮| 日韩欧美 国产精品| 舔av片在线| 免费黄网站久久成人精品| 精品福利观看| 国产真实伦视频高清在线观看 | 麻豆成人午夜福利视频| 精品午夜福利视频在线观看一区| 看黄色毛片网站| 日本成人三级电影网站| 亚洲自拍偷在线| 麻豆国产97在线/欧美| 在现免费观看毛片| 精品久久久久久成人av| 中国美白少妇内射xxxbb| 国产午夜福利久久久久久| 亚洲av.av天堂| 岛国在线免费视频观看| 99热6这里只有精品| 国产精品一区二区免费欧美| 国产精品1区2区在线观看.| 美女免费视频网站| 亚洲午夜理论影院| 黄色丝袜av网址大全| 国产成人一区二区在线| 国内少妇人妻偷人精品xxx网站| 丰满的人妻完整版| 少妇裸体淫交视频免费看高清| 三级毛片av免费| 日韩大尺度精品在线看网址| 成年版毛片免费区| 日韩一本色道免费dvd| 深夜精品福利| a在线观看视频网站| 女的被弄到高潮叫床怎么办 | 波多野结衣高清作品| 久久久久九九精品影院| 在线播放国产精品三级| 两性午夜刺激爽爽歪歪视频在线观看| 国产三级在线视频| 十八禁网站免费在线| 少妇猛男粗大的猛烈进出视频 | 国产高清不卡午夜福利| 国产真实伦视频高清在线观看 | 在线观看舔阴道视频| 我的老师免费观看完整版| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩卡通动漫| 欧美色欧美亚洲另类二区| 国产在线精品亚洲第一网站| 在线天堂最新版资源| 天堂动漫精品| 高清在线国产一区| 亚洲国产精品sss在线观看| 国产亚洲av嫩草精品影院| 国产伦在线观看视频一区| 他把我摸到了高潮在线观看| 韩国av在线不卡| 波野结衣二区三区在线| 嫩草影院精品99| 在线播放国产精品三级| 人人妻人人澡欧美一区二区| 国产亚洲av嫩草精品影院| 日本熟妇午夜| 国产男人的电影天堂91| 欧美国产日韩亚洲一区| 午夜福利18| 俺也久久电影网| 亚洲性久久影院| 亚州av有码| 免费看日本二区| av在线蜜桃| 丝袜美腿在线中文| 国产亚洲精品av在线| 一进一出好大好爽视频| 亚洲avbb在线观看| 波多野结衣高清无吗| 天天躁日日操中文字幕| 韩国av一区二区三区四区| 亚洲综合色惰| 老司机福利观看| 国内少妇人妻偷人精品xxx网站| 免费搜索国产男女视频| 日本熟妇午夜| 亚洲性夜色夜夜综合| 最近最新中文字幕大全电影3| 亚洲国产精品sss在线观看| 欧美色视频一区免费| 欧美不卡视频在线免费观看| 国产黄色小视频在线观看| 大又大粗又爽又黄少妇毛片口| eeuss影院久久| 成人av在线播放网站| 成人国产综合亚洲| 欧美一区二区亚洲| 成人欧美大片| АⅤ资源中文在线天堂| 精品久久久久久久久久久久久| 精品久久久久久,| 国内精品宾馆在线| 亚洲经典国产精华液单| 两个人视频免费观看高清| 天堂av国产一区二区熟女人妻| 自拍偷自拍亚洲精品老妇| 亚洲内射少妇av| 两个人的视频大全免费| 男人舔女人下体高潮全视频| 久久精品综合一区二区三区| 直男gayav资源| 午夜影院日韩av| 亚洲欧美日韩高清专用| 亚洲人成网站在线播放欧美日韩| 岛国在线免费视频观看| 欧美绝顶高潮抽搐喷水| 亚洲国产欧洲综合997久久,| 欧美xxxx性猛交bbbb| 亚洲无线在线观看| 日本一二三区视频观看| 精品久久久久久久久久久久久| 国产免费一级a男人的天堂| av在线天堂中文字幕| 久久精品国产亚洲网站| 精品一区二区三区人妻视频| 国产乱人视频| 最新中文字幕久久久久| 久久6这里有精品| 亚洲第一电影网av| 一进一出抽搐gif免费好疼| 午夜福利在线观看吧| 搞女人的毛片| 最近视频中文字幕2019在线8| 欧美一级a爱片免费观看看| 3wmmmm亚洲av在线观看| 欧美中文日本在线观看视频| 男人舔奶头视频| 桃红色精品国产亚洲av| 女同久久另类99精品国产91| 色播亚洲综合网| 免费av观看视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产欧洲综合997久久,| 国产综合懂色| 高清日韩中文字幕在线| 久久久久久久久久久丰满 | 能在线免费观看的黄片| 免费在线观看影片大全网站| 国产精品嫩草影院av在线观看 | 国国产精品蜜臀av免费| 无人区码免费观看不卡| 日韩在线高清观看一区二区三区 | 欧美精品国产亚洲| 国产精品一区www在线观看 | 精品无人区乱码1区二区| 夜夜看夜夜爽夜夜摸| 99久久无色码亚洲精品果冻| 伦理电影大哥的女人| 嫩草影院新地址| 成人美女网站在线观看视频| 男人舔奶头视频| 欧美日本视频| 亚洲乱码一区二区免费版| 日本爱情动作片www.在线观看 | 一边摸一边抽搐一进一小说| 一卡2卡三卡四卡精品乱码亚洲| 国产白丝娇喘喷水9色精品| 久久中文看片网| 一边摸一边抽搐一进一小说| 在线播放国产精品三级| 三级国产精品欧美在线观看| 超碰av人人做人人爽久久| 老司机福利观看| 一个人看的www免费观看视频| 国产高清激情床上av| 日本欧美国产在线视频| 国产亚洲精品综合一区在线观看| 亚洲乱码一区二区免费版| 精品久久久久久久末码| 国产aⅴ精品一区二区三区波| 91久久精品电影网| 嫩草影院精品99| 老司机福利观看| av在线蜜桃| 亚洲精品国产成人久久av| 成熟少妇高潮喷水视频| 亚洲国产高清在线一区二区三| 波多野结衣高清作品| 国产中年淑女户外野战色| 不卡视频在线观看欧美| 久久亚洲精品不卡| av女优亚洲男人天堂| 亚洲欧美日韩东京热| 色吧在线观看| 日本黄色片子视频| 日韩中文字幕欧美一区二区| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 国产精品久久久久久av不卡| 在线观看免费视频日本深夜| 亚洲欧美清纯卡通| 日本黄色视频三级网站网址| 国产精品嫩草影院av在线观看 | 久久亚洲真实| 色综合色国产| 色播亚洲综合网| 色综合色国产| 在线观看午夜福利视频| 成人特级黄色片久久久久久久| 成人午夜高清在线视频| 欧美精品啪啪一区二区三区| 久久6这里有精品| 免费观看精品视频网站| 久久这里只有精品中国| 久久精品国产鲁丝片午夜精品 | 干丝袜人妻中文字幕| 国产极品精品免费视频能看的| 禁无遮挡网站| 身体一侧抽搐| 精品人妻熟女av久视频| 国产精品综合久久久久久久免费| 国产精品一区www在线观看 | 99久久成人亚洲精品观看| 黄色女人牲交| 观看免费一级毛片| 欧美精品国产亚洲| 免费观看精品视频网站| 丰满乱子伦码专区| 熟妇人妻久久中文字幕3abv| 天美传媒精品一区二区| 成年版毛片免费区| 日韩欧美精品免费久久| 国内精品美女久久久久久| 99在线视频只有这里精品首页| 三级国产精品欧美在线观看| 亚洲人与动物交配视频| 2021天堂中文幕一二区在线观| 国产美女午夜福利| 自拍偷自拍亚洲精品老妇| 色视频www国产| 亚洲自拍偷在线| 亚洲精品456在线播放app | 最近最新中文字幕大全电影3| 美女黄网站色视频| 亚洲三级黄色毛片| 久久久久久大精品| 欧美激情国产日韩精品一区| 99国产精品一区二区蜜桃av| av福利片在线观看| 麻豆国产av国片精品| 身体一侧抽搐| 国产日本99.免费观看| 午夜福利在线在线| 国产v大片淫在线免费观看| 成人国产麻豆网| 国产av一区在线观看免费| 亚州av有码| 精品一区二区三区av网在线观看| 亚洲精品色激情综合| 99热这里只有是精品50| 亚洲熟妇熟女久久| 午夜免费男女啪啪视频观看 | 蜜桃亚洲精品一区二区三区| 亚洲美女视频黄频| 观看免费一级毛片| 国产一区二区在线av高清观看| 看片在线看免费视频| 欧美性猛交黑人性爽| 日韩欧美在线乱码| 日本爱情动作片www.在线观看 | www.色视频.com| 黄色丝袜av网址大全| 日韩欧美三级三区| 99久久精品一区二区三区| 狂野欧美激情性xxxx在线观看| 老司机深夜福利视频在线观看| 日本爱情动作片www.在线观看 | 精品乱码久久久久久99久播| 国产欧美日韩精品一区二区| 99九九线精品视频在线观看视频| 免费看光身美女| 日本爱情动作片www.在线观看 | 亚洲图色成人| 三级男女做爰猛烈吃奶摸视频| 色吧在线观看| 99精品久久久久人妻精品| 欧美日本视频| 久久久成人免费电影| 很黄的视频免费| 国产免费一级a男人的天堂| 国模一区二区三区四区视频| 一个人看视频在线观看www免费| 久久久久久久午夜电影| 成人美女网站在线观看视频| 国产精品av视频在线免费观看| 国产精品久久久久久久久免| 一级毛片久久久久久久久女| 91在线观看av| 狂野欧美激情性xxxx在线观看| 特级一级黄色大片| 中国美女看黄片| 老熟妇仑乱视频hdxx| 国产视频内射| 九九久久精品国产亚洲av麻豆| 亚洲在线自拍视频| 欧美一区二区亚洲| 婷婷精品国产亚洲av| 亚洲欧美日韩高清专用| 精品久久久久久,| 久久精品国产亚洲av天美| 国产欧美日韩精品亚洲av| 日韩在线高清观看一区二区三区 | 成人精品一区二区免费| 极品教师在线免费播放| h日本视频在线播放| 久久午夜亚洲精品久久| 午夜免费成人在线视频| 国产单亲对白刺激| 自拍偷自拍亚洲精品老妇| 亚洲国产精品久久男人天堂| 亚洲av美国av| 亚洲性久久影院| videossex国产| 精品人妻偷拍中文字幕| 特大巨黑吊av在线直播| 日本精品一区二区三区蜜桃| 亚洲,欧美,日韩| 欧美xxxx黑人xx丫x性爽| 久久久久久九九精品二区国产| 免费看av在线观看网站| 他把我摸到了高潮在线观看| a级一级毛片免费在线观看| 欧美日韩乱码在线| 亚洲人成网站在线播放欧美日韩| 少妇人妻精品综合一区二区 | 欧美最新免费一区二区三区| 男人的好看免费观看在线视频| 一区二区三区高清视频在线| 一级av片app| 久久精品夜夜夜夜夜久久蜜豆| 老女人水多毛片| 乱人视频在线观看| 在线观看66精品国产| 乱系列少妇在线播放| 男人舔奶头视频| 亚洲精品乱码久久久v下载方式| 免费不卡的大黄色大毛片视频在线观看 | 天天一区二区日本电影三级| 春色校园在线视频观看| 波多野结衣高清无吗| 亚洲久久久久久中文字幕| 欧美精品国产亚洲| 成人永久免费在线观看视频| 一区二区三区免费毛片| 嫩草影视91久久| 中文字幕免费在线视频6| 亚洲熟妇中文字幕五十中出| 国产精品乱码一区二三区的特点| 亚洲最大成人手机在线| 97人妻精品一区二区三区麻豆| 一本精品99久久精品77| 淫秽高清视频在线观看| 国产老妇女一区| 熟女人妻精品中文字幕| 日本 av在线| 午夜视频国产福利| 一个人观看的视频www高清免费观看| ponron亚洲| 午夜福利成人在线免费观看| 国产精品三级大全| 少妇熟女aⅴ在线视频| 亚洲欧美清纯卡通| 国产三级在线视频| av在线老鸭窝| a级毛片免费高清观看在线播放| 婷婷精品国产亚洲av| 国产精品美女特级片免费视频播放器| 国产精品久久久久久亚洲av鲁大| 直男gayav资源| 亚洲无线在线观看| 欧美区成人在线视频| 成年免费大片在线观看| 美女黄网站色视频| 一个人看视频在线观看www免费| 国产精华一区二区三区| 亚洲成人中文字幕在线播放| 日韩大尺度精品在线看网址| 少妇裸体淫交视频免费看高清| 永久网站在线| 蜜桃亚洲精品一区二区三区| 变态另类丝袜制服| 看免费成人av毛片| 简卡轻食公司| 国产免费av片在线观看野外av| 日日摸夜夜添夜夜添小说| 欧美成人a在线观看| 亚洲精品在线观看二区| 一本精品99久久精品77| 天美传媒精品一区二区| 国产高清视频在线播放一区| 欧美极品一区二区三区四区| a级毛片免费高清观看在线播放| 黄色日韩在线| 成年免费大片在线观看| 深夜精品福利| 婷婷亚洲欧美| 欧美色视频一区免费| 精品国产三级普通话版| 中文字幕av在线有码专区| 日本爱情动作片www.在线观看 | 国内少妇人妻偷人精品xxx网站| 成年女人毛片免费观看观看9| 热99在线观看视频| 黄色欧美视频在线观看| 我要搜黄色片| 成人毛片a级毛片在线播放| 国产精品永久免费网站| 国产精品自产拍在线观看55亚洲| 永久网站在线| 久久草成人影院| 天堂动漫精品| 少妇高潮的动态图| 亚洲欧美激情综合另类| 免费不卡的大黄色大毛片视频在线观看 | 日本熟妇午夜| 成人av在线播放网站| 免费av毛片视频| 国产精品亚洲美女久久久| 久久久色成人| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕av在线有码专区| a在线观看视频网站| 精品人妻1区二区| 国产大屁股一区二区在线视频| 少妇裸体淫交视频免费看高清| 欧美一区二区精品小视频在线| 国内久久婷婷六月综合欲色啪| 久久久久久久久久黄片| 免费人成在线观看视频色| 制服丝袜大香蕉在线| 一边摸一边抽搐一进一小说| 搡老熟女国产l中国老女人| 我要看日韩黄色一级片| 国产午夜精品论理片| 欧美在线一区亚洲| 成人高潮视频无遮挡免费网站| 性插视频无遮挡在线免费观看|