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

    哈爾濱西泉眼水庫夏季浮游植物群落結(jié)構(gòu)動態(tài)特征*

    2015-09-25 10:18:56王建國于洪賢馬成學姚允龍汪義杰馬金龍唐紅亮劉壯添
    湖泊科學 2015年4期
    關(guān)鍵詞:泉眼綠藻藍藻

    王建國,于洪賢,馬成學,姚允龍,汪義杰,馬金龍,唐紅亮,劉壯添,李 麗

    (1:珠江水利委員會珠江水利科學研究院,廣州510611)

    (2:東北林業(yè)大學野生動物資源學院,哈爾濱150040)

    水庫及其流域水系是內(nèi)陸水體的重要組成部分,也是人類生產(chǎn)、生活和社會發(fā)展需要的重要基礎.隨著水庫流域經(jīng)濟的發(fā)展,其接受的污染負荷日益增加,從而導致其水環(huán)境質(zhì)量惡化,水體富營養(yǎng)化加速,已成為制約我國社會經(jīng)濟可持續(xù)發(fā)展的限制性因素[1].據(jù)2009年《中國環(huán)境狀況公報》數(shù)據(jù),在我國26個國控重點湖泊(水庫)中,IV類及以下水質(zhì)有20個.盡管大部分湖庫已經(jīng)嚴重富營養(yǎng)化并部分暴發(fā)藍藻水華,但仍不得不作為飲用水源地[2],從而對居民飲水安全埋下隱患.浮游植物是生態(tài)系統(tǒng)中的初級生產(chǎn)者[3],其種類組成和群落結(jié)構(gòu)的動態(tài)變化將直接影響水生態(tài)系統(tǒng)的結(jié)構(gòu)與功能,對維持淡水生態(tài)系統(tǒng)平衡起到至關(guān)重要的作用[4-5].已有研究發(fā)現(xiàn),營養(yǎng)鹽、水溫、光照強度及浮游動物的捕食是影響浮游植物群落結(jié)構(gòu)的主要因素[6-7].氮或者磷通常成為淡水、河口和沿海生態(tài)系統(tǒng)中生物初級生產(chǎn)力的限制因子[8],水溫對太湖浮游植物組成起到重要作用[2],但目前對我國東北地區(qū)的水庫及河流藻類研究較少,尤其是哈爾濱這一北部高寒區(qū)與其他地區(qū)顯著不同,如溫度低、降水少、營養(yǎng)鹽遠低于我國南方平原庫區(qū).2005年松花江水污染事件后,西泉眼水庫作為哈爾濱市第二飲用水源地,對其研究不斷增多.如王越等[9]、遲晉峰等[10]對西泉眼水庫的營養(yǎng)鹽來源、水環(huán)境現(xiàn)狀進行初步分析,陳翠翠等[11]對西泉眼水庫消落帶的水生植物群落結(jié)構(gòu)進行研究,高孜娟等[12]對西泉眼水庫魚類多樣性與群落結(jié)構(gòu)進行初步分析,呂東珂[13]對西泉眼水庫這一泥炭型水庫水-氣界面CO2通量進行研究.然而,已有研究往往缺少對浮游植物種群動態(tài)與水環(huán)境因子更深層次的分析與探討,且未對水庫和來源入庫河流水體進行對比性分析研究.因此,本文對西泉眼水庫夏季浮游植物群落特征與水環(huán)境因子的主要驅(qū)動因子動態(tài)變化關(guān)系進行研究,可為我國東北地區(qū)重要水源地的保護和生態(tài)治理提供基礎依據(jù).

    1 材料與方法

    1.1 調(diào)查區(qū)域概況

    西泉眼水庫(45°12'~45°20'N,127°17'~127°25'E)位于哈爾濱市境內(nèi),于 1996 年竣工,大(Ⅱ)型水庫,常水位庫容為4.78×108m3,平均水深約為12 m.該水庫是攔蓄松花江一級支流阿什河而成,并另有二道河和黃泥河匯入.西泉眼水庫控制流域面積為1146 km2,占全流域面積的32%;流域內(nèi)生長著大片次生林,土壤為黑土,周邊森林覆蓋率高.該區(qū)域?qū)儆诖箨懶詼貛Ъ撅L氣候,夏季高溫、多雨,多年平均降雨量為563 mm,其中6--9月份降水占全年降水量的70%.

    1.2 調(diào)查方法

    1.2.1 樣點位置 根據(jù)西泉眼水庫庫盆形狀及流域地理特征,本研究共布設13個采樣點(圖1),具體如下:1#、2#和3#(阿什河入庫口左、中、右);4#(黃泥河入庫口);5#、6#和7#(庫心左、中、右);8#(雙龜山);9#(大壩出水口);10#(平山鐵路橋);11#(平山西橋);12#(平山南橋);13#(阿什河二道河匯流處).同時,采樣點用Garmin公司生產(chǎn)的GPS12型全球衛(wèi)星定位系統(tǒng)定位.采樣時間為2010年7、8月.

    1.2.2 浮游植物樣品采集、保存及鑒定 定性樣品采用25#浮游生物網(wǎng)在水面下0.5 m處進行“∞”字形撈取.定量樣品采用5 L采水器根據(jù)水層深度分表層(水下0.5 m)、中層、底層(距底0.5 m)采集混合水樣,所有樣品均用魯哥試劑固定,靜置24 h后虹吸,再靜置48 h后濃縮至30 ml.吸取0.1 ml水樣置于0.1 ml浮游生物計數(shù)框內(nèi)用Motic BA400顯微鏡計數(shù)、分析、鑒定[14].

    圖1 西泉眼水庫浮游植物采樣點Fig.1 Sampling sites of phytoplankton in Xiquanyan Reservoir

    1.2.3 水環(huán)境理化指標測定 水溫、pH、溶解氧(DO)、電導率(COND)等指標用美國YSI-6200多功能水質(zhì)分析儀于現(xiàn)場測定;透明度用塞氏盤法測定;總氮(TN)、總磷(TP)、硝態(tài)氮(-N)、氨態(tài)氮(NH3-N)、葉綠素a(Chl.a)、高錳酸鹽指數(shù)(CODMn)和Cl-均參照文獻[15]進行測定.

    1.3 數(shù)據(jù)分析

    浮游植物豐度、生物量的計算參照《淡水浮游生物研究方法》[16];按照Lampitt等提出的優(yōu)勢度指數(shù)(y),以優(yōu)勢度 y>0.02來確定優(yōu)勢種[17];并采用 Shannon-Wiener多樣性指數(shù)(H')[18]對浮游植物生物多樣性進行分析.典范對應分析(CCA)是基于種類數(shù)據(jù)和環(huán)境數(shù)據(jù)集合的排序方法[19-20],本文在對環(huán)境因子進行正態(tài)分布檢驗后,應用CCA對西泉眼水庫浮游植物群落與水環(huán)境因子進行相關(guān)分析.

    數(shù)據(jù)統(tǒng)計分析及制圖運用OriginPro 8.1軟件.利用CANOCO 4.5軟件包進行CCA分析,采用浮游植物豐度指標來反映浮游植物種群與水環(huán)境的關(guān)系,為浮游植物豐度的數(shù)據(jù)獲得正態(tài)分布,對浮游植物豐度值進行l(wèi)g(x+1)轉(zhuǎn)換.

    2 結(jié)果與分析

    2.1 西泉眼水庫水體理化參數(shù)

    研究期間,西泉眼水庫水深在0.3~16.5 m之間,庫區(qū)顯著高于入庫河流區(qū);水溫在18.65~26.25℃之間,庫區(qū)略高于入庫河流區(qū);透明度在0.30~1.35 m之間,濁度在5.1~34.4 FTU之間,與透明度呈相反的變化趨勢;pH值在7.4~9.92之間,庫區(qū)略高于入庫河流區(qū);DO濃度在4.75~7.53 mg/L之間,入庫河流區(qū)明顯高于庫區(qū)水體;COND在0.002~0.204 μS/cm之間,入庫河流區(qū)明顯高于庫區(qū)水體;Cl-濃度在2.7~13.5 mg/L 之間;-N濃度在0.21~1.67 mg/L之間,入庫河流區(qū)明顯高于庫區(qū)水體;CODMn在3.0~6.4 mg/L之間,庫區(qū)略高于入庫河流(表1).總體來看,西泉眼水庫庫區(qū)與入庫河流區(qū)的理化指標之間差異顯著(P<0.05,n=10),夏季水溫的逐漸上升也比較適于大部分浮游植物種類的快速生長,尤其是綠藻門和藍藻門種類.

    表1 西泉眼水庫夏季水環(huán)境因子時空分布特征Tab.1 The spatial-temporal disrtibution of environmental variables in the summer in Xiquanyan Reservoir

    同時,根據(jù)GB 3838 -2002評價標準,西泉眼水庫的NH3-N總體優(yōu)于Ⅲ類水質(zhì)標準;TP和TN濃度在大多數(shù)監(jiān)測點位超出了Ⅲ類水質(zhì)標準,尤其在7月(除2#、9#等個別點位外),均明顯高于8月份(圖2),表明夏季水庫庫區(qū)及入庫河流區(qū)的營養(yǎng)鹽濃度較高,為浮游植物的快速生長提供了營養(yǎng)源.此外,根據(jù)Chl.a富營養(yǎng)化評價標準[21],8月Chl.a濃度明顯高于7月,既表明水體已呈現(xiàn)富營養(yǎng)化趨勢,也表明浮游植物豐度對營養(yǎng)物質(zhì)的響應呈一定的滯后性.總體來看,7月和8月的水環(huán)境因子中Chl.a濃度差異顯著(P<0.05,n=10),其余各指標整體差異不顯著(P >0.05,n=10).

    2.2 浮游植物種類組成

    經(jīng)鏡檢,共鑒定出浮游植物7門169種(含變種及變型),其中綠藻門最多,83種,占總種數(shù)的49.11%;硅藻門次之,47種,占總種數(shù)的27.81%;藍藻門15種,占總種數(shù)的8.88%,而裸藻門(12種,7.10%)和金藻門(6種,3.6%)相對較少,甲藻門(4種,2.4%)和隱藻門(2種,1.2%)則最低.結(jié)果表明,西泉眼水庫夏季浮游植物群落以綠藻-硅藻-藍藻型為主.

    水庫庫區(qū)的浮游植物群落結(jié)構(gòu)為綠藻-硅藻-藍藻型,而入庫河流則為硅藻-綠藻-藍藻型;庫區(qū)綠藻和藍藻種數(shù)分別為入庫河流區(qū)的3倍多,而入庫河流區(qū)的硅藻種數(shù)則是庫區(qū)的近2倍(圖3).對比研究發(fā)現(xiàn),西泉眼水庫庫區(qū)與入庫河流區(qū)浮游植物群落結(jié)構(gòu)的差異顯著(P<0.05,n=13).分析認為,由于水庫綠藻和藍藻種數(shù)遠多于入庫河流區(qū),說明水庫庫區(qū)暴發(fā)綠藻、藍藻水華的潛在風險較大,尤其8月份庫區(qū)各采樣點的Chl.a濃度均超出10 μg/L;而入庫河流區(qū)的天然流水生境比庫區(qū)緩流水體更適宜硅藻生存.

    根據(jù)浮游植物在各采樣點出現(xiàn)的頻率與豐度,夏季7月和8月共發(fā)現(xiàn)優(yōu)勢種10種,形成以硅藻門、綠藻門和藍藻門占優(yōu)勢地位的群落結(jié)構(gòu)(表2).而且,從浮游植物優(yōu)勢種的污染指示等級來看,西泉眼水庫已處于中度富營養(yǎng)水平,并呈現(xiàn)出由中度富營養(yǎng)型水體向重度富營養(yǎng)型過渡的趨勢.

    圖2 西泉眼水庫主要營養(yǎng)鹽及Chl.a濃度時空變化特征Fig.2 The spatial-temporal variations of the main nutrients and chlorophyll-a concentration in Xiquanyan Reservoir

    另外,值得注意的是,在夏季7月,藍藻門中的阿氏席藻豐度低于肘狀針桿藻和尖針桿藻,為第三優(yōu)勢種;在夏季8月雖然綠藻門優(yōu)勢種類相對較少,但是綠藻門中的游絲藻豐度卻最高,為第一優(yōu)勢種,阿氏席藻仍為第三優(yōu)勢種.研究表明,西泉眼水庫藍藻種群具備良好的種質(zhì)資源,一旦水溫適宜、營養(yǎng)物質(zhì)充足,藍藻種群的某些種類可能會急劇生長,甚至暴發(fā)水華.

    圖3 西泉眼水庫夏季庫內(nèi)與入庫河流區(qū)浮游植物定量種類組成對比Fig.3 The species of phytoplankton in reservoir area compared with inflow-river area in Xiquanyan Reservoir

    2.3 浮游植物豐度、生物量及群落多樣性

    西泉眼水庫庫區(qū)浮游植物的豐度和生物量從阿什河入庫口(1#、2#和3#)和黃泥河入庫口(4#)向庫心處(5#、6#和7#)、大壩口處(8#和9#)呈遞減趨勢;入庫河流區(qū)各樣點浮游植物豐度和生物量變化波動不大,差異不明顯;同時,庫區(qū)內(nèi)各采樣點浮游植物的豐度和生物量均大于入庫河流水體各采樣點;而且,阿什河入庫口(1#、2#和3#)和黃泥河入庫口(4#)浮游植物的豐度和生物量遠高于其他采樣點(圖4).此外,8月5#采樣點的Chl.a濃度略低于6#采樣點和4#采樣點,但其生物量最高;分析認為,這與8月5#采樣點浮游植物定量樣品中采集到數(shù)量較多且個體較大的尖針桿藻和某種游絲藻有關(guān),而其個體體積對生物量的貢獻量高過其對 Chl.a濃度的影響.

    表2 西泉眼水庫夏季浮游植物優(yōu)勢種及其豐度、污染指示等級*Tab.2 Dominant species and abundance,environmental indicator class of phytoplankton in summer in Xiquanyan Reservoir

    總體上看,西泉眼水庫庫區(qū)浮游植物豐度和生物量在7月和8月變化趨勢高度擬合.分析認為,在水庫庫區(qū)這一緩流水體中更適于綠藻和藍藻種類的生長,而入庫河流區(qū)的天然流水則對硅藻尤其對著生硅藻更為適宜;而河流與水庫交匯處,則由于營養(yǎng)鹽充足且水體流態(tài)多樣,從而為藻類大量繁衍提供了良好的生境條件.

    西泉眼水庫及其入庫河流浮游植物Shannon-Wiener指數(shù)的空間異質(zhì)性分布變化特征見圖5.根據(jù)評價標準[21],7月僅有5#采樣點的H'為2.73,低于中污染水平基準線,其余各采樣點位均高于此標準;8月5#~9#采樣點位的H'均低于中污染基準線,其余采樣點位高于此標準.同時,入庫河流浮游植物的Shannon-Wiener指數(shù)均高于庫區(qū).結(jié)果表明,入庫天然河流的水質(zhì)整體上要優(yōu)于水庫緩流水體,且水庫庫區(qū)的水質(zhì)從河流入庫口處到庫心區(qū),再到大壩口處,水質(zhì)呈一定的惡化趨勢.

    圖4 西泉眼水庫及入庫河流夏季浮游植物豐度(a)和生物量(b)的空間分布Fig.4 The spatial distribution of phytoplankton's abundance(a)and biomass(b)in Xiquanyan Reservoir and the inflow-rivers

    2.4 浮游植物與環(huán)境因子的關(guān)系

    對西泉眼水庫庫區(qū)6個典型樣點和入庫河流區(qū)的4個樣點,按出現(xiàn)頻度及其豐度進行選擇,并將夏季7、8月環(huán)境因子數(shù)據(jù)和浮游植物豐度合并分析:庫區(qū)種類:頻次≥4;入庫河流區(qū)種類:頻次≥3;在豐度上保證所選種類豐度占整個水庫浮游植物豐度的95%以上,以確保所選主要浮游植物種類構(gòu)成及豐度與西泉眼水庫整體狀況基本一致(表 3).

    參照董旭輝等[22]的研究成果,對包括所有環(huán)境變量的數(shù)據(jù)庫進行初步CCA分析,對膨脹因子(VIF)大于20的環(huán)境變量進行選擇性刪除.通過CCA表明,環(huán)境變量中的水深、水溫、透明度、pH值、濁度5個變量的VIF均大于20,予以刪除.

    西泉眼水庫浮游植物群落CCA分析的統(tǒng)計信息(表4)可以看出,水庫庫區(qū)的總變量為0.3761,可解釋擬合變量占85.64%(可解釋變量占76.77%);排列測試結(jié)果顯示所有軸pseudo-F=3.2,P=0.002;入庫河流區(qū)的總變量為0.3071,可解釋擬合變量占81.60%(可解釋變量占74.25%);排列測試結(jié)果:所有軸pseudo-F=4.3,P=0.005.總體來看,西泉眼水庫庫區(qū)和入庫河流區(qū)軸1的特征值均明顯大于其它各軸的特征值,能夠更好地描述物種與環(huán)境因子的關(guān)系;同時,軸1貢獻的可解釋變量的累計百分比、可解釋擬合變量的累計百分比均比其余各軸高.因此CCA分析可以很好地反映浮游植物群落與環(huán)境因子之間的相關(guān)關(guān)系.

    CCA分析結(jié)果(圖6A,B)表明,浮游植物物種與環(huán)境因子的排序,大體反映了該物種的綜合環(huán)境生態(tài)位和生態(tài)學特征.即排序圖中空間位置聚集在一起的物種,其對水環(huán)境因子等生境要求較為相近;而空間位置分散的物種,說明物種之間對生境的要求差異較大.

    CCA排序圖6A和圖6B顯示,在西泉眼水庫庫區(qū)和入庫河流區(qū)均明顯地將夏季7月和8月劃分開來,其中庫區(qū)以軸2為界,入庫河流區(qū)以軸1為界,表明相同區(qū)域相同季節(jié)內(nèi)不同月份也呈現(xiàn)明顯的時間異質(zhì)性.同時,還可以看出,庫區(qū)樣點河流入庫口處(2#和4#)、庫心(6#和8#)、大壩口處(9#)空間異質(zhì)性分布規(guī)律明顯;入庫河流樣點上游(10#和11#)、中游(12#)和下游(13#)空間異質(zhì)性分布顯著.

    此外,從庫區(qū)的浮游植物種類組成來看(圖6A),藍藻門和部分綠藻門、硅藻門的種類主要分布在軸2右側(cè),而部分綠藻門、金藻門、裸藻門和隱藻門的藻類主要分布在軸1左側(cè),尤其藍藻門種類全部集中分布在軸2右側(cè),而左側(cè)的綠藻門種類明顯多于右側(cè).同時,軸1與-N、NH3-N、TN濃度和CODMn呈顯著正相關(guān),表明庫區(qū)氮類營養(yǎng)鹽和CODMn是藍藻門種類和部分喜氮綠藻門的主要環(huán)境變量;而與DO、電導率則呈顯著負相關(guān),也表明庫區(qū)充足的DO、電導率也是大部分綠藻門和金藻門種類的主要環(huán)境變量.從入庫河流區(qū)的浮游植物種類組成來看(圖6B),金藻門和大部分硅藻門種類主要分布在軸1上方,表明入庫河流中充足的DO是金藻和大部分硅藻的主要環(huán)境變量;而絕大多數(shù)綠藻門種類、部分硅藻門種類,以及某些裸藻門種類主要分布在軸1下方,表明入庫河流中氮類營養(yǎng)鹽和CODMn是其主要環(huán)境變量.總體來看,庫區(qū)浮游植物種類組成及豐度主要與NH3-N等氮類營養(yǎng)鹽濃度顯著相關(guān),卻幾乎不受TP濃度的影響;而入庫河流區(qū)浮游植物種類組成及豐度則受氮、磷等營養(yǎng)鹽物的雙重影響.

    圖5 西泉眼水庫及入庫河流浮游植物Shannon-Wiener指數(shù)變化趨勢Fig.5 The distribution trend of Shannon-Wiener index of phytoplankton in Xiquanyan Reservoir and the inflow-rivers

    表3 西泉眼水庫及入庫河流夏季浮游植物主要種類組成及編碼Tab.3 Phytoplankton species composition and codes of Xiquanyan Reservoir and the inflow-rivers in summer

    表4 CCA排序圖各軸的特征值及種類與環(huán)境排序軸間的統(tǒng)計量Tab.4 Eigen values and correlation coefficients of species and environmental factors on the axes of CCA

    圖6 西泉眼水庫夏季浮游植物豐度與環(huán)境因子的CCA三維排序圖(A:庫區(qū);B:入庫河流)(① Δ:物種;О:樣點;→:環(huán)境因子.②浮游植物編碼見表3.③采樣斷面的編碼為“樣點編號+采樣月份”,如“2#-7”表示“2#斷面7月份樣品”)Fig.6 CCA triplots between phytoplankton abundance and environmental factors in summer in Xiquanyan Reservoir(A:reservoir area;B:inflow-rivers)

    3 討論

    3.1 浮游植物種類組成

    西泉眼水庫作為大型(Ⅱ)水利樞紐工程,為半天然半人工的水體,近年來水環(huán)境質(zhì)量日益惡化,導致水體富營養(yǎng)化加速;而入庫河流為天然水體,水環(huán)境質(zhì)量相對較好.為此本研究鑒定出浮游植物種數(shù)要明顯高于單一類型水體(圖3).同時,研究表明,西泉眼水庫庫區(qū)以綠藻門、硅藻門和藍藻門為主.然而,庫區(qū)藍藻門種數(shù)要明顯高于入庫河流區(qū),一方面表明西泉眼水庫庫區(qū)水質(zhì)由于流動性較低而比入庫河流區(qū)差,即從入庫河流區(qū)到庫區(qū)水質(zhì)有惡化的趨勢;另一方面也表明水庫這種靜止水體藍、綠藻水華暴發(fā)的潛在幾率要高于天然流動水體.西泉眼水庫已有研究[23]也表明,庫區(qū)浮游植物種類以綠藻門、藍藻門和硅藻門為主,而且夏季豐水期種類最多.鏡泊湖[24]、海浪河[25]浮游植物的調(diào)查研究發(fā)現(xiàn),其浮游植物種類組成同樣以綠藻門、硅藻門和藍藻門為主.對比來看,即使同處于高緯度地區(qū),藻類組成相似,但在天然湖泊與人工湖泊(水庫)之間、天然河流之間的浮游植物種數(shù)具有區(qū)域性的差異.

    西泉眼水庫的優(yōu)勢種群主要是硅藻門、綠藻門和藍藻門的種類,金藻門和隱藻門的個別種類也為優(yōu)勢種.其中夏季8月的游絲藻豐度、生物量和優(yōu)勢度均最高,屬于第一優(yōu)勢種,而取樣時TP濃度在30~110 mg/m3之間,均值為60 mg/m3.Jensen等對比了丹麥境內(nèi)不同營養(yǎng)水平的湖泊,總結(jié)出當湖泊中TP濃度超出1000 mg/m3時,藍藻通常被綠藻所取代[26].而在我國太湖,當TP濃度在150~200 mg/m3時,絲狀綠藻中的某種游絲藻在藻類暴發(fā)時占優(yōu)勢[2].Scheffer研究認為,這種現(xiàn)象有可能是依賴于湖泊的自然環(huán)境、緯度以及浮游植物種群的“內(nèi)穩(wěn)態(tài)”機制[27].然而,本研究中游絲藻急劇增加時,其TP濃度并不高;通過CCA排序圖(圖6A)也表明,TP濃度與游絲藻無顯著相關(guān)性.分析認為,某種游絲藻占明顯優(yōu)勢,不僅要分析TP濃度的影響,還應將其納入到浮游植物群落層次進行探討.

    3.2 浮游植物豐度與生物量的變化

    西泉眼水庫夏季浮游植物豐度和生物量的時空分布特征變化趨勢較為一致(圖4).在溫帶地區(qū),水溫是引起水體浮游生物季節(jié)性變化的主要因素之一[28].當夏季7月水溫最高時,浮游植物藍藻的種類數(shù)量、豐度、生物量達到峰值,尤其在水庫庫區(qū)各樣點.同時,夏季浮游植物豐度和生物量大量增加,與營養(yǎng)鹽的增加有密切聯(lián)系.而本文選取夏季作為典型季節(jié),其各月的水體溫差相對較小,但在CCA排序圖中仍將7月(均值為24.16℃)和8月(均值為21.88℃)的采樣斷面明顯分開,表明水溫對浮游植物群落結(jié)構(gòu)在同一季節(jié)內(nèi)仍有明顯的影響.同時,西泉眼水庫7月水溫高于8月,其藍藻豐度和生物量也明顯高于8月,在一定程度上表明,水溫與藍藻豐度及其生物量存在顯著正相關(guān)關(guān)系.分析認為,一方面由于采樣月份的水溫均較高,另一方面,采樣月份西泉眼水庫屬于豐水期,營養(yǎng)鹽來源相對較多,面源污染物可能有所增加.

    水環(huán)境質(zhì)量決定生物種群和群落結(jié)構(gòu)特征,而生物種群和群落也對水環(huán)境質(zhì)量的變化及時作出響應,可以客觀反映水體質(zhì)量的變化.實際應用中,Shannon-Wiener指數(shù)是使用最多的多樣性指數(shù)[29].生物多樣性分析表明,西泉眼水庫庫區(qū)夏季屬于中度富營養(yǎng)型水體,入庫河流屬于輕度富營養(yǎng)型水體.分析認為,入庫河流區(qū)流經(jīng)低山丘陵區(qū),森林覆蓋率高達80%,加之常年流動,水質(zhì)相對較好;而庫區(qū)水體由于流動性差,加上建庫時淹沒的大量農(nóng)田和林地,庫底存在大量的營養(yǎng)鹽類沉積物,在水溫較高和豐水期沖刷的干擾下,沉積物中營養(yǎng)鹽上泛,從而導致其水質(zhì)下降.

    3.3 浮游植物與環(huán)境因子的關(guān)系

    研究表明,浮游植物與水環(huán)境因子之間的關(guān)系非常密切[30].浮游植物的生長需要在一定的光照和水溫下進行,在溫帶和寒帶地區(qū),水溫通常是限制因子,影響著浮游植物的生長[3,31].Komarkova等通過CCA分析發(fā)現(xiàn),微囊藻與較高的水溫和高Na+濃度有較強的相關(guān)性,而脆桿藻和甲藻則與夏季水溫上升有關(guān)[32].Chen等對太湖的研究表明,短時期浮游植物消長基本上受水溫、風和濁度的影響,而長期浮游植物生物量的動態(tài)變化則受到營養(yǎng)水平的影響[2].在亞熱帶地區(qū),水溫和營養(yǎng)鹽負荷是影響大鏡山水庫浮游植物群落結(jié)構(gòu)動態(tài)變化的主要驅(qū)動因子[33],而降雨、水溫及水力滯留時間則是影響富春江水庫浮游植物群落結(jié)構(gòu)及密度變化的主要驅(qū)動因子[34].本研究表明,西泉眼水庫屬于溫帶地區(qū),夏季水溫顯著增高,其中7月水溫均值為24.16℃、8月為21.88℃,滿足浮游植物藍藻門和綠藻門種類快速繁殖的溫度條件,水溫與藻類種數(shù)、豐度呈一定的正相關(guān).此外,研究發(fā)現(xiàn),西泉眼水庫庫區(qū)多數(shù)綠藻門和藍藻門種類與NH3-N等氮類營養(yǎng)鹽濃度呈顯著的正相關(guān)關(guān)系,而入庫河流區(qū)的多數(shù)綠藻門、隱藻門和硅藻門種類與NH3-N等氮類營養(yǎng)鹽、TP濃度呈顯著正相關(guān)關(guān)系,表明NH3-N、-N、TN和TP濃度是影響西泉眼水庫浮游植物群落結(jié)構(gòu)及豐度的主要驅(qū)動因子.

    同時,Smith研究表明,當?shù)粗蠺N/TP<29時,藍藻就出現(xiàn)急劇暴發(fā)水華的趨勢;而當TN/TP>29時,藍、綠藻豐度就出現(xiàn)減少的趨勢[35].本研究表明,西泉眼水庫夏季TN/TP平均為25.37(其中,7月為26.84、8月為24.25),加上充足的氮源和磷源供給,更有利于藍藻和綠藻門浮游植物的急劇生長,在一定程度上導致藍藻和綠藻豐度急劇升高.而且,根據(jù)浮游植物CCA排序圖監(jiān)測斷面分析(圖6),總體上將西泉眼水庫分為庫區(qū)生態(tài)系統(tǒng)和入庫河流生態(tài)系統(tǒng)兩個部分,又進一步將7月和8月的時間差異性凸顯,符合西泉眼水庫生態(tài)環(huán)境現(xiàn)狀,又能夠更好地解釋兩種生態(tài)系統(tǒng)各自的特征.

    4 結(jié)論及建議

    1)經(jīng)調(diào)查分析,西泉眼水庫夏季浮游植物的豐度和生物量相對較高.豐度在10.92×105~88.44×105cells/L間波動,以夏季8月的3#采樣點最高;生物量在2.61~22.44 mg/L之間波動,以夏季8月5#采樣點最高.Shannon-Wiener指數(shù)在2.46~4.69之間,以夏季7月12#采樣點最高.夏季以梅尼小環(huán)藻、肘狀針桿藻、阿氏席藻和游絲藻為主要優(yōu)勢藻類,卵形隱藻、棕鞭藻和固氮魚腥藻為次優(yōu)勢藻.

    2)運用CCA分析得出,NH3-N、-N、TP和TN等營養(yǎng)鹽濃度是影響浮游植物的主要環(huán)境變量,其次是水溫、水深、透明度和濁度等因子.同時,總體上將西泉眼水庫分為水庫生態(tài)系統(tǒng)和河流生態(tài)系統(tǒng)兩部分,以及7月和8月兩個采樣時間段,符合實際狀況.而且,在庫區(qū)樣點河流入庫口處、庫心、大壩口處空間異質(zhì)性分布規(guī)律明顯;入庫河流樣點來源、中游和下游空間異質(zhì)性分布顯著.

    3)西泉眼水庫庫區(qū)水體處于中度富營養(yǎng)型向重度富營養(yǎng)型水體過渡階段.為此,建議水庫管理部門建立水生態(tài)系統(tǒng)健康監(jiān)測與評估預警體系和相應的水生態(tài)管理制度,并采取適宜的生態(tài)修復技術(shù)等生態(tài)管理措施,以確保哈爾濱市供水水源安全.

    致謝:感謝任磊、鞠永富、梁雄偉、費滕、高孜娟、李佳娟、陳翠翠等同學和西泉眼水庫管理處在采樣及制圖方面提供的幫助.

    [1]劉建康.東湖生態(tài)學研究(二).北京:科學出版社,1995.

    [2]Chen YW,Qin BQ,Teubner K et al.Long-term dynamics of phytoplankton assemblages:Microcystis-domination in Lake Taihu,a large shallow lake in China.Journal of Plankton Research,2003,25(4):445-453.

    [3]Reynolds T,Cliff N.An interactive preference ordering model and its monte carlo evaluation.Psychometrika,1984,49(2):247-255.

    [4]Irigoien X,Harris RP,Verheye HM et al.Copepod hatching success in marine ecosystems with high diatom concentrations.Nature,2002,419(6905):387-389.

    [5]Cardinale BJ,Palmer MA,Collins SL.Species diversity enhances ecosystem functioning through interspecific facilitation.Nature,2002,415(6870):426-429.

    [6]Figueredo CC,Giani A.Seasonal variation in the diversity and species richness of phytoplankton in a tropical eutrophic reservoir.Hydrobiologia,2001,445(1/2/3):165-174.

    [7]Negro AI,DE Hoyos C,Vega JC.Phytoplankton structure and dynamics in Lake Sanabria and Valparaíso reservoir(NW Spain).Hydrobiologia,2000,424(1/2/3):25-37.

    [8]Downing JA.Marine nitrogen:phosphorus stoichiometry and the global N:P cycle.Biogeochemistry,1997,37(3):237-252.

    [9]王 越,孫 強,吳 昊.西泉眼水庫富營養(yǎng)化的源分析.北方環(huán)境,2005,30(2):29-30,6.

    [10]遲晉峰,劉玉黛,遲晉旭.哈爾濱市西泉眼水庫水環(huán)境現(xiàn)狀及水源保護對策.黑龍江水利科技,2010,38(1):174-175.

    [11]陳翠翠,于洪賢,姚允龍等.西泉眼水庫消落帶水生植物群落及分布特點.東北林業(yè)大學學報,2012,40(1):110-113.

    [12]高孜娟,于洪賢,馬成學等.西泉眼水庫魚類多樣性與群落結(jié)構(gòu)初步研究.水產(chǎn)學雜志,2012,25(1):19-24.

    [13]呂東珂.哈爾濱周邊泥炭型水庫水-氣界面CO2通量研究[學位論文].哈爾濱:東北林業(yè)大學,2013.

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

    [15]國家環(huán)境保護總局《水和廢水監(jiān)測分析方法》編委會.水和廢水監(jiān)測分析方法:第4版.北京:中國環(huán)境科學出版社,2002.

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

    [17]Lampitt S,Wishner E,Turley M et al.Marine snow studies in the Northeast Atlantic Ocean:distribution,composition and roles as a food source for migrating plankton.Marine Biology,1993,116:680-702.

    [18]Shannon CE.A mathematical theory of communication.ACM SIGMOBILE Mobile Computing and Communications Review,2001,5(1):3-55.

    [19]胡勝華,張世羊,賀 鋒等.武漢月湖水體主要元素與浮游植物的PCA與CCA分析.生態(tài)環(huán)境學報,2010,19(1):17-22.

    [20]TER Braak CJF.CANOCO -an extension of DECORANA to analyze species-environment relationships.Hydrobiologia,1989,184(3):169-170.

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

    [22]董旭輝,羊向東,王 榮等.長江中下游地區(qū)湖泊硅藻-總磷轉(zhuǎn)換函數(shù).湖泊科學,2006,18(1):1-12.

    [23]張 穎,白羽軍,白 昕等.西泉眼水庫富營養(yǎng)化狀況及浮游藻類群落分析.哈爾濱商業(yè)大學學報:自然科學版,2007,23(5):523-525,567.

    [24]馬成學,劉曼紅,黃璞祎等.鏡泊湖枯水期和豐水期浮游植物群落結(jié)構(gòu).東北林業(yè)大學學報,2010,38(2):35-37.

    [25]張麗娜.海浪河浮游植物群落結(jié)構(gòu)特征和水質(zhì)營養(yǎng)狀態(tài)評價[學位論文].哈爾濱:東北林業(yè)大學,2010.

    [26]Jensen J,Jeppesen E,Olrik K et al.Impact of nutrients and physical factors on the shift from cyanobacterial to chlorophyte dominance in shallow Danish lakes.Canadian Journal of Fisheries and Aquatic Sciences,1994,51(8):1692-1699.

    [27]Scheffer M.Ecology of shallow lakes.New York:Springer,2004.

    [28]Reynolds C.The ecology of freshwater phytoplankton.Cambridge:Cambridge University Press,1984.

    [29]Sournia A.Phytoplankton manual:Monographs on oceanographic methodology.UNESCO,1971,6(1):337.

    [30]Scala S,Carels N,F(xiàn)alciatore A et al.Genome properties of the diatom Phaeodactylum tricornutum.Plant Physiology,2002,129(3):993.

    [31]Sommer U.The impact of light intensity and daylength on silicate and nitrate competition among marine phytoplankton.Limnology and Oceanography,1994,39(7):1680-1688.

    [32]Komarkova J,Tavera R.Steady state of phytoplankton assemblage in the tropical Lake Catemaco(Mexico).Hydrobiologia,2003,502(1):187-196.

    [33]李秋華,韓博平.基于CCA的典型調(diào)水水庫浮游植物群落動態(tài)特征分析.生態(tài)學報,2007,27(6):2355-2364.

    [34]盛海燕,虞左明,韓軼才等.亞熱帶大型河流型水庫——富春江水庫浮游植物群落及其與環(huán)境因子的關(guān)系.湖泊科學,2010,22(2):235-243.

    [35]Smith VH.Low nitrogen to phosphorus ratios favor dominance by blue-green algae in lake phytoplankton.Science,1983,221(4611):669.

    猜你喜歡
    泉眼綠藻藍藻
    大地的取水者
    遼河(2022年1期)2022-02-14 05:15:04
    一起走過的日子
    南美白對蝦養(yǎng)殖池塘藍藻水華處理舉措
    香榧綠藻的生物學特性及物種鑒定
    針對八月高溫藍藻爆發(fā)的有效處理方案
    無綠藻的微生物學特性和菌種鑒定
    多變的小泉眼
    全球近10年無綠藻病例報道文獻的回顧
    以成長為主題解讀《窗燈》
    可怕的藍藻
    777米奇影视久久| 韩国高清视频一区二区三区| 亚洲在久久综合| 丰满少妇做爰视频| 午夜福利网站1000一区二区三区| 成人免费观看视频高清| 搞女人的毛片| 午夜免费观看性视频| 欧美 日韩 精品 国产| 亚洲高清免费不卡视频| 国产精品偷伦视频观看了| 亚洲国产日韩一区二区| 亚洲精品中文字幕在线视频 | 97在线视频观看| 日本色播在线视频| 亚洲av国产av综合av卡| 国产永久视频网站| av在线蜜桃| 狠狠精品人妻久久久久久综合| 久久久久网色| 日日摸夜夜添夜夜爱| 天美传媒精品一区二区| 女人久久www免费人成看片| 国产一区有黄有色的免费视频| 婷婷色av中文字幕| 国产精品三级大全| 亚洲欧美成人综合另类久久久| freevideosex欧美| 女人久久www免费人成看片| 日韩在线高清观看一区二区三区| 中文在线观看免费www的网站| 极品教师在线视频| 国产一区二区三区av在线| 高清视频免费观看一区二区| 亚洲不卡免费看| 少妇熟女欧美另类| 一区二区三区四区激情视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品国产露脸久久av麻豆| 国产精品一区二区三区四区免费观看| 亚洲av在线观看美女高潮| 下体分泌物呈黄色| 少妇猛男粗大的猛烈进出视频 | 亚洲精品456在线播放app| 午夜老司机福利剧场| 中文精品一卡2卡3卡4更新| 神马国产精品三级电影在线观看| 免费观看av网站的网址| 成人亚洲精品一区在线观看 | 少妇人妻久久综合中文| 亚洲欧美精品自产自拍| 成人高潮视频无遮挡免费网站| 秋霞伦理黄片| 国产v大片淫在线免费观看| 能在线免费看毛片的网站| 91久久精品国产一区二区三区| 五月玫瑰六月丁香| 日韩欧美精品v在线| 在线a可以看的网站| 久久久久久久久久人人人人人人| 亚洲国产高清在线一区二区三| 99精国产麻豆久久婷婷| 亚洲三级黄色毛片| 免费看日本二区| 精品一区二区三区视频在线| 啦啦啦中文免费视频观看日本| 精品久久国产蜜桃| 久久99精品国语久久久| 欧美日韩综合久久久久久| 99re6热这里在线精品视频| 97超碰精品成人国产| a级一级毛片免费在线观看| 久久久a久久爽久久v久久| 婷婷色av中文字幕| 熟女人妻精品中文字幕| 97精品久久久久久久久久精品| 久久久久精品性色| a级毛色黄片| 免费黄色在线免费观看| 国产黄色视频一区二区在线观看| 国产免费一区二区三区四区乱码| 欧美日韩视频高清一区二区三区二| 欧美日韩视频高清一区二区三区二| 噜噜噜噜噜久久久久久91| 亚洲av男天堂| 精品国产三级普通话版| 亚洲精品自拍成人| 国内少妇人妻偷人精品xxx网站| 中文资源天堂在线| 老师上课跳d突然被开到最大视频| 成人亚洲精品av一区二区| 免费在线观看成人毛片| 男人爽女人下面视频在线观看| 婷婷色综合www| 婷婷色综合www| 18禁裸乳无遮挡免费网站照片| 亚洲国产最新在线播放| 日本欧美国产在线视频| 最近中文字幕高清免费大全6| 国产国拍精品亚洲av在线观看| 一级毛片 在线播放| eeuss影院久久| 精品人妻偷拍中文字幕| 午夜福利网站1000一区二区三区| 成人高潮视频无遮挡免费网站| 欧美97在线视频| 日韩一区二区视频免费看| 免费人成在线观看视频色| 亚洲在线观看片| 少妇的逼水好多| 又黄又爽又刺激的免费视频.| 久久久久久久大尺度免费视频| 最近手机中文字幕大全| 日韩欧美精品免费久久| 日本爱情动作片www.在线观看| 我要看日韩黄色一级片| 亚洲欧美日韩卡通动漫| 久久综合国产亚洲精品| 男女国产视频网站| 国国产精品蜜臀av免费| 麻豆国产97在线/欧美| 国产视频首页在线观看| 在线观看国产h片| 99热网站在线观看| 五月伊人婷婷丁香| 三级国产精品欧美在线观看| 你懂的网址亚洲精品在线观看| 成人亚洲精品一区在线观看 | 亚洲,一卡二卡三卡| 欧美激情在线99| 水蜜桃什么品种好| 婷婷色综合www| 亚洲va在线va天堂va国产| 亚洲高清免费不卡视频| 色视频在线一区二区三区| 日韩强制内射视频| 国产一区二区亚洲精品在线观看| 亚洲一级一片aⅴ在线观看| 99热国产这里只有精品6| 精品人妻视频免费看| 日产精品乱码卡一卡2卡三| 91久久精品电影网| 晚上一个人看的免费电影| 国产黄色视频一区二区在线观看| 国产精品嫩草影院av在线观看| 午夜老司机福利剧场| 免费观看无遮挡的男女| 视频中文字幕在线观看| 国国产精品蜜臀av免费| 国产黄a三级三级三级人| 国产欧美另类精品又又久久亚洲欧美| 搡女人真爽免费视频火全软件| 亚洲国产精品专区欧美| 国产一区二区亚洲精品在线观看| 一区二区三区四区激情视频| 亚洲国产精品成人综合色| 久久99热这里只频精品6学生| 国产综合精华液| 欧美日韩视频高清一区二区三区二| 噜噜噜噜噜久久久久久91| 深爱激情五月婷婷| 少妇猛男粗大的猛烈进出视频 | 在线播放无遮挡| 亚洲性久久影院| 免费大片黄手机在线观看| av天堂中文字幕网| 亚洲激情五月婷婷啪啪| 亚洲真实伦在线观看| 日韩一本色道免费dvd| 亚洲自拍偷在线| 自拍偷自拍亚洲精品老妇| 中文天堂在线官网| 亚洲真实伦在线观看| 91aial.com中文字幕在线观看| 国产午夜精品久久久久久一区二区三区| 久久精品国产自在天天线| 久久精品国产自在天天线| 在线观看一区二区三区| www.色视频.com| 国产男女内射视频| 日本猛色少妇xxxxx猛交久久| 亚洲精品国产成人久久av| 久久影院123| 欧美激情在线99| 丝袜美腿在线中文| 色播亚洲综合网| 国产精品女同一区二区软件| 精品久久久久久久久av| 久久6这里有精品| 欧美成人a在线观看| 日韩 亚洲 欧美在线| 国产v大片淫在线免费观看| 国产色爽女视频免费观看| a级一级毛片免费在线观看| 欧美日本视频| 国产av码专区亚洲av| 大话2 男鬼变身卡| 亚洲av中文av极速乱| 亚洲熟女精品中文字幕| 亚洲欧美成人综合另类久久久| 亚洲伊人久久精品综合| 国产精品人妻久久久影院| 久久ye,这里只有精品| 色视频在线一区二区三区| 亚洲欧美清纯卡通| 熟女人妻精品中文字幕| 中文字幕制服av| 日本一本二区三区精品| 国产伦理片在线播放av一区| 国产精品一区二区在线观看99| 1000部很黄的大片| 成人国产麻豆网| 亚洲伊人久久精品综合| 亚洲最大成人av| 国产成人福利小说| 亚洲最大成人av| 偷拍熟女少妇极品色| videossex国产| 99热6这里只有精品| 国模一区二区三区四区视频| 国产精品久久久久久av不卡| 欧美国产精品一级二级三级 | 欧美 日韩 精品 国产| 内地一区二区视频在线| 国产日韩欧美亚洲二区| 一级毛片aaaaaa免费看小| 亚洲精华国产精华液的使用体验| 九草在线视频观看| 成人亚洲欧美一区二区av| 神马国产精品三级电影在线观看| 亚洲av.av天堂| 91在线精品国自产拍蜜月| 2021少妇久久久久久久久久久| 欧美 日韩 精品 国产| 成人美女网站在线观看视频| 欧美+日韩+精品| 亚洲欧美精品自产自拍| 日本午夜av视频| 日本猛色少妇xxxxx猛交久久| 91精品伊人久久大香线蕉| 久久久色成人| 久久99热6这里只有精品| 欧美激情在线99| 偷拍熟女少妇极品色| 性色av一级| 亚洲欧美成人精品一区二区| 国产成人精品久久久久久| 免费av毛片视频| 午夜亚洲福利在线播放| 国产亚洲精品久久久com| 免费黄频网站在线观看国产| 久久精品国产亚洲av涩爱| 国产免费一级a男人的天堂| 亚洲av成人精品一二三区| 欧美日韩国产mv在线观看视频 | 国产精品蜜桃在线观看| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品国产成人久久av| 欧美成人a在线观看| 日韩一区二区三区影片| 亚洲精品中文字幕在线视频 | 国产毛片在线视频| 国产毛片在线视频| 国产精品久久久久久久电影| 国产乱人视频| 国产精品无大码| 日韩精品有码人妻一区| 亚州av有码| 性插视频无遮挡在线免费观看| 秋霞在线观看毛片| av在线老鸭窝| 欧美激情久久久久久爽电影| 日韩一本色道免费dvd| 精品久久久久久电影网| 少妇的逼好多水| 亚洲高清免费不卡视频| 精品人妻熟女av久视频| 欧美极品一区二区三区四区| 国产色婷婷99| 国产成人精品福利久久| 国产黄频视频在线观看| 又粗又硬又长又爽又黄的视频| 秋霞伦理黄片| 国产成人一区二区在线| 寂寞人妻少妇视频99o| 成年av动漫网址| 国产爽快片一区二区三区| 精品久久久久久久人妻蜜臀av| 三级男女做爰猛烈吃奶摸视频| 亚洲天堂av无毛| 日韩av免费高清视频| 午夜老司机福利剧场| 男人狂女人下面高潮的视频| 国产成人freesex在线| 成年女人看的毛片在线观看| 91aial.com中文字幕在线观看| 国产成年人精品一区二区| 国产精品精品国产色婷婷| 日韩免费高清中文字幕av| 成人二区视频| 联通29元200g的流量卡| 精品国产露脸久久av麻豆| 日韩不卡一区二区三区视频在线| 国国产精品蜜臀av免费| 免费黄色在线免费观看| 日韩电影二区| 国产一区二区三区av在线| 久久久精品94久久精品| 国产在视频线精品| 国产综合懂色| 日韩一本色道免费dvd| 久久精品国产自在天天线| 青青草视频在线视频观看| 国产精品久久久久久精品古装| 中文字幕人妻熟人妻熟丝袜美| av天堂中文字幕网| 五月伊人婷婷丁香| 久久综合国产亚洲精品| 简卡轻食公司| 欧美97在线视频| 亚洲第一区二区三区不卡| 日本wwww免费看| 国精品久久久久久国模美| 如何舔出高潮| 免费看光身美女| 午夜福利在线在线| 国产精品一及| 91精品伊人久久大香线蕉| 国产精品久久久久久久久免| 欧美激情国产日韩精品一区| 麻豆久久精品国产亚洲av| 亚洲av成人精品一二三区| 久久ye,这里只有精品| 毛片女人毛片| 九九爱精品视频在线观看| 不卡视频在线观看欧美| 亚洲精品乱码久久久久久按摩| 国产精品一及| 91精品伊人久久大香线蕉| 热re99久久精品国产66热6| 国产午夜精品一二区理论片| 日本免费在线观看一区| 亚洲av成人精品一二三区| 嫩草影院精品99| 亚洲欧美一区二区三区国产| 亚洲精品一二三| 日本与韩国留学比较| 综合色av麻豆| 五月天丁香电影| 亚洲精品一区蜜桃| 国产 精品1| 老女人水多毛片| av福利片在线观看| 男女下面进入的视频免费午夜| 午夜免费鲁丝| 久久久久性生活片| 国产精品不卡视频一区二区| 国产永久视频网站| 欧美性猛交╳xxx乱大交人| 久久人人爽人人片av| 黄色欧美视频在线观看| 欧美xxⅹ黑人| 一级毛片 在线播放| 一个人观看的视频www高清免费观看| 久久精品夜色国产| 久久久精品免费免费高清| 日韩av在线免费看完整版不卡| 91久久精品国产一区二区三区| 性色av一级| 尾随美女入室| 国产精品福利在线免费观看| 久久综合国产亚洲精品| 91精品国产九色| 日韩,欧美,国产一区二区三区| 成人黄色视频免费在线看| 亚洲精品久久午夜乱码| 蜜桃久久精品国产亚洲av| 伦理电影大哥的女人| 日韩精品有码人妻一区| 一个人观看的视频www高清免费观看| 97精品久久久久久久久久精品| 小蜜桃在线观看免费完整版高清| 亚洲av二区三区四区| 日韩制服骚丝袜av| 性色avwww在线观看| 国产综合懂色| 在线天堂最新版资源| 精品久久久久久电影网| av在线播放精品| 亚洲不卡免费看| 成年女人看的毛片在线观看| 国产成人福利小说| 国产视频内射| 久久久久国产网址| 亚洲精品乱久久久久久| 在线免费十八禁| 亚洲三级黄色毛片| 男人舔奶头视频| 狂野欧美激情性xxxx在线观看| 久久久久国产精品人妻一区二区| 日韩一区二区视频免费看| 日韩欧美 国产精品| 欧美另类一区| 一级av片app| 国产黄片美女视频| 97超视频在线观看视频| 国内精品宾馆在线| 青春草视频在线免费观看| 中文字幕av成人在线电影| 在线观看国产h片| 久热久热在线精品观看| 亚洲成色77777| 自拍欧美九色日韩亚洲蝌蚪91 | 老师上课跳d突然被开到最大视频| 99精国产麻豆久久婷婷| 亚洲精品色激情综合| 狂野欧美白嫩少妇大欣赏| 中文资源天堂在线| 日产精品乱码卡一卡2卡三| 国产黄片视频在线免费观看| 又爽又黄无遮挡网站| 小蜜桃在线观看免费完整版高清| 国产黄频视频在线观看| 综合色丁香网| 中文字幕亚洲精品专区| 欧美日韩综合久久久久久| 日韩强制内射视频| 国产精品久久久久久久久免| 亚洲经典国产精华液单| 亚洲自偷自拍三级| 成人鲁丝片一二三区免费| 麻豆成人av视频| 精品久久久久久久久亚洲| 国产成人a∨麻豆精品| 国产高清三级在线| 天堂中文最新版在线下载 | 我要看日韩黄色一级片| 日本黄色片子视频| 亚洲国产日韩一区二区| 日本午夜av视频| 婷婷色综合www| 日韩在线高清观看一区二区三区| 国产一区二区在线观看日韩| 99久国产av精品国产电影| 日韩欧美一区视频在线观看 | 下体分泌物呈黄色| 日日摸夜夜添夜夜添av毛片| 亚洲欧美成人综合另类久久久| 午夜老司机福利剧场| 久久这里有精品视频免费| 插逼视频在线观看| 久久久久久伊人网av| 亚洲精品成人av观看孕妇| 99久久九九国产精品国产免费| 啦啦啦在线观看免费高清www| 日韩一区二区三区影片| 最近中文字幕2019免费版| 亚洲av电影在线观看一区二区三区 | 高清毛片免费看| 22中文网久久字幕| 性插视频无遮挡在线免费观看| 日韩电影二区| 黄片wwwwww| 美女cb高潮喷水在线观看| 免费看不卡的av| 五月伊人婷婷丁香| 色婷婷久久久亚洲欧美| 久久这里有精品视频免费| 99久久精品热视频| 色视频www国产| 成人亚洲精品一区在线观看 | 国产黄片美女视频| 一级毛片我不卡| 久久精品国产亚洲av天美| 亚洲精品aⅴ在线观看| av在线播放精品| 亚洲人成网站在线观看播放| 国产精品秋霞免费鲁丝片| 国产69精品久久久久777片| 三级国产精品欧美在线观看| 69人妻影院| 成人综合一区亚洲| 少妇被粗大猛烈的视频| 五月天丁香电影| 亚洲av二区三区四区| 精品久久国产蜜桃| 国产精品蜜桃在线观看| 亚洲欧美一区二区三区黑人 | 亚洲国产精品成人久久小说| 白带黄色成豆腐渣| 在线观看美女被高潮喷水网站| 精品久久久噜噜| 大码成人一级视频| 亚洲图色成人| 3wmmmm亚洲av在线观看| 少妇 在线观看| 麻豆久久精品国产亚洲av| 亚洲四区av| h日本视频在线播放| 99九九线精品视频在线观看视频| 日韩视频在线欧美| 国产伦精品一区二区三区视频9| 亚洲欧美日韩另类电影网站 | 性色av一级| 我的老师免费观看完整版| a级毛片免费高清观看在线播放| 乱系列少妇在线播放| 夫妻性生交免费视频一级片| 99热国产这里只有精品6| 精品久久久久久久末码| 欧美+日韩+精品| 国产精品三级大全| 老女人水多毛片| 成人国产麻豆网| 国产精品伦人一区二区| 91午夜精品亚洲一区二区三区| 国精品久久久久久国模美| av专区在线播放| 国产精品国产三级国产av玫瑰| 99热这里只有精品一区| 国产精品偷伦视频观看了| 午夜爱爱视频在线播放| 国产男女超爽视频在线观看| 黄色怎么调成土黄色| 午夜精品一区二区三区免费看| 久久99热这里只有精品18| 欧美成人精品欧美一级黄| freevideosex欧美| 春色校园在线视频观看| 国产黄片美女视频| 亚洲精品久久午夜乱码| 一个人看视频在线观看www免费| 一级a做视频免费观看| 伊人久久国产一区二区| 免费观看a级毛片全部| 久久99热6这里只有精品| 永久免费av网站大全| 日本欧美国产在线视频| 国产精品一区二区在线观看99| 啦啦啦中文免费视频观看日本| 国产精品人妻久久久久久| 成人国产av品久久久| 精品人妻熟女av久视频| 制服丝袜香蕉在线| 能在线免费看毛片的网站| 亚洲国产日韩一区二区| 中国国产av一级| 成人一区二区视频在线观看| 中国三级夫妇交换| 亚洲av在线观看美女高潮| 欧美日韩精品成人综合77777| 久久久久性生活片| 亚洲av欧美aⅴ国产| 久久久久九九精品影院| 网址你懂的国产日韩在线| 久久99热这里只频精品6学生| 人人妻人人看人人澡| 亚洲内射少妇av| 成人美女网站在线观看视频| 久久久久久久久大av| 九九爱精品视频在线观看| 日韩伦理黄色片| 欧美成人精品欧美一级黄| 亚洲精品,欧美精品| 亚洲成人一二三区av| 日产精品乱码卡一卡2卡三| 久久精品国产亚洲网站| 汤姆久久久久久久影院中文字幕| 乱系列少妇在线播放| 两个人的视频大全免费| av女优亚洲男人天堂| 精品熟女少妇av免费看| 一级毛片我不卡| 大陆偷拍与自拍| 性色avwww在线观看| 国产男女超爽视频在线观看| 在线观看免费高清a一片| 少妇高潮的动态图| av在线播放精品| 国产精品一区二区三区四区免费观看| 日日摸夜夜添夜夜爱| 大香蕉久久网| 少妇猛男粗大的猛烈进出视频 | a级毛片免费高清观看在线播放| 波多野结衣巨乳人妻| 一级毛片电影观看| 看黄色毛片网站| 中国国产av一级| 国产伦理片在线播放av一区| 免费大片18禁| 伊人久久精品亚洲午夜| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲av涩爱| 校园人妻丝袜中文字幕| 亚洲最大成人中文| 国精品久久久久久国模美| 国产成人精品福利久久| 少妇猛男粗大的猛烈进出视频 | 我要看日韩黄色一级片| 免费看不卡的av| 久久久精品欧美日韩精品| 中文资源天堂在线| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放| 日日摸夜夜添夜夜爱| 成人亚洲精品一区在线观看 | 亚洲精品日韩av片在线观看| 蜜桃亚洲精品一区二区三区| 久久精品久久久久久噜噜老黄| 亚洲av在线观看美女高潮| 我的女老师完整版在线观看| 午夜精品一区二区三区免费看| 欧美精品国产亚洲| 午夜免费观看性视频| 一级二级三级毛片免费看|