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

    三峽水庫運行期間原生動物群落的時空異質(zhì)性

    2015-02-07 06:51:05鄭金秀池仕運汪紅軍胡菊香
    生態(tài)學報 2015年11期
    關鍵詞:原生動物纖毛蟲水區(qū)

    鄭金秀, 池仕運, 李 聃, 汪紅軍, 胡菊香

    水利部中國科學院水工程生態(tài)研究所, 水利部水工程生態(tài)效應與生態(tài)修復重點實驗室, 武漢 430079

    三峽水庫運行期間原生動物群落的時空異質(zhì)性

    鄭金秀, 池仕運, 李 聃, 汪紅軍, 胡菊香*

    水利部中國科學院水工程生態(tài)研究所, 水利部水工程生態(tài)效應與生態(tài)修復重點實驗室, 武漢 430079

    于2010年10月—2011年6月三峽水庫正常運行周期內(nèi)對庫區(qū)干流原生動物進行調(diào)查,研究其空間分布及水庫周期排蓄期間的變化。共檢測到原生動物99種,蓄水后纖毛蟲有增多的趨勢。水庫運行的不同階段優(yōu)勢種不同,大致演變?yōu)椋荷皻だw毛蟲(蓄水期)—非砂殼類纖毛蟲(高位運行期間)—有殼肉足蟲(低位運行期間)。不同水域優(yōu)勢種也存在差異,從上游到下游特征指示種變化為:有殼肉足蟲(變動回水區(qū))—纖毛蟲(湖泊區(qū))。結果表明,三峽水庫原生動物的分布具有明顯的時空異質(zhì)性(P<0.05),影響原生動物時空分布的主要因素有透明度、溫度、電導率和葉綠素a。原生動物平均密度為952.19 個/L,平均生物量為8.14 μg/L。蓄水期上游現(xiàn)存量高于下游,低位運行期間則低于下游。原生動物Marglef和Shannonn-Weiver多樣性指數(shù)平均值分別為3.78和2.18,1月份最低,6月份最高。蓄水175 m后上游變動回水區(qū)原生動物具有較高的豐度。水庫冬蓄夏排的運行模式模糊了河流本身的季節(jié)變化,使原生動物的種類和現(xiàn)存量更多的受水庫水動力學的影響,使水體理化因子和水文特征呈現(xiàn)明顯的時空差異,最終形成原生動物種群的時空異質(zhì)性分布。

    三峽; 原生動物; 運行; 時空異質(zhì)性

    三峽水庫2008年進行175 m試驗性蓄水,并于2010年10月26日首次蓄水到175 m設計目標,至此,三峽水庫年水位落差達30 m。在175 m蓄水水位時,三峽干流回水末端到達重慶。變動回水區(qū)從黃草峽至江津紅花磧,變動回水區(qū)長約146 km。變動回水區(qū)的水體滯留時間進一步延長,水體逐漸從河流型轉變?yōu)檫^渡型。三峽水庫“冬蓄夏排”的周期性蓄水模式,不僅使干流不同區(qū)域形成不同類型水體,而且模糊了河流本身的季節(jié)變化,使三峽水庫浮游生物在季節(jié)演變上呈現(xiàn)出與湖泊不同的規(guī)律[1]。三峽水庫蓄水后,眾多學者對水庫的生態(tài)環(huán)境特征[2- 4]、浮游生物群落結構[5- 7]、富營養(yǎng)化問題[8- 9]等開展了研究,但是對三峽水庫原生動物的研究較少,龔循舉等人在蓄水前對三峽地區(qū)的原生動物區(qū)系做了研究[10],蓄水后筆者在三峽庫區(qū)做了初步調(diào)查[11]。此外,蓄水后僅見三峽庫區(qū)小江支流的報道[12]。

    本文調(diào)查期間正處于三峽水庫正常運行初期,水生態(tài)系統(tǒng)仍處于不穩(wěn)定階段。變動回水區(qū)范圍改變,水位波動幅度增大,都將對水庫水生態(tài)系統(tǒng)造成影響。原生動物作為水生態(tài)系統(tǒng)食物鏈中的重要環(huán)節(jié),對環(huán)境的細微變化敏感,對水生態(tài)系統(tǒng)變化響應迅速。通過對原生動物空間分布及水庫周期排蓄期間的變化研究,可以反映三峽水庫水生態(tài)系統(tǒng)的狀況及其發(fā)展趨勢。本文研究不僅旨在探討三峽水庫運行對原生動物種群結構的影響,了解原生動物的演變過程,為闡明三峽水庫獨特的水生態(tài)系統(tǒng)的變化提供有力支持,也希望能為三峽庫區(qū)的原生動物分布狀況積累基礎數(shù)據(jù),為生態(tài)環(huán)境保護提供參考。

    1 材料與方法

    1.1 樣品采集

    從上游永川起至下游秭歸,即從175 m回水區(qū)末端至壩前,共設置13個監(jiān)測斷面,依水面寬度及實際采樣點狀況設左右岸采樣點。三峽水庫正常蓄水后,枯水期175 m水位回水末端位于江津附近的紅花磧,汛期145 m水位回水末端位于涪陵附近。依照水庫不同蓄水水位所達到的回水區(qū)域,可分為變動回水區(qū)和湖泊區(qū)兩組,即從水位145 m回水區(qū)末端至175 m回水區(qū)末端區(qū)域為變動回水區(qū),該區(qū)域由于水位周期變動使生境處于河流-湖泊交替變換狀態(tài),包括以下采樣點:朱沱(ZT)、魚洞(YD)、寸灘(CT)、扇沱(ST)、清溪場(QX);湖泊區(qū)為水庫永久的淹沒區(qū),不隨水庫調(diào)度周期變化,主要包括:蘇家左(SJL)、蘇家右(SJR)、桐園左(TYL)、桐園右(TYR)、苦草沱左(KCTL)、苦草沱右(KCTR)、白帝城左(BDCL)、白帝城右(BDCR)、培石(PS)、黃蠟石(HLS)、屈原祠(QYC)、銀杏沱(YXT)(圖1)。

    根據(jù)水庫蓄水運行周期,設置3次調(diào)查時間,分別為2010年10月(蓄水期),水庫蓄水達到175 m,調(diào)查期間水位波動在172—174 m;2011年1月(高位運行期),水庫在高水位運行,水位維持在174 m左右;2011年6月(低位運行期),水位降為148 m。

    圖1 三峽采樣點示意圖

    1.2 樣品處理

    原生動物定性樣品用25號浮游生物網(wǎng)(孔徑為64 μm)大范圍拖撈獲得,取50 mL用5%的福爾馬林固定。河流定量樣品用2.5 L采水器采水10 L,混合后取2 L水立即用15%的魯哥氏液固定。三峽水庫庫區(qū),定量樣品用2.5 L采水器于5個深度(0、5、10、20、40 m)各采水5 L,混合后取2 L水用15%的魯哥氏液固定。所獲樣品均帶回實驗室進行分類、鑒定[13- 14]。

    用YSI6600水質(zhì)分析儀同步監(jiān)測水質(zhì)參數(shù),包括:水溫(T)、pH值、溶解氧(DO)、電導率(Cond)、氧化還原電位(ORP);用賽氏盤測透明度(SD);采集水樣,抽濾保存,回實驗室測定Chl a等。

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

    原生動物密度用0.1 mL計數(shù)框,視野法計算;生物量按體積估算。具體計算參照章宗涉和黃祥飛的方法[15]。

    原生動物多樣性分析采用Marglef指數(shù)(d)和Shannonn-Wiever指數(shù)(H′)計算。計算公式如下:

    d= (S-1) / lnN

    式中,S為原生動物種類數(shù),N為原生動物個體數(shù)。

    H′ =-∑Pi×log2Pi,Pi=ni/N

    式中,Pi為樣品中第i種(屬)個體所占比例;N為采集樣品中原生動物的總個體數(shù);ni為第i種個體數(shù)。

    采用除趨勢對應分析(DCA)對原生動物時空分布進行排序,并采用檢驗組間差異性的非參數(shù)程序多響應置換過程分析(MRPP)法進行不同條件之間的差異性檢驗。再將DCA軸與環(huán)境因子做相關性分析,探討影響原生動物分布的主要環(huán)境因子。

    指示種分析采用指示種分析(ISA)檢測和描述站點原生動物對不同條件下的指示作用,并用蒙特卡羅檢驗來決定指示值(IV)的統(tǒng)計顯著性。

    采用單因素方差分析(One-way ANOVA)對變動回水區(qū)和湖泊區(qū)原生動物的分布進行差異性分析,多重比較(LSD)分析法檢驗數(shù)據(jù)的差異性。

    DCA分析、MRPP分析、ISA分析皆用PC-ORD軟件進行。One-way ANOVA分析和DCA軸與環(huán)境因子的相關性分析用SPSS17.0軟件進行。

    2 結果與分析

    2.1 種類組成

    調(diào)查期間共檢出原生動物99種,其中纖毛蟲50種。10月份蓄水期檢出49種,纖毛蟲占49.98%;1月份高位運行期間檢出50種,纖毛蟲占56.86%;6月份低位運行期間檢出66種,纖毛蟲占47.06%。總體來說,肉足蟲與纖毛蟲種數(shù)相當,在不同運行期間比例略有變化,蓄水后纖毛蟲有增多的趨勢,因而1月份高位運行期間纖毛蟲種類多于另兩個運行階段。

    原生動物在水庫運行的不同階段優(yōu)勢種組成不同,蓄水期常見種主要有表殼蟲(Arcella)、砂殼蟲(Difflugia)等有殼肉足蟲、旋回俠盜蟲(Strobilidiumgyrans)和筒殼蟲(Tintinnidium)、似鈴殼蟲(Tintinnopsis)等砂殼纖毛蟲,優(yōu)勢種為旋回俠盜蟲和小筒殼蟲(T.pusillum),高位運行期間其它種類纖毛蟲增多,膜袋蟲(Cyclidium)和鐘蟲(Vorticella)成為常見種,其優(yōu)勢種主要有旋回俠盜蟲、小筒殼蟲和膜袋蟲;低位運行期間則以有殼肉足蟲、砂殼纖毛蟲和旋回俠盜蟲為常見種,優(yōu)勢種是旋回俠盜蟲、表殼蟲、匣殼蟲(Centropyxis)和圓殼蟲(Cyclopyxis)。旋回俠盜蟲是全年的優(yōu)勢種。為辨別3個時期種類組成上的差異,用PC-ORD軟件計算分析3個時期的指示種,取IV(indicator value)值較大且達到顯著水平(P<0.05)的為該時期的特征指示種(表1)??梢钥闯霾煌瑫r期原生動物特征種類的演替,大致為:砂殼纖毛蟲(蓄水期)—非砂殼類纖毛蟲(高位運行期間)—有殼肉足蟲(低位運行期間)。

    表1 不同時期原生動物優(yōu)勢種

    原生動物在不同水域種類組成有較大差異,計算分析變動回水區(qū)和湖泊區(qū)的主要特征指示種(表2)。可以看出變動回水區(qū)原生動物特征指示種主要是有殼肉足蟲,下游湖泊區(qū)纖毛蟲種類增多,并以巧刺日蟲(Raphidiophryselegans)和陀螺俠盜蟲(Strobilidiumvelox)為該水域的特征指示種。

    表2 不同水域原生動物優(yōu)勢種

    2.2 密度和生物量

    干流原生動物平均密度為952.19 個/L,10月份密度754.41 個/L,1月份密度1138.25 個/L,6月份密度963.92 個/L;干流平均生物量為8.14 μg/L,10月份生物量為 6.65 μg/L,1月份為7.95 μg/L,6月份為9.80 μg/L。不同站點原生動物現(xiàn)存量如圖2和圖3。蓄水期上游變動回水區(qū)原生動物現(xiàn)存量明顯高于下游湖泊區(qū),但低位運行期間湖泊區(qū)的現(xiàn)存量要比變動回水區(qū)高許多。高位運行期,變動回水區(qū)密度為612 個/L,遠低于下游湖泊區(qū)1288 個/L,但生物量卻比湖泊區(qū)高,分別為12.43 μg/L和6.28 μg/L。這是由于高位運行期間湖泊區(qū)小個體纖毛蟲數(shù)量增多,如睥睨蟲(體長30—35 μm)、膜袋蟲(體長22—30 μm)、蚤中縊蟲(Mesodiniumpulex)(體長20—25 μm)、櫛毛蟲(Didinium)(體長25—60 μm)等大量出現(xiàn),但變動回水區(qū)仍以個體稍大的原生動物為主,如表殼蟲(殼長70—150 μm)、睫杵蟲(Ophryoglena)(體長120—160 μm)、長圓靴纖蟲(Cothurniaoblonga)(鞘長72—80 μm)等。通過One-way ANOVA分析,變動回水區(qū)和湖泊區(qū)原生動物密度和生物量都達到差異極顯著水平(P<0.01)。

    圖2 原生動物密度變化

    圖3 原生動物生物量變化

    2.3 理化因子變化及與原生動物的相關性

    三峽庫區(qū)水體理化因子檢測結果如表3所示。從表中可看出,除了pH值變化較小外,其余水質(zhì)參數(shù)波動較大,從上游到下游也有很明顯的趨勢走向。其中透明度、電導率、溫度、Chl a呈上升趨勢,而溶解氧、氧化還原電位呈下降趨勢。各指標在清溪場和蘇家之間有明顯躍動界限。清溪場以上為變動回水區(qū),蘇家以下則是湖泊區(qū),可見兩個區(qū)域的水質(zhì)特征有明顯的差異。通過One-way ANOVA分析證明,這兩個區(qū)域的差異達到極顯著水平(P<0.01)。葉綠素a含量是表征藻類現(xiàn)存量的重要指標,本文檢測出葉綠素a除了苦草沱斷面在夏季>10 μg/L,為富營養(yǎng)狀態(tài)之外,其余斷面都處于貧營養(yǎng)-中營養(yǎng)水平,且大部分為貧營養(yǎng)狀態(tài)。

    表3 三峽水庫水體理化因子

    三峽庫區(qū)水體理化因子與原生動物的相關性如表4。種類數(shù)和多樣性指數(shù)受透明度、電導率、溫度共同影響,呈顯著相關(P<0.05,P<0.01)。且在湖泊區(qū),種類數(shù)與透明度、溫度和電導率的相關性更為顯著(P<0.01)。纖毛蟲比例與透明度成正相關(P<0.01),與溶解氧成負相關(P<0.05)。而原生動物現(xiàn)存量與理化因子關系不明顯。葉綠素a對原生動物的影響不明顯,只在過渡區(qū)與種類數(shù)成負相關(P<0.05)。

    表4 水體理化因子與原生動物的相關性

    (1)表示過渡區(qū);(2)表示湖泊區(qū)

    2.4 原生動物多樣性

    計算Marglef指數(shù)(d)和Shannonnn-wiever指數(shù)(H′),平均值分別為3.78和2.18。兩種多樣性指數(shù)變化趨勢相同,如圖4。Marglef指數(shù)以1月份高位運行期指數(shù)值最低,為3.46,6月份低位運行期最高,為4.22。指數(shù)在不同水域不同季節(jié)變化較明顯,10月份和6月份之間達到差異極顯著水平(P<0.01)。蓄水期和高位運行期變動回水區(qū)的指數(shù)值較高,但低位運行期湖泊區(qū)的marglef指數(shù)遠高于變動回水區(qū)。Shannonn-Wiener指數(shù)(H′)也是1月份最低,6月份最高,但季節(jié)之間差異不大(P>0.05)。蓄水期不同水域之間H′指數(shù)也差別不大(P>0.05)。

    2.5 原生動物在水庫不同運行時期的分布

    根據(jù)原生動物檢出情況,做DCA分析,結果如圖5。進行MRPP分析,統(tǒng)計量T=-9.765,A=0.059,P<0.001,達到極顯著水平,組間差異分析結果如表5,P值皆小于0.05,表明水庫不同運行時期的原生動物分布差異顯著。

    圖4 不同水域多樣性指數(shù)季節(jié)變化

    表5 原生動物季節(jié)變化的MRPP分析結果

    圖5 原生動物季節(jié)變化的DCA分析

    2.6 原生動物的空間分布

    根據(jù)不同水域原生動物的種類組成,做不同站點的DCA分析,結果如圖6。再進行MRPP分析檢驗差異性,統(tǒng)計量T=-6.1642561,A=0.02614683,P=0.00001544,達到極顯著水平。表明三峽水庫變動回水區(qū)與湖泊區(qū)之間原生動物種類組成和分布差異顯著。

    圖6 原生動物空間分布的DCA分析

    2.7 理化因子與DCA軸的相關性

    將理化參數(shù)與DCA軸數(shù)值做Spearman相關性檢驗,相關系數(shù)如表6??梢?,影響原生動物時空分布的主要因素主要有透明度、溫度、電導率和Chl a。軸一主要反映透明度的變化;軸二與Chl a顯著相關;軸三與溫度和電導率二者相關性顯著。 總方差4.3493,第三軸累積貢獻率為:(0.3194150031+0.2796075940+0.2058645636)/4.3493=18.5%

    表6 理化參數(shù)與DCA軸的相關性

    3 討論

    3.1 三峽水庫原生動物時空分布的異質(zhì)性

    水庫特有的形態(tài)結構和運行特性導致水庫在物理、化學上均存在一個縱向分區(qū),在生境上表現(xiàn)出由激流環(huán)境到靜水環(huán)境的過渡[16],也使得浮游動物在分布上具有顯著差異。

    隨著三峽水庫的運行,水文節(jié)律發(fā)生改變,原生動物種群結構呈現(xiàn)明顯的季節(jié)性差異。蓄水期水位上升,直至接近175 m,淹沒區(qū)擴大,上游江段形成新的回水區(qū),水流明顯減緩。蓄水后水體滯留時間延長,原生動物密度和生物量上升較快。此時的生境還不穩(wěn)定,原生動物主要以有殼肉足蟲、砂殼纖毛蟲和旋回俠盜蟲為主。水庫蓄水175 m后,在該水位運行數(shù)月,至第2年1至2月。第2次采樣正處于水庫高位運行期間,水庫生態(tài)系統(tǒng)較蓄水期穩(wěn)定,原生動物類群發(fā)生明顯變化,纖毛蟲比例上升,并成為優(yōu)勢種,密度和生物量也較蓄水期增多。6月份為三峽水庫低位運行期間,水位148 m,此時回水區(qū)至萬州附近,萬州上游江段重新處于河流狀態(tài)。從物種組成可以看出,此時原生動物表殼科和砂殼科肉足蟲明顯增多,并成為優(yōu)勢種。低位運行期間原生動物密度有所下降,但生物量繼續(xù)上升,可見大個體的種類比例上升而小型種類減少。

    從空間分布上看,變動回水區(qū)和湖泊區(qū)之間原生動物的種群分布差異較大,三峽水庫正常運行后,上游從朱沱至萬州由原來的河流區(qū)轉變?yōu)檫^渡區(qū),隨著水庫運行階段的不同,其生境不斷改變,使該水域原生動物種群也處在動態(tài)變化過程,并與下游湖泊區(qū)呈明顯差異。

    3.2 環(huán)境因子對原生動物種群結構的影響

    環(huán)境因子(如水溫、透明度等)對浮游動物群落演替具有重要作用[17]。本研究中,溫度、透明度、電導率和葉綠素a是影響原生動物時空異質(zhì)性的主要因素。

    水溫是三峽水庫環(huán)境特征的控制性因子一[2]。水溫以各種形式影響水庫水動力、理化和生物過程[18]。原生動物最適溫度范圍多在10—25 ℃之間[19]。在一定溫度范圍內(nèi),原生動物種群增長率r與水溫呈正相關[20]。三峽水庫湖泊區(qū)原生動物種類數(shù)與水溫相關性極顯著(P<0.01),且由于湖泊區(qū)范圍遠大于變動回水區(qū)范圍,使得水溫的影響在原生動物結構形成過程中起了重要作用。

    除水溫外,透明度、電導率等理化條件也呈現(xiàn)縱向梯度,共同影響了三峽水庫的生態(tài)過程[21- 23]。電導率在不同水域間的差異可能和流域的面積有關,電導率一般隨流域面積的增大而上升[24]。在同一水域電導率的時間變化則可能是溫度及其它環(huán)境因子決定的。三峽水庫6月夏季水溫較高,相應的溶解度大,而可溶性的離子越多,電導率越大。但在上游過渡區(qū),電導率并沒有隨溫度升高而上升,可能和上游蓄水,環(huán)境發(fā)生改變有關。營養(yǎng)鹽主要來自徑流的輸送,蓄水期上游輸送的營養(yǎng)鹽及淹沒土地釋放出的營養(yǎng)鹽,都會增大水體的電導率。多樣性指數(shù)是群落演替的一個重要指標,群落多樣性的變化是生物生存環(huán)境發(fā)生變化的重要標志。三峽原生動物多樣性指數(shù)受透明度、電導率、溫度共同影響,成顯著相關(P<0.05)。

    葉綠素a大體上能反映水體的營養(yǎng)水平。有研究表明,原生動物群落演變與水體營養(yǎng)水平關系密切[25- 27]。本文研究,原生動物的種類分布受葉綠素a影響,但密度和生物量隨葉綠素a的消長不明顯,其種群結構應是受生物和非生物因素的共同作用的結果。

    在成庫初期,水力條件變動是直接影響浮游生物分布的主導因素[28]。三峽水庫上游過渡區(qū)的水力、水環(huán)境條件尚處于變動期,水力條件波動及其引起的水環(huán)境改變對該水域的原生動物分布影響較大;而下游庫區(qū)已成庫多年,原生動物經(jīng)過一定時間的適應,基本上已形成自身的分布和演變規(guī)律,受環(huán)境因子和生物因子的影響更為明顯。因此,是生物因素和非生物因素綜合作用,共同影響了三峽水庫原生動物的時空分布。

    3.3 水庫運行對原生動物種群結構的影響

    一般來說,水庫成庫后要經(jīng)過4—5年的波動后水生生物群落結構才趨于穩(wěn)定[29]。水庫運行的特點之一就是對流量的調(diào)節(jié)。受徑流的影響,河流區(qū)原生動物的種類和數(shù)量相對較少。由河流向水庫過渡后,由于水力滯留時間的延長,流速減緩,環(huán)境趨于穩(wěn)定,給浮游動物提供了良好的棲息環(huán)境,因此,庫區(qū)內(nèi)浮游動物的種類及現(xiàn)存量上要高于上游原河流區(qū)(即過渡區(qū)朱沱至清溪江段)。然而,水庫過渡區(qū)浮游動物卻具有較高豐度[30]。從圖2,圖3可知,蓄水期,河流區(qū)轉變?yōu)檫^渡區(qū)后原生動物豐度明顯高于中下游庫區(qū)。這可能是由蓄水期入庫水流攜帶的大量有機物質(zhì)造成的。當水大量流入水庫時,由于水動力學原理以及懸浮物的濃度的增加,會引起浮游動物種群的波動[31- 32]。比較而言,原生動物中的有殼肉足蟲比纖毛蟲對懸浮物更為敏感[33],因此,過渡區(qū)原生動物的種類組成中,有殼肉足蟲占有更大的比重。此外,外來的懸浮物還帶來了大量碳源,對過渡區(qū)的食物網(wǎng)造成影響,進而影響了原生動物的種群結構[30]。黃海強[34]對蓄水后三峽庫區(qū)浮游生物與泥沙沉降的關系進行研究,認為泥沙沉降有利于浮游生物的生長。Kirt的研究[35]表明,懸浮泥沙可以改變浮游動物的群落結構。

    除了形成空間上的不同分區(qū)外,水庫流量調(diào)節(jié)在洪水季蓄水削峰和枯水季補充水量也對原生動物造成了較大影響。通過調(diào)節(jié),三峽水庫可以對營養(yǎng)物質(zhì)、藻類組成和現(xiàn)存量等進行調(diào)控[36],從而影響原生動物的結構變化。與丹江口水庫[37],陸水水庫[38]等供水水庫相比,三峽水庫浮游動物在季節(jié)變化上并沒有遵循夏秋季高峰,冬季較少的規(guī)律。三峽水庫原生動物在1月份枯水季的現(xiàn)存量高于10月份平水季節(jié),1月份原生動物的密度甚至比6月份豐水季節(jié)的密度高??梢?,三峽水庫運行,對水量的調(diào)節(jié)使水體在較短時期內(nèi)發(fā)生波動,模糊了河流本身的季節(jié)變化。吳生桂等人也認為,水庫調(diào)度蓄洪補枯影響了原生動物的季節(jié)變化規(guī)律[39]。水庫運行,通過改變水動力學條件,進而影響水體物理化學過程,使水體理化因子和水文特征呈現(xiàn)明顯的時空差異,最終形成浮游生物群落的時空異質(zhì)性分布。

    [1] 王英才, 邱光勝, 陳水松, 胡圣. 三峽庫區(qū)試驗性蓄水期間浮游生物群落特點研究. 人民長江, 2012, 43(12): 4- 9.

    [2] 羅專溪, 張遠, 鄭丙輝, 富國, 陸兆華. 三峽水庫蓄水初期水生態(tài)環(huán)境特征分析. 長江流域科學與環(huán)境, 2005, 14(6): 781- 785.

    [3] 胡征宇, 蔡慶華. 三峽水庫蓄水前后水生態(tài)系統(tǒng)動態(tài)的初步研究. 水生生物學報, 2006, 30(1): 1- 6.

    [4] 胡念三, 劉德富, 紀道斌, 楊正健. 三峽水庫干流倒灌對支流庫灣營養(yǎng)鹽分布的影響. 環(huán)境科學與技術, 2012, 35(10): 6- 11.

    [5] 姚建良. 三峽水庫蓄水后浮游甲殼動物生態(tài)學研究 [D]. 烏魯木齊: 新疆農(nóng)業(yè)大學, 2007.

    [6] 馮坤, 萬成炎, 彭建華, 張志永, 鄭志偉, 胡蓮, 鄒曦. 三峽庫區(qū)26條支流浮游甲殼動物的群落結構. 水生態(tài)學雜志, 2012, 33(4): 40- 48.

    [7] 朱愛民, 喬曄, 梁友光, 鄭志偉, 周蓮鳳, 鄒清, 汪紅軍. 三峽水庫支流童莊河擬多甲藻水華的監(jiān)測. 水生態(tài)學雜志, 2012, 33(4): 49- 53.

    [8] 李斌, 盧偉, 閆海. 長江三峽庫區(qū)湖北段水華發(fā)生情況及成因分析. 環(huán)境科學與技術, 2012, 35(6I): 214- 216.

    [9] 鄧春光, 龔玲. 三峽庫區(qū)富營養(yǎng)化發(fā)展趨勢研究. 農(nóng)業(yè)環(huán)境科學學報, 2007, 26(增刊): 279- 282.

    [10] 龔循矩, 肖化忠, 沈韞芬. 長江三峽地區(qū)的原生動物區(qū)系研究. 水生生物學報, 1990, 14(4): 289- 297.

    [11] 鄭金秀, 胡菊香, 周連鳳, 馬沛明, 王翔, 吳生桂. 長江上游原生動物的群落生態(tài)學研究. 水生態(tài)學雜志, 2009, 2(2): 88- 93.

    [12] 陳小娟, 潘曉潔, 馮坤, 鄭志偉, 彭建華, 萬成炎, 趙秀蘭. 小江回水區(qū)原生動物和輪蟲群落結構特征研究. 水生態(tài)學雜志, 2012, 33(5): 31- 35.

    [13] 沈韞芬,章宗涉,龔循矩,顧曼如,施之新,魏印心. 微型生物監(jiān)測新技術. 北京: 中國建筑工業(yè)出版社, 1990.

    [14] 蔣燮治,沈韞芬,龔循矩.西藏水生無脊椎動物. 北京: 科學出版社, 1983.

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

    [16] 林秋奇, 韓博平. 水庫生態(tài)系統(tǒng)特征研究及其在水庫水質(zhì)管理中的應用. 生態(tài)學報, 2001, 21(6): 1034- 1040.

    [17] Sommer U. Plankton Ecology: Succession in Plankton Communities. Berlin: Spring-Verlag, 1989: 253- 296.

    [18] 蔡慶華, 孫志禹. 三峽水庫水環(huán)境與水生態(tài)研究的進展與展望. 湖泊科學, 2012, 24(2): 169- 177.

    [19] 沈韞芬. 原生動物學. 北京: 科學出版社, 1999.

    [20] 徐潤林, 沈韞芬, 顧曼如. 武漢東湖浮游原生動物的生產(chǎn)量 // 第四次中國海洋湖沼科學會議論文集. 北京: 科學出版社, 1991: 164- 173.

    [21] 田佳玉. 原生動物生長和分布的影響因子. 哈爾濱師范大學自然科學學報, 2012, 28(4): 61- 63, 70- 70.

    [22] Patterson D J, Hedley S. Free-Living Freshwater Protozoa: A Colour Guide. London: Manson Publishing Ltd, 1996: 7- 221.

    [23] Pratt J R, Cairns J Jr. Functional groups in the protozoa: roles in differing ecosystems. The Journal of Protozoology, 1985, 32(3): 415- 423.

    [24] Kalff J. 湖沼學——內(nèi)陸水生態(tài)系統(tǒng). 古濱河, 劉正文, 李寬意, 等譯. 北京: 高等教育出版社, 2011.

    [25] Graham J M, Kent A D, Lauster G H, Yannarell A C, Graham L E, Triplett E W. Seasonal dynamics of phytoplankton and planktonic protozoan communities in a northern temperate humic lake: Diversity in a dinoflagellate dominated system. Microbial Ecology, 2004, 48(4): 528- 540.

    [26] Li F C, Kang X J, Yang W B, Guan Y Q, Zhang X H, Liu W W, Shen G M, Li J L, Wang H W. Protozoan community character in relation to trophic level in the Beijing section of the Juma River. Biodiversity Science, 2006, 14(4): 327- 332.

    [27] 胡菊香, 吳生桂, 唐會元, 鄒清, 胡望斌, 陳金生. 巢湖的原生動物及其對富營養(yǎng)化的響應. 水利漁業(yè), 2007, 27(1): 76- 79.

    [28] 張晟. 三峽庫區(qū)水體中營養(yǎng)鹽與浮游生物量分布特征 [D]. 重慶: 西南農(nóng)業(yè)大學, 2005.

    [29] Ioriya T, Inoue S, Haga M, Yogot N. Change of chemical and biological water environment at a newly constructed reservoir. Water Science and Technology, 1998, 37(2): 187- 194.

    [30] Nogueira M G. Zooplankton composition, dominance and abundance as indicators of environmental compartmentalization in Jurumirim Reservoir (Paranapanema River), S?o Paulo, Brazil. Hydrobiologia, 2001, 455(1/3): 1- 18.

    [31] Arruda J A, Marzolf G R, Faulk R T. The role of suspended sediments in the nutrition of zooplankton in turbid reservoirs. Ecology, 1983, 64(5): 1225- 1235.

    [32] McCabe G D, O′Brien W J. The effects of suspended silt on feeding and reproduction of Daphnia pulex. American Midland Naturalist, 1983, 110(2): 324- 337.

    [33] Jack J D, Gilbert J J. The effect of suspended clay on ciliate population growth rates. Freshwater Biology, 1993, 29(3): 385- 394.

    [34] 黃海強. 三峽庫區(qū)泥沙沉降對浮游生物的影響研究 [D]. 重慶: 西南大學, 2008.

    [35] Kirk K L, Gilbert J J. Suspended clay and the population dynamics of planktonic rotifers and cladocerans. Ecology, 1990, 71(5): 1741- 1755.

    [36] 韓德舉, 胡菊香, 高少波, 常秀林, 彭建華, 吳生桂, 胡傳林. 三峽水庫135 m蓄水過程壩前水域浮游生物變化的研究. 水利漁業(yè), 2005, 25(5): 55- 58.

    [37] 韓德舉, 彭建華, 簡東, 鄒清, 鄔紅娟. 丹江口水庫的餌料生物資源及水體營養(yǎng)狀態(tài)評價. 湖泊科學, 1997, 9(1): 57- 62.

    [38] 韓德舉, 吳生桂, 鄒清, 彭建華, 胡菊香. 陸水水庫的浮游生物及營養(yǎng)類型. 湖泊科學, 1996, 8(4): 351- 358.

    [39] 吳生桂, 簡東, 曾強, 張曉敏. 東江水庫對耒水中下游原生動物的影響. 長江流域資源與環(huán)境, 2000, 9(1): 124- 128.

    Spatio-temporal heterogeneity of protozoan communities during the period of Three Gorges Reservoir operation

    ZHENG Jinxiu, CHI Shiyun, LI Dan, WANG Hongjun, HU Juxiang*

    KeyLaboratoryofEcologicalImpactsofHydraulic-ProjectsandRestorationofAquaticEcosystemofMinistryofWaterResources,InstituteofHydroecology,MinistryofWaterResourcesandChineseAcademyofSciences.Wuhan430079,China

    Three Gorges Reservoir was impounded to 175 m for the first time in October 2010, of which the backwater area arrived in Chongqing, and the upstream water retention time extended further and the reaches changed to transitional waters. The heterogenrity of protozoan communities was studied in the main stem of Three Gorges Reservoir during the period of Three Gorges Reservoir normal operation (October 2010—June 2011). The spatial variations in community structure and species distribution were observed during reservoir drainage and storage. A total of 99 protozoan species were identified. Of these taxa, ciliates represented an increasing tendency after impoundment. The dominant species differed in operating stages, with the indicator species generally evolved from tintinnids (impounding period), ciliates (high water level period) to testate amoebae (low water level period). The dominant species varied in different reservoir areas, while the indicator species testate amoebae (transition zone) were changed into ciliates (lake zone) from upstream to downstream.Detrended correspondence analysis (DCA) and multi-response permutation procedure (MRPP) analysis showed that protozoan community structures and species distributions showed a clear spatio-temporal heterogeneity in Three Gorges Reservoir. The temporal heterogeneity peaked during impounding stage (P<0.05). The statistical analysis showed that the spatio-temporal variations in community pattern and species distribution were significantly correlated to environmental variables, especially transparency, temperature, conductivity and chlorophyll a. The average density and biomall of protozoa were 952.19 ind./L and 8.14 g/L, respectively. The standing crops of protozoa in upstream were higher than that in downstream during impounding stage, but lower than that during low water level period. The mean values of Marglef′s and Shannon-Weiner diversity indices were 3.78 and 2.18, respectively, with minimum value in January and maximum in June. The abundance of protozoa in transition zone was high after 175 m impoundment. It may be due to large amount of organic matter brought by water flow during the first impounding stage. Reservoir operation pattern of “drainage in winter and storage in summer” “blurs” seasonal variation of the river. Protozoan species composition and standing crops are mainly affected by reservoir hydrodynamic. Physic-chemical factors and hydrological characteristics showed different spatio-temporal variations, and resulted in the heterogeneity in community pattern and species distribution of protozoa on both temporal and spatial scales.

    Three Gorges Reservoir; protozoa; operation; spatio-temporal heterogeneity

    國家自然科學基金面上項目(51279112); 水利部公益性行業(yè)科研專項(201001001- 03)

    2013- 07- 24;

    2014- 06- 12

    10.5846/stxb201307241942

    *通訊作者Corresponding author.E-mail: hujuxiang2005@163.com

    鄭金秀, 池仕運, 李聃, 汪紅軍, 胡菊香.三峽水庫運行期間原生動物群落的時空異質(zhì)性.生態(tài)學報,2015,35(11):3569- 3579.

    Zheng J X, Chi S Y, Li D, Wang H J, Hu J X.Spatio-temporal heterogeneity of protozoan communities during the period of Three Gorges Reservoir operation.Acta Ecologica Sinica,2015,35(11):3569- 3579.

    猜你喜歡
    原生動物纖毛蟲水區(qū)
    拉薩河中下游原生動物群落物種分布格局及其驅(qū)動因素
    渤海灣浮游纖毛蟲豐度和生物量的周年變化
    瓶瓶罐罐大世界之沙水區(qū)玩具
    DNA條形碼在原生動物系統(tǒng)分類學中的應用
    動物學實驗教學中“原生動物的結構觀察”伸縮泡觀察方法的優(yōu)化
    英國城市更新研究——以克萊德濱水區(qū)為例
    渤海灣近岸海域浮游纖毛蟲豐度和生物量的季節(jié)變化*
    后申遺時代杭州段運河濱水區(qū)可持續(xù)保護利用思路探討
    海洋浮游纖毛蟲攝食研究綜述*
    海洋與湖沼(2016年1期)2016-01-15 03:52:09
    拉魯濕地夏、秋季纖毛蟲物種多樣性研究
    西藏科技(2015年6期)2015-09-26 12:12:08
    人成视频在线观看免费观看| av女优亚洲男人天堂| 晚上一个人看的免费电影| 秋霞伦理黄片| 美女大奶头黄色视频| 亚洲欧美成人综合另类久久久| 黄片播放在线免费| 在线免费观看不下载黄p国产| 亚洲av成人精品一区久久| 一区在线观看完整版| 黑人高潮一二区| 啦啦啦视频在线资源免费观看| 在线观看www视频免费| 两个人免费观看高清视频| 成年美女黄网站色视频大全免费 | 51国产日韩欧美| 日日啪夜夜爽| 99国产综合亚洲精品| 久久久久精品久久久久真实原创| av在线播放精品| 欧美丝袜亚洲另类| 看免费成人av毛片| 免费播放大片免费观看视频在线观看| 日本爱情动作片www.在线观看| 秋霞在线观看毛片| 蜜桃在线观看..| 天天躁夜夜躁狠狠久久av| xxx大片免费视频| 最近手机中文字幕大全| 国产亚洲av片在线观看秒播厂| videosex国产| 只有这里有精品99| 人人妻人人添人人爽欧美一区卜| 特大巨黑吊av在线直播| 黑丝袜美女国产一区| 精品国产乱码久久久久久小说| 欧美日韩视频高清一区二区三区二| 一级黄片播放器| 男人添女人高潮全过程视频| 国产永久视频网站| 久久午夜综合久久蜜桃| 婷婷成人精品国产| 成人免费观看视频高清| 国产伦精品一区二区三区视频9| 热re99久久国产66热| 成年av动漫网址| 欧美激情极品国产一区二区三区 | 色视频在线一区二区三区| 日韩一区二区三区影片| 欧美日本中文国产一区发布| 人人妻人人澡人人看| 一区二区三区四区激情视频| 欧美成人午夜免费资源| 精品熟女少妇av免费看| 九色成人免费人妻av| 免费日韩欧美在线观看| 日韩一区二区视频免费看| 人妻夜夜爽99麻豆av| 国产男人的电影天堂91| 亚洲欧美一区二区三区国产| 亚洲性久久影院| 国产精品一国产av| 国产av国产精品国产| 美女视频免费永久观看网站| 看免费成人av毛片| 亚洲国产精品专区欧美| 亚洲国产色片| 国产精品不卡视频一区二区| 精品久久久久久久久亚洲| 亚洲精品久久午夜乱码| 伊人久久国产一区二区| 我的老师免费观看完整版| 在线观看三级黄色| 韩国av在线不卡| 18+在线观看网站| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 免费黄色在线免费观看| 国产探花极品一区二区| 久久精品国产自在天天线| 亚洲欧美日韩卡通动漫| 在线观看人妻少妇| 久久热精品热| 国产高清有码在线观看视频| 精品亚洲成a人片在线观看| 99精国产麻豆久久婷婷| 精品国产一区二区久久| 丰满少妇做爰视频| 一个人免费看片子| 午夜福利在线观看免费完整高清在| 777米奇影视久久| 国产老妇伦熟女老妇高清| 国产淫语在线视频| av播播在线观看一区| 欧美日韩视频高清一区二区三区二| 亚洲精品456在线播放app| 久久鲁丝午夜福利片| 亚洲av成人精品一二三区| 热re99久久国产66热| 国产亚洲av片在线观看秒播厂| 亚洲欧美一区二区三区黑人 | 美女国产视频在线观看| 最近手机中文字幕大全| 国产不卡av网站在线观看| av在线观看视频网站免费| 搡老乐熟女国产| 国产69精品久久久久777片| 十八禁网站网址无遮挡| 26uuu在线亚洲综合色| av卡一久久| 精品亚洲成国产av| 国产精品 国内视频| 久久午夜综合久久蜜桃| 97在线人人人人妻| 日韩一本色道免费dvd| 久久久久人妻精品一区果冻| 亚洲美女搞黄在线观看| 边亲边吃奶的免费视频| 久久97久久精品| 黑人巨大精品欧美一区二区蜜桃 | 亚州av有码| 男人爽女人下面视频在线观看| 日韩欧美一区视频在线观看| 3wmmmm亚洲av在线观看| 中文乱码字字幕精品一区二区三区| 午夜免费男女啪啪视频观看| 视频区图区小说| videossex国产| 亚洲欧美色中文字幕在线| 91午夜精品亚洲一区二区三区| 国产精品女同一区二区软件| 免费黄频网站在线观看国产| 中文天堂在线官网| 国精品久久久久久国模美| 乱码一卡2卡4卡精品| 纵有疾风起免费观看全集完整版| 伦理电影免费视频| 亚洲精品中文字幕在线视频| 中文字幕人妻丝袜制服| 美女xxoo啪啪120秒动态图| 亚洲国产成人一精品久久久| 我的老师免费观看完整版| 久久影院123| 18禁动态无遮挡网站| 97超视频在线观看视频| 久久精品国产亚洲网站| 久久久国产精品麻豆| 一本—道久久a久久精品蜜桃钙片| 一本一本综合久久| 国产欧美亚洲国产| 成人国产av品久久久| 永久免费av网站大全| 狂野欧美激情性xxxx在线观看| 亚洲av综合色区一区| videosex国产| 久久久久人妻精品一区果冻| 精品99又大又爽又粗少妇毛片| 夫妻性生交免费视频一级片| 成人影院久久| 中文字幕亚洲精品专区| av网站免费在线观看视频| 9色porny在线观看| 日本91视频免费播放| 免费av不卡在线播放| 欧美精品国产亚洲| 欧美3d第一页| 人成视频在线观看免费观看| 精品少妇内射三级| kizo精华| 18禁裸乳无遮挡动漫免费视频| 亚洲精品第二区| 久久鲁丝午夜福利片| 亚洲精品乱久久久久久| 久久精品国产亚洲网站| 国产成人freesex在线| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说| 精品少妇黑人巨大在线播放| 各种免费的搞黄视频| 久久人人爽人人片av| 国产成人精品福利久久| 一级片'在线观看视频| 国产乱来视频区| 国产成人精品在线电影| 搡老乐熟女国产| 亚洲精品一区蜜桃| 只有这里有精品99| 美女国产高潮福利片在线看| videosex国产| 久久人人爽av亚洲精品天堂| 黄色欧美视频在线观看| 又大又黄又爽视频免费| 久久人人爽人人片av| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图 | 免费观看a级毛片全部| 国语对白做爰xxxⅹ性视频网站| 成人午夜精彩视频在线观看| 亚洲欧洲国产日韩| 妹子高潮喷水视频| 伊人亚洲综合成人网| 91在线精品国自产拍蜜月| 丝袜喷水一区| 国产一区二区在线观看日韩| 亚洲欧美一区二区三区国产| 亚洲经典国产精华液单| 欧美性感艳星| 精品视频人人做人人爽| 亚洲怡红院男人天堂| 国产 精品1| 成年人免费黄色播放视频| 自拍欧美九色日韩亚洲蝌蚪91| 91精品国产九色| 成年美女黄网站色视频大全免费 | 熟女电影av网| 黄色毛片三级朝国网站| 赤兔流量卡办理| www.色视频.com| 久久人人爽人人片av| 美女脱内裤让男人舔精品视频| 99九九线精品视频在线观看视频| h视频一区二区三区| 日本午夜av视频| 曰老女人黄片| 人体艺术视频欧美日本| 美女xxoo啪啪120秒动态图| 国产高清三级在线| 精品久久久久久电影网| 国产精品久久久久久精品古装| 免费观看av网站的网址| 国产免费视频播放在线视频| 日韩成人伦理影院| 日本猛色少妇xxxxx猛交久久| 特大巨黑吊av在线直播| 国产黄频视频在线观看| 国产 精品1| 老司机影院毛片| 搡老乐熟女国产| 又黄又爽又刺激的免费视频.| 日本色播在线视频| 99热6这里只有精品| 全区人妻精品视频| 一区二区三区四区激情视频| 欧美激情 高清一区二区三区| 一本色道久久久久久精品综合| 丰满饥渴人妻一区二区三| 免费观看在线日韩| 热re99久久国产66热| 欧美老熟妇乱子伦牲交| 国产男女内射视频| 丝袜喷水一区| 美女国产视频在线观看| 九九爱精品视频在线观看| 国产女主播在线喷水免费视频网站| 少妇猛男粗大的猛烈进出视频| 午夜免费观看性视频| 亚洲av男天堂| 熟妇人妻不卡中文字幕| 日本av免费视频播放| 老司机影院成人| 亚洲成色77777| 亚洲情色 制服丝袜| 国产免费又黄又爽又色| 99热网站在线观看| 久久久久久久国产电影| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 天堂中文最新版在线下载| 久久精品久久精品一区二区三区| 国产精品人妻久久久久久| 黑人高潮一二区| 内地一区二区视频在线| 一区在线观看完整版| 狂野欧美白嫩少妇大欣赏| 你懂的网址亚洲精品在线观看| 人人妻人人添人人爽欧美一区卜| 成人午夜精彩视频在线观看| 日本猛色少妇xxxxx猛交久久| 国产日韩欧美亚洲二区| 视频区图区小说| 国产69精品久久久久777片| 国产女主播在线喷水免费视频网站| 好男人视频免费观看在线| 91精品三级在线观看| 免费av中文字幕在线| 乱码一卡2卡4卡精品| 日韩伦理黄色片| 一边摸一边做爽爽视频免费| 男女啪啪激烈高潮av片| 国产精品久久久久久精品电影小说| 人人妻人人澡人人爽人人夜夜| 成人漫画全彩无遮挡| 美女中出高潮动态图| 99久久中文字幕三级久久日本| 欧美亚洲日本最大视频资源| 亚洲av中文av极速乱| 精品一区二区三区视频在线| 午夜福利在线观看免费完整高清在| 嫩草影院入口| 久久久久国产网址| 9色porny在线观看| 亚洲性久久影院| 精品一区在线观看国产| 99久久精品一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 国产色爽女视频免费观看| 少妇人妻 视频| 热99久久久久精品小说推荐| 九色成人免费人妻av| 久久久a久久爽久久v久久| 草草在线视频免费看| 国产精品久久久久久精品古装| 亚洲精品日韩在线中文字幕| 国产成人aa在线观看| 亚洲第一区二区三区不卡| 毛片一级片免费看久久久久| 五月开心婷婷网| 9色porny在线观看| 国产精品99久久久久久久久| 久久综合国产亚洲精品| 色5月婷婷丁香| 多毛熟女@视频| 中文字幕久久专区| xxx大片免费视频| 丝袜喷水一区| 日韩精品有码人妻一区| 高清毛片免费看| 亚洲天堂av无毛| 秋霞在线观看毛片| 免费人妻精品一区二区三区视频| 韩国av在线不卡| 日韩av不卡免费在线播放| 久久久久精品性色| 91久久精品电影网| 简卡轻食公司| 男女啪啪激烈高潮av片| 久久精品国产亚洲网站| 国产精品成人在线| 青春草国产在线视频| 伊人亚洲综合成人网| 中文乱码字字幕精品一区二区三区| 免费观看在线日韩| 久久久久网色| 亚洲av欧美aⅴ国产| av不卡在线播放| 精品国产一区二区三区久久久樱花| 99久久综合免费| 欧美激情极品国产一区二区三区 | 亚洲人成网站在线播| 国产 一区精品| 在线观看www视频免费| 免费人妻精品一区二区三区视频| 免费看av在线观看网站| 欧美日韩视频精品一区| 欧美日韩亚洲高清精品| av福利片在线| 熟妇人妻不卡中文字幕| 亚洲av男天堂| 国产日韩欧美在线精品| 色哟哟·www| 久久狼人影院| 自拍欧美九色日韩亚洲蝌蚪91| 又黄又爽又刺激的免费视频.| 成人国语在线视频| 青春草视频在线免费观看| 国产精品人妻久久久影院| 国产成人精品一,二区| 日韩大片免费观看网站| 欧美精品高潮呻吟av久久| 精品久久蜜臀av无| 国产高清三级在线| 一区在线观看完整版| 最近2019中文字幕mv第一页| 欧美激情极品国产一区二区三区 | 丰满乱子伦码专区| 久久久久久久国产电影| 午夜日本视频在线| 最近手机中文字幕大全| av有码第一页| 黄片播放在线免费| 久久久精品免费免费高清| 我的女老师完整版在线观看| 如何舔出高潮| 男女边摸边吃奶| 99九九在线精品视频| 精品少妇黑人巨大在线播放| 久久影院123| 免费人成在线观看视频色| 亚洲综合色惰| 晚上一个人看的免费电影| 久久热精品热| 99久久精品国产国产毛片| 亚洲精品美女久久av网站| 一级毛片电影观看| 伊人久久国产一区二区| 久久精品国产自在天天线| 日本黄色片子视频| 久久久亚洲精品成人影院| 久久午夜福利片| 国产成人freesex在线| 成人18禁高潮啪啪吃奶动态图 | 久久久午夜欧美精品| 亚洲欧美中文字幕日韩二区| 中文字幕人妻丝袜制服| 美女内射精品一级片tv| 久久精品国产亚洲av涩爱| 国国产精品蜜臀av免费| 久久久国产欧美日韩av| 又大又黄又爽视频免费| 18禁动态无遮挡网站| 在线观看免费日韩欧美大片 | 久久99蜜桃精品久久| 久久国产精品男人的天堂亚洲 | 丰满饥渴人妻一区二区三| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图 | 边亲边吃奶的免费视频| 色视频在线一区二区三区| 精品人妻在线不人妻| 亚州av有码| av线在线观看网站| 日本午夜av视频| 韩国av在线不卡| 亚洲精品aⅴ在线观看| 99久久综合免费| 国产成人freesex在线| 91久久精品国产一区二区三区| 国产亚洲精品第一综合不卡 | 五月玫瑰六月丁香| 亚洲人成网站在线播| 日本av免费视频播放| 亚洲欧美日韩另类电影网站| 欧美xxⅹ黑人| 性色avwww在线观看| 在线天堂最新版资源| 国内精品宾馆在线| 一级,二级,三级黄色视频| 国产视频内射| 久久精品国产自在天天线| 国产日韩欧美亚洲二区| 国产精品三级大全| 伊人久久精品亚洲午夜| 一区二区三区免费毛片| 亚洲av欧美aⅴ国产| 青春草国产在线视频| 国产黄频视频在线观看| av线在线观看网站| 九色成人免费人妻av| 日韩成人伦理影院| 国产成人av激情在线播放 | 一级毛片电影观看| 国产探花极品一区二区| 国语对白做爰xxxⅹ性视频网站| 午夜福利网站1000一区二区三区| 久久精品国产亚洲网站| 日韩不卡一区二区三区视频在线| 久久国产精品男人的天堂亚洲 | 丰满饥渴人妻一区二区三| 久久久久精品久久久久真实原创| 精品久久久噜噜| 欧美另类一区| 亚洲在久久综合| 毛片一级片免费看久久久久| 两个人的视频大全免费| 蜜桃国产av成人99| 十分钟在线观看高清视频www| 国产亚洲精品第一综合不卡 | 精品国产国语对白av| 建设人人有责人人尽责人人享有的| 国产精品无大码| 亚洲精品一区蜜桃| 久久人人爽人人爽人人片va| 亚洲,欧美,日韩| 99九九在线精品视频| 中文天堂在线官网| 国产精品久久久久久久久免| 爱豆传媒免费全集在线观看| 中国美白少妇内射xxxbb| 91久久精品国产一区二区成人| 免费观看av网站的网址| 久久99热6这里只有精品| 99热国产这里只有精品6| 乱码一卡2卡4卡精品| 国产精品欧美亚洲77777| 国产在线视频一区二区| av电影中文网址| 日韩av在线免费看完整版不卡| 国产精品久久久久久精品古装| 久久久久久久久大av| 两个人免费观看高清视频| 国产高清有码在线观看视频| 亚洲欧美一区二区三区黑人 | 狂野欧美激情性bbbbbb| 99热国产这里只有精品6| 亚洲欧美色中文字幕在线| 亚洲av电影在线观看一区二区三区| 欧美日本中文国产一区发布| 亚洲av免费高清在线观看| 婷婷色麻豆天堂久久| 特大巨黑吊av在线直播| 国产成人精品婷婷| 伦精品一区二区三区| 一区二区三区乱码不卡18| 伊人亚洲综合成人网| 亚洲国产成人一精品久久久| 国产成人精品婷婷| 亚洲国产日韩一区二区| 五月天丁香电影| 午夜av观看不卡| av在线老鸭窝| 国产日韩一区二区三区精品不卡 | 色视频在线一区二区三区| 国产熟女欧美一区二区| 建设人人有责人人尽责人人享有的| 日本wwww免费看| 80岁老熟妇乱子伦牲交| 久久 成人 亚洲| 亚洲欧美一区二区三区黑人 | 亚洲精品成人av观看孕妇| 日韩伦理黄色片| 久久精品熟女亚洲av麻豆精品| 成人18禁高潮啪啪吃奶动态图 | 亚洲,一卡二卡三卡| 满18在线观看网站| 久久久久久久大尺度免费视频| 免费不卡的大黄色大毛片视频在线观看| 少妇被粗大的猛进出69影院 | 女性被躁到高潮视频| 国产黄频视频在线观看| tube8黄色片| 黑人欧美特级aaaaaa片| 欧美三级亚洲精品| 999精品在线视频| 99久久中文字幕三级久久日本| 亚洲欧美成人精品一区二区| 插阴视频在线观看视频| 超色免费av| 高清在线视频一区二区三区| 国产精品欧美亚洲77777| 看十八女毛片水多多多| 高清欧美精品videossex| 欧美精品一区二区免费开放| 日本91视频免费播放| 久久人人爽av亚洲精品天堂| 成年美女黄网站色视频大全免费 | 各种免费的搞黄视频| 国产精品99久久久久久久久| 色婷婷久久久亚洲欧美| 久久精品国产自在天天线| h视频一区二区三区| av卡一久久| 婷婷色麻豆天堂久久| 精品亚洲成a人片在线观看| 欧美性感艳星| 寂寞人妻少妇视频99o| 亚洲精品中文字幕在线视频| 久久人人爽人人爽人人片va| 少妇高潮的动态图| 欧美日韩av久久| 亚洲怡红院男人天堂| 在线播放无遮挡| 高清欧美精品videossex| 国产 一区精品| 国产淫语在线视频| 丝袜美足系列| 久久久午夜欧美精品| 久久久a久久爽久久v久久| 黄色怎么调成土黄色| 国产女主播在线喷水免费视频网站| 亚洲人成网站在线播| 91久久精品国产一区二区三区| 精品国产一区二区三区久久久樱花| 欧美丝袜亚洲另类| 国产国拍精品亚洲av在线观看| 日韩熟女老妇一区二区性免费视频| 精品一区二区三区视频在线| 搡女人真爽免费视频火全软件| 97超碰精品成人国产| 一区二区三区乱码不卡18| 制服诱惑二区| 大又大粗又爽又黄少妇毛片口| 男男h啪啪无遮挡| 这个男人来自地球电影免费观看 | 9色porny在线观看| 满18在线观看网站| 亚洲精品乱码久久久v下载方式| 丰满乱子伦码专区| 一级毛片黄色毛片免费观看视频| 欧美3d第一页| 精品久久久久久久久av| 国产白丝娇喘喷水9色精品| 国产日韩欧美视频二区| 日韩三级伦理在线观看| 精品酒店卫生间| 午夜激情福利司机影院| av免费在线看不卡| 亚洲在久久综合| 国产成人aa在线观看| 日韩大片免费观看网站| 看免费成人av毛片| 秋霞在线观看毛片| xxx大片免费视频| videosex国产| 性色av一级| 国产成人精品婷婷| 久久精品人人爽人人爽视色| 久久久久国产精品人妻一区二区| 国产精品国产三级国产av玫瑰| 精品人妻熟女毛片av久久网站| 最近2019中文字幕mv第一页| 日韩大片免费观看网站| 国产精品一区二区三区四区免费观看| 久久精品国产亚洲网站| 大片电影免费在线观看免费| 久久精品人人爽人人爽视色|