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

    呼倫湖及其周圍水域春季浮游植物功能群及其影響因子

    2020-06-28 08:21:02李星醇于洪賢竇華山潘海峰馬成學(xué)
    水產(chǎn)學(xué)雜志 2020年3期
    關(guān)鍵詞:呼倫湖需氧量水域

    李星醇,于洪賢,竇華山,潘海峰,馬成學(xué)

    (1.東北林業(yè)大學(xué)野生動(dòng)物與自然保護(hù)地學(xué)院,黑龍江 哈爾濱 150040;2.內(nèi)蒙古達(dá)賚湖(呼倫湖)國家級(jí)自然保護(hù)區(qū)管理局,內(nèi)蒙古 海拉爾 021008)

    浮游植物是水域生態(tài)系統(tǒng)中重要的初級(jí)生產(chǎn)者,其群落結(jié)構(gòu)和物種變化、演替與環(huán)境因子關(guān)系極為密切,可直接而靈敏地反映水環(huán)境的動(dòng)態(tài)變化[1,2]。因此,對(duì)浮游植物的研究有助于評(píng)價(jià)水體水質(zhì)的營養(yǎng)狀態(tài)[3]。傳統(tǒng)的浮游植物群落研究以種屬為基本分類單元[4],忽略了浮游植物群體的環(huán)境生態(tài)特征,在具體應(yīng)用中暴露出很多不足[5]。2002 年Reynolds 等[5]提出浮游藻類功能群理論,2009 年P(guān)adisák 等[6]在Reynolds 的基礎(chǔ)上進(jìn)行了補(bǔ)充和完善。迄今這種生態(tài)功能分類方法共劃分和描述了39組浮游藻類功能群。這種方法有利于更加準(zhǔn)確地描述浮游植物的生境特征,預(yù)測(cè)特定生境下浮游植物的群落分布狀況[7,8]。

    近年來,對(duì)熱帶、亞熱帶、溫帶地區(qū)各種不同類型水域的各個(gè)浮游植物功能群進(jìn)行了大量研究[9-12],但對(duì)寒冷地區(qū)浮游植物功能群的相關(guān)研究較少[13,14]。呼倫湖作為北方寒冷堿性水體的代表,目前相關(guān)研究還主要集中在氣候[15]、湖泊面積[16]以及環(huán)境變化等[17],運(yùn)用浮游植物功能群來響應(yīng)呼倫湖環(huán)境變化的研究很少,且重點(diǎn)都只放到了呼倫湖湖域[18]。近年來,呼倫湖水位逐年下降,水體營養(yǎng)鹽程度急劇升高,整個(gè)水體環(huán)境遭到顯著破壞,周圍河流流域也受到相應(yīng)污染。要從源頭控制污染,必須將呼倫湖的周圍水域也納入研究之中。本文通過探究呼倫湖及其周圍水域浮游植物功能群及其與環(huán)境因子之間的關(guān)系,以期為寒區(qū)高堿度水體浮游植物功能群的研究及呼倫湖的保護(hù)工作提供重要理論依據(jù)。

    1 材料與方法

    1.1 研究地概況

    呼倫湖(48°33'~49°20'N,116°58'~117°48'E)又名達(dá)賚湖,位于內(nèi)蒙古呼倫貝爾市新巴爾虎左旗、新巴爾虎右旗和滿洲里市之間[19],是內(nèi)蒙古自治區(qū)第一大湖,中國第五大湖,素有“草原明珠”的美譽(yù)[20]。呼倫湖湖面呈不規(guī)則斜長方形,長軸為西南至東北方向,湖長93km,最大寬度為41km,湖周長447km,總面積為2339km2。湖面水位最高時(shí)海拔545.59m,最大水深為10m,平均水深為5.7m,蓄水量為138.5億m3。呼倫湖流域?qū)儆谥袦貛О敫珊挡菰瓪夂?,氣候寒冷干燥,冬長夏短,春季干旱多風(fēng)沙,秋季霜凍早,降溫快[21],封凍期達(dá)170~180d,最大冰厚為1.3m。呼倫湖西岸為山巒和懸崖陡壁,湖底為泥底;東岸和南岸地勢(shì)平坦開闊,湖底多為細(xì)砂礫層[20]。湖泊主要的補(bǔ)給來源除了降水和地下水外,還來自于周圍的4 條河流,分別為烏爾遜河、克魯倫河、海拉爾河和新開河。

    烏爾遜河發(fā)源于貝爾湖北岸,北流注入呼倫湖,長223km。上游為額爾古納河水系的哈拉哈河,該河位于呼倫貝爾市新巴爾虎右旗和新巴爾虎左旗交界處。河流兩岸地形平坦,河寬60~70m。河漫灘上葦柳叢生,濕地較多,河兩岸牧業(yè)較發(fā)達(dá)。

    克魯倫河發(fā)源于蒙古國的肯特山東麓。流經(jīng)呼倫貝爾盟新巴爾虎右旗,東流注入呼倫湖。流域面積7153km2,兩岸為半荒漠的低山圍繞,河谷寬約35km,河寬60~70m。兩岸沼澤濕地多,較高的階地上生長著優(yōu)良牧草,牧業(yè)發(fā)達(dá)。11 月到次年4 月結(jié)冰。上游用于灌溉,流送木材。沿岸牧草豐富,自古為重要農(nóng)牧業(yè)地帶。

    海拉爾河位于中國內(nèi)蒙古自治區(qū)呼倫貝爾市境內(nèi),是額爾古納河的上游河段。歷史上曾是呼倫湖補(bǔ)充流量的重要水源之一。海拉爾河發(fā)源于大興安嶺西側(cè),呈東至西流向,干流全長714.9km,流域面積5.45 萬km2。年中洪峰有兩次:5 月融雪期和8月夏雨期[22]。2010 年8 月,由于氣候干旱等原因,呼倫湖水位持續(xù)下降,使呼倫湖內(nèi)濕地萎縮,湖區(qū)生態(tài)環(huán)境退化。國家環(huán)境保護(hù)開展“引河濟(jì)湖”工程。從海拉爾河引水,恢復(fù)呼倫湖水位。自此,海拉爾河通過引河濟(jì)湖工程成了呼倫湖的重要水源之一。

    新開河位于呼倫湖東北部,河長25km,是1965年為保護(hù)扎賚諾爾煤礦安全,人工修建的運(yùn)河。新開河的流向受呼倫湖和海拉爾河水位變化的制約,順逆不定,是呼倫湖的吞吐河流。

    呼倫湖與周圍4 條河流特殊的連通關(guān)系,使得這片由五大水系組成的水域飽含神秘色彩。它遠(yuǎn)離海洋,冷暖氣流交匯頻繁,在調(diào)節(jié)氣候、涵養(yǎng)水源、防止荒漠化、保護(hù)生物多樣性乃至我國北方的生態(tài)安全中發(fā)揮著不可替代的作用。同時(shí),作為響應(yīng)全球氣候變化最敏感的地區(qū)之一,它在呼倫貝爾市的生態(tài)保護(hù)中更發(fā)揮著相當(dāng)重要的作用。

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

    2018 年5 月在呼倫湖及其周圍水域采集水樣及浮游植物樣品??紤]到呼倫湖與周圍4 條河流特殊的連通關(guān)系,將整個(gè)調(diào)查水域劃分為5 大水系(呼倫湖體、克魯倫河、烏爾遜河、新開河、海拉爾河)。為了能夠較為系統(tǒng)地表征呼倫湖及其周圍水域浮游植物的分布情況,遵循均勻布點(diǎn),使采集的樣品遍布全區(qū),具有全區(qū)代表性的原則,共設(shè)置28個(gè)采樣點(diǎn)。其中,在呼倫湖水體的中心區(qū)(M22-湖心3、M23-湖心2、M24-湖心1)、主要進(jìn)出水口附近(M7-新開河入(出)湖口、M13-克魯倫河入湖口、M19-烏爾遜河入湖口、M21-海拉爾河入湖口)和沿岸區(qū)(西岸M8-M12;南岸M15;東岸M16-M18、M20)設(shè)代表性的采樣點(diǎn)共計(jì)17 個(gè);在克魯倫河水域設(shè)置了1 個(gè)采樣點(diǎn)(M14);在烏爾遜河水域設(shè)置4 個(gè)采樣點(diǎn)(M25、M26、M27、M28),在新開河設(shè)置1個(gè)采樣點(diǎn)(M6),在海拉爾河水域設(shè)置5 個(gè)采樣點(diǎn)(M1-M5)(圖1、表1)。

    圖1 春季呼倫湖及其周圍水域浮游植物采樣點(diǎn)分布圖Fig.1 Sampling site in Hulun Lake and surrounding waters in spring

    表1 春季呼倫湖及其周圍水域浮游植物采樣點(diǎn)布設(shè)情況Tab.1 Distribution of phytoplankton sampling sites in Hulun Lake and surrounding waters in spring

    1.3 浮游植物樣品的采集與鑒定

    浮游植物樣品用1L 采水器采集上、中、下水層混合水樣1L,用15%魯哥氏液現(xiàn)場(chǎng)向樣品中加入15mL 以固定,帶回實(shí)驗(yàn)室,沉淀48h,濃縮至30mL,用浮游植物計(jì)數(shù)框在光學(xué)顯微鏡下進(jìn)行觀察和統(tǒng)計(jì)。浮游植物種類根據(jù)文獻(xiàn)[23,24]進(jìn)行鑒定,并根據(jù)不同浮游植物的體積相關(guān)公式計(jì)算生物量[25]。

    1.4 水環(huán)境因子的測(cè)定

    水溫(WT)用表面水溫計(jì)測(cè)定,pH、溶解氧(DO)、電導(dǎo)率(Cond)、氨氮(NH4+-N)等理化指標(biāo)使用YSI-6600 多功能水質(zhì)分析儀現(xiàn)場(chǎng)測(cè)定;總氮(TN)、總磷(TP)、化學(xué)需氧量(CODCr)等理化指標(biāo)根據(jù)《水和廢水監(jiān)測(cè)分析方法》[26]測(cè)定。

    1.5 數(shù)據(jù)處理與分析

    1.5.1 浮游植物優(yōu)勢(shì)度指數(shù)(y)

    y=fi×Pi

    式中,fi為第i 種的出現(xiàn)頻率,Pi為第i 種個(gè)體數(shù)占所有點(diǎn)位中總個(gè)體數(shù)的比例,y 為優(yōu)勢(shì)度。當(dāng)y>0.02 時(shí),為浮游植物優(yōu)勢(shì)種。

    1.5.2 Pearson 分析

    運(yùn)用SPSS19.0 軟件,通過Pearson 線性相關(guān)分析研究浮游植物重要功能群同水環(huán)境因子間的相關(guān)性。

    1.5.3 RDA(冗余分析)分析

    對(duì)浮游植物功能群的生物量數(shù)據(jù)和環(huán)境數(shù)據(jù)進(jìn)行Lg(X+1)的量綱化處理后,采用Canoco for Windows 4.5 軟件對(duì)浮游植物功能群數(shù)據(jù)進(jìn)行DCA 分析。結(jié)果表明,物種SD(梯度)的單峰響應(yīng)值大于2而小于3。蒙特卡洛置換試驗(yàn)后,采用RDA 進(jìn)一步分析,研究功能群特征與環(huán)境因子之間的相關(guān)性。

    2 結(jié)果和分析

    2.1 春季呼倫湖及其周圍水域浮游植物組成

    經(jīng)鏡檢,在呼倫湖及其周圍水域中共檢測(cè)出浮游植物7 門80 種(含變種及變型),其中呼倫湖6門62 種,綠藻門最多,共鑒定出32 種,硅藻門次之,共鑒定出19 種,藍(lán)藻門共鑒定出7 種,種類數(shù)最低的是裸藻門、金藻門和隱藻門,分別鑒定出2種、1 種和1 種。而呼倫湖周圍河流流域共鑒定出7門72 種,硅藻門28 種,綠藻門26 種,藍(lán)藻門9 種,種類數(shù)最少的依然是金藻門、裸藻門、隱藻門和甲藻門,分別為3 種、3 種、2 種、1 種。各門浮游植物的生物量是:河流系統(tǒng)主要以硅藻門為主,而湖泊系統(tǒng)則以綠藻門為主(圖2)。

    圖2 春季呼倫湖及其周圍水域浮游植物生物量Fig.2 Phytoplankton biomass in Hulun Lake and its surrounding waters in spring

    根據(jù)不同種類浮游植物在各采樣點(diǎn)出現(xiàn)的頻率與豐度,春季在呼倫湖及其周圍水域共發(fā)現(xiàn)優(yōu)勢(shì)種3 門7 種(表2),其中以綠藻門為主,共鑒定出4種,硅藻門2 種,藍(lán)藻門1 種,硅藻門的梅尼小環(huán)藻Cyclotella meneghiniana 優(yōu)勢(shì)度最高,達(dá)到0.239,頻率為1,即在每一個(gè)采樣點(diǎn)中都有分布,其次為綠藻門的四足十字藻Crucigenia tetrapedia,優(yōu)勢(shì)度為0.161。

    表2 春季呼倫湖及其周圍水域浮游植物優(yōu)勢(shì)種Tab.2 Phytoplankton dominant species in Hulun Lake and its surrounding waters in spring

    2.2 春季呼倫湖及其周圍水域浮游植物功能群生物量的水平分布

    呼倫湖及周圍水域的浮游植物共劃分出20 個(gè)功能群,分別為L0、S1、H1、MP、M、P、TB、D、C、X3、Y、X2、W1、W2、X1、J、F、G、N 和Ws(表3),將各采樣點(diǎn)中浮游植物生物量大于10%的功能群定義為重要功能群,則有4 個(gè)重要功能群,分別為C、X2、J、D,其余功能群劃分為其他??傮w來看,春季呼倫湖及其周圍水域浮游植物功能群水平分布差異顯著,各采樣點(diǎn)間各功能群的生物量不同(圖3)。春季呼倫湖及其周圍水域浮游植物的生物量在0.70~11.73mg/L 之間波動(dòng)。最高點(diǎn)為東烏珠爾(M4),為11.73mg/L,次高為嵯崗(M5),為11.17mg/L,這兩個(gè)采樣點(diǎn)都位于海拉爾河下游。最低為湖的西岸M11,為0.70mg/L。從整體趨勢(shì)上看,M4、M5 生物量尤其突出,其他采樣點(diǎn)生物量均穩(wěn)定在一定水平。

    圖3 春季呼倫湖及其周圍水域浮游植物重要功能群組成及其生物量變化Fig.3 Composition and biomass changes of phytoplankton important functional groups in Hulun Lake and its surrounding waters in spring

    由于河流系統(tǒng)和湖泊系統(tǒng)存在差異,因此浮游植物功能群的分布也不同。呼倫湖及其周圍水域浮游植物重要功能群的分布情況見圖4。結(jié)合圖3 可以發(fā)現(xiàn),河流流域浮游植物生物量遠(yuǎn)大于湖內(nèi),各重要功能群的分布也存在顯著差異。在河流生態(tài)系統(tǒng)中,功能群C 分布最廣,在海拉爾河(M4、M5)分布最多,僅是在這兩個(gè)樣點(diǎn)的分布即可達(dá)到總分布量的59%;而在克魯倫河(M14)所占比重極大,可高達(dá)該樣點(diǎn)總生物量的75.2%。在湖泊生態(tài)系統(tǒng)中,功能群J 分布最廣,在湖內(nèi)各樣點(diǎn)的分布較為平均,平均生物量可達(dá)0.47mg/L,遠(yuǎn)超河流流域(0.10mg/L)。其次,功能群D 在河流生態(tài)系統(tǒng)中分布也較廣,僅在海拉爾河(M4、M5)的分布量可達(dá)總分布量的70.6%;對(duì)于X2 功能群來說,在湖泊系統(tǒng)分布的生物量要略大于河流系統(tǒng),海拉爾河(M4)和貝爾湖(M28)兩采樣點(diǎn)共占據(jù)了高達(dá)39%的X2 功能群的生物量。

    表3 春季呼倫湖及其周圍水域浮游植物重要功能群組成Tab.3 Composition of important phytoplankton functional groups in Hulun Lake and its surrounding waters in spring

    圖4 春季呼倫湖及其周圍水域不同水體類型中浮游植物重要功能群的分布Fig.4 Distribution of phytoplankton important functional groups in Hulun Lake and its surrounding waters in spring

    2.3 呼倫湖及其周圍水域水體理化參數(shù)

    春季呼倫湖及其周圍水域水環(huán)境因子中電導(dǎo)率、總磷、總氮含量和化學(xué)需氧量水平分布差異顯著,而水溫、pH、溶解氧和氨氮含量水平分布差異不顯著(表4)。電導(dǎo)率在193~1786μs/cm 之間,湖區(qū)(平均1646.88μs/cm)明顯高于入湖河流(平均410.55μs/cm)??偭撞▌?dòng)于0.04~0.48mg/L,東西岸(M9~M13、M16)和湖中心(M22~M24)顯著高于其他采樣點(diǎn),均顯著高于國家V 類水標(biāo)準(zhǔn)(0.2mg/L)。呼和諾爾(M3)和烏蘭諾爾(M27)也超過國家V 類水標(biāo)準(zhǔn)(0.2mg/L);總氮波動(dòng)于0.3~3.4mg/L,呼倫湖西岸(M8~M9)、東岸(M16、M18、M19)以及嵯崗M5均超過國家V 類水標(biāo)準(zhǔn)(2mg/L),化學(xué)需氧量波動(dòng)于7~72mg/L,湖內(nèi)(平均54.35mg/L)明顯高于入湖河流(平均28.45mg/L),湖內(nèi)(M7~M13,M15~M24)17 個(gè)采樣點(diǎn)均超過國家V 類水標(biāo)準(zhǔn)(40mg/L),呼和諾爾(M3)和烏蘭諾爾(M27)也均超過國家V 類水標(biāo)準(zhǔn)(40mg/L)??傮w來看,湖區(qū)與入湖河流的理化指標(biāo)之間差異顯著。呼倫湖內(nèi)(M9、M16)水質(zhì)極差,總氮、總磷和化學(xué)耗氧量均超過國家V 類水標(biāo)準(zhǔn)(2mg/L、0.2mg/L、40mg/L),呼倫湖內(nèi)西岸(M10~M13)、湖中心(M22~M24)以及周圍水域的呼和諾爾(M3)、烏蘭諾爾(M27)水質(zhì)也不容樂觀,總磷和化學(xué)耗氧量污染極為嚴(yán)重。

    表4 春季(5 月)呼倫湖及其周圍水域各采樣點(diǎn)水化學(xué)因子Tab.4 Hydrochemical factors of Hulun Lake and its surrounding waters in spring(May)

    2.4 浮游植物功能群同水環(huán)境因子的關(guān)系

    2.4.1 重要功能群同水環(huán)境因子的Pearson 相關(guān)分析

    重要功能群的生物量同水環(huán)境因子的相關(guān)性見表5。功能群C 的生物量與電導(dǎo)率、化學(xué)需氧量呈極顯著負(fù)相關(guān);功能群X2 的生物量與溫度、溶解氧呈極顯著正相關(guān),與電導(dǎo)率呈顯著負(fù)相關(guān);功能群J的生物量與pH、電導(dǎo)率和化學(xué)需氧量呈極顯著正相關(guān);功能群D 的生物量與電導(dǎo)率和化學(xué)需氧量呈極顯著負(fù)相關(guān)。

    2.4.2 重要功能群同水環(huán)境因子的RDA 分析

    對(duì)物種數(shù)據(jù)進(jìn)行除趨勢(shì)(DCA)結(jié)果顯示,四個(gè)軸梯度長度最大值為2.102<3,故對(duì)數(shù)據(jù)進(jìn)行冗余分析(RDA)。選擇4 個(gè)重要功能群,環(huán)境因子包括溫度(WT)、酸堿度(pH)、溶解氧(DO)、電導(dǎo)率(Cond)、氨氮(NH4+-N)、總磷(TP)、總氮(TN)和化學(xué)需氧量(CODCr)共8 個(gè)。RDA 分析結(jié)果如表6 和圖5 所示。軸1、軸2 共解釋了物種數(shù)據(jù)累計(jì)變化率的75.8%,表明浮游植物重要功能群與環(huán)境變量之間存在顯著相關(guān)性。從軸1 來看,最主要的負(fù)相關(guān)因子為電導(dǎo)率(Cond),相關(guān)系數(shù)達(dá)-0.7313,其次為化學(xué)需氧量(CODCr),相關(guān)系數(shù)達(dá)-0.6850;從軸2看,最主要的正相關(guān)因子是溫度(WT),相關(guān)系數(shù)達(dá)0.6176,其次為pH,相關(guān)系數(shù)達(dá)0.5749。浮游植物重要功能群X2、C、D 與溶解氧(DO)、溫度(WT)呈正相關(guān),同電導(dǎo)率(Cond)、化學(xué)需氧量(CODCr)、氨氮(NH4+-N)等呈負(fù)相關(guān)。重要功能群J 則相反,同溫度(WT)、溶解氧(DO)呈負(fù)相關(guān),與化學(xué)需氧量(CODCr)、電導(dǎo)率(Cond)等呈正相關(guān)??偟═N)對(duì)浮游植物功能群生長影響較小,不是限制因子。

    圖5 重要功能群同水環(huán)境因子的RDA 分析Fig.5 RDA analysis of functional groups and environment factors

    表5 重要功群生物量同水環(huán)境因子的Pearson 相關(guān)分析Tab.5 Significant Pearson correlation between the biomass of the main phytoplankton functional groups and environmental variables

    表6 排序軸特征值、種類與環(huán)境因子排序軸的相關(guān)系數(shù)Tab.6 Redundancy analysis results for phytoplankton functional groups

    3 討論

    3.1 春季呼倫湖及其周圍水域浮游植物功能群的特點(diǎn)

    本研究發(fā)現(xiàn),春季呼倫湖及其周圍水域浮游植物的優(yōu)勢(shì)功能群為C、X2、J、D,這與優(yōu)勢(shì)種所對(duì)應(yīng)的功能群相似,但并不完全一致。優(yōu)勢(shì)種尖針桿藻Synedra acus、梅尼小環(huán)藻Cyclotella meneghiniana、四足十字藻Crucigenia tetrapedia 對(duì)應(yīng)的功能群與重要功能群一致,但其他優(yōu)勢(shì)種對(duì)應(yīng)的功能群并不是重要功能群。春季呼倫湖水域河流生態(tài)系統(tǒng)主要以功能群C 和D 為主,春季呼倫湖水域湖泊生態(tài)系統(tǒng)主要以功能群J 為主。功能群X2 在河流和湖泊系統(tǒng)中均有分布,但在湖泊的分布量略多。功能群C以華麗星桿藻Asterionella formosa 和梅尼小環(huán)藻為代表,功能群D 以尖針桿藻、肘狀針桿藻Synedra ulna 和平片針桿藻Synedra tabulata 為代表。藻功能群C 和功能群D 對(duì)低光照都有一定的耐受性[27],二者均為硅藻。溫度是影響浮游植物功能群分布的重要因素之一,它通過控制細(xì)胞內(nèi)化學(xué)反應(yīng)強(qiáng)度來影響藻類的生長[28]。一般而言,藍(lán)藻和綠藻適宜生長在較高的水溫中,硅藻則恰好相反,更適宜生活在較低的水溫中[29]。春季呼倫湖及其周圍水域體溫較低,透明度較低,濁度高,光照弱,功能群C 和功能群D 耐低溫和耐低光照的特點(diǎn)使其成為春季的重要功能群。除溫度外,水動(dòng)力學(xué)也是影響浮游植物功能群分布的主要因素之一。河流區(qū)水流更為湍急,懸浮顆粒更多[30,31],因此更適宜于硅藻的生長[32],功能群J 是以十字藻和柵藻為代表的綠藻,適合生長在混合型的富營養(yǎng)水體中,其適宜生長的環(huán)境與pH 有正相關(guān)性[33]。呼倫湖內(nèi),尤其是呼倫湖西岸和湖中心營養(yǎng)鹽顯著高于周圍河流水域,總磷和化學(xué)需氧量顯著超標(biāo),水體富營養(yǎng)嚴(yán)重,因此,功能群J在呼倫湖內(nèi)占顯著優(yōu)勢(shì)。另一方面,呼倫湖作為北方典型偏堿性水體,對(duì)比湖泊和周圍河流流域,呼倫湖內(nèi)水體pH 平均值可高達(dá)8.92,高于周圍河流流域的pH(8.56),也有利于功能群J 的生長。由卵形衣藻和球衣藻為主構(gòu)成的功能群X2,在春季也占一定優(yōu)勢(shì),X2 功能群適合生長在中營養(yǎng)型到高度營養(yǎng)型的易受擾動(dòng)的淺水水體中,對(duì)濾食十分敏感[14]。呼倫湖及周圍水域生活污水和畜牧糞便的大量排放,導(dǎo)致水體出現(xiàn)富營養(yǎng)化,使水中大型濾食性浮游動(dòng)物顯著減少,浮游動(dòng)物的濾食功能顯著降低,更適宜功能群X2 的生長,這也與韓博平等[34]的研究結(jié)果一致。

    3.2 春季呼倫湖及其周圍水域浮游植物功能群同環(huán)境因子之間的關(guān)系

    溫度、光照、水體動(dòng)力學(xué)特征、營養(yǎng)鹽狀況、浮游動(dòng)物牧食等是影響浮游植物出現(xiàn)或消失的重要因素,這些因素能在特定生境中影響浮游植物功能群的組成[5,35]。無論是南方熱帶、亞熱帶水域還是北方寒冷地區(qū),營養(yǎng)鹽通常是最主要的驅(qū)動(dòng)因子[36-38]。同樣,呼倫湖作為北方內(nèi)陸典型寒冷地區(qū)高堿性水體,營養(yǎng)鹽仍是夏季浮游植物生長的最主要驅(qū)動(dòng)因子[13],但春季則明顯不同。本研究中Pearson 相關(guān)分析和RDA 分析均表明,化學(xué)需氧量和電導(dǎo)率是影響呼倫湖及其周圍水域春季浮游植物功能群組成及演替的主要環(huán)境因子。春季呼倫湖周圍水域和呼倫湖體水域化學(xué)需氧量和電導(dǎo)率差異顯著,呼倫湖湖體水域化學(xué)需氧量和電導(dǎo)率均顯著高于呼倫湖周圍水域。

    化學(xué)需氧量主要是水樣中所具有的還原性物質(zhì)被強(qiáng)氧化劑(如重鉻酸鉀)氧化的量[39]。它能科學(xué)測(cè)量水樣中有機(jī)物的含量[40],直接反映水體受還原性物質(zhì)污染的程度[41]。在許多關(guān)于春季浮游植物功能群和環(huán)境因子關(guān)系的研究中,CODCr都是最主要的影響因子[42,43]。在河流和湖泊的水質(zhì)監(jiān)測(cè)中,若化學(xué)需氧量的數(shù)值較高,也就說明該地區(qū)水質(zhì)含有較多的有機(jī)污染物。本研究發(fā)現(xiàn):呼倫湖內(nèi)的化學(xué)需氧量嚴(yán)重超標(biāo)(平均54.35mg/L),遠(yuǎn)高于周圍河流流域(平均28.45mg/L),這也與文獻(xiàn)[44]的研究相符。正常情況下,水體中的有機(jī)物通常來自于兩個(gè)方面:一是外界向水體中排放的如農(nóng)藥、化肥以及城市污水、工業(yè)廢水一類的有機(jī)物;二是水體底泥釋放及水體中的生物產(chǎn)生。呼倫湖地區(qū)以放牧為主,畜牧業(yè)十分發(fā)達(dá),這就導(dǎo)致湖中有機(jī)物質(zhì)大多以懸浮的植物殘?bào)w存在,在春季溫度較低且供氧充足的條件下,有機(jī)物的生物分解作用緩慢,使湖水中化學(xué)需氧量增加[45]。湖上游有熟皮廠和硝礦,主要污染物是氨氮、農(nóng)藥,毒物污染以氟化物、砷、汞、酚為主[46]。從5 月開始,呼倫湖水域降雨量開始增加,雨量的增加一方面增大了地表沖刷,使攜帶的枯腐植物、牲畜糞便等物質(zhì)增多,另一方面,湖水上漲,淹沒周圍的一些湖灘草地,這也導(dǎo)致其中的一些有機(jī)物隨之溶入湖中,使得呼倫湖內(nèi)有較高有機(jī)物含量?,F(xiàn)場(chǎng)勘測(cè)也發(fā)現(xiàn),呼倫湖內(nèi)水體流速更低,較為渾濁,水體中的泥沙也是測(cè)得較高的CODCr的原因之一。D 功能群耐受擾動(dòng),而水體流動(dòng)會(huì)使有機(jī)質(zhì)含量減少,使水體CODCr降低,因此功能群D 與化學(xué)需氧量(CODCr)顯著負(fù)相關(guān)。更有研究表明,CODCr和水體氮、磷含量呈明顯正相關(guān)[47],這也與本研究結(jié)果相符。高CODCr含量,高氮、磷含量,使水極度富營養(yǎng)化,高濃度的營養(yǎng)鹽水平使得功能群J 得以在湖域大量繁殖,因此功能群J 與化學(xué)需氧量(CODCr)呈極顯著正相關(guān)。

    電導(dǎo)率是與水體中浮游植物生長密切相關(guān)的另一重要環(huán)境變量。測(cè)定電導(dǎo)率可以間接推測(cè)水中離子成分的總濃度[48],判斷溶解性總固體物的濃度,判斷水質(zhì)狀況[49]。一般而言,水體中的電導(dǎo)率同溶解于水中的離子濃度、溫度、碳含量有明顯相關(guān)關(guān)系。呼倫湖內(nèi)接收大量上游及兩岸排放的污染物,水體的離子濃度自然居高不下。經(jīng)過冬季漫長冰封期的積累作用之后,水體滯留時(shí)間較長,因此在湖內(nèi)聚集了大量的有機(jī)質(zhì)。再加上春季隨著氣溫的逐漸升高,水體中光合和代謝作用不斷增強(qiáng),水體中無機(jī)碳含量逐漸升高。整個(gè)湖內(nèi)碳含量升高,因此測(cè)得較高的電導(dǎo)率。本研究中功能群C 大量分布在河流流域,一是因?yàn)楹恿魈己肯鄬?duì)低,功能群C 有耐低碳的特點(diǎn);二也證實(shí)了功能群C 與電導(dǎo)率呈明顯負(fù)相關(guān)關(guān)系。有研究表明,部分水生植物可以降低水體電導(dǎo)率,提高透明度[50]。春季呼倫湖水溫較低,含鹽量較高,這樣的水體特點(diǎn)決定了湖體水草稀少[46],較高的電導(dǎo)率也得不到抑制。通常情況下,水體純凈度越高,電導(dǎo)率越低[51]。研究期間,呼倫湖湖體電導(dǎo)率均值為1646.9μs/cm,周圍河流流域的電導(dǎo)率均值為410.5μs/cm。湖內(nèi)顯著高于河流流域也在一定程度上說明了湖內(nèi)水中溶解鹽含量高,水質(zhì)更差。電導(dǎo)率與水中氮、磷營養(yǎng)鹽的含量密切相關(guān)。通常電導(dǎo)率越大,水體中氮、磷營養(yǎng)鹽含量越高[52]。這也與本研究中的電導(dǎo)率與總氮(0.5233)、總磷(0.4337)呈正相關(guān)以及電導(dǎo)率與功能群J 呈顯著正相關(guān)的結(jié)論相吻合。除了電導(dǎo)率和化學(xué)需氧量外,對(duì)于寒冷地區(qū)的水域來說,溫度變化范圍廣、溫差大,水溫顯著變化影響對(duì)浮游植物功能群的分布。本研究中X2 功能群與溫度就表現(xiàn)出極顯著的正相關(guān)性。RDA 分析結(jié)果中軸2 最主要的正相關(guān)因子分別是溫度和pH,功能群J 對(duì)于pH也表現(xiàn)出了顯著的正相關(guān)性,這也與Sommer 等[33]的研究結(jié)果一致。呼倫湖比其他水體受周圍牧場(chǎng)和人為干擾更多,呼倫湖內(nèi)營養(yǎng)鹽含量高,總氮含量也明顯偏高,所以整體來說氮元素并不缺乏,總氮對(duì)浮游植物功能群影響很小,這與尹子龍[18]的研究結(jié)果相一致。

    3.3 呼倫湖及其周圍水域目前現(xiàn)存問題及相應(yīng)保護(hù)措施

    近些年來,呼倫湖水域污染嚴(yán)重,水質(zhì)持續(xù)惡化,呼倫湖內(nèi)的西岸及湖中心富營養(yǎng)化尤為嚴(yán)重。河流流域的M3(呼和諾爾)、M27(烏蘭諾爾)水質(zhì)極差。湖體的富營養(yǎng)程度要比周圍河流流域嚴(yán)重。地表徑流和干草的大量入湖,湖泊水體中總氮、總磷含量不斷升高,伴隨著蒸發(fā)量的不斷增大、湖內(nèi)水位持續(xù)下降、湖水不能有效流通,湖內(nèi)營養(yǎng)物質(zhì)、鹽分、有機(jī)污染物不斷濃縮,富營養(yǎng)化程度逐漸加深。而周圍河流流域由于旅游業(yè)興起、人為活動(dòng)、工業(yè)畜牧業(yè)污染所導(dǎo)致的水質(zhì)惡化加重了整體水域的富營養(yǎng)情況。這是目前呼倫湖流域所面對(duì)的最嚴(yán)峻生態(tài)問題之一。水體嚴(yán)重的富營養(yǎng)化不僅使得魚類資源大大減少,生物多樣性降低,更引發(fā)了浮游動(dòng)物小型化、藍(lán)藻爆發(fā)、大面積湖濱蘆葦濕地喪失、草場(chǎng)退化、有毒藻類泛濫、漁業(yè)資源瀕臨枯竭、大量珍稀鳥類遷移等一系列問題。

    究其原因,主要有以下兩個(gè)方面。一是自然原因。隨著全球氣候變暖,湖體蒸發(fā)量逐漸增大,入湖徑流的減少,使得湖體蓄水量減少,水體更新速度變慢,自凈能力下降。二是人為原因。湖體的污染來源自于入湖河流和湖體水域。污染不僅局限于水域周邊發(fā)達(dá)的畜牧業(yè),更與人為活動(dòng)不可分割。隨著資源的不斷開發(fā)利用,旅游業(yè)的興起,居民用水量的增加,致使水位下降;生活污染物的大量排放,更使得水質(zhì)大大降低。因此,要想從根本上改善呼倫湖水域的水質(zhì),需從多方面著手。首先,要減少氮磷排放量,降低水體富營養(yǎng)化程度,保護(hù)區(qū)內(nèi)應(yīng)禁止放牧,防止畜牧糞便匯集到水域中使富營養(yǎng)情況加重。其次,周邊應(yīng)禁止開設(shè)旅游區(qū),避免人為因素的干擾。第三,各管理部門應(yīng)加強(qiáng)對(duì)呼倫湖及周圍水域的實(shí)時(shí)監(jiān)測(cè),及時(shí)掌握水域動(dòng)態(tài)。最后,嘗試用不同修復(fù)手段控制水域內(nèi)富營養(yǎng)。如利用工程性措施修復(fù),機(jī)械清淤,生態(tài)方法如植被修復(fù)法、水生動(dòng)物修復(fù)法、微生物修復(fù)法等。

    3.4 結(jié)論

    (1)呼倫湖及其周圍水域共鑒定出7 門80 種浮游植物。呼倫湖湖區(qū)浮游植物群落以綠藻門為主,而周圍河流流域則以硅藻門為主。浮游植物功能群可以劃分為20 個(gè),分別為L0、S1、H1、MP、M、P、TB、D、C、X3、Y、X2、W1、W2、X1、J、F、G、N、Ws,其中重要功能群為C、X2、J、D。

    (2)化學(xué)需氧量(CODCr)和電導(dǎo)率(Cond)是影響呼倫湖及其周圍水域浮游植物功能群最主要的限制因子。

    (3)功能群C 和D 同化學(xué)需氧量和電導(dǎo)率顯著負(fù)相關(guān),功能群J 同化學(xué)需氧量和電導(dǎo)率顯著正相關(guān),而功能群X2 受水溫和溶解氧影響較大,同水溫和溶解氧顯著正相關(guān)。

    (4)GB3838-2002 評(píng)價(jià)結(jié)果表明,呼倫湖內(nèi)(M9、M16)水質(zhì)極差,總氮、總磷和化學(xué)耗氧量均超過國家V 類水標(biāo)準(zhǔn)(2mg/L、0.2mg/L、40mg/L);呼倫湖內(nèi)西岸(M10~M13)、湖中心(M22~M24)以及周圍水域的呼和諾爾(M3)、烏蘭諾爾(M27)水中總磷和化學(xué)需氧量污染極為嚴(yán)重??傮w來看,周圍河流水質(zhì)明顯優(yōu)于呼倫湖內(nèi)水質(zhì)。

    (5)呼倫湖水域水體的富營養(yǎng)化是目前最嚴(yán)重的問題之一。要想從源頭控制污染,采取禁止放牧、禁止開設(shè)旅游區(qū)等措施,嘗試?yán)霉こ绦源胧┬迯?fù)清淤、生態(tài)方法植被修復(fù)、水生動(dòng)物修復(fù)法、微生物修復(fù)法來修復(fù)水體環(huán)境。

    猜你喜歡
    呼倫湖需氧量水域
    坐在呼倫湖畔看星星
    遼河(2022年1期)2022-02-14 19:45:39
    坐在呼倫湖畔看星星
    遼河(2022年1期)2022-02-14 05:15:04
    進(jìn)博會(huì)水域環(huán)境保障研究及展望
    水上消防(2020年5期)2020-12-14 07:16:14
    黃河口附近海域化學(xué)需氧量和石油烴分布及其關(guān)鍵控制環(huán)境因子分析
    柳江水域疍民的歷史往事
    文史春秋(2019年11期)2020-01-15 00:44:36
    新巴爾虎寫意(組詩)
    北極光(2019年5期)2019-11-13 04:54:18
    城市水域生態(tài)景觀設(shè)計(jì)探討
    智能城市(2018年7期)2018-07-10 08:30:00
    《水質(zhì)化學(xué)需氧量的測(cè)定》新舊標(biāo)準(zhǔn)區(qū)別探討
    連續(xù)流動(dòng)分光光度法測(cè)定水中化學(xué)需氧量
    治淮(2013年1期)2013-03-11 20:05:15
    我國水域?qū)垂δ芏ㄎ环诸惐Wo(hù)
    日本免费a在线| 国产高清视频在线播放一区| 一本精品99久久精品77| 黄片大片在线免费观看| 国产伦人伦偷精品视频| 18禁在线播放成人免费| 90打野战视频偷拍视频| 神马国产精品三级电影在线观看| 国产精品一区二区三区四区免费观看 | 天天躁日日操中文字幕| 国产精品久久视频播放| 在线看三级毛片| 欧美大码av| 午夜日韩欧美国产| 欧美激情久久久久久爽电影| 悠悠久久av| 久久久久久人人人人人| 在线天堂最新版资源| 在线观看免费视频日本深夜| 国产精品亚洲美女久久久| 9191精品国产免费久久| 免费av观看视频| 国产午夜精品论理片| 久久精品亚洲精品国产色婷小说| 两个人看的免费小视频| 最近最新中文字幕大全电影3| 在线播放国产精品三级| 少妇人妻精品综合一区二区 | 内地一区二区视频在线| 亚洲欧美激情综合另类| 18禁黄网站禁片午夜丰满| 少妇人妻精品综合一区二区 | 亚洲精品亚洲一区二区| 欧美又色又爽又黄视频| 白带黄色成豆腐渣| 亚洲国产日韩欧美精品在线观看 | 两性午夜刺激爽爽歪歪视频在线观看| 99久久99久久久精品蜜桃| 黑人欧美特级aaaaaa片| www.熟女人妻精品国产| 亚洲不卡免费看| 欧美丝袜亚洲另类 | 国产男靠女视频免费网站| 国产真人三级小视频在线观看| 97超视频在线观看视频| 性色av乱码一区二区三区2| 亚洲欧美日韩高清专用| 久久伊人香网站| 真人做人爱边吃奶动态| 手机成人av网站| 国产野战对白在线观看| 欧美日韩中文字幕国产精品一区二区三区| 91久久精品国产一区二区成人 | 欧美色视频一区免费| 欧美成人性av电影在线观看| 嫩草影院精品99| 久久亚洲精品不卡| 人人妻人人澡欧美一区二区| 日韩人妻高清精品专区| 国内少妇人妻偷人精品xxx网站| 国产99白浆流出| 亚洲人成网站在线播放欧美日韩| 亚洲av二区三区四区| 亚洲最大成人中文| 国产精品99久久久久久久久| 久久久久国产精品人妻aⅴ院| 亚洲国产欧美网| 亚洲午夜理论影院| 国产午夜精品论理片| 久久久久久人人人人人| 亚洲熟妇中文字幕五十中出| 中亚洲国语对白在线视频| 欧美黄色淫秽网站| 国产亚洲精品一区二区www| 亚洲一区二区三区色噜噜| www日本黄色视频网| 有码 亚洲区| 国模一区二区三区四区视频| 91麻豆av在线| 成人永久免费在线观看视频| 神马国产精品三级电影在线观看| 最近最新免费中文字幕在线| 国产成人aa在线观看| 日韩精品青青久久久久久| 亚洲精品乱码久久久v下载方式 | 最后的刺客免费高清国语| 人人妻人人澡欧美一区二区| 国产熟女xx| 亚洲男人的天堂狠狠| 久久午夜亚洲精品久久| 亚洲av免费在线观看| 国产午夜福利久久久久久| 中文字幕久久专区| 欧美黄色淫秽网站| 男女视频在线观看网站免费| 男女做爰动态图高潮gif福利片| 男女下面进入的视频免费午夜| 国产高潮美女av| 欧美日本视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利在线观看吧| 亚洲国产欧洲综合997久久,| 最近视频中文字幕2019在线8| 国产欧美日韩一区二区三| 日本成人三级电影网站| 国产精品 欧美亚洲| 俺也久久电影网| 乱人视频在线观看| 久久人人精品亚洲av| 久久精品亚洲精品国产色婷小说| 欧美一级a爱片免费观看看| 国产中年淑女户外野战色| 欧美激情在线99| 99精品欧美一区二区三区四区| 日韩成人在线观看一区二区三区| 欧美区成人在线视频| 国产亚洲精品一区二区www| 中文在线观看免费www的网站| xxxwww97欧美| 最近最新免费中文字幕在线| 亚洲av日韩精品久久久久久密| 欧美bdsm另类| 亚洲自拍偷在线| 国内精品美女久久久久久| 国产精品99久久久久久久久| 美女被艹到高潮喷水动态| 国产色爽女视频免费观看| 少妇人妻精品综合一区二区 | 老司机福利观看| 97碰自拍视频| 国产亚洲精品综合一区在线观看| 亚洲熟妇熟女久久| 老熟妇仑乱视频hdxx| 国产99白浆流出| 黄色成人免费大全| 欧美日韩瑟瑟在线播放| 日本熟妇午夜| 日本熟妇午夜| 亚洲av五月六月丁香网| 国产精品99久久99久久久不卡| 精品福利观看| 国产一级毛片七仙女欲春2| 欧美最黄视频在线播放免费| 18+在线观看网站| 伊人久久精品亚洲午夜| 色老头精品视频在线观看| 女同久久另类99精品国产91| 国内精品美女久久久久久| 伊人久久大香线蕉亚洲五| 色在线成人网| 一区福利在线观看| 天堂动漫精品| 亚洲av日韩精品久久久久久密| 久久久精品欧美日韩精品| 亚洲性夜色夜夜综合| av视频在线观看入口| 高清毛片免费观看视频网站| 亚洲国产色片| 欧美黄色淫秽网站| 欧洲精品卡2卡3卡4卡5卡区| 久久国产精品人妻蜜桃| 可以在线观看毛片的网站| 人妻丰满熟妇av一区二区三区| 日日摸夜夜添夜夜添小说| 91九色精品人成在线观看| 伊人久久精品亚洲午夜| 婷婷亚洲欧美| 九九在线视频观看精品| 日韩欧美 国产精品| 国产成人系列免费观看| 国产精品久久久久久精品电影| 国产国拍精品亚洲av在线观看 | 亚洲精品久久国产高清桃花| 免费观看的影片在线观看| 无人区码免费观看不卡| 色综合欧美亚洲国产小说| 国产aⅴ精品一区二区三区波| 亚洲欧美日韩高清在线视频| 国产熟女xx| av福利片在线观看| 91久久精品电影网| 亚洲精品日韩av片在线观看 | 中文字幕熟女人妻在线| 母亲3免费完整高清在线观看| 午夜免费激情av| 国产麻豆成人av免费视频| 特大巨黑吊av在线直播| 黄色丝袜av网址大全| 窝窝影院91人妻| 久久精品国产亚洲av涩爱 | 国产精品女同一区二区软件 | 国产伦一二天堂av在线观看| 欧美激情久久久久久爽电影| 亚洲成人免费电影在线观看| 老熟妇仑乱视频hdxx| 亚洲精品成人久久久久久| 亚洲精品色激情综合| 亚洲av日韩精品久久久久久密| 色播亚洲综合网| 亚洲成人久久爱视频| 亚洲第一欧美日韩一区二区三区| 精品一区二区三区av网在线观看| 在线播放无遮挡| 国产极品精品免费视频能看的| 少妇人妻一区二区三区视频| 九九久久精品国产亚洲av麻豆| 99在线人妻在线中文字幕| av视频在线观看入口| 国产精品美女特级片免费视频播放器| 欧美一级毛片孕妇| 一本精品99久久精品77| 好看av亚洲va欧美ⅴa在| 一个人看的www免费观看视频| 麻豆久久精品国产亚洲av| 午夜影院日韩av| 中文字幕人妻丝袜一区二区| 极品教师在线免费播放| www.www免费av| 女人被狂操c到高潮| 亚洲人成网站高清观看| 欧美日韩国产亚洲二区| 夜夜看夜夜爽夜夜摸| 少妇人妻一区二区三区视频| 国产成年人精品一区二区| 18禁美女被吸乳视频| 久久久色成人| 又爽又黄无遮挡网站| 看免费av毛片| 白带黄色成豆腐渣| 亚洲 国产 在线| 九九久久精品国产亚洲av麻豆| 午夜福利免费观看在线| 国产精品 国内视频| 可以在线观看的亚洲视频| 欧美中文综合在线视频| 日日夜夜操网爽| 亚洲精品一卡2卡三卡4卡5卡| 好男人电影高清在线观看| 99久久精品热视频| 成人特级av手机在线观看| 嫁个100分男人电影在线观看| 蜜桃亚洲精品一区二区三区| 中文字幕高清在线视频| 成年女人看的毛片在线观看| 午夜福利欧美成人| 国产一区二区三区在线臀色熟女| 给我免费播放毛片高清在线观看| 国内精品久久久久久久电影| 精品久久久久久久久久免费视频| 国产精品亚洲美女久久久| 两个人的视频大全免费| 草草在线视频免费看| 欧美日韩乱码在线| 国产淫片久久久久久久久 | 99国产精品一区二区蜜桃av| 亚洲成人久久爱视频| 亚洲天堂国产精品一区在线| 免费看光身美女| 91字幕亚洲| 亚洲内射少妇av| 亚洲五月婷婷丁香| 午夜福利成人在线免费观看| 可以在线观看毛片的网站| 日本成人三级电影网站| 国产午夜精品久久久久久一区二区三区 | 中文字幕人成人乱码亚洲影| 国产欧美日韩精品亚洲av| 最后的刺客免费高清国语| 熟女少妇亚洲综合色aaa.| 两性午夜刺激爽爽歪歪视频在线观看| 91久久精品电影网| 热99在线观看视频| 少妇裸体淫交视频免费看高清| 亚洲成人久久爱视频| 国产在线精品亚洲第一网站| 国产精品99久久99久久久不卡| 亚洲成人精品中文字幕电影| 嫩草影院精品99| 中文字幕精品亚洲无线码一区| 天堂影院成人在线观看| 黄色日韩在线| av欧美777| 69人妻影院| 桃色一区二区三区在线观看| 精品久久久久久久人妻蜜臀av| 国产一区二区三区视频了| 亚洲激情在线av| 日韩欧美精品免费久久 | 久久久国产成人精品二区| 一级毛片高清免费大全| 亚洲一区二区三区色噜噜| 99久久精品热视频| 男女视频在线观看网站免费| 亚洲欧美精品综合久久99| 18+在线观看网站| 国内精品美女久久久久久| 免费电影在线观看免费观看| 国模一区二区三区四区视频| 国产精品久久电影中文字幕| 内射极品少妇av片p| 亚洲片人在线观看| а√天堂www在线а√下载| 夜夜夜夜夜久久久久| 美女被艹到高潮喷水动态| 亚洲人成网站在线播放欧美日韩| 久久这里只有精品中国| 亚洲七黄色美女视频| 91麻豆av在线| 国产精品乱码一区二三区的特点| 99久国产av精品| netflix在线观看网站| 午夜福利在线在线| 舔av片在线| 成年免费大片在线观看| 又黄又粗又硬又大视频| 3wmmmm亚洲av在线观看| 夜夜爽天天搞| 欧美极品一区二区三区四区| 欧美精品啪啪一区二区三区| 在线视频色国产色| 99国产精品一区二区蜜桃av| 亚洲av不卡在线观看| 日韩国内少妇激情av| 日韩成人在线观看一区二区三区| 免费观看人在逋| 国产成人欧美在线观看| 国产探花极品一区二区| 午夜福利18| 国内毛片毛片毛片毛片毛片| 狂野欧美激情性xxxx| 欧美激情久久久久久爽电影| 五月伊人婷婷丁香| 91久久精品国产一区二区成人 | 国产激情偷乱视频一区二区| 国产av在哪里看| 国产伦人伦偷精品视频| 亚洲av电影不卡..在线观看| 手机成人av网站| 高清日韩中文字幕在线| 久久久久九九精品影院| 国产又黄又爽又无遮挡在线| 三级国产精品欧美在线观看| 国产精品久久久久久精品电影| 久久精品91无色码中文字幕| 在线播放国产精品三级| 国语自产精品视频在线第100页| 3wmmmm亚洲av在线观看| 国产精品 国内视频| www.熟女人妻精品国产| 国产精品国产高清国产av| 久久精品国产亚洲av香蕉五月| a在线观看视频网站| 亚洲天堂国产精品一区在线| av天堂中文字幕网| 亚洲国产欧洲综合997久久,| 午夜精品久久久久久毛片777| netflix在线观看网站| 97人妻精品一区二区三区麻豆| 嫩草影院入口| 真人做人爱边吃奶动态| 国产主播在线观看一区二区| 一本久久中文字幕| 啦啦啦观看免费观看视频高清| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 精品久久久久久久久久久久久| 小说图片视频综合网站| 99热只有精品国产| 麻豆成人午夜福利视频| 在线国产一区二区在线| 午夜免费成人在线视频| 18禁美女被吸乳视频| svipshipincom国产片| 夜夜看夜夜爽夜夜摸| 九九热线精品视视频播放| 精品午夜福利视频在线观看一区| 小说图片视频综合网站| 哪里可以看免费的av片| 免费看日本二区| 欧美最新免费一区二区三区 | 国产精品98久久久久久宅男小说| 岛国视频午夜一区免费看| x7x7x7水蜜桃| 网址你懂的国产日韩在线| 亚洲熟妇中文字幕五十中出| 美女高潮的动态| 在线观看免费视频日本深夜| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成a人片在线一区二区| 国产成+人综合+亚洲专区| 亚洲美女视频黄频| 久久久国产成人免费| 国产伦精品一区二区三区视频9 | 亚洲成a人片在线一区二区| 国产精品99久久99久久久不卡| 国产高清视频在线观看网站| 少妇人妻精品综合一区二区 | 99精品欧美一区二区三区四区| 少妇人妻精品综合一区二区 | 在线免费观看的www视频| 成人av在线播放网站| 熟妇人妻久久中文字幕3abv| 国产精品一及| 99热6这里只有精品| 亚洲色图av天堂| 欧美黄色片欧美黄色片| 亚洲人与动物交配视频| 亚洲国产精品合色在线| 波多野结衣高清作品| e午夜精品久久久久久久| 网址你懂的国产日韩在线| 国产成人福利小说| 欧美一级毛片孕妇| 日本黄色视频三级网站网址| 欧美一区二区亚洲| 黄色丝袜av网址大全| 成人国产一区最新在线观看| 亚洲欧美日韩卡通动漫| 日本一本二区三区精品| 波多野结衣高清无吗| 欧美xxxx黑人xx丫x性爽| 在线国产一区二区在线| 搡女人真爽免费视频火全软件 | 亚洲精品国产精品久久久不卡| 少妇的丰满在线观看| 中文字幕人妻熟人妻熟丝袜美 | 日本 欧美在线| 国产伦精品一区二区三区四那| 午夜福利18| 高清毛片免费观看视频网站| 日韩欧美在线乱码| 日本 av在线| 欧美日韩福利视频一区二区| 天天一区二区日本电影三级| 一二三四社区在线视频社区8| 亚洲精品456在线播放app | 亚洲aⅴ乱码一区二区在线播放| 中文字幕人妻熟人妻熟丝袜美 | 国产亚洲精品久久久com| 欧美高清成人免费视频www| 日本 av在线| 熟女人妻精品中文字幕| 欧美绝顶高潮抽搐喷水| 欧美午夜高清在线| 久久久国产成人免费| 亚洲一区二区三区不卡视频| 国产精品,欧美在线| 亚洲一区高清亚洲精品| 日本黄大片高清| 欧美成人a在线观看| 日韩欧美在线乱码| 国产精品久久久久久精品电影| av片东京热男人的天堂| 午夜福利在线观看免费完整高清在 | 很黄的视频免费| 亚洲 欧美 日韩 在线 免费| 欧美zozozo另类| 色综合婷婷激情| 亚洲av免费在线观看| 免费看日本二区| 亚洲av成人不卡在线观看播放网| 国模一区二区三区四区视频| 欧美乱码精品一区二区三区| 国产欧美日韩精品亚洲av| 亚洲av电影不卡..在线观看| 亚洲精品亚洲一区二区| 久久天躁狠狠躁夜夜2o2o| 国产av麻豆久久久久久久| 人妻夜夜爽99麻豆av| 嫩草影视91久久| 国产三级在线视频| 欧美bdsm另类| av中文乱码字幕在线| 色老头精品视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 搡女人真爽免费视频火全软件 | 久久香蕉国产精品| 精品久久久久久成人av| 嫩草影视91久久| 久久久色成人| 小说图片视频综合网站| 舔av片在线| 国内精品久久久久精免费| 久久99热这里只有精品18| 欧美丝袜亚洲另类 | 欧美色视频一区免费| 麻豆成人av在线观看| 叶爱在线成人免费视频播放| 美女被艹到高潮喷水动态| 久久99热这里只有精品18| 最新在线观看一区二区三区| 亚洲一区二区三区色噜噜| 日日干狠狠操夜夜爽| 天堂av国产一区二区熟女人妻| 18禁国产床啪视频网站| 色综合欧美亚洲国产小说| av中文乱码字幕在线| 亚洲成av人片免费观看| 欧美日韩瑟瑟在线播放| 听说在线观看完整版免费高清| 亚洲欧美日韩高清在线视频| 久久这里只有精品中国| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看| 欧美成人一区二区免费高清观看| 成人永久免费在线观看视频| 国产成人福利小说| 亚洲第一电影网av| 午夜福利18| 国产成+人综合+亚洲专区| 99在线人妻在线中文字幕| 国产在视频线在精品| 亚洲精品成人久久久久久| 操出白浆在线播放| 亚洲一区二区三区色噜噜| 嫩草影院入口| 香蕉久久夜色| 国产精品自产拍在线观看55亚洲| 怎么达到女性高潮| 亚洲av美国av| 精品日产1卡2卡| 国产三级在线视频| 国产成人aa在线观看| 黄色日韩在线| 露出奶头的视频| 日日摸夜夜添夜夜添小说| 国产精品久久久久久精品电影| 狂野欧美激情性xxxx| 亚洲美女黄片视频| 免费在线观看亚洲国产| 身体一侧抽搐| 亚洲精品久久国产高清桃花| 日韩精品中文字幕看吧| 国产伦人伦偷精品视频| 国产亚洲精品综合一区在线观看| 亚洲国产欧洲综合997久久,| 级片在线观看| 热99在线观看视频| 99久久综合精品五月天人人| 亚洲自拍偷在线| 色视频www国产| 男人舔奶头视频| 别揉我奶头~嗯~啊~动态视频| 免费在线观看影片大全网站| 成年人黄色毛片网站| 日本成人三级电影网站| 中亚洲国语对白在线视频| 亚洲va日本ⅴa欧美va伊人久久| 18禁在线播放成人免费| 9191精品国产免费久久| a级毛片a级免费在线| 欧美高清成人免费视频www| 国产一区二区在线av高清观看| 高潮久久久久久久久久久不卡| 久久精品影院6| 免费在线观看日本一区| 长腿黑丝高跟| 观看免费一级毛片| 精品乱码久久久久久99久播| 男女那种视频在线观看| 久久亚洲精品不卡| 精品人妻一区二区三区麻豆 | 岛国在线观看网站| 国产伦人伦偷精品视频| 757午夜福利合集在线观看| 日日干狠狠操夜夜爽| 国产精品99久久99久久久不卡| 欧美黑人巨大hd| 婷婷亚洲欧美| 精品久久久久久久末码| 网址你懂的国产日韩在线| 欧美成人性av电影在线观看| 一本综合久久免费| 国产精品女同一区二区软件 | 亚洲av熟女| 白带黄色成豆腐渣| 午夜福利在线观看吧| 国产视频一区二区在线看| www.www免费av| 日本黄色视频三级网站网址| a在线观看视频网站| 欧美日韩综合久久久久久 | 国产在线精品亚洲第一网站| 岛国视频午夜一区免费看| 淫妇啪啪啪对白视频| 国产国拍精品亚洲av在线观看 | 成人国产一区最新在线观看| 国内精品久久久久久久电影| 色在线成人网| 成年女人毛片免费观看观看9| www日本黄色视频网| 国产v大片淫在线免费观看| 热99在线观看视频| 午夜久久久久精精品| 亚洲av成人av| 色综合婷婷激情| 舔av片在线| 无限看片的www在线观看| 欧美激情在线99| 一级黄色大片毛片| 黄片大片在线免费观看| ponron亚洲| 一级黄片播放器| 色在线成人网| 欧美日韩亚洲国产一区二区在线观看| 黄色女人牲交| 国产成人av激情在线播放| 亚洲国产欧美人成| 香蕉久久夜色| 夜夜爽天天搞| 亚洲成人免费电影在线观看| 亚洲成av人片免费观看| 国产爱豆传媒在线观看| e午夜精品久久久久久久| 国产不卡一卡二|