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

    西藏尼洋河水生生物群落時(shí)空動(dòng)態(tài)及與環(huán)境因子的關(guān)系:1.浮游植物*

    2013-05-28 07:36:16劉海平葉少文楊雪峰張良松鐘國(guó)輝李鐘杰
    湖泊科學(xué) 2013年5期
    關(guān)鍵詞:尼洋河底質(zhì)綠藻

    劉海平,葉少文,楊雪峰,張良松,鐘國(guó)輝,李鐘杰**

    (1:西藏大學(xué)農(nóng)牧學(xué)院,林芝 860000)

    (2:中國(guó)科學(xué)院水生生物研究所,淡水生態(tài)與生物技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,武漢 430072)

    (3:中國(guó)科學(xué)院大學(xué),北京 100049)

    (4:錫林浩特市第六中學(xué),錫林浩特 026000)

    (5:福建省水產(chǎn)技術(shù)推廣總站,福州 350003)

    由于青藏高原的隆升,對(duì)高原河流湖泊水系的演化起到了積極的推動(dòng)作用[1],因而其獨(dú)特的生態(tài)位以及科研價(jià)值愈發(fā)凸顯.一方面高原河流湖泊對(duì)全球變暖有著積極的響應(yīng)[2],擁有特殊的生態(tài)系統(tǒng)服務(wù)功能[3],另外,由于高原湖泊沉積物中儲(chǔ)存的信息折射了礦物學(xué)[4]、同位素地球化學(xué)[5]、生物學(xué)[6]等方面內(nèi)容,得到了科研工作者的推崇.不容忽視的是,青藏高原氣溫在逐年上升[7-8],預(yù)計(jì)到2100年青藏高原氣溫將上升2.0 ~3.6℃[9],加之人類干擾程度的不斷加劇[10],對(duì)高原河流生態(tài)系統(tǒng)產(chǎn)生了巨大的沖擊,水資源脆弱性發(fā)生了改變[11],地表徑流減少,湖泊萎縮,濕地退化.同時(shí),由于高原河流蘊(yùn)藏著巨大的水能資源,水電站的開發(fā)和建設(shè)不可避免地踴躍出現(xiàn),河流的生態(tài)系統(tǒng)完整性迎來(lái)了挑戰(zhàn)[12].生態(tài)系統(tǒng)完整性不僅要求生態(tài)系統(tǒng)結(jié)構(gòu)的完整以及結(jié)構(gòu)合理,也要求生態(tài)系統(tǒng)功能健全以及功能正常發(fā)揮,比較而言,高原河流水資源較為脆弱,水資源承載能力與區(qū)域經(jīng)濟(jì)、社會(huì)發(fā)展和生態(tài)環(huán)境保護(hù)方面容易發(fā)生錯(cuò)位[13],促進(jìn)高原河流流域的可持續(xù)發(fā)展則顯得尤為重要[14-15].

    西藏是我國(guó)河流較多的省區(qū)之一,流域面積大于10000 km2的河流有20 余條,該區(qū)水量充沛,僅次于四川,居全國(guó)各省區(qū)的第二位[16].隨著社會(huì)和經(jīng)濟(jì)的發(fā)展、旅游業(yè)的興起、水電站開發(fā)、現(xiàn)代農(nóng)牧業(yè)的推廣等,逐漸演變?yōu)槟嵫蠛又苓叺貐^(qū)發(fā)展趨勢(shì),加之尼洋河魚類區(qū)系結(jié)構(gòu)較為簡(jiǎn)單[17],易受到環(huán)境條件的波動(dòng).本系列論文選擇尼洋河水生生物作為研究對(duì)象,分析各類水生生物群落結(jié)構(gòu)的時(shí)空動(dòng)態(tài),判斷尼洋河水質(zhì)理化因子與各類水生生物群落之間的相互關(guān)系,通過(guò)篩選關(guān)鍵水質(zhì)理化因子,推斷和演繹它們與水生生物群落結(jié)構(gòu)的作用機(jī)制,從而為尼洋河水域生態(tài)的可持續(xù)發(fā)展提供理論基礎(chǔ).

    浮游植物是一個(gè)生態(tài)學(xué)單位,它包括所有生活在水中營(yíng)浮游方式生活的微小植物,通常就指浮游藻類,而不包括細(xì)菌和其他植物[18].浮游植物可作為漁產(chǎn)潛力估算重要的依據(jù)[19],在水域生態(tài)系統(tǒng)的上行效應(yīng)和下行效應(yīng)中發(fā)揮著重要的作用,很有必要對(duì)浮游植物進(jìn)行功能類群劃分[20],以便更好地理解浮游植物群落,控制水域生態(tài)系統(tǒng)的能量流和物質(zhì)流,從而保障水域生態(tài)系統(tǒng)的可持續(xù)發(fā)展.研究表明,氣候變暖促進(jìn)了浮游植物春季的繁殖,浮游植物最大生物量以及個(gè)體大小呈現(xiàn)減小的趨勢(shì)[21],浮游植物群落結(jié)構(gòu)的演替從耐溫性上表現(xiàn)出冷水性物種減少而暖水性物種增多的趨勢(shì),而在營(yíng)養(yǎng)結(jié)構(gòu)上表現(xiàn)出嗜氮性物種增多的演替趨勢(shì)[22].由于海拔和環(huán)境的影響,高原河流浮游植物群落比較穩(wěn)定,種類分布較為均勻,表現(xiàn)為貧營(yíng)養(yǎng)型河流的特征,水域環(huán)境表現(xiàn)良好[23].在全球氣候變暖和人類活動(dòng)對(duì)自然影響的框架之下,尼洋河水域生態(tài)安全屏障愈發(fā)突出,通過(guò)有效檢測(cè)浮游植物群落,積極推動(dòng)尼洋河水域生態(tài)系統(tǒng)健康發(fā)展則顯得尤為重要.

    1 實(shí)驗(yàn)設(shè)計(jì)和數(shù)據(jù)處理

    1.1 研究區(qū)域概況

    尼洋河流域位于西藏林芝地區(qū)西部(29°26' ~ 29°55'N,92°22' ~94°27'E).全長(zhǎng) 286 km,流域面積為17535 km2,河流大體由西偏北往東偏南流,其中,上游大致為東西流向,到中下游尼西附近轉(zhuǎn)為南東流向.尼洋河全干流的落差為2080 m,平均坡降為7.27‰[16].

    1.2 采樣點(diǎn)設(shè)置

    根據(jù)河流形態(tài)特征設(shè)置4 個(gè)采樣點(diǎn),每個(gè)采樣點(diǎn)采集3 次平行樣進(jìn)行分析.4 個(gè)采樣點(diǎn)分別位于尼洋河中上游(樣點(diǎn)Ⅰ;29°53'N,93°16'E;海拔 3391 m)、尼洋河中游(樣點(diǎn)Ⅱ;30°00'N,93°54'E;海拔 3225 m)、尼洋河中下游(樣點(diǎn)Ⅲ;29°31'N,94°27'E;海拔2946 m)以及尼洋河與雅魯藏布江交匯處(樣點(diǎn)Ⅳ;29°26'N,94°28'E;海拔 2919 m),詳見圖 1,采樣范圍為 29°26' ~30°00'N,93°16' ~94°28'E,海拔高度從 3391 m 降至2919 m,落差接近500 m.

    圖1 尼洋河采樣點(diǎn)分布Fig.1 Distribution of sampling sites in Niyang River

    1.3 采樣時(shí)間

    從2008年11月開始采集水樣,用于理化分析,除2月份水樣因部分路面結(jié)冰無(wú)法采樣外,其他月份水樣均已按照月份進(jìn)行采集和分析;各類水生生物樣本按照每個(gè)樣點(diǎn)每個(gè)季節(jié)采集1 次,即:春季(3月份)、夏季(6月份)、秋季(9月份)、冬季(12月份),共 16 個(gè)樣本.

    1.4 環(huán)境因子的采集和分析

    為了探析各類水生生物與環(huán)境因子之間的關(guān)系,初步選擇11 項(xiàng)環(huán)境因子作為研究對(duì)象,包括:1 個(gè)分組變量,即:河流底質(zhì)類型(砂石、泥土);10 個(gè)連續(xù)變量:采樣點(diǎn)海拔、表層水溫(WT)、表層pH、硬度(Hardness)、礦化度(TDS)、表層溶解氧(DO)、總氮(TN)、銨態(tài)氮(NH4+-N)、總磷(TP)和總堿度(Alkalinity).

    表層水溫和pH 用德國(guó)CX-401 便攜式多參數(shù)水質(zhì)測(cè)定儀現(xiàn)場(chǎng)測(cè)定.硬度、礦化度、溶解氧、總氮、銨態(tài)氮、總磷和總堿度參照文獻(xiàn)[24]開展野外采集和保存工作,帶回實(shí)驗(yàn)室分析測(cè)定.以當(dāng)量濃度(meq/L)計(jì)Ca2+和Mg2+濃度,二者之和乘以系數(shù)2.804,即為水體的硬度,單位為°DH(德國(guó)度).Mg2+和Ca2+均參照GB/T 11905-1989(火焰原子吸收法)測(cè)定;水樣在103 ~105℃下烘干,所得為殘?jiān)?,在沒有污染的水體中,殘?jiān)康韧瑸榈V化度[24];表層溶解氧參照GB/T 7489-1987(碘量法)測(cè)定;總氮參照GB/T 11894-1989(堿性過(guò)硫酸鉀消解紫外分光光度法)測(cè)定;銨態(tài)氮參照GB/T 7479-1987(鈉氏試劑光度法)測(cè)定;總磷參照GB/T 11893-1989(鉬銻抗分光光度法)測(cè)定;總堿度參照酸堿指示劑滴定法測(cè)定[24].

    1.5 數(shù)據(jù)統(tǒng)計(jì)和分析

    用主成分分析方法(PCA)分析尼洋河浮游植物時(shí)空特征,用Duncan 法檢驗(yàn)各采樣點(diǎn)以及季節(jié)之間各類水生生物相關(guān)參數(shù)的差異性,用分類回歸樹模型(Classification and Regression Tree Model,CART)和典范對(duì)應(yīng)分析方法(Canonical Correlation Analysis,CCA)分析尼洋河各類水生生物與環(huán)境因子之間的關(guān)系.分析前,須對(duì)按照月份采集的水體理化指標(biāo)根據(jù)季節(jié)歸類整理,即:春季(1、3月份)、夏季(4-6月份)、秋季(7-9月份)和冬季(10-12月份).在CCA 分析前對(duì)各類水生生物數(shù)據(jù)和環(huán)境因子數(shù)據(jù)進(jìn)行l(wèi)g(x+1)對(duì)數(shù)轉(zhuǎn)化的標(biāo)準(zhǔn)化處理,用R 語(yǔ)言軟件(2.14.1)對(duì)數(shù)據(jù)進(jìn)行分析.

    1.6 浮游植物的采集和分析

    用300 目的浮游生物網(wǎng)過(guò)濾20 L 水,得到的浮游植物定量樣品加入魯哥試劑保存,帶回實(shí)驗(yàn)室,靜置48 h后濃縮.計(jì)數(shù)時(shí),將濃縮樣充分搖勻后吸取0.1 ml 置于計(jì)數(shù)框內(nèi),在10×40 倍顯微鏡下觀察,每片有46行,每片計(jì)數(shù)框計(jì)數(shù)N 行,每個(gè)樣品計(jì)數(shù)2 次,取其平均值.2 次計(jì)數(shù)結(jié)果與其平均值之差小于10%說(shuō)明結(jié)果有效,否則須計(jì)數(shù)第3 片,鑒定到種或?qū)賉25-26].

    按照如下公式對(duì)浮游植物進(jìn)行定量分析:

    為了進(jìn)一步反映尼洋河浮游植物群落的時(shí)空特征,對(duì)Shannon-Wiener 多樣性指數(shù)(H')、Pielou 均勻度指數(shù)(J')、物種豐富度和總豐度等進(jìn)行分析.

    Shannon-Wiener 多樣性指數(shù)[27]為:

    Pielou 均勻度指數(shù)[28]為:

    式中,S 為浮游植物種類數(shù),pi為浮游植物i 種占所有浮游植物的比例(%),H'max為理論最大多樣性指數(shù).

    物種豐富度指某一采樣點(diǎn)或者季節(jié)出現(xiàn)的物種個(gè)數(shù).

    總豐度指某一采樣點(diǎn)或者季節(jié)單位體積出現(xiàn)的浮游植物個(gè)體數(shù)量(cells/L).

    出現(xiàn)頻率指某一浮游植物出現(xiàn)的次數(shù)占所有調(diào)查樣點(diǎn)數(shù)的百分比.

    相對(duì)豐度指某一浮游植物個(gè)體數(shù)占所有物種個(gè)體數(shù)的百分比.

    2 結(jié)果與分析

    2.1 尼洋河浮游植物種類以及優(yōu)勢(shì)種

    尼洋河共有浮游植物7 門29 科48 屬,包括硅藻門、藍(lán)藻門、綠藻門、隱藻門、黃藻門、甲藻門、裸藻門,其中硅藻門8 科20 屬,藍(lán)藻門6 科8 屬,綠藻門10 科13 屬,甲藻門2 科2 屬,裸藻門1 科3 屬,隱藻門和黃藻門各1 科1 屬(表1).

    將浮游植物出現(xiàn)頻率在90%以上的定義為優(yōu)勢(shì)種,結(jié)果表明,尼洋河浮游植物優(yōu)勢(shì)類群以硅藻為主,包括脆桿藻科的脆桿藻和等片藻,舟形藻科的舟形藻,橋彎藻科的橋彎藻,異極藻科的異極藻,曲殼藻科的曲殼藻,菱板藻科的菱形藻,以上7 屬硅藻的相對(duì)豐度之和達(dá)79.25%,其中以橋彎藻的相對(duì)豐度最大,達(dá)19.70%,其次是曲殼藻和舟形藻,分別達(dá)18.54%和10.62%(表1).

    2.2 尼洋河浮游植物時(shí)空變化特征

    尼洋河浮游植物主要包括硅藻、綠藻和藍(lán)藻,對(duì)其總豐度自然對(duì)數(shù)化,縮小其數(shù)量級(jí)差距,分析尼洋河3大藻類的時(shí)空演替特征(圖2).在尼洋河沿程以及季節(jié)方面,硅藻的豐度占有絕對(duì)優(yōu)勢(shì),為尼洋河的優(yōu)勢(shì)藻類.就尼洋河沿程方面而言,綠藻總豐度在尼洋河下游出現(xiàn)一個(gè)峰值,以此為界,尼洋河其他河段綠藻總豐度處于減少趨勢(shì),硅藻總豐度與綠藻的變化趨勢(shì)一致,藍(lán)藻總豐度變化趨勢(shì)與綠藻有所不同,尼洋河中游藍(lán)藻總豐度最低,低于綠藻的總豐度,在下游略有上升,但也未及綠藻的總豐度.

    在季節(jié)方面,硅藻和綠藻總豐度在夏季均出現(xiàn)最大值,冬季則不同,硅藻總豐度達(dá)到了最低值,而綠藻則有所回升,但仍比夏季低,藍(lán)藻總豐度在秋季出現(xiàn)最低值,以此為界,尼洋河其他河段的藍(lán)藻總豐度呈增加趨勢(shì),最大值出現(xiàn)在冬季.

    表1 尼洋河浮游植物種類出現(xiàn)頻率和相對(duì)豐度*Tab.1 Occurrence frequency and relative abundance of phytoplankton species in Niyang River

    圖2 尼洋河浮游植物硅藻、綠藻、藍(lán)藻時(shí)空演替特征Fig.2 Spatio-temporal succession characteristics for phytoplankton about Bacillariophyta,Cyanophyta and Chlorophyta in Niyang River

    圖3 基于PCA 分析尼洋河浮游植物時(shí)空特征(圖中數(shù)字為樣點(diǎn)標(biāo)記,其中1 ~4分別代表采樣點(diǎn)Ⅰ的春、夏、秋、冬季,5 ~8 代表采樣點(diǎn)Ⅱ的四季,9 ~12 代表采樣點(diǎn)Ⅲ的四季,13 ~16 代表采樣點(diǎn)Ⅳ的四季;S1 代表采樣點(diǎn)Ⅰ,S2 代表采樣點(diǎn)Ⅱ,S3 代表采樣點(diǎn)Ⅲ,S4 代表采樣點(diǎn)Ⅳ.第一主成分解釋率為18.5%(圖中左上圖第一條黑色柱所示),第二主成分解釋率為12.5%(圖中左上圖第二條黑色柱所示)Fig.3 Spatio-temporal characteristics for phytoplankton in Niyang River based on PCA

    用PCA 方法分析了尼洋河沿程、季節(jié)以及河道底質(zhì)方面浮游植物群落結(jié)構(gòu)特征(圖3)發(fā)現(xiàn),尼洋河下游浮游植物群落較其他河段囊括的范圍大,尼洋河中游浮游植物群落囊括的范圍最小.尼洋河夏季的浮游植物群落結(jié)構(gòu)與其他季節(jié)的分離,說(shuō)明夏季的浮游植物群落較為獨(dú)特,秋季的浮游植物群落囊括范圍較春季和冬季大,春季囊括的游植物群落范圍最小.尼洋河底質(zhì)為黏土的河段浮游植物群落囊括的范圍遠(yuǎn)遠(yuǎn)大于底質(zhì)為砂石的河段,可見底質(zhì)為黏土的河道對(duì)浮游植物群落結(jié)構(gòu)有著很重要的影響.

    尼洋河沿程浮游植物群落4 個(gè)生物多樣性指標(biāo)不存在顯著差異(P >0.05),物種豐富度、總豐度以及Shannon-Wiener 多樣性指數(shù)呈現(xiàn)一致的趨勢(shì),即:尼洋河下游(采樣點(diǎn)Ⅲ)最高,中上游以及交匯處呈現(xiàn)下降趨勢(shì),物種豐富度平均值在20 左右,總豐度在2×105cells/L 左右,Shannon-Wiener 多樣性指數(shù)在2.2 左右.Pielou 均勻度指數(shù)則在尼洋河中游出現(xiàn)峰值,其他河段呈下降趨勢(shì),Pielou 均勻度指數(shù)平均值在0.75 左右,交匯處的Shannon-Wiener 多樣性指數(shù)以及Pielou 均勻度指數(shù)較其他河段均低(圖4).

    圖4 尼洋河浮游植物時(shí)空變化特征(用Duncan 法檢驗(yàn)各采樣點(diǎn)以及季節(jié)之間浮游植物生物多樣性差異性,不同字母表示差異達(dá)顯著水平(P <0.05))Fig.4 Spatio-temporal dynamic characteristics for phytoplankton in Niyang River

    圖5 基于CCA 方法分析尼洋河硅藻(a)、藍(lán)藻(b)、綠藻(c)與環(huán)境因子、樣點(diǎn)之間的關(guān)系Fig.5 Relationships among phytoplankton(a:Bacillariophyta;b:Cyanophyta;c:Chlorophyta),environment factors and sampling sites in Niyang River based on CCA

    4 個(gè)多樣性指標(biāo)在季節(jié)之間則有較大的變化幅度.夏季的浮游植物物種豐富度與其他季節(jié)存在顯著差異(P <0.05),這個(gè)季節(jié)的總豐度與秋、冬季存在顯著差異(P <0.05),物種豐富度與總豐度在該季節(jié)均為最大值,分別約為25 和3×105cells/L,物種豐富度最小值出現(xiàn)在春季,剛剛超過(guò)15,總豐度最小值出現(xiàn)在冬季,在105cells/L 以下.冬季的浮游植物Shannon-Wiener 多樣性指數(shù)在顯著差異(P <0.05),冬季 Pielou 均勻度指數(shù)與其他 3 個(gè)季節(jié)的存在顯著差異(P <0.05),Shannon-Wiener 多樣性指數(shù)與Pielou 均勻度指數(shù)在該季節(jié)均為最大值,Shannon-Wiener 多樣性指數(shù)接近 2.4,Pielou 均勻度指數(shù)剛愈0.85,Shannon-Wiener 多樣性指數(shù)最低值出現(xiàn)在春季,在2.0 以下,Pielou 均勻度指數(shù)最低值出現(xiàn)在夏季,在0.70 以上.

    2.3 尼洋河浮游植物與環(huán)境因子的關(guān)系及關(guān)鍵預(yù)測(cè)因子

    為了更加清晰地表達(dá)尼洋河浮游植物與環(huán)境因子之間的關(guān)系,將浮游植物主要的3 個(gè)類群與9項(xiàng)環(huán)境指標(biāo)進(jìn)行了CCA 分析,浮游植物群落包括硅藻、藍(lán)藻和綠藻,9 項(xiàng)環(huán)境指標(biāo)包括表層水溫、表層pH,硬度、礦化度、表層溶解氧、總氮、銨態(tài)氮、總磷、總堿度.結(jié)果顯示,硅藻門部分藻類與環(huán)境因子存在較大的關(guān)聯(lián),如:舟形藻科的美壁藻以及羽紋藻、雙菱藻科的雙菱藻,這些藻類的密度與理化因子銨態(tài)氮、表層pH 值、表層水溫有著很大的關(guān)聯(lián);曲殼藻科的卵形藻的密度與表層溶解氧以及總堿度有著很大的關(guān)聯(lián).藍(lán)藻門的部分藻類與環(huán)境因子存在著較大的關(guān)聯(lián),如:色球藻科的微囊藻以及顫藻科的束藻的密度與總堿度、總磷相關(guān)聯(lián),管孢藻科的管孢藻的密度與銨態(tài)氮有著較大的關(guān)聯(lián),念珠藻科的念珠藻的密度與表層pH 值以及表層水溫有著較大的關(guān)聯(lián).綠藻門僅有卵囊藻科的卵囊藻的密度與銨態(tài)氮存在較大的關(guān)聯(lián),鼓藻科的新月藻的密度與表層水溫存在較大的關(guān)聯(lián)(圖5).

    選擇浮游植物總豐度、Shannon-Wiener 多樣性指數(shù)、Pielou 均勻度指數(shù)3 項(xiàng)指標(biāo),采用CART 模型預(yù)測(cè)浮游植物群落與環(huán)境因子的相互作用關(guān)系,結(jié)果顯示,尼洋河浮游植物總豐度受到pH 值、季節(jié)以及河道底質(zhì)等因素的影響(圖6).具體來(lái)講,當(dāng)pH 值≥8.0 且為夏、冬季時(shí),浮游植物總豐度為102541.66 cells/L;當(dāng)pH 值≥8.0 且為秋季和春季時(shí),浮游植物總豐度為145666.65 cells/L;當(dāng)pH 值<8.0 且河道底質(zhì)為黏土?xí)r,浮游植物總豐度為37950 cells/L;當(dāng)pH 值<8.0、河道底質(zhì)為砂石且為秋季時(shí),浮游植物總豐度為57499.98 cells/L;當(dāng)pH 值<8.0、河道底質(zhì)為砂石且季節(jié)為夏季和冬季時(shí),浮游植物總豐度為43700.03 cells/L;當(dāng)pH 值<8.0、河道底質(zhì)為砂石且季節(jié)為春季時(shí),浮游植物總豐度為161766.67 cells/L.綜上所述,尼洋河浮游植物總豐度受到pH 值和底質(zhì)的影響較大.

    圖6 基于分類回歸樹分析尼洋河浮游植物總豐度(a)、Shannon-Wiener 多樣性指數(shù)(b)、Pielou 均勻度指數(shù)(c)與環(huán)境因子之間的關(guān)系(分支左側(cè)代表?xiàng)l件成立,右側(cè)代表相反條件成立)Fig.6 Relationship between total abundance(a),Shannon-Wiener diversity index(b),Pielou evenness index(c)of phytoplankton and environment factors in Niyang River based on CART

    尼洋河浮游植物Shannon-Wiener 多樣性指數(shù)(H')受到總堿度、pH 值、海拔以及底質(zhì)的影響,具體來(lái)講,當(dāng)總堿度≥42.22 mg/L 且海拔≥3308 m 時(shí),H'為2.056;當(dāng)總堿度≥42.22 mg/L 且海拔 <3308 m 時(shí),H'為1.949;當(dāng)總堿度 <42.22 mg/L 且 pH <7.8 時(shí),H'為 1.912;當(dāng)總堿度 <42.22 mg/L、pH≥7.8,同時(shí)海拔 <2932 m 時(shí),H'為 1.732;當(dāng)總堿度 <42.22 mg/L、pH≥7.8、海拔≥2932 m,同時(shí)河道底質(zhì)為砂石時(shí),H'為2.245.當(dāng)總堿度<42.22 mg/L、pH≥7.8、海拔≥2932 m,同時(shí)河道底質(zhì)為黏土?xí)r,H'為2.334.綜上所述,尼洋河浮游植物H'受到河道底質(zhì)的影響較大,底質(zhì)為黏土的水域浮游植物H'較底質(zhì)為砂石的大.

    尼洋河浮游植物J'受到總氮、海拔、河道底質(zhì)、季節(jié)以及pH 值的影響,具體來(lái)說(shuō),當(dāng)總氮<5.82 mg/L 且海拔 <3086 m 時(shí),J'為 0.700;當(dāng)總氮 <5.82 mg/L 且海拔≥3086 m 時(shí),J'為 0.698;當(dāng)總氮≥5.82 mg/L 且河道底質(zhì)為黏土?xí)r,J'為0.706;當(dāng)總氮≥5.82 mg/L、河道底質(zhì)為砂石且季節(jié)為秋季時(shí),J'為0.656;當(dāng)總氮≥5.82 mg/L、河道底質(zhì)為砂石、季節(jié)為春、夏、冬季,同時(shí)pH≥8.0 時(shí),J'為0.658;當(dāng)總氮≥5.82 mg/L、河道底質(zhì)為砂石、季節(jié)為春、夏、冬季,同時(shí)pH <8.0 時(shí),J'為0.721.綜上所述,尼洋河浮游植物J'受到pH 值影響較大,pH 值低于8.0 的水域浮游植物J'比pH 值高于8.0 的大.

    3 討論

    演替是一個(gè)特定水體內(nèi)的物理(如光、溫度)、化學(xué)(營(yíng)養(yǎng)物、水質(zhì)、毒素)和生物(競(jìng)爭(zhēng)、捕食)因素的改變所引致的種類變化[18].浮游植物群落存在季節(jié)演替[29],外界環(huán)境條件的改變較之群落初始物種組成對(duì)浮游植物群落的演替順序更具決定性作用[30],溫度、氮磷比和水體穩(wěn)定性是影響浮游植物群落演替的主要因子[31].

    海拔決定了一個(gè)地區(qū)的溫度和光照等環(huán)境因素的變化,屬于宏觀尺度的環(huán)境因子,通常情況下,由于海拔的升高,水域溫度降低,冰凍期延長(zhǎng)[32-33],物種的豐富度也隨之降低[34],浮游植物的總豐度呈現(xiàn)下降趨勢(shì)[35],李芳[36]指出,尼洋河浮游植物總豐度與海拔呈負(fù)相關(guān).海拔升高,水溫降低,一般認(rèn)為,水溫降到浮游植物的最適溫度以下時(shí),浮游植物的生長(zhǎng)逐漸減緩直至停止生長(zhǎng)、甚至消失,其光合作用的強(qiáng)度也隨之逐漸降低[37].

    一方面,水溫是影響浮游植物種類組成及生物量的關(guān)鍵因子[38].隨著水溫的升高,大量浮游植物得以快速生長(zhǎng)和繁殖[39],也有資料表明水溫與浮游植物豐度呈負(fù)相關(guān)[40],原因在于研究區(qū)域(長(zhǎng)江近海岸)的固定優(yōu)勢(shì)種(具槽帕拉藻Paralia sulcata)易聚集于透明度<10 m、溫度<15℃的水體之中,合適的較低水溫使得優(yōu)勢(shì)種在水體上部得以生長(zhǎng)和繁殖.由于浮游植物物種豐富度、總豐度以及多樣性指數(shù)所闡述的生物學(xué)意義有所不同,對(duì)水溫的響應(yīng)有所區(qū)別,吳衛(wèi)菊等[23]指出,水溫變化與浮游植物生物量呈極顯著正相關(guān)關(guān)系,與浮游植物物種豐富度呈極顯著負(fù)相關(guān)關(guān)系,與浮游植物生物多樣性呈顯著負(fù)相關(guān)關(guān)系.

    另一方面,隨著氣候變暖,浮游植物群落的演替表現(xiàn)出物種對(duì)溫度的生態(tài)適應(yīng)性,這種演替對(duì)水域生態(tài)系統(tǒng)結(jié)構(gòu)和功能有著深遠(yuǎn)的意義[41].水溫高于14℃時(shí),微囊藻生物量迅速增加;束絲藻生物量隨溫度增加而增加,高于22℃時(shí)迅速減少[38].平均溫度低于20℃時(shí),藍(lán)藻沒有大規(guī)模生長(zhǎng),硅藻門、綠藻門生物量急劇降低,溫度超過(guò)25℃時(shí),藍(lán)藻迅速增長(zhǎng)并很快成為絕對(duì)優(yōu)勢(shì)種[42].

    中間高度膨脹(mid-altitude bulge)理論指出,往往在中間海拔區(qū)域內(nèi),物種的多樣性較高,在低海拔或者高海拔區(qū)域內(nèi),物種的多樣性較低[43].本研究中,尼洋河浮游植物多樣性(Shannon-Wiener 多樣性指數(shù)和Pielou 均勻度指數(shù))符合中間高度膨脹理論,表現(xiàn)為尼洋河沿程浮游植物多樣性變化特點(diǎn)基本一致,即:尼洋河中下游(中游)浮游植物Shannon-Wiener 多樣性指數(shù)(Pielou 均勻度指數(shù))最大,其他河段Shannon-Wiener多樣性指數(shù)呈下降趨勢(shì).也有資料表明,海拔與浮游植物Shannon-Wiener 多樣性指數(shù)和Pielou 均勻度指數(shù)呈正相關(guān)[44],初步推斷出現(xiàn)這一矛盾在于海拔的梯度是否達(dá)到浮游植物多樣性峰值轉(zhuǎn)折點(diǎn).

    不容忽視的是,營(yíng)養(yǎng)鹽在影響浮游植物生物量的同時(shí),對(duì)浮游植物的群落結(jié)構(gòu)也產(chǎn)生影響,不同浮游植物類群對(duì)營(yíng)養(yǎng)鹽的敏感性不同[45].值得強(qiáng)調(diào)的是,N、P 含量及其組成與浮游植物關(guān)系較為密切,是影響浮游植物演替的重要因子.不同氮磷比對(duì)浮游植物群落的多樣性指數(shù)、物種組成及演替過(guò)程均有顯著影響[46].另外,在同一季節(jié),硅藻和甲藻利用不同層次的氮鹽,能有效減少兩者間的競(jìng)爭(zhēng)壓力[47].不同的硅藻種類對(duì)環(huán)境因子敏感性有所差異[48],如運(yùn)動(dòng)型(motile ecological guild)的硅藻對(duì)溫度較為敏感,而喜靜水型(low profile guild)以及喜流水型(high profile guild)的硅藻則對(duì)總氮、總磷等理化因子較為敏感[49].

    實(shí)際上,浮游植物群落結(jié)構(gòu)往往是多維環(huán)境因子綜合作用的結(jié)果,如本研究中,硅藻門舟形藻科的美壁藻以及羽紋藻、雙菱藻科的雙菱藻,這些藻類密度與理化因子銨態(tài)氮、表層pH 值、表層水溫有著很大的關(guān)聯(lián),曲殼藻科的卵形藻密度與表層溶解氧以及總堿度有著很大的關(guān)聯(lián);藍(lán)藻門色球藻科的微囊藻以及顫藻科的束藻密度與總堿度、總磷相關(guān)聯(lián),管孢藻科的管孢藻密度與銨態(tài)氮有著較大的關(guān)聯(lián),念珠藻科的念珠藻密度與表層pH 值以及表層水溫有著較大的關(guān)聯(lián);綠藻僅有卵囊藻科的卵囊藻密度與銨態(tài)氮存在較大的關(guān)聯(lián),鼓藻科的新月藻密度與表層水溫存在較大的關(guān)聯(lián).

    研究結(jié)果表明,尼洋河浮游植物類群以硅藻占有絕對(duì)優(yōu)勢(shì),綠藻和藍(lán)藻伴隨存在,這主要由于尼洋河一年四季均以冷水為主,較為適合喜冷水的硅藻生長(zhǎng)與繁殖[50].在尼洋河沿程方面,這3 類藻類演替受到河流的形態(tài)特征的影響,尼洋河下游較其它河段平緩,形成了獨(dú)特的“尼洋河河谷”風(fēng)景區(qū),其河道底質(zhì)以泥沙和黏土為主,具有較大的比表面積,更容易為浮游植物的生長(zhǎng)和繁殖提供各類營(yíng)養(yǎng)鹽[51],同時(shí)由于河谷河段水流較其他河段平緩,減少了因水流給浮游植物帶來(lái)的生存威脅[49],與海拔因素相比,尼洋河浮游植物物種豐富度或者總豐度對(duì)底質(zhì)的響應(yīng)更多一些,即以泥沙和黏土為底質(zhì)的河道浮游植物物物種豐富度以及總豐度較砂石為底質(zhì)的河道更大.在季節(jié)方面,這3 類藻類演替主要受到外源性補(bǔ)給水的影響,包括冰雪融水和降水,這些氣候因素是水體營(yíng)養(yǎng)鹽,如硅磷比,發(fā)生改變的主要驅(qū)動(dòng)因子[52],導(dǎo)致在夏季硅藻和綠藻總豐度上升為最大值.同時(shí),由于受到環(huán)境同質(zhì)性的影響,藍(lán)藻呈下降趨勢(shì),在冬季,由于外源性補(bǔ)給水的減少,尼洋河水體受外界干擾作用減弱,浮游植物多樣性,包括Shannon-Wiener 多樣性指數(shù)和Pielou 均勻度指數(shù),均達(dá)到了最大值.

    [1]秦大河,陳宜瑜,李學(xué)勇等.中國(guó)氣候與環(huán)境演變(上卷):中國(guó)氣候與環(huán)境的演變與預(yù)測(cè).北京:科學(xué)出版社,2005:389-390.

    [2]Ju JT,Zhu LP,F(xiàn)eng JL et al.Hydrodynamic process of Tibetan Plateau lake revealed by grain size:Case study of Pumayum Co.Chinese Sci Bull,2012,57(19):2433-2441.

    [3]賀桂芹.西藏高寒濕地生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估及濕地保護(hù)對(duì)策研究[學(xué)位論文].楊凌:西北農(nóng)林科技大學(xué),2007.

    [4]Wu WH,Xu SJ,Lu HY et al.Mineralogy,major and trace element geochemistry of riverbed sediments in the headwaters of the Yangtze,Tongtian River and Jinsha River.J Asian Earth Sci,2011,40(2):611-621.

    [5]Wu WH,Zheng HB,Xu SJ et al.Geochemistry and provenance of bed sediments of the large rivers in the Tibetan Plateau and Himalayan region.Int J Earth Sci,2012,101(5):1357-1370.

    [6]Bai JH,Cui BS,Chen B et al.Spatial distribution and ecological risk assessment of heavy metals in surface sediments from a typical plateau lake wetland,China.Ecol Model,2009,222(2):301-306.

    [7]謝 虹.青藏高原蒸散發(fā)及其對(duì)氣候變化的響應(yīng)(1970-2010)[學(xué)位論文].蘭州:蘭州大學(xué),2012.

    [8]張?jiān)萍t,許長(zhǎng)軍,亓 青等.青藏高原氣候變化及其生態(tài)效應(yīng)分析.青海大學(xué)學(xué)報(bào):自然科學(xué)版,2011,29(4):18-22.

    [9]沈永平,王根緒,吳青柏等.長(zhǎng)江-黃河源區(qū)未來(lái)氣候情景下的生態(tài)環(huán)境變化.冰川凍土,2002,24(3):308-314.

    [10]黃 琦.基于GIS 的三江源地區(qū)生態(tài)環(huán)境變化與人類活動(dòng)影響研究[學(xué)位論文].北京:中央民族大學(xué),2012.

    [11]秦大河,陳宜瑜,李學(xué)勇等.中國(guó)氣候與環(huán)境演變(下卷):氣候與環(huán)境變化的影響與適應(yīng)、減緩對(duì)策.北京:科學(xué)出版社,2005:98-109.

    [12]翟紅娟.縱向嶺谷區(qū)水電工程脅迫對(duì)河流生態(tài)完整性影響的研究[學(xué)位論文].北京:北京師范大學(xué),2009.

    [13]段順瓊,王 靜,馮少輝等.云南高原湖泊地區(qū)水資源脆弱性評(píng)價(jià)研究.中國(guó)農(nóng)村水利水電,2011,9:55-59.

    [14]胡元林,鄭 文.高原湖泊流域可持續(xù)發(fā)展研究.生態(tài)經(jīng)濟(jì),2011,236(3):168-171,183.

    [15]洛?!れ`智多杰.青藏高原水資源的保護(hù)與利用.資源科學(xué),2005,27(2):23-27.

    [16]關(guān)志華,陳傳友.西藏河流與湖泊.北京:科學(xué)出版社,1984.

    [17]沈紅保,郭 麗.西藏尼洋河魚類組成調(diào)查與分析.河北漁業(yè),2008,173(5):51-54,60.

    [18]劉建康.高級(jí)水生生物學(xué).北京:科學(xué)出版社,1999:176-198.

    [19]胡 蓮,潘曉潔,鄒 曦等.三道河水庫(kù)浮游植物群落結(jié)構(gòu)特征及其漁產(chǎn)潛力分析.水生態(tài)學(xué)志,2012,33(4):90-95.

    [20]劉足根,張 柱,張 萌等.贛江流域浮游植物群落結(jié)構(gòu)與功能類群劃分.長(zhǎng)江流域資源與環(huán)境,2012,21(3):375-384.

    [21]Sommer U,Aberle N,Lengfellner K et al.The Baltic Sea spring phytoplankton bloom in a changing climate:an experimental approach.Mar Biol,2012,159(11):2479-2490.

    [22]劉東艷.膠州灣浮游植物與沉積物中硅藻群落結(jié)構(gòu)演替的研究[學(xué)位論文].青島:中國(guó)海洋大學(xué),2004.

    [23]吳衛(wèi)菊,楊 凱,汪志聰?shù)?云貴高原漁洞水庫(kù)浮游植物群落結(jié)構(gòu)及季節(jié)演替.水生態(tài)學(xué)雜志,2012,33(2):69-75.

    [24]國(guó)家環(huán)境保護(hù)總局《水和廢水監(jiān)測(cè)分析方法》編委會(huì).水和廢水監(jiān)測(cè)分析方法:第4 版.北京:中國(guó)環(huán)境科學(xué)出版社,2002:45-46,109,121-124.

    [25]朱蕙忠,陳嘉佑.中國(guó)西藏硅藻.北京:科學(xué)出版社,2000.

    [26]胡鴻鈞,李堯英,魏心印等.中國(guó)淡水藻類.上海:上海科學(xué)技術(shù)出版社,1980.

    [27]Shannon EE,Weaver W.The mathematical theory of communication.London:University Illinois Press,1949:125.

    [28]Pielou EC.Ecological diversity.New York:John Wiley,1975:165.

    [29]Romanov RE,Kirillov VV.Analysis of the seasonal dynamics of river phytoplankton based on succession rate indices for key event identification.Hydrobiologia,2012,695(1):293-304.

    [30]黃 偉,朱旭宇,曾江寧等.氮磷比對(duì)浙江近岸浮游植物群落結(jié)構(gòu)影響的實(shí)驗(yàn)研究.海洋學(xué)報(bào),2012,34(5):128-138.

    [31]岳 強(qiáng),黃 成,史元康等.廣東南水水庫(kù)富營(yíng)養(yǎng)化與浮游植物群落動(dòng)態(tài).環(huán)境科學(xué)與技術(shù),2012,34(8):112-116.

    [32]Rundle SD,Jenkins A,Ormerod SJ.Macroinvertebrate communities in streams in Himalaya,Nepal.Freshwater Biol,1993,30(1):169-180.

    [33]Ao M,Alfred JRB,Gupta A.Studies on some lotic systems in the north-eastern hill regions of India.Limnologica,1984,15(1):135-141.

    [34]Suren AM.Macroinvertebrate communities of streams in western Nepal:effects of altitude and land use.Freshwater Biol,1994,32(2):323-336.

    [35]張軍燕.瑪曲至湖口段黃河干流浮游生物群落結(jié)構(gòu)特征研究[學(xué)位論文].西安:西北大學(xué),2009.

    [36]李 芳.西藏尼洋河流域水生生物研究及水電工程對(duì)其影響的預(yù)測(cè)評(píng)價(jià)[學(xué)位論文].西安:西北大學(xué),2009.

    [37]Talling JF.The relative growth rate of three plankton diatoms in relation to underwater radiation and temperature.Ann Bot,1995,19(2):329-341.

    [38]代龔圓,李 杰,李 林等.滇池北部湖區(qū)浮游植物時(shí)空格局及相關(guān)環(huán)境因子.水生生物學(xué)報(bào),2012,36(5):946-956.

    [39]張才學(xué),龔玉艷,孫省利.湛江港灣潛在赤潮生物的時(shí)空分布及其影響因素.生態(tài)學(xué)雜志,2012,31(7):1763-1770.

    [40]郭術(shù)津,孫 軍,戴民漢等.2009年冬季東海浮游植物群集.生態(tài)學(xué)報(bào),2012,32(10):3266-3278.

    [41]Winder M,Sommer U.Phytoplankton response to a changing climate.Hydrobiologia,2012,698(1):5-16.

    [42]蔡琳琳,朱廣偉,朱夢(mèng)圓等.太湖梅梁灣湖岸帶浮游植物群落演替及其與水華形成的關(guān)系.生態(tài)科學(xué),2012,31(4):345-351.

    [43]Zhao CM,Chen WL,Tian ZQ et al.Altitudinal pattern of plant species diversity in Shennongjia Mountains,central China.J Integr Plant Biol,2005,47(12):1431-1449.

    [44]劉 艷.額爾齊斯河及鄰近內(nèi)陸河流域浮游植物生態(tài)學(xué)研究[學(xué)位論文].上海:上海海洋大學(xué),2011.

    [45]賴俊翔,俞志明,宋秀賢等.利用特征色素研究長(zhǎng)江口海域浮游植物對(duì)營(yíng)養(yǎng)鹽加富的響應(yīng).海洋科學(xué),2012,36(5):42-52.

    [46]Huang W,Zhu XY,Zeng JN et al.Responses in growth and succession of the phytoplankton community to different N/P ratios near Dongtou Island in the East China Sea.J Exp Mar Biol Ecol,2012,434:102-109.

    [47]袁 騏,王云龍,沈新強(qiáng).N 和P 對(duì)東海中北部浮游植物的影響研究.海洋環(huán)境科學(xué),2005,24(4):5-8.

    [48]宋書群.黃、東海浮游植物功能群研究[學(xué)位論文].青島:中國(guó)科學(xué)院海洋研究所,2010.

    [49]Stenger KC,Lengyel E,Crossetti LO et al.Diatom ecological guilds as indicators of temporally changing stressors and disturbances in the small Torna-stream,Hungary.Ecol Indic,2013,24:138-147.

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

    [51]張 智,張顯忠,楊駿驊.植物凈化床對(duì)雙龍湖水體有機(jī)污染物的去除效果分析.生態(tài)環(huán)境,2006,15(4):708-713.

    [52]Morabito G,Oggioni A,Austoni M.Resource ratio and human impact:how diatom assemblages in Lake Maggiore responded to oligotrophication and climatic variability.Hydrobiologia,2012,698(1):47-60.

    猜你喜歡
    尼洋河底質(zhì)綠藻
    不同發(fā)育階段中華絨螯蟹對(duì)底質(zhì)的喜好性研究
    一種ReliefF和隨機(jī)森林模型組合的多波束海底底質(zhì)分類方法
    用于海底目標(biāo)識(shí)別與底質(zhì)分類的多波束水體波形預(yù)處理
    尼洋河:風(fēng)光旖旎的天然畫卷
    桃花節(jié)赴林芝
    香榧綠藻的生物學(xué)特性及物種鑒定
    無(wú)綠藻的微生物學(xué)特性和菌種鑒定
    文蛤的底質(zhì)選擇性及潛沙能力研究
    全球近10年無(wú)綠藻病例報(bào)道文獻(xiàn)的回顧
    以成長(zhǎng)為主題解讀《窗燈》
    国产淫片久久久久久久久| 久久精品熟女亚洲av麻豆精品| 成人国产麻豆网| 熟女av电影| 一区二区三区免费毛片| 一本色道久久久久久精品综合| 亚洲成色77777| 成人鲁丝片一二三区免费| 伊人久久精品亚洲午夜| 一区二区三区四区激情视频| 亚洲综合色惰| 亚洲经典国产精华液单| 久久久久久久大尺度免费视频| 最近最新中文字幕免费大全7| 岛国毛片在线播放| 国产精品伦人一区二区| 亚洲色图综合在线观看| 少妇人妻久久综合中文| 伊人久久精品亚洲午夜| 成年人午夜在线观看视频| 中文精品一卡2卡3卡4更新| 日产精品乱码卡一卡2卡三| 国产成人午夜福利电影在线观看| 日韩欧美精品免费久久| 亚洲欧洲国产日韩| 欧美日韩综合久久久久久| av在线老鸭窝| 麻豆乱淫一区二区| 色播亚洲综合网| 国产av不卡久久| 亚洲成人av在线免费| 欧美极品一区二区三区四区| 91久久精品国产一区二区成人| 国产精品一区二区性色av| 国产色爽女视频免费观看| 亚洲国产色片| 听说在线观看完整版免费高清| 精华霜和精华液先用哪个| 免费观看av网站的网址| 免费看光身美女| 国产精品国产三级专区第一集| 日本黄色片子视频| 日韩av不卡免费在线播放| 一级黄片播放器| 国产成人91sexporn| 亚洲欧美日韩东京热| 真实男女啪啪啪动态图| 毛片女人毛片| 久久ye,这里只有精品| 97人妻精品一区二区三区麻豆| 熟女电影av网| 国产淫语在线视频| 熟女人妻精品中文字幕| 亚洲精品,欧美精品| 欧美潮喷喷水| 日本欧美国产在线视频| 国产男人的电影天堂91| 日本免费在线观看一区| 亚洲人成网站高清观看| 久久久久九九精品影院| 一级毛片电影观看| 亚洲天堂av无毛| 国产成人一区二区在线| 日本三级黄在线观看| 午夜老司机福利剧场| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久色成人| 中国国产av一级| 有码 亚洲区| 日本与韩国留学比较| 中文资源天堂在线| 亚洲精品乱码久久久v下载方式| 亚洲精品日本国产第一区| 永久网站在线| 久久精品久久久久久久性| 日本爱情动作片www.在线观看| 精品久久久久久久久av| 在线观看av片永久免费下载| 一级爰片在线观看| 免费看av在线观看网站| 一本一本综合久久| 国产精品女同一区二区软件| 汤姆久久久久久久影院中文字幕| 2021少妇久久久久久久久久久| 亚洲精品成人久久久久久| 国产伦理片在线播放av一区| 麻豆乱淫一区二区| 99九九线精品视频在线观看视频| av免费在线看不卡| 男女啪啪激烈高潮av片| 夜夜看夜夜爽夜夜摸| 天天一区二区日本电影三级| 99热这里只有是精品50| 狂野欧美激情性bbbbbb| 一级a做视频免费观看| 男人舔奶头视频| 国产精品一区二区在线观看99| 又爽又黄a免费视频| 欧美一区二区亚洲| 日韩欧美精品v在线| 三级经典国产精品| 日韩一区二区视频免费看| 91精品国产九色| 久久久久网色| 麻豆国产97在线/欧美| 精品人妻视频免费看| 伦精品一区二区三区| 欧美高清性xxxxhd video| 天天躁日日操中文字幕| 亚洲国产欧美人成| 国产一区二区在线观看日韩| 一区二区三区精品91| 亚洲欧美精品自产自拍| 国产精品一区二区在线观看99| 观看免费一级毛片| 狂野欧美白嫩少妇大欣赏| 嫩草影院入口| 久久综合国产亚洲精品| 日产精品乱码卡一卡2卡三| 午夜福利网站1000一区二区三区| 高清视频免费观看一区二区| 另类亚洲欧美激情| 国产成人精品一,二区| 成人亚洲精品av一区二区| 国产69精品久久久久777片| 在线免费十八禁| 能在线免费看毛片的网站| 菩萨蛮人人尽说江南好唐韦庄| 欧美三级亚洲精品| 欧美一级a爱片免费观看看| 老司机影院成人| 97超视频在线观看视频| 中文字幕av成人在线电影| 国产91av在线免费观看| 我要看日韩黄色一级片| 草草在线视频免费看| 国产精品一区二区三区四区免费观看| 人人妻人人看人人澡| 在线亚洲精品国产二区图片欧美 | 热re99久久精品国产66热6| 成人漫画全彩无遮挡| 一边亲一边摸免费视频| 97超视频在线观看视频| 午夜福利在线在线| 精品熟女少妇av免费看| 看十八女毛片水多多多| 久久久欧美国产精品| 国产色婷婷99| 大香蕉97超碰在线| 国产成人免费无遮挡视频| 97在线视频观看| 一级毛片 在线播放| 黄片wwwwww| 午夜日本视频在线| 在线观看一区二区三区激情| 日韩成人av中文字幕在线观看| 人人妻人人澡人人爽人人夜夜| 又爽又黄a免费视频| 丝瓜视频免费看黄片| 九九爱精品视频在线观看| 国产v大片淫在线免费观看| 亚洲精品国产av成人精品| 青春草亚洲视频在线观看| 嫩草影院精品99| 国产国拍精品亚洲av在线观看| av在线观看视频网站免费| 1000部很黄的大片| 在线a可以看的网站| 国产高清不卡午夜福利| 中文字幕久久专区| 亚洲欧美一区二区三区国产| 国产亚洲精品久久久com| 欧美zozozo另类| 18禁在线播放成人免费| 欧美一区二区亚洲| 国产精品成人在线| 亚洲欧美一区二区三区黑人 | 毛片一级片免费看久久久久| 日韩人妻高清精品专区| 国产成人精品福利久久| 国产大屁股一区二区在线视频| 国产精品久久久久久精品电影小说 | 国产精品一及| 自拍欧美九色日韩亚洲蝌蚪91 | 天堂中文最新版在线下载 | 日本一本二区三区精品| 欧美97在线视频| 韩国高清视频一区二区三区| 亚洲欧美中文字幕日韩二区| 亚洲va在线va天堂va国产| 日本爱情动作片www.在线观看| 99re6热这里在线精品视频| 成人毛片60女人毛片免费| 欧美潮喷喷水| 国产91av在线免费观看| 国产淫片久久久久久久久| 一本一本综合久久| 午夜老司机福利剧场| 亚洲精品国产av成人精品| 黄色视频在线播放观看不卡| 亚洲最大成人手机在线| 六月丁香七月| 久久97久久精品| 亚洲精品国产成人久久av| 丰满人妻一区二区三区视频av| 最新中文字幕久久久久| 亚洲三级黄色毛片| 高清在线视频一区二区三区| 国产 一区 欧美 日韩| 最近最新中文字幕大全电影3| 又黄又爽又刺激的免费视频.| 成人一区二区视频在线观看| 人妻 亚洲 视频| a级毛色黄片| 少妇 在线观看| 国产老妇伦熟女老妇高清| 久久久久精品性色| 国产 精品1| 国产黄色视频一区二区在线观看| 久久久色成人| 欧美xxxx性猛交bbbb| 亚洲精品成人av观看孕妇| 麻豆精品久久久久久蜜桃| 一级av片app| 国产亚洲一区二区精品| 免费看av在线观看网站| 免费黄网站久久成人精品| av在线老鸭窝| 美女视频免费永久观看网站| 黄色视频在线播放观看不卡| 在线免费十八禁| 国产高清有码在线观看视频| 国产一区二区亚洲精品在线观看| 777米奇影视久久| 汤姆久久久久久久影院中文字幕| 欧美区成人在线视频| 肉色欧美久久久久久久蜜桃 | 边亲边吃奶的免费视频| 人妻夜夜爽99麻豆av| 天堂中文最新版在线下载 | 国产精品99久久久久久久久| 国产毛片a区久久久久| 少妇高潮的动态图| 国产精品久久久久久久久免| 青春草国产在线视频| 成年免费大片在线观看| 男人爽女人下面视频在线观看| 综合色av麻豆| 99热网站在线观看| 99热网站在线观看| 国产男人的电影天堂91| 男女国产视频网站| 噜噜噜噜噜久久久久久91| 91久久精品国产一区二区三区| 日韩成人av中文字幕在线观看| 国产一区二区亚洲精品在线观看| 日产精品乱码卡一卡2卡三| 国产女主播在线喷水免费视频网站| 尾随美女入室| 国产精品国产av在线观看| 九九爱精品视频在线观看| 免费黄网站久久成人精品| 国产一区有黄有色的免费视频| 蜜桃亚洲精品一区二区三区| 国产精品蜜桃在线观看| 国产色爽女视频免费观看| 别揉我奶头 嗯啊视频| 黄色视频在线播放观看不卡| 日本猛色少妇xxxxx猛交久久| 午夜亚洲福利在线播放| 一区二区三区乱码不卡18| 日韩亚洲欧美综合| 国产精品嫩草影院av在线观看| 性插视频无遮挡在线免费观看| videossex国产| 我要看日韩黄色一级片| 搡女人真爽免费视频火全软件| 在线播放无遮挡| 国产大屁股一区二区在线视频| 免费少妇av软件| 日韩在线高清观看一区二区三区| 国产一区二区三区综合在线观看 | 亚洲av电影在线观看一区二区三区 | 在线观看一区二区三区激情| 久久久久久久国产电影| 男女无遮挡免费网站观看| 内地一区二区视频在线| 免费黄网站久久成人精品| 在线观看一区二区三区激情| 国产精品人妻久久久久久| 日本熟妇午夜| 肉色欧美久久久久久久蜜桃 | 亚洲熟女精品中文字幕| 精品国产三级普通话版| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产成人一精品久久久| 热re99久久精品国产66热6| a级毛色黄片| 熟妇人妻不卡中文字幕| 国产色爽女视频免费观看| 国产免费又黄又爽又色| 在线观看人妻少妇| 久久久精品免费免费高清| 成年女人看的毛片在线观看| 黄色一级大片看看| 内射极品少妇av片p| 国产毛片a区久久久久| 只有这里有精品99| 少妇人妻久久综合中文| 精品久久久久久电影网| 久久ye,这里只有精品| 国产免费视频播放在线视频| 最近的中文字幕免费完整| 内地一区二区视频在线| 亚洲三级黄色毛片| 色综合色国产| 亚洲真实伦在线观看| 色吧在线观看| 欧美性感艳星| 嫩草影院入口| 丝袜脚勾引网站| 欧美极品一区二区三区四区| 欧美变态另类bdsm刘玥| 成人欧美大片| 国产白丝娇喘喷水9色精品| 91午夜精品亚洲一区二区三区| 日韩中字成人| 国产成人福利小说| 毛片一级片免费看久久久久| 久久精品夜色国产| 久久久久久久久久久丰满| 中文资源天堂在线| 免费观看av网站的网址| 国内揄拍国产精品人妻在线| 99久久精品热视频| 有码 亚洲区| 欧美少妇被猛烈插入视频| 一级二级三级毛片免费看| 国产综合精华液| 久久精品久久精品一区二区三区| 亚洲精品国产av成人精品| 五月天丁香电影| 97精品久久久久久久久久精品| 国产 一区 欧美 日韩| 人妻一区二区av| 国产白丝娇喘喷水9色精品| 自拍偷自拍亚洲精品老妇| 国产视频首页在线观看| 最近的中文字幕免费完整| 久久亚洲国产成人精品v| 亚洲av中文av极速乱| 男女边吃奶边做爰视频| 国产美女午夜福利| 欧美变态另类bdsm刘玥| 两个人的视频大全免费| 少妇人妻久久综合中文| 黄色怎么调成土黄色| 亚洲精品影视一区二区三区av| 婷婷色av中文字幕| 亚洲人成网站在线观看播放| 制服丝袜香蕉在线| 免费电影在线观看免费观看| 六月丁香七月| 女的被弄到高潮叫床怎么办| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久久久久免费av| 最近中文字幕2019免费版| 亚洲成人精品中文字幕电影| 男女国产视频网站| 亚洲四区av| 熟妇人妻不卡中文字幕| 国产精品不卡视频一区二区| 永久网站在线| 少妇丰满av| av在线天堂中文字幕| 精品人妻熟女av久视频| 国产一区有黄有色的免费视频| 成年版毛片免费区| 日日啪夜夜撸| 夫妻午夜视频| 亚洲精品国产av成人精品| 99热这里只有精品一区| 高清视频免费观看一区二区| 亚洲av中文av极速乱| 亚洲欧洲日产国产| 欧美日韩在线观看h| av免费在线看不卡| 亚洲aⅴ乱码一区二区在线播放| 国产老妇伦熟女老妇高清| 亚洲人成网站在线播| 久久久久久久久久成人| 午夜福利在线在线| 精品午夜福利在线看| 亚洲高清免费不卡视频| 欧美日韩亚洲高清精品| 久久精品国产亚洲av涩爱| 久久99精品国语久久久| 精品久久国产蜜桃| 一区二区三区免费毛片| 欧美日韩亚洲高清精品| 日日摸夜夜添夜夜添av毛片| 国产高潮美女av| 啦啦啦在线观看免费高清www| 亚洲国产最新在线播放| 久久精品国产亚洲网站| 极品教师在线视频| 一级爰片在线观看| 亚洲国产日韩一区二区| 九九在线视频观看精品| 99热国产这里只有精品6| 啦啦啦在线观看免费高清www| 狂野欧美激情性xxxx在线观看| 天天躁日日操中文字幕| 韩国av在线不卡| 日韩成人伦理影院| 亚洲精品乱码久久久久久按摩| 日韩不卡一区二区三区视频在线| 丰满少妇做爰视频| 又黄又爽又刺激的免费视频.| 欧美日韩综合久久久久久| 免费看a级黄色片| 一级二级三级毛片免费看| 久久精品国产亚洲网站| 久久亚洲国产成人精品v| 欧美精品一区二区大全| 2018国产大陆天天弄谢| 精品一区二区三区视频在线| 亚洲美女搞黄在线观看| 免费大片黄手机在线观看| 日本色播在线视频| 免费观看性生交大片5| 国产精品伦人一区二区| 特级一级黄色大片| 蜜臀久久99精品久久宅男| av一本久久久久| 精品久久久精品久久久| 丝袜喷水一区| 成人毛片60女人毛片免费| 日韩三级伦理在线观看| 国产老妇伦熟女老妇高清| 成人一区二区视频在线观看| 在线观看av片永久免费下载| 国产亚洲午夜精品一区二区久久 | 国产久久久一区二区三区| 18+在线观看网站| 免费看a级黄色片| av在线天堂中文字幕| 韩国av在线不卡| 免费看光身美女| 国产精品偷伦视频观看了| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av天美| 欧美日韩一区二区视频在线观看视频在线 | 日本与韩国留学比较| 久久久久久久久久久免费av| 精华霜和精华液先用哪个| 日韩强制内射视频| 97精品久久久久久久久久精品| 夜夜爽夜夜爽视频| 在线看a的网站| 日韩大片免费观看网站| 欧美xxxx性猛交bbbb| 亚洲无线观看免费| 99热国产这里只有精品6| 六月丁香七月| 亚洲av成人精品一区久久| 两个人的视频大全免费| 精品久久久精品久久久| 特级一级黄色大片| 免费电影在线观看免费观看| 综合色av麻豆| 69av精品久久久久久| 精品一区二区三卡| 国产一区二区三区av在线| 最近手机中文字幕大全| 亚洲精华国产精华液的使用体验| 又大又黄又爽视频免费| 国产淫片久久久久久久久| 精品少妇黑人巨大在线播放| 国产在视频线精品| 少妇的逼水好多| 亚洲av男天堂| 超碰97精品在线观看| 免费观看av网站的网址| 少妇裸体淫交视频免费看高清| 日韩视频在线欧美| 久久精品久久久久久噜噜老黄| 国产黄色视频一区二区在线观看| 在线看a的网站| 国产成人a∨麻豆精品| 久热这里只有精品99| 国产精品av视频在线免费观看| 日韩一区二区三区影片| 99久久九九国产精品国产免费| 高清av免费在线| videos熟女内射| 亚洲欧美精品专区久久| av女优亚洲男人天堂| 在线a可以看的网站| 国产亚洲最大av| 丝袜喷水一区| 草草在线视频免费看| 亚洲成人久久爱视频| 看十八女毛片水多多多| 我要看日韩黄色一级片| 精品人妻一区二区三区麻豆| 亚洲成人精品中文字幕电影| 国产亚洲精品久久久com| 日本熟妇午夜| 日韩强制内射视频| 国产精品麻豆人妻色哟哟久久| 少妇被粗大猛烈的视频| 99热这里只有是精品在线观看| 免费观看av网站的网址| 青青草视频在线视频观看| 可以在线观看毛片的网站| 成人国产av品久久久| 国产亚洲91精品色在线| 国产爽快片一区二区三区| 深夜a级毛片| 97精品久久久久久久久久精品| 大香蕉久久网| 国产男女内射视频| 国产久久久一区二区三区| 22中文网久久字幕| 国语对白做爰xxxⅹ性视频网站| 国产成人午夜福利电影在线观看| 国产伦精品一区二区三区视频9| 99热全是精品| 国产熟女欧美一区二区| 色视频在线一区二区三区| 国产成人精品一,二区| 少妇人妻一区二区三区视频| 久久久久精品性色| 国产黄片美女视频| 91在线精品国自产拍蜜月| kizo精华| 亚洲欧洲日产国产| 男人添女人高潮全过程视频| 亚洲国产av新网站| 特大巨黑吊av在线直播| 国产黄频视频在线观看| 黄色视频在线播放观看不卡| 肉色欧美久久久久久久蜜桃 | 18+在线观看网站| 2018国产大陆天天弄谢| 国产精品麻豆人妻色哟哟久久| 久久99热这里只有精品18| 晚上一个人看的免费电影| 街头女战士在线观看网站| 水蜜桃什么品种好| 午夜福利网站1000一区二区三区| 久久国内精品自在自线图片| 三级男女做爰猛烈吃奶摸视频| 女人十人毛片免费观看3o分钟| 国产毛片在线视频| av一本久久久久| 国产探花在线观看一区二区| 亚洲精品色激情综合| 免费黄频网站在线观看国产| 3wmmmm亚洲av在线观看| 三级国产精品片| 国产黄色免费在线视频| 老司机影院成人| 制服丝袜香蕉在线| 插阴视频在线观看视频| 大陆偷拍与自拍| 丰满少妇做爰视频| 老司机影院毛片| 国产免费一区二区三区四区乱码| h日本视频在线播放| 性插视频无遮挡在线免费观看| 成人免费观看视频高清| 丝袜美腿在线中文| 成人免费观看视频高清| 久久97久久精品| 久久久久国产网址| 久久99热这里只频精品6学生| 91狼人影院| 亚洲精品久久午夜乱码| 在线精品无人区一区二区三 | 自拍欧美九色日韩亚洲蝌蚪91 | 3wmmmm亚洲av在线观看| 久热久热在线精品观看| 亚洲成人精品中文字幕电影| 欧美丝袜亚洲另类| 国产亚洲5aaaaa淫片| a级毛色黄片| 只有这里有精品99| 国产视频内射| 一区二区三区精品91| a级毛色黄片| 日韩欧美精品v在线| 如何舔出高潮| 黄片无遮挡物在线观看| 免费看av在线观看网站| 成人黄色视频免费在线看| 国产一级毛片在线| 精品人妻视频免费看| 特大巨黑吊av在线直播| 精品久久久久久久久亚洲| 亚洲av不卡在线观看| 3wmmmm亚洲av在线观看| 熟妇人妻不卡中文字幕| 国产伦精品一区二区三区视频9| 亚洲美女搞黄在线观看| 菩萨蛮人人尽说江南好唐韦庄| 伊人久久国产一区二区| 七月丁香在线播放| 国产淫语在线视频| 舔av片在线| 亚洲国产日韩一区二区| 嫩草影院入口| 全区人妻精品视频| 亚洲av免费高清在线观看|