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

    新豐江水庫浮游動物群落結構特征

    2023-04-29 00:00:00秦云霞周雯喬永民陳瑞楊洪允
    水生態(tài)學雜志 2023年4期

    摘要:分析華南大型飲用水源水庫浮游動物群落結構特征及主要環(huán)境影響因子,為水庫生態(tài)環(huán)境保護提供基礎資料和參考依據(jù)。于2019年8月至2020年6月在新豐江水庫設置 12個采樣站點進行每2個月1次的浮游動物監(jiān)測,設置采樣站點12個。結果顯示,調查期間共采集到浮游動物19屬30種,其中輪蟲8屬17種,枝角類5屬7種,橈足類6屬6種,夏季出現(xiàn)的種類數(shù)最多(25種),春季最少(14種);各站點浮游動物密度和生物量月均值分別為11.52~161.52 個/L和0.92~17.65 mg/L,具有極顯著的季節(jié)差異(Plt;0.01),表現(xiàn)為夏季高、冬季低的特點;空間分布上呈現(xiàn)為上游和下游河道型庫區(qū)高,中心敞水區(qū)低的分布趨勢;浮游動物Shannon-Wiener多樣性指數(shù)為0.87~1.19,Margalef豐富度指數(shù)為0.97~1.53,Pielou均勻度指數(shù)為0.54~0.93。浮游動物密度與環(huán)境因子的相關分析和冗余分析(RDA)結果表明,水溫、葉綠素a和pH是影響新豐江水庫浮游動物群落結構的主要環(huán)境因子。

    關鍵詞:浮游動物;密度;生物量;環(huán)境因子;冗余分析;新豐江水庫

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

    浮游動物是淡水生態(tài)系統(tǒng)中的重要生物類群,既是浮游植物的攝食者,又是魚類等高營養(yǎng)級捕食者的開口餌料(楊宇峰和黃祥飛,1992;Lomartire et al,2021),在物質轉化、能量流動、信息傳遞等生態(tài)過程中起著承上啟下的作用(Jia et al,2016)。浮游動物作為中間營養(yǎng)級,受上行效應和下行效應機制影響顯著(Zhou et al,2020),對水體環(huán)境理化因子以及魚類、浮游植物等生物因子的變化反應敏感,其種類組成、密度、生物量變化及其生物多樣性特征被視為反映水環(huán)境質量的重要指征(Shen et al,2019)。因此探究浮游動物群落與水體環(huán)境因子的關系,對于深入了解水生態(tài)系統(tǒng)結構特點,科學評估其環(huán)境健康狀態(tài),確保水質與水生態(tài)安全具有重要意義(楊蘇文等,2015)。

    新豐江水庫是廣東省第一大水庫,庫容1.39×1010 m3,是珠三角地區(qū)重要的飲用水源地和東江水量調節(jié)樞紐,是支撐粵港澳大灣區(qū)社會經(jīng)濟發(fā)展的戰(zhàn)略資源(文曉慧和鐘標城, 2016)。新豐江水庫主要水質指標常年保持在國家地表水Ⅰ類標準,綜合營養(yǎng)狀態(tài)指數(shù)(TLI)為18~27,水庫整體處于貧營養(yǎng)狀態(tài),具有獨特的生態(tài)系統(tǒng)結構(張輝等, 2022)。1995年,白慶笙等(1995)對新豐江水庫浮游動物群落特點進行了初步研究,討論了新豐江水庫水質狀況與浮游動物群落結構的關系;2007年,趙帥營和韓博平(2007)研究了貧營養(yǎng)環(huán)境狀態(tài)下浮游動物的群落特征,探討了銀魚對浮游動物群落結構的影響。此后對新豐江水庫浮游動物的相關研究鮮有報道。近年來,隨著人為干擾的加劇,新豐江水庫水質和水生態(tài)發(fā)生了較為明顯的變化,如偶發(fā)水質異常、藻類短期內快速增殖和水色變化等現(xiàn)象(張輝等,2022)。因此,從水質和水生態(tài)保護的角度對該水庫進行深入的生態(tài)學研究,可為水庫生態(tài)環(huán)境保護提供基礎資料和參考依據(jù)。

    本研究在優(yōu)化調查點位和頻次的基礎上,對新豐江水庫浮游動物進行了2個月1次為期1年的野外定量采樣調查,分析了新豐江水庫浮游動物的群落結構特征及其時空變化;結合水質分析數(shù)據(jù),并通過相關分析和冗余分析(RDA)研究了新豐江水庫浮游動物群落結構與環(huán)境因子的相關關系;探討了貧營養(yǎng)水庫中影響浮游動物群落結構的關鍵環(huán)境因子,以期為新豐江水庫生態(tài)系統(tǒng)長期變化研究積累資料。

    1 " 材料與方法

    1.1 " 研究區(qū)域概況

    新豐江水庫位于廣東省河源市東源縣(114°15’~ 114°50’E,23°40’~23°10’N),是全國第七大水庫和廣東省最大飲用水源保護地,水域面積364 km2,平均年入庫水量0.6×1010 m3,多年平均水深28.7 m,最大水深93 m(周梅等,2018)。新豐江水庫所在區(qū)域屬亞熱帶季風氣候,年平均氣溫為20.7℃,氣候溫和。年均降水量1 793.2 mm,年最大降水量2 732 mm,最小1 050.9 mm(李可見等,2018)。

    1.2 " 樣點設置

    根據(jù)《湖泊調查技術規(guī)程》(中國科學院南京地理與湖泊研究所,2015)、《全國淡水生物物種資源調查技術規(guī)定(試行)》(環(huán)境保護部生態(tài)司,2010)中關于水庫浮游動物定量采樣相關樣點布設標準及新豐江水庫的水文特點,在調查水庫中共設置12個采樣站點(圖1),于2019年8月至2020年6月進行了浮游動物定量采樣調查,采樣頻次為2個月1次。

    1.3 " 樣品采集與處理

    浮游動物定性樣品采集使用13號浮游生物網(wǎng),網(wǎng)口在水面下50 cm作“∞”狀多次拖網(wǎng)采集,水樣收集到50 mL聚乙烯瓶中,加入4%福爾馬林溶液進行現(xiàn)場固定。定量樣品采用采水器分別在0.5、5、10、15 m處采集水樣置于水桶中混合均勻,經(jīng)25號浮游生物網(wǎng)過濾濃縮,加入4%福爾馬林溶液現(xiàn)場固定,靜置24 h后吸取上清液濃縮定容至30 mL,在光學顯微鏡下進行物種鑒定與計數(shù)。浮游動物的種類鑒定參考《中國淡水輪蟲志》(王家楫,1961)、《中國動物志?淡水枝角類》(蔣燮治和堵南山,1979)和《中國動物志?淡水橈足類》(中國科學院動物研究所甲殼動物研究組,1979)。浮游動物計數(shù)以及密度和生物量的計算方法參照《淡水浮游生物研究方法》(章宗涉和黃祥飛,1991)。

    浮游動物采樣的同時,現(xiàn)場測定部分水體理化指標。透明度(SD)用塞氏透明度盤測定,水溫(WT)、pH、溶解氧(DO)和電導率(EC)使用YSI水質分析儀(Hash, DS5X)測定??偭祝═P)、總氮(TN)、高錳酸鹽指數(shù)(CODMn)、葉綠素a(Chl-a)采集水樣帶回實驗室檢測,分析方法參照《國家地表水環(huán)境質量監(jiān)測網(wǎng)監(jiān)測任務作業(yè)指導書(試行)》(環(huán)境保護部環(huán)境監(jiān)測司,2017)。

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

    1.4.1 " 優(yōu)勢度 " 浮游動物優(yōu)勢物種依據(jù)其優(yōu)勢度值(Y)確定,當Ygt;0.02時為優(yōu)勢種。

    [Y=niN·f][i]

    式中:[ni]為第i種在各采樣點的總個體數(shù);N為各采樣點中浮游動物的總個體數(shù);[fi]為第i種在各采樣點出現(xiàn)的頻率。

    1.4.2 " 多樣性指數(shù) " Margalef豐富度指數(shù)(D)、Shannon-Wiener多樣性指數(shù)(H')和Pielou均勻度指數(shù)(J)計算公式如下:

    D = (S-1)/lnN " " " " " " " "①

    式中:D為物種豐富度指數(shù);S為浮游動物種類數(shù);N 為浮游動物總個體數(shù)。

    H' =-[(Ni/N)·ln(Ni/N)] ②

    式中:H'為生物多樣性指數(shù);[Ni]為第i種的個體數(shù);N為浮游動物總個體數(shù)。

    J = H'/lnS " " " " " " " " " " "③

    式中:H'為Shannon-Wiener多樣性指數(shù);S為浮游動物種類數(shù)。

    1.5 " 統(tǒng)計分析

    采用Canoco for Windows5軟件對浮游動物密度數(shù)據(jù)進行去趨勢分析(DCA),結果顯示,排序軸梯度長度(LGA)最大長度為2.3,因此對浮游動物密度與環(huán)境因子進行冗余分析(RDA)。其中以該物種在各樣點出現(xiàn)的頻率≥25%且相對密度≥1%的條件進行篩選(Lopes et al,2005)。采用SPSS17.0、Excel2013、ArcGIS10.6和Origin pro 8.0軟件對數(shù)據(jù)進行處理和統(tǒng)計繪圖。

    2 " 結果與分析

    2.1 " 水質理化指標

    調查期間新豐江水庫的平均水溫為(25.4±0.5)℃,變化范圍在20.7~30.9℃,具有典型的亞熱帶水域特點;溶解氧濃度為7.08~8.43 mg/L,均值為(7.84±0.08)mg/L;pH平均值為7.56±0.05,變化范圍為7.16~8.31,水體呈弱堿性;葉綠素a濃度為0.4~3.3 μg/L,均值為(1.4±0.5)μg/L。TP濃度為0.003~0.021 mg/L,均值為(0.009±0.001)mg/L,TP含量在春夏季高于秋冬季;TN濃度0.4~0.99 mg/L,均值為(0.54±0.04)mg/L;高錳酸鹽指數(shù)的變化范圍為0.9~2.4 mg/L,均值為(1.52±0.07)mg/L。研究期間,除總氮外,新豐江水庫其他主要水質指標濃度均值均達到《地表水環(huán)境質量標準》(GB3838-2002)Ⅰ類水質標準。

    2.2 " 浮游動物種類組成

    調查期間共鑒定出浮游動物19屬30種,以輪蟲種類數(shù)最多,為8屬17種,占總種類數(shù)的56.7%;枝角類次之,為5屬7種,占總種類數(shù)的23.3%;橈足類為6屬6種,占20.0%(表1)。其中夏季浮游動物種類數(shù)最多(25種),春季最少(14種)。舌狀葉鏢水蚤(Phyllodiaptomus tunguidus)、右突新鏢水蚤(Neodiaptomus schmackeri)、廣布中劍水蚤(Mesocyclops leuckarti)、臺灣溫劍水蚤(Thermocyclops taihokuensis)和長額象鼻溞(Bosmina longirostris)在各月份均有出現(xiàn)。

    調查期間共出現(xiàn)浮游動物優(yōu)勢種8種(Y≥ 0.02),其中橈足類4種:舌狀葉鏢水蚤、右突新鏢水蚤、廣布中劍水蚤和臺灣溫劍水蚤;枝角類4種:短尾秀體溞(Diaphanosoma brachyurum)、僧帽溞(Daphnia cucullata)、長額象鼻溞和微型裸腹溞(Moina micrura)。從時間分布特征來看,僅有廣布中劍水蚤密度在各調查月份均處于優(yōu)勢地位,其次為長額象鼻溞,其密度優(yōu)勢在調查期間也非常明顯,僅在2020年4月的優(yōu)勢度低于0.02(表2)。

    2.3 " 浮游動物密度和生物量的時間分布

    新豐江水庫浮游動物密度范圍是0.20 ~ 880.00 個/L,平均值為(60.13±14.47)個/L。浮游動物密度具顯著季節(jié)性變化(Plt;0.01),月均密度最高值出現(xiàn)在2019年8月,為(171.66±72.79)個/L,此后大幅下降,2020年2月出現(xiàn)最低值(0.75±0.19)個/L。由春季至夏季,浮游動物密度緩慢上升,至2020年6月達到(81.31±29.65)個/L。從密度組成特征來看,橈足類和枝角類占絕對優(yōu)勢(圖2a)。

    生物量在調查期間變化范圍為0.02 ~ 98.96 mg/L,均值為(5.00±1.48)mg/L。生物量季節(jié)變化呈現(xiàn)為夏季高于冬季。從月度變化看,浮游動物月均生物量在2019年8月最高,為(16.99±7.99)mg/L,此后,生物量大幅下降,至2020年2月出現(xiàn)最低值(0.09±0.02)mg/L,與密度變化相似。浮游動物生物量在種類組成上均以橈足類的生物量占比最大(圖2b)。

    2.4 " 浮游動物密度和生物量的空間分布

    研究期間,各站點月均密度水平為11.52~161.52 個/L,表現(xiàn)為水庫中心敞水區(qū)低,上游和下游的河道型庫區(qū)高的特點。其中,最高值出現(xiàn)在下游S12站點,為(161.52±143.76)個/L,次高值出現(xiàn)在上游的S1站點,為(124±66.50)個/L,最低值出現(xiàn)在水庫中心敞水區(qū)的S10站點,為(11.52±2.75)個/L(圖3a);生物量在各站點變化范圍為0.92 ~ 17.65 mg/L,其空間分布特征與密度相似,表現(xiàn)為開闊的中心區(qū)域低,河道型庫區(qū)的上游和下游高的特點,其中最高值出現(xiàn)在S12采樣點,在S10采樣點出現(xiàn)最低值(圖3b)。

    2.5 " 浮游動物多樣性時空變化

    浮游動物Margalef豐富度指數(shù)(D)均值為1.21±1.32,變化范圍為0.97~1.53。最高值出現(xiàn)在2019年12月(1.53±0.80),2019年10月和2020年2月的豐富度指數(shù)略低于2019年12月,表現(xiàn)為秋冬季偏高的特點(圖4a);在各站點的變化范圍為0.41~1.95,上游峽谷型庫區(qū)較高,然后降低,至水庫開闊的中心區(qū)域再次出現(xiàn)高值,最高值出現(xiàn)在S9站點(圖4b)。Shannon-Wiener多樣性指數(shù)(H')均值為1.03±0.33,變化范圍在0.87~1.19之間,其時間變化與D值呈相反趨勢,最低值出現(xiàn)在2019年12月,最高值出現(xiàn)在2020年6月(圖4a);從空間分布來看,由上游至下游呈下降趨勢,并在S12站點出現(xiàn)極低值(圖4b)。Pielou均勻度指數(shù)(J)均值為0.66±0.17,變化范圍在0.54~0.93之間,其時間變化有別于D值和H'值,在2020年2月份出現(xiàn)最高值,最低值出現(xiàn)在2019年12月(圖4a);在空間分布上無顯著變化,變化范圍為0.59~0.71(圖4b)。

    2.6 " 浮游動物與環(huán)境因子相關性分析

    浮游動物密度和生物量與環(huán)境因子進行Pearson相關性分析,結果表明浮游動物的密度和生物量與溫度、pH、葉綠素a之間呈顯著的正相關關系(Plt;0.05),表明水溫、pH和以葉綠素a為指征的浮游植物對浮游動物的生長與繁殖起著重要的調節(jié)作用;密度和生物量與透明度呈顯著負相關,表明浮游動物在透明度高的水域,密度和生物量較低;密度和生物量與總氮、總磷以及溶解氧的相關性不顯著。

    2.7 nbsp; 浮游動物與環(huán)境因子多元分析

    從采集到的浮游動物種類中篩選出共10種與環(huán)境因子進行RDA分析,物種及代碼見圖5和表1。結果表明2個環(huán)境排序軸的相關系數(shù)為0,Monte Carlo置換檢驗所有排序軸均達到顯著水平(Plt;0.05),排序結果可靠(表4)。第一排序軸和第二排序軸累計解釋了88.9%的物種信息量。

    由表5和圖5可知,溫度、pH、透明度和葉綠素a與浮游動物群落結構密切相關。其中水溫、pH和葉綠素a與第一排序軸為正相關,相關系數(shù)分別為0.6779、0.6054和0.6497;總氮與第二排序軸呈正相關,相關系數(shù)為0.4909。就物種與環(huán)境因子的關系而言,舌狀葉鏢水蚤、右突新鏢水蚤、廣布中劍水蚤、臺灣溫劍水蚤、長額象鼻溞和頸溝基合溞與pH、葉綠素a和水溫呈正相關,與透明度呈負相關;微型裸腹溞、僧帽溞和無節(jié)幼體與總氮和總磷呈正相關關系,螺形龜甲輪蟲主要受透明度影響。

    3 " 討論

    3.1 " 浮游動物群落結構特征

    1995年在新豐江水庫進行的浮游動物調查發(fā)現(xiàn)輪蟲14種,枝角類2種,橈足類9種,優(yōu)勢種為螺形龜甲輪蟲和針簇多肢輪蟲(Polyarthra trigla)(白慶笙等,1995);2002年的調查發(fā)現(xiàn)輪蟲10種,枝角類7種,橈足類3種,在豐水期以輪蟲居多,枯水期以橈足類為主(陳麗芬等,2002);2007年的調查發(fā)現(xiàn)輪蟲32種,枝角類7種,橈足類7種,浮游動物在數(shù)量上以橈足類占優(yōu)勢(趙帥營和韓博平,2007)。而本研究的調查,鑒定發(fā)現(xiàn)輪蟲17種,枝角類7種,橈足類6種,浮游動物的優(yōu)勢種和密度均以橈足類和枝角類為主,研究結果與前期的調查報道稍有差異。綜上可知,新豐江水庫浮游動物在群落組成上由以輪蟲、橈足類為主的群落特征向以橈足類和枝角類占優(yōu)勢的格局轉變,在組成個體上趨向大型化。該現(xiàn)象產(chǎn)生的原因,一方面取決于調查范圍和時間的差異;另一方面還受到下行效應的影響。浮游動物作為魚類的餌料,魚類捕食常常會導致浮游動物群落結構中的個體大小發(fā)生變化(Chaparro et al,2014)。根據(jù)大小效率假說(Brooks amp; Dodson,1965),魚類在攝食選擇上,通常會優(yōu)先選擇個體較大的浮游動物。因此,在魚類捕食壓力增加的情況下,個體較小的種類(如輪蟲、小型枝角類)因能夠逃避捕食而占據(jù)優(yōu)勢(Yang et al,2005),反之則以大型枝角類和橈足類數(shù)量占優(yōu)勢(Jemberg et al,2017)。因此魚類捕食強度是調控浮游動物重要因素之一,近年來,新豐江水庫魚類捕撈強度有增無減。實際調查顯示,漁獲量與10年前相比大幅下降,部分種類生物量已不足過去的1%,捕食壓力的降低,為大型枝角類增殖和發(fā)育創(chuàng)造了有利的條件,其密度達到26.46 個/L,約為之前調查結果的10倍。Korponai等(2003)對Balaton湖浮游動物群落特征的研究也發(fā)現(xiàn)了類似現(xiàn)象。橈足類和枝角類作為輪蟲食物競爭的有力對手,它們的大量存在會使輪蟲的生長受到抑制(溫新利等,2006),這是新豐江水庫輪蟲密度相對較低的主要原因之一。

    一般認為,浮游動物密度與水體營養(yǎng)水平呈正相關(Pinto-Coelho et al,2011;魯敏和謝平,2002)。本次調查發(fā)現(xiàn)新豐江水庫浮游動物的平均密度為60.13 個/L,遠低于處于富營養(yǎng)化水平的湖泊,如太湖梅梁灣浮游動物的平均密度為1 774.33 個/L(楊佳等,2020),滇池為6 478 個/L(楊蘇文等,2015)。浮游動物密度較低的特點與新豐江水庫所處的貧營養(yǎng)狀態(tài)有關。潘繼征等(2009)對撫仙湖浮游動物的群落研究發(fā)現(xiàn),湖泊的貧營養(yǎng)狀態(tài)是浮游動物密度低的主要原因之一。營養(yǎng)鹽作為浮游植物生長的必需元素,通過上行效應來影響浮游動物的分布(Li et al,2019)。當水體營養(yǎng)鹽不足,浮游植物的種群數(shù)量會處于較低水平。新豐江水庫浮游植物在數(shù)量上屬于比較低的水平(張輝等,2022),作為浮游動物食物的浮游植物在水體中的數(shù)量較少必然會使浮游動物的數(shù)量受到限制。有研究表明,藻類密度的降低會導致橈足類體型變小和產(chǎn)卵量降低(Makino amp; Ban,2000)。

    新豐江水庫浮游動物密度分布的顯著特點是水庫上游和下游的河道型庫區(qū)高于中心敞水區(qū)。水庫上游有入庫支流惡馬坑、南溪坑、大寨水和小水河的匯入,營養(yǎng)鹽相對豐富(TP濃度為0.15 mg/L,TN濃度為2.2 mg/L),增加了浮游動物的種類和密度,如微型裸腹溞、短尾秀體溞等主要在這一區(qū)域出現(xiàn);而中心敞水區(qū)是連接大壩的開闊水域,此區(qū)深度和透明度較高,營養(yǎng)鹽含量相對較低(TP:0.001~0.009 mg/L;TN:為0.40~0.62 mg/L),不利于浮游動物的生長繁殖(陳佳琪等,2020)。

    3.2 " 浮游動物群落結構與環(huán)境因子的關系

    不同的水體中環(huán)境因素對浮游動物群落結構的影響存在差異(Cai et al,2020)。相關分析和RDA分析顯示新豐江水庫浮游動物群落結構與溫度、葉綠素a和pH呈極顯著正相關關系。其中水溫是影響浮游動物生長、繁殖發(fā)育、群落組成和數(shù)量變化等極為重要的環(huán)境因子(Zhao et al,2020;Kubo et al,2013),也是影響浮游動物季節(jié)演替的關鍵因子(Record et al,2010)。新豐江水庫地處亞熱帶地區(qū),夏秋季節(jié)具有良好的光照和水溫條件,而春冬季溫度相對較低,水溫條件的季節(jié)差異為浮游動物群落結構的季節(jié)性變化創(chuàng)造了條件。在本次調查中,浮游動物密度呈現(xiàn)出夏秋季節(jié)高,春冬季節(jié)低的特點。溫度的升高可以使浮游動物的繁殖速度加快,發(fā)育期縮短(劉寶貴等,2016),進而使浮游動物豐度在夏季(8月)達到峰值。金瓊貝等(1991)研究發(fā)現(xiàn),橈足類、枝角類和無節(jié)幼體的個體數(shù)隨著溫度的升高呈上升趨勢,并在30℃時達到最高值,其中廣布中劍水蚤受溫度影響最為顯著,其個體數(shù)占橈足類總數(shù)的63.93%。RDA分析也表明,在新豐江水庫中,常見的浮游動物種和優(yōu)勢種受水溫影響,與之呈密切的正相關性。

    浮游動物的季節(jié)性變化與食物濃度也有著密切的關系(Rinke amp; Vijverberg,2005)。浮游植物作為水生態(tài)系統(tǒng)的初級生產(chǎn)者,是浮游動物重要的食物來源,而葉綠素a作為浮游植物生物量的代表參數(shù)(James et al,1988),能夠直接影響浮游動物的分布。本研究中,葉綠素a季節(jié)性變化與浮游動物一致,表現(xiàn)為葉綠素a含量在夏季明顯高于冬季。其原因主要是浮游植物由于夏季高溫大量繁殖,導致葉綠素a含量明顯上升,為浮游動物提供了充足的食物來源。研究表明,橈足類傾向于攝食藻類和有生命的高營養(yǎng)顆粒(Russo et al,2016),在整個調查期間,新豐江水庫浮游植物生物量以小型的小球藻屬(Chlorella)和小環(huán)藻屬(Cyclotella)為主,是浮游動物主要攝食類群,而枝角類傾向于攝食直徑為3 ~ 15 [μ]m的藻類(如Chlorella),100 [μ]m以上的則不能利用(Sommer et al,2003)。本研究中,RDA分析也顯示了浮游動物與葉綠素a呈極顯著的正相關關系,同時表明了以小球藻屬和小環(huán)藻屬為主的浮游植物為浮游動物繁殖發(fā)育提供了一定的餌料資源。

    浮游動物的種類、數(shù)量變化與水體pH值密切相關(楊麗麗等,2011)。研究表明輪蟲和橈足類適宜生活在偏堿性水體中,枝角類則在酸性、中性及微堿性水體中生活得更好(高原等,2013)。此次調查期間,新豐江水庫pH為7~9,屬于典型的弱堿性水體。其中夏季pH明顯高于其他季節(jié)(Plt;0.05),其主要原因之一在于夏季藻類增殖和光合作用增強,消耗大量溶解態(tài)CO2,使水體酸堿平衡向堿性偏移,受其影響新豐江水庫夏季浮游動物密度也明顯高于其他季節(jié),相關分析也顯示浮游動物的密度和生物量與pH呈極顯著的正相關關系。類似的現(xiàn)象也出現(xiàn)于百花水庫(胡藝等,2020)、鄱陽湖(陳佳琪等,2020)、西涼湖(吳利等,2011)等水體。

    4 " 結論

    (1)調查期間在新豐江水庫共鑒定出浮游動物19屬30種,其中優(yōu)勢種為舌狀葉鏢水蚤、右突新鏢水蚤、廣布中劍水蚤、臺灣溫劍水蚤、短尾秀體溞、僧帽溞、長額象鼻溞和微型裸腹溞。魚類的長期過度捕撈是新豐江水庫浮游動物群落以大型橈足類和枝角類為優(yōu)勢的主要原因。

    (2)新豐江水庫浮游動物密度為0.20~880.00個/L,月均密度為(60.13±14.47)個/L,其季節(jié)表現(xiàn)為夏季gt;秋季gt;春季gt;冬季;在空間上河道型庫區(qū)浮游動物較為豐富,中心敞水區(qū)浮游動物較少,營養(yǎng)鹽時空分異是該現(xiàn)象產(chǎn)生的主要調控因素;生物量時空變化與密度相似,其變化范圍為0.02~98.96 mg/L,均值為(5.00±1.48)mg/L。

    (3)相關分析和RDA分析結果表明,水溫是影響新豐江水庫浮游動物季節(jié)變化的主導因素,其次為葉綠素a和pH值。

    參考文獻

    白慶笙,徐潤林,余日清,等,1995. 新豐江水庫浮游動物群落結構研究[J].中山大學學報叢論,(1): 135-141.

    陳佳琪,趙坤,曹玥,2020. 鄱陽湖浮游動物群落結構及其與環(huán)境因子的關系[J].生態(tài)學報,40(18): 1-14.

    陳麗芬,林秋奇,胡韌,等,2002. 亞熱帶大型水庫——新豐江水庫的浮游生物群落特征[J].生態(tài)學報,21(2): 104-107.

    胡藝,李秋華,何應,等,2020. 貴州高原水庫浮游動物分布特征及影響因子——以阿哈水庫為例[J].中國環(huán)境科學,40(1): 227-236.

    高原,賴子尼,李捷,等,2013. 連江浮游動物多樣性的空間分布[J].生物多樣性,21(6): 699-708.

    環(huán)境保護部環(huán)境監(jiān)測司,2017.國家地表水環(huán)境質量監(jiān)測網(wǎng)監(jiān)測任務作業(yè)指導書(試行)[M]. 北京:中國環(huán)境出版社.

    環(huán)境保護部生態(tài)司,2010. 全國淡水生物物種資源調查技術規(guī)定(試行)[M].北京:環(huán)境保護部.

    蔣燮治,堵南山,1979. 中國動物志?淡水枝角類[M].北京:科學出版社.

    金瓊貝,盛連喜,張然,1991. 溫度對浮游動物群落的影響[J].東北師大學報(自然科學版), (4):103-111.

    李可見,黃芳芳,張衛(wèi)強,等,2018. 新豐江水庫及庫區(qū)地表水水質分析與評價[J].林業(yè)與環(huán)境科學,34(3): 6-13.

    劉寶貴,劉霞,吳瑤,等,2016. 鄱陽湖浮游甲殼動物群落結構特征[J].生態(tài)學報,36(24): 8205-8213.

    魯敏,謝平,2002. 武漢東湖不同湖區(qū)浮游甲殼動物[J].海洋與湖沼,33(2): 174-181.

    潘繼征,熊飛,李文朝,等,2009. 撫仙湖浮游甲殼動物群落結構與空間分布[J].湖泊科學,21(3): 408-414.

    王家楫,1961. 中國淡水輪蟲志[M].北京:科學出版社.

    溫新利,席貽龍,張雷,等,2006. 蕪湖市鏡湖輪蟲群落結構分析及水質的生態(tài)學評價[J].水生生物學報,30(2): 152-158.

    文曉慧,鐘標城,2016. 新豐江水庫水質現(xiàn)狀研究[J].環(huán)境與發(fā)展防治與處理,6(12): 55-58.

    吳利,馮偉松,張?zhí)昧?,等?011. 湖北省西涼湖浮游動物群落周年動態(tài)變化及其與環(huán)境因子的關系[J].湖泊科學,23(4): 619-625.

    楊佳,周健,秦伯強,等,2020. 太湖梅梁灣浮游動物群落結構長期變化特征(1997-2017年)[J].環(huán)境科學,41(3): 1246-1255.

    楊麗麗,周小玉,劉其根,等,2011. 新疆烏倫古湖浮游甲殼動物的季節(jié)演替及與環(huán)境因子的關系[J].生態(tài)學報,31(15): 4468-4476.

    楊蘇文,金位棟,何佳,等,2015. 滇池、洱海浮游動植物環(huán)境圖譜[M].北京:科學出版社.

    楊宇峰,黃祥飛,1992. 鰱鳙對浮游甲殼動物群落結構的影響[J].湖泊科學,4(3): 78-86.

    張輝,彭宇瓊,鄒賢妮,等,2021. 華南大型深水水庫浮游植物群落特征及營養(yǎng)狀態(tài)評價[J].環(huán)境科學與技術,44(4): 186-194.

    張輝,彭宇瓊,鄒賢妮,等,2022. 新豐江水庫浮游植物功能分組特征及其與環(huán)境因子的關系[J].中國環(huán)境科學,42(1):380-392.

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

    趙帥營,韓博平,2007. 大型深水貧營養(yǎng)水庫——新豐江水庫浮游動物群落分析[J].湖泊科學,19(3): 305-314.

    中國科學院動物研究所甲殼動物研究組,1979. 中國動物志?淡水橈足類[M].北京:科學出版社.

    中國科學院南京地理與湖泊研究所,2015. 湖泊調查技術規(guī)程[M].北京:科學出版社.

    周梅,葉麗菲,張超,等,2018. 廣東新豐江水庫表層水體CO2分壓及其影響因素[J].湖泊科學,30(3): 770-781.

    Brooks J L,Dodson S I,1965. Predation,Body size,and Composition of Plankton[J]. Science,150(3692): 28-35.

    Cai W W,Xia J H,Yang M Z,2020. Cross-basin analysis of freshwater ecosystem health based on a zooplankton-based index of biotic integrity: Models and application[J].Ecological Indicators,114:106333.

    Chaparro G,F(xiàn)ontanarrosa M S,Cataldo D,et al,2014. Hydrology driven factors might weaken fish predation effects on zooplankton structure in a vegetated warm temperate floodplain lake[J]. Hydrobiologia,752(1): 187-202.

    James J,Elser,Monica,et al,1988. Zooplankton-mediated transitions between N- and P- limited algal growth[J]. Limnology amp; Oceanography,33(1): 1-14.

    Jemberg S,Lehtiniemi M,Uusitalo L,2017. Evaluating zooplankton indicators using signal detection theory[J]. Ecological Indicators,77: 14-22.

    Jia Z N,Swadling K M,Meiners K M,et al,2016. The zooplankton food web under East Antarctic pack ice-A stable isotope study[J]. Deep-See Research Part Ⅱ,131: 189-202.

    Korponai J,Paulovits G,Mátyás K,et al,2003. Long-term changes of cladoceran community in a shallow hypertrophic reservoir in Hungary[J]. Hydrobiologia,504(1/3): 193-201.

    Kubo T,Asai T,Okumura Y,2013. Environmental study of the seasonal succession of mesozooplankton in a brackish water Odra River estuary during 2003-2005[J]. Journal of Ecological Engineering,14(2): 12-25.

    Li Y,Liu L,Cui S Z,et al,2019. Long-term effects of nutrient changes on rotifer communities in a subtropical lake[J]. Limnology,20(1): 191-201.

    Lomartire S,Marques J C,Goncalves A M,2021. The key role of zooplankton in ecosystem services: A perspective of interaction between zooplankton and fish recruitment[J]. Ecological Indicators,129: 107867.

    Lopes M,Bicudo C,F(xiàn)erragut M C,2005. Short term spatial and temporal variation of phytoplankton in a shallow tropical oligotrophic reservoir,southeast Brazil[J]. Hydrobiologia,542(1): 235-247.

    Makino W,Ban S,2000. Response of life history traits to food conditions in a Cyclopoid copepod from an oligotrophic environment[J]. Limnology amp; Oceanography,45(2): 396-407.

    Pinto-Coelho R,Pinel-Alloul B,Méthot G,et al,2011. Crustacean zooplankton in lakes and reservoirs of temperate and tropical regions: Variation with trophic status[J]. Canadian Journal of Fisheries and Aquatic Sciences,62(2): 348-361.

    Record N R,Pershing A J,Jossi J W,2010. Biodiversity as a dynamic variable in the Gulf of Maine continuous plankton recorder transect[J]. Journal of Plankton Research,32(12): 1675-1684.

    Rinke K,Vijverberg J,2005. A model approach to evaluate the effect of temperature and food concentration on individual life-history and population dynamics of Daphnia[J]. Ecological Modeling,186(3): 326-344.

    Russo E,F(xiàn)ranke K,Hager H,et al,2016. Modifying the functional diversity in the zooplankton assemblage of an oligotrophic lake differentially affects pelagic community structure and biomass[J]. Food Webs,(8): 23-31.

    Shen J C,Zhang Z H,Liu R,et al,2019. Ecological restoration of eroded karst utilizing pioneer moss and vascular plant species with selection based on vegetation diversity and underlying soil chemistry[J]. International Journal of Phytoremediation,20(7): 1-11.

    Sommer U,Sommer F,Santer B,et al,2003. Daphnia versus copepod impact on summer phytoplankton: functional compensation at both trophic levels[J]. Oecologia,135: 639-647.

    Yang Y F,Huang X F,Liu J K,et al,2005. Effects of fish stocking on the zooplankton community structure in a shallow lake in China[J].Fisheries Management and Ecology,12(2): 81-89.

    Zhao Q Y,Liu S T,Niu X L,2020. Effect of water temperature on the dynamic behavior of phytoplankton-zooplankton model[J]. Applied Mathematics and Computation,378: 125211.

    Zhou J,Qin B Q,Zhu G W,et al,2020. Long-term variation of zooplankton communities in a large,heterogenous lake: Implications for future environmental change scenarios[J].Enviromental Research,187(12): 109704.

    (責任編輯 " 鄭金秀)

    Zooplankton Community Structure in Xinfengjiang Reservoir

    QIN Yun‐xia1, ZHOU Wen2, QIAO Yong‐min1, CHEN Rui1, YANG Hong‐yun1

    (1. Institute of Hydrobiology,Jinan University,Guangzhou " 510632, P.R. China;

    2. South China Institute of Environmental Science,MEE,Guangzhou " 510655, P.R. China)

    Abstract:Xinfengjiang reservoir is the largest reservoir in Guangdong Province, serving both as a drinking water source for the Pearl River Delta region and as the water volume regulation hub for Dongjiang River. In this study, we explored the community structure and spatiotemporal variation of zooplankton in Xinfengjiang reservoir based on zooplankton monitoring at 12 sampling sites once every two months from August 2019 to June 2020. Then, combined with the results of water quality monitoring, we analyzed the relationship of the zooplankton community with water environmental parameters and discussed the primary environmental factors in large drinking water source reservoirs of South China. The aim was to accumulate baseline data for a study of long-term ecosystem variation in Xinfengjiang reservoir and provide a reference to support the protection of water resources and aquatic ecology. A total of 30 zooplankton species from 19 genera were observed including 17 rotifer species from 8 genera,7 cladoceran species from 5 genera and 6 copepod species from 6 genera. There were 8 dominant species (Y≥0.02) identified during the investigation, consisting of 4 cladoceran species and 4 copepod species. The number of zooplankton species was highest (25 species) in summer and lowest (14 species) in spring. The ranges of zooplankton density and biomass were 11.52-161.52 ind/L and 0.92-17.65 mg/L, with extreme seasonal variation (Plt;0.01), and both were generally higher in summer and lower in winter. Spatially, zooplankton density and biomass were higher in the riverine sections, upstream and downstream, and lower in the central lacustrine area. The ranges of the Shannon-Wiener diversity index, Margalef richness index and Pielou evenness index were, respectively, 0.87-1.19,0.97-1.53 and 0.54-0.93. Correlation and redundancy analyses of zooplankton density and water environmental parameters show that water temperature, Chl-a and pH were the primary environmental factors affecting zooplankton community structure.

    Key words:zooplankton; density; biomass; environmental factors; redundancy analysis; Xinfengjiang reservoir

    1024香蕉在线观看| 在线看a的网站| 一区二区三区激情视频| 欧美久久黑人一区二区| 国产欧美日韩综合在线一区二区| 欧美成人精品欧美一级黄| 看免费av毛片| 人成视频在线观看免费观看| 精品人妻在线不人妻| 午夜激情av网站| 久久综合国产亚洲精品| 亚洲一卡2卡3卡4卡5卡精品中文| 久久亚洲精品不卡| 老司机在亚洲福利影院| 久久国产精品人妻蜜桃| 99香蕉大伊视频| 亚洲图色成人| 久久人人爽av亚洲精品天堂| 精品一区二区三区四区五区乱码 | netflix在线观看网站| 日韩av免费高清视频| 在线观看免费日韩欧美大片| 国产一区二区激情短视频 | 精品一区二区三区四区五区乱码 | h视频一区二区三区| 欧美性长视频在线观看| 久久久久国产一级毛片高清牌| 亚洲一码二码三码区别大吗| 久热这里只有精品99| 中文字幕高清在线视频| 国产亚洲欧美精品永久| 欧美日韩视频精品一区| 久久精品人人爽人人爽视色| 国产午夜精品一二区理论片| 亚洲精品国产av蜜桃| 下体分泌物呈黄色| 亚洲欧洲精品一区二区精品久久久| 欧美精品一区二区免费开放| 性高湖久久久久久久久免费观看| 三上悠亚av全集在线观看| 91九色精品人成在线观看| 国产男女内射视频| 国产精品偷伦视频观看了| 丝袜喷水一区| 欧美日韩精品网址| 亚洲精品久久午夜乱码| 日本vs欧美在线观看视频| 电影成人av| 亚洲成色77777| 国产淫语在线视频| 国产视频一区二区在线看| 一级片'在线观看视频| 99re6热这里在线精品视频| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品九九99| 久久人人97超碰香蕉20202| 免费高清在线观看日韩| 精品亚洲乱码少妇综合久久| 真人做人爱边吃奶动态| 少妇人妻久久综合中文| 国产有黄有色有爽视频| 亚洲国产精品999| 亚洲激情五月婷婷啪啪| 欧美人与善性xxx| 好男人电影高清在线观看| 超碰成人久久| 高清不卡的av网站| 一级片'在线观看视频| 最新的欧美精品一区二区| www.999成人在线观看| 欧美+亚洲+日韩+国产| 黄色视频在线播放观看不卡| 我的亚洲天堂| 亚洲成国产人片在线观看| 欧美精品一区二区免费开放| 国产有黄有色有爽视频| 国产精品二区激情视频| 日韩视频在线欧美| 国产成人欧美| 欧美人与性动交α欧美软件| 国产片内射在线| 十分钟在线观看高清视频www| 国产精品偷伦视频观看了| 国产伦人伦偷精品视频| 国产日韩欧美视频二区| 国产1区2区3区精品| 午夜91福利影院| 中文字幕精品免费在线观看视频| 我要看黄色一级片免费的| av不卡在线播放| 黄色视频在线播放观看不卡| 免费看av在线观看网站| 亚洲国产成人一精品久久久| 亚洲伊人色综图| 久热这里只有精品99| 岛国毛片在线播放| 亚洲成色77777| 久久99一区二区三区| 亚洲五月婷婷丁香| 别揉我奶头~嗯~啊~动态视频 | 最新的欧美精品一区二区| 欧美成狂野欧美在线观看| 菩萨蛮人人尽说江南好唐韦庄| √禁漫天堂资源中文www| 50天的宝宝边吃奶边哭怎么回事| 精品少妇久久久久久888优播| 日日摸夜夜添夜夜爱| 无限看片的www在线观看| 9热在线视频观看99| 婷婷成人精品国产| 久久久久久久大尺度免费视频| 亚洲av成人不卡在线观看播放网 | 十八禁高潮呻吟视频| 日韩电影二区| 80岁老熟妇乱子伦牲交| 男女床上黄色一级片免费看| 人体艺术视频欧美日本| 一级片免费观看大全| 久久免费观看电影| 青草久久国产| 亚洲成国产人片在线观看| 亚洲欧洲日产国产| 久久99热这里只频精品6学生| 中文字幕精品免费在线观看视频| 精品国产超薄肉色丝袜足j| 18禁观看日本| 丁香六月天网| 午夜福利一区二区在线看| 亚洲精品美女久久久久99蜜臀 | 男男h啪啪无遮挡| 2021少妇久久久久久久久久久| 中文乱码字字幕精品一区二区三区| 亚洲五月色婷婷综合| 午夜精品国产一区二区电影| 成年动漫av网址| 亚洲久久久国产精品| 国产av精品麻豆| 亚洲一区二区三区欧美精品| 久久精品久久久久久噜噜老黄| 首页视频小说图片口味搜索 | 国产av一区二区精品久久| 90打野战视频偷拍视频| 在线天堂中文资源库| 亚洲精品美女久久久久99蜜臀 | 久久久久久亚洲精品国产蜜桃av| 青青草视频在线视频观看| 国产精品久久久久成人av| 免费看av在线观看网站| 国产一区有黄有色的免费视频| 亚洲精品美女久久av网站| 黄频高清免费视频| 日本av免费视频播放| 色综合欧美亚洲国产小说| 日本av手机在线免费观看| 高清欧美精品videossex| 日韩中文字幕欧美一区二区 | 香蕉丝袜av| 51午夜福利影视在线观看| 成人手机av| 国产成人一区二区三区免费视频网站 | 欧美中文综合在线视频| 欧美日韩视频高清一区二区三区二| netflix在线观看网站| 国产精品久久久久久精品古装| 涩涩av久久男人的天堂| 国产成人一区二区在线| 成在线人永久免费视频| 男人操女人黄网站| 日韩av免费高清视频| 久久精品国产综合久久久| 超碰97精品在线观看| 青春草视频在线免费观看| 国产成人精品在线电影| 国产黄色免费在线视频| 日韩制服骚丝袜av| 国产人伦9x9x在线观看| 色播在线永久视频| 久久久久国产精品人妻一区二区| 自线自在国产av| 1024香蕉在线观看| 一级毛片黄色毛片免费观看视频| 两性夫妻黄色片| 日韩av免费高清视频| 少妇猛男粗大的猛烈进出视频| 久久久久国产精品人妻一区二区| 女人被躁到高潮嗷嗷叫费观| 日韩欧美一区视频在线观看| 夫妻性生交免费视频一级片| 这个男人来自地球电影免费观看| 久久九九热精品免费| 午夜老司机福利片| 成人影院久久| 亚洲午夜精品一区,二区,三区| 国产在线一区二区三区精| 中文字幕人妻丝袜一区二区| 男女免费视频国产| 99精品久久久久人妻精品| 伦理电影免费视频| 久久精品亚洲熟妇少妇任你| 国产欧美日韩一区二区三区在线| 亚洲人成网站在线观看播放| 久久人妻福利社区极品人妻图片 | 欧美 亚洲 国产 日韩一| 欧美日韩亚洲高清精品| kizo精华| 亚洲国产精品一区二区三区在线| 亚洲五月婷婷丁香| 亚洲欧洲日产国产| 在现免费观看毛片| 丁香六月天网| 日本黄色日本黄色录像| av网站免费在线观看视频| 十分钟在线观看高清视频www| 国产成人av激情在线播放| 久久精品亚洲熟妇少妇任你| 亚洲图色成人| 国产一卡二卡三卡精品| 叶爱在线成人免费视频播放| 国产色视频综合| 久久ye,这里只有精品| 午夜91福利影院| 午夜精品国产一区二区电影| 人妻一区二区av| 色网站视频免费| 国产免费福利视频在线观看| 午夜福利,免费看| 天堂8中文在线网| 国产精品一区二区精品视频观看| 午夜影院在线不卡| 香蕉国产在线看| 亚洲av国产av综合av卡| 午夜av观看不卡| 一边亲一边摸免费视频| 丝袜美足系列| 18禁观看日本| 久久精品国产亚洲av涩爱| 1024视频免费在线观看| 精品一区二区三区四区五区乱码 | 一级a爱视频在线免费观看| 少妇精品久久久久久久| 国产av精品麻豆| 日本vs欧美在线观看视频| 日日爽夜夜爽网站| 蜜桃国产av成人99| 热99久久久久精品小说推荐| 51午夜福利影视在线观看| 日本vs欧美在线观看视频| 精品高清国产在线一区| 亚洲精品国产区一区二| 另类亚洲欧美激情| 久久人妻熟女aⅴ| 日本vs欧美在线观看视频| 中文字幕高清在线视频| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区| 国产精品亚洲av一区麻豆| 欧美日韩国产mv在线观看视频| 婷婷成人精品国产| 97人妻天天添夜夜摸| 91国产中文字幕| 亚洲精品一二三| 一区二区三区乱码不卡18| 亚洲国产av新网站| 国产伦人伦偷精品视频| 欧美在线黄色| 欧美成人午夜精品| 搡老乐熟女国产| 菩萨蛮人人尽说江南好唐韦庄| 美女脱内裤让男人舔精品视频| 男人添女人高潮全过程视频| 日韩一区二区三区影片| 丁香六月欧美| 又粗又硬又长又爽又黄的视频| 只有这里有精品99| 亚洲精品美女久久av网站| 在线观看一区二区三区激情| 亚洲专区国产一区二区| 男男h啪啪无遮挡| 国产亚洲午夜精品一区二区久久| 午夜av观看不卡| 最近最新中文字幕大全免费视频 | 少妇人妻久久综合中文| 这个男人来自地球电影免费观看| 日韩欧美一区视频在线观看| 校园人妻丝袜中文字幕| 下体分泌物呈黄色| 一区二区三区精品91| 99国产综合亚洲精品| 国产男女超爽视频在线观看| bbb黄色大片| 男男h啪啪无遮挡| 中文字幕精品免费在线观看视频| 人妻人人澡人人爽人人| 国产精品一区二区免费欧美 | 欧美人与性动交α欧美精品济南到| 久久免费观看电影| 成年美女黄网站色视频大全免费| www.999成人在线观看| 久久久精品免费免费高清| 中文乱码字字幕精品一区二区三区| 热re99久久国产66热| 亚洲欧美日韩高清在线视频 | 亚洲专区国产一区二区| 中文乱码字字幕精品一区二区三区| av国产久精品久网站免费入址| 亚洲,一卡二卡三卡| 亚洲熟女毛片儿| kizo精华| 男女免费视频国产| 国产精品国产三级专区第一集| 日本wwww免费看| 亚洲综合色网址| 国产高清国产精品国产三级| 国产1区2区3区精品| 男人爽女人下面视频在线观看| 在线观看人妻少妇| 免费黄频网站在线观看国产| 午夜免费成人在线视频| 好男人电影高清在线观看| 亚洲精品国产av蜜桃| 美女视频免费永久观看网站| 国产亚洲一区二区精品| 黄网站色视频无遮挡免费观看| 日韩 亚洲 欧美在线| 亚洲欧美中文字幕日韩二区| 欧美人与性动交α欧美软件| 夜夜骑夜夜射夜夜干| 一级毛片我不卡| 精品少妇内射三级| 久久国产亚洲av麻豆专区| 亚洲精品av麻豆狂野| 天堂8中文在线网| 中文字幕最新亚洲高清| 国产视频首页在线观看| 国产亚洲午夜精品一区二区久久| 熟女少妇亚洲综合色aaa.| 女人精品久久久久毛片| 香蕉国产在线看| 亚洲精品自拍成人| 欧美精品人与动牲交sv欧美| 久久久亚洲精品成人影院| 国产福利在线免费观看视频| 51午夜福利影视在线观看| 国产熟女欧美一区二区| 叶爱在线成人免费视频播放| 天天添夜夜摸| 夫妻性生交免费视频一级片| 亚洲成色77777| 亚洲av电影在线进入| 亚洲男人天堂网一区| 亚洲久久久国产精品| 久久亚洲国产成人精品v| 91精品三级在线观看| 久久鲁丝午夜福利片| 午夜福利免费观看在线| 国产男女超爽视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美激情高清一区二区三区| 亚洲九九香蕉| 少妇裸体淫交视频免费看高清 | 国产福利在线免费观看视频| a级毛片在线看网站| 最新的欧美精品一区二区| 视频在线观看一区二区三区| 永久免费av网站大全| 国产伦人伦偷精品视频| 欧美日韩亚洲综合一区二区三区_| 成人午夜精彩视频在线观看| 在线观看人妻少妇| 午夜免费观看性视频| 欧美人与善性xxx| 免费高清在线观看日韩| 人成视频在线观看免费观看| 青草久久国产| 久久久久久亚洲精品国产蜜桃av| av网站在线播放免费| 久久热在线av| 99久久综合免费| 天天躁夜夜躁狠狠躁躁| 国产在线观看jvid| 18禁国产床啪视频网站| 亚洲av电影在线进入| a级毛片在线看网站| 日韩伦理黄色片| av又黄又爽大尺度在线免费看| 国产精品一区二区在线不卡| 男男h啪啪无遮挡| 国产精品一区二区精品视频观看| 激情视频va一区二区三区| 精品亚洲成a人片在线观看| 1024香蕉在线观看| 美女脱内裤让男人舔精品视频| 色精品久久人妻99蜜桃| 国产精品秋霞免费鲁丝片| 国产精品 欧美亚洲| 啦啦啦在线免费观看视频4| 男人爽女人下面视频在线观看| 少妇粗大呻吟视频| 久久精品熟女亚洲av麻豆精品| 精品少妇黑人巨大在线播放| 国产免费福利视频在线观看| 波多野结衣一区麻豆| av网站在线播放免费| 久久国产亚洲av麻豆专区| 777久久人妻少妇嫩草av网站| 飞空精品影院首页| 成年动漫av网址| 国产精品.久久久| 亚洲精品自拍成人| av片东京热男人的天堂| 9色porny在线观看| 女性被躁到高潮视频| 国产高清视频在线播放一区 | 天堂俺去俺来也www色官网| 91成人精品电影| 高清视频免费观看一区二区| 国产伦人伦偷精品视频| 别揉我奶头~嗯~啊~动态视频 | 久久九九热精品免费| 可以免费在线观看a视频的电影网站| 9色porny在线观看| 国产片特级美女逼逼视频| 三上悠亚av全集在线观看| 国产又色又爽无遮挡免| 国产极品粉嫩免费观看在线| 天堂俺去俺来也www色官网| 老司机影院毛片| 久久久久视频综合| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 菩萨蛮人人尽说江南好唐韦庄| 在现免费观看毛片| 赤兔流量卡办理| 国产精品久久久久久人妻精品电影 | videos熟女内射| 女性生殖器流出的白浆| 国产黄频视频在线观看| 男女高潮啪啪啪动态图| 免费在线观看视频国产中文字幕亚洲 | 观看av在线不卡| 青春草亚洲视频在线观看| 国产亚洲欧美精品永久| 色视频在线一区二区三区| 亚洲欧美一区二区三区国产| 丝袜喷水一区| 90打野战视频偷拍视频| 观看av在线不卡| 50天的宝宝边吃奶边哭怎么回事| 久久久久国产精品人妻一区二区| 搡老岳熟女国产| 中文字幕人妻丝袜制服| 少妇人妻 视频| 久久国产亚洲av麻豆专区| 欧美亚洲日本最大视频资源| 国产精品国产av在线观看| 久久99一区二区三区| 久久久久久亚洲精品国产蜜桃av| 亚洲中文日韩欧美视频| 老司机靠b影院| 国产精品二区激情视频| 亚洲精品日韩在线中文字幕| 中文精品一卡2卡3卡4更新| 婷婷色综合www| 精品亚洲乱码少妇综合久久| 成人免费观看视频高清| 美女主播在线视频| av国产精品久久久久影院| 建设人人有责人人尽责人人享有的| 亚洲人成网站在线观看播放| 国产成人精品无人区| 亚洲综合色网址| 国产片特级美女逼逼视频| 久久国产精品男人的天堂亚洲| 色网站视频免费| 亚洲免费av在线视频| 王馨瑶露胸无遮挡在线观看| 国产老妇伦熟女老妇高清| 在线观看一区二区三区激情| 男女之事视频高清在线观看 | 中文字幕色久视频| 无限看片的www在线观看| 国产成人系列免费观看| 精品国产一区二区三区久久久樱花| 国产成人91sexporn| 午夜日韩欧美国产| 亚洲图色成人| 国产高清不卡午夜福利| 日本vs欧美在线观看视频| 国语对白做爰xxxⅹ性视频网站| 亚洲国产精品国产精品| 久久久精品免费免费高清| 最新的欧美精品一区二区| 9色porny在线观看| 欧美成人精品欧美一级黄| 欧美激情极品国产一区二区三区| 国产极品粉嫩免费观看在线| 午夜日韩欧美国产| 欧美精品高潮呻吟av久久| 永久免费av网站大全| 久久久久精品国产欧美久久久 | 国产精品一国产av| 国精品久久久久久国模美| 乱人伦中国视频| 人体艺术视频欧美日本| 午夜激情久久久久久久| 下体分泌物呈黄色| 宅男免费午夜| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 久久国产精品男人的天堂亚洲| 日韩一本色道免费dvd| 国产精品一区二区在线不卡| 少妇粗大呻吟视频| 叶爱在线成人免费视频播放| 久久久国产一区二区| 天堂俺去俺来也www色官网| 九草在线视频观看| 日韩大片免费观看网站| 亚洲欧美中文字幕日韩二区| 尾随美女入室| 最黄视频免费看| 一本综合久久免费| avwww免费| 久久精品亚洲av国产电影网| 免费在线观看日本一区| 精品亚洲成a人片在线观看| 亚洲av欧美aⅴ国产| 中国国产av一级| 日韩制服丝袜自拍偷拍| 国产成人精品无人区| 青青草视频在线视频观看| 亚洲精品一区蜜桃| 一级,二级,三级黄色视频| 欧美日韩视频高清一区二区三区二| 纯流量卡能插随身wifi吗| 最近中文字幕2019免费版| 大片电影免费在线观看免费| 啦啦啦啦在线视频资源| 久久女婷五月综合色啪小说| 婷婷丁香在线五月| 亚洲三区欧美一区| 99国产综合亚洲精品| 国产视频首页在线观看| 嫩草影视91久久| 满18在线观看网站| 50天的宝宝边吃奶边哭怎么回事| 91成人精品电影| 91国产中文字幕| 国产在视频线精品| 亚洲男人天堂网一区| 最新在线观看一区二区三区 | 精品国产一区二区三区久久久樱花| 又黄又粗又硬又大视频| 亚洲黑人精品在线| 狂野欧美激情性bbbbbb| 久久人妻熟女aⅴ| 亚洲欧美日韩高清在线视频 | 欧美精品一区二区大全| 黄色毛片三级朝国网站| 午夜免费观看性视频| 少妇人妻 视频| 欧美成人午夜精品| 丝袜美腿诱惑在线| 成人影院久久| av在线老鸭窝| 亚洲一卡2卡3卡4卡5卡精品中文| 五月天丁香电影| 久久精品国产亚洲av涩爱| 又粗又硬又长又爽又黄的视频| 妹子高潮喷水视频| 亚洲国产精品一区二区三区在线| 精品一区二区三区av网在线观看 | 欧美97在线视频| bbb黄色大片| 国产精品国产三级国产专区5o| 肉色欧美久久久久久久蜜桃| 久久午夜综合久久蜜桃| 国产精品麻豆人妻色哟哟久久| 亚洲av国产av综合av卡| 中文字幕精品免费在线观看视频| 欧美日本中文国产一区发布| 大片电影免费在线观看免费| 亚洲图色成人| 久久人妻福利社区极品人妻图片 | 欧美黑人精品巨大| 亚洲情色 制服丝袜| 国产精品.久久久| 午夜福利影视在线免费观看| 免费观看人在逋| 99国产精品一区二区三区| 国产成人精品久久久久久| 欧美 亚洲 国产 日韩一| 肉色欧美久久久久久久蜜桃| 国产精品久久久久成人av| 国产成人精品久久久久久| 老司机亚洲免费影院| 国产免费福利视频在线观看| 久久青草综合色| 国产成人一区二区在线| 免费久久久久久久精品成人欧美视频| 在线观看人妻少妇| 久久这里只有精品19| 国产高清视频在线播放一区 | 色网站视频免费| 久久国产亚洲av麻豆专区| 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲日产国产| 精品国产一区二区三区四区第35| 操美女的视频在线观看| 午夜视频精品福利| 男人操女人黄网站| 在线观看www视频免费| 这个男人来自地球电影免费观看| 曰老女人黄片| 亚洲一码二码三码区别大吗|