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

    烏梁素海冰封期浮游植物群落特征與環(huán)境因子CCA分析

    2018-05-03 06:04:57何婷婷勾芒芒
    關(guān)鍵詞:冰體烏梁素海甲藻

    李 興,何婷婷,勾芒芒

    (1.內(nèi)蒙古師范大學(xué)內(nèi)蒙古節(jié)水農(nóng)業(yè)工程研究中心,呼和浩特 010022;2.內(nèi)蒙古機(jī)電職業(yè)技術(shù)學(xué)院水利與土木建筑工程系,呼和浩特 010070)

    寒旱區(qū)湖泊具有獨(dú)特區(qū)域和氣候特征,冰封期長、冰層厚、流量小、污染重,為國內(nèi)外研究熱點(diǎn)。湖泊、水庫、池塘等水體在冰封時期水質(zhì)變化較大,冰體生長過程中,對水質(zhì)指標(biāo)有排斥作用,水體內(nèi)物質(zhì)濃度不同程度增加,影響水產(chǎn)健康養(yǎng)殖[1]。因此,揭示冰封期湖庫水體和冰體水質(zhì)與浮游植物關(guān)系及變化過程具有重要意義。目前,關(guān)于湖庫、池塘等水域浮游植物與環(huán)境因子響應(yīng)關(guān)系方面研究較多,主要集中在溫暖濕潤地區(qū)水庫、海灣[2-3]、湖泊[4-6]和河流[7]等區(qū)域,研究區(qū)域浮游植物與環(huán)境因子響應(yīng)關(guān)系基本明確。針對寒冷干旱地區(qū)冰封期湖泊水體和冰體中環(huán)境因子與浮游植物間響應(yīng)關(guān)系研究鮮見報道。冰封期湖泊水體、冰體內(nèi)環(huán)境因子與浮游植物響應(yīng)關(guān)系與非冰封期不同,冰封期冰體阻礙水體與大氣復(fù)氧,冰封期積雪覆蓋影響太陽凈輻射等,使冰封期水體與冰體內(nèi)浮游植物群落結(jié)構(gòu)特征有一定典型性和特殊性。內(nèi)蒙古烏梁素海典型氣候特征和地域環(huán)境決定該區(qū)域湖泊水體中污染物遷移轉(zhuǎn)化過程及污染分布特性與其他區(qū)域湖泊不同,冰封期冰-水介質(zhì)中藻類及污染物濃度及分布特性與湖泊水華現(xiàn)象關(guān)系密切。因此,本文以寒旱區(qū)湖泊烏梁素海為研究對象,揭示寒旱區(qū)湖泊冰封期藻類及環(huán)境因子在冰-水介質(zhì)中分布特性及響應(yīng)關(guān)系,對治理和恢復(fù)寒旱區(qū)湖泊生態(tài)系統(tǒng)、發(fā)展水產(chǎn)健康養(yǎng)殖具有重要意義。

    烏梁素海位于內(nèi)蒙古自治區(qū)烏拉特前旗,地理坐標(biāo)為 40°36′~41°03′N,108°43′~108°57′E,湖泊南北長、東西窄,湖岸線長為130 km。湖泊蓄水量2.5~3億m3,水深0.5~3.0 m,平均水深1.5 m,水域面積273.32 km2,為我國第八大淡水湖,是黃河中上游重要保水、蓄水和調(diào)水場所[8]。湖泊流域內(nèi)降雨量少蒸發(fā)量大,每年11月初開始結(jié)冰,翌年3月份解凍,冰封期長,在蒙新高原湖區(qū),烏梁素海氣候區(qū)域特征具有代表性。烏梁素海承接河套地區(qū)大量農(nóng)田退水,有凈化水質(zhì)、調(diào)節(jié)水量、防止河套灌區(qū)土壤鹽堿化、減緩水體富營養(yǎng)化進(jìn)程等多種生態(tài)功能[9],對河套地區(qū)乃至黃河流域生態(tài)環(huán)境有重要調(diào)節(jié)作用,在我國濕地生態(tài)系統(tǒng)中居重要地位。

    1 材料與方法

    1.1 采樣點(diǎn)布置

    按照我國水環(huán)境及濕地生態(tài)系統(tǒng)調(diào)查規(guī)范,根據(jù)烏梁素海污染源、蘆葦、水動力及環(huán)境特征,將烏梁素海在空間上以2 km×2 km正方形網(wǎng)格剖分,在網(wǎng)格交點(diǎn)以梅花形布置取樣監(jiān)測點(diǎn)10處(見圖1)[10]。每次取樣點(diǎn)位置GPS定位。取樣時間為2016年1月15日烏梁素海冰層穩(wěn)定期。監(jiān)測點(diǎn)包括進(jìn)水區(qū)域(J11)、I12、出水區(qū)域河口(V3)、北部明水區(qū)域(L15)、蘆葦區(qū)域(N13)、旅游區(qū)域(Q10)、西大灘(L11)、大卜洞(O10)、Q8、二點(diǎn)。

    1.2 樣品采集與處理

    2016年1月15日在烏梁素海冰層穩(wěn)定期使用冰鉆采集器采樣,破冰后,用自制壓力采水器采集1 000 mL冰下水樣待測,同時將冰樣按厚度每10 cm分層,因不同采樣點(diǎn)冰厚不同,分3~5層,置于2 000 mL塑料瓶內(nèi),黑暗低溫條件下自然融化,立即測定冰體中總氮(過硫酸鉀氧化紫外分光光度法)、溶解性總磷(鉬銻抗分光光度法)、總磷(鉬銻抗分光光度法)、硝態(tài)氮(酚二磺酸光度法)、亞硝態(tài)氮(重氮偶合分光光度法)、氨態(tài)氮(納氏試劑分光光度法)和化學(xué)需氧量(重鉻酸鉀法)水體指標(biāo),以上室內(nèi)水質(zhì)指標(biāo)測定方法參照文獻(xiàn)[11]。野外監(jiān)測指標(biāo)主要包括冰上氣象站監(jiān)測指標(biāo)(環(huán)境溫度、風(fēng)速、風(fēng)向等)和冰下水體環(huán)境指標(biāo)(水溫、pH、溶解氧、氧化還原電位),冰下水體野外現(xiàn)場測試環(huán)境指標(biāo)用ODEON Multy 8330便攜式多參數(shù)水質(zhì)分析儀測定(購自PONSEL公司)。

    圖1 烏梁素海采樣點(diǎn)分布Fig.1 Sampling point in the Wuliangsuhai Lake

    浮游植物定量樣品用采水器在現(xiàn)場取1 000 mL水面下0.5 m處水樣,加入15 mL魯哥試液固定;于水平試驗(yàn)臺上靜置沉淀30 h,虹吸法去掉上清液,余下20~25 mL沉淀物轉(zhuǎn)入30 mL定量瓶中,定容至30 mL。將室溫條件下融化冰塊作為冰體中浮游植物定量樣品,取1 000 mL,處理步驟同上。鏡檢方法、浮游植物種類鑒定等參照文獻(xiàn)[12]。

    藻類密度參照文獻(xiàn)[13]中藻類計數(shù)計算公式:

    式中,N為藻類密度(104cell·mL-1);A為浮游生物計數(shù)框面積(mm2);Ac為計數(shù)面積,即視野面積×視野數(shù)(mm2);Vw為l L水樣經(jīng)沉淀濃縮后樣品體積(mL);V為計數(shù)框體積(mL);n為計數(shù)所得藻類數(shù)量。

    以浮游植物優(yōu)勢度Y分析烏梁素海冰封期浮游植物群落結(jié)構(gòu)特征,將優(yōu)勢度Y>0.02藻類定為優(yōu)勢種。各參數(shù)計算公式如下:

    式中,ni為第i種浮游植物數(shù)量;N為浮游植物總量;fi為第i個浮游植物在所采樣品中頻率。浮游藻類分類鑒定依據(jù)《中國淡水藻類-系統(tǒng):生態(tài)及分類》。其中浮游植物優(yōu)勢種藻類鑒定到種,其他種類鑒定到屬和種。

    用Canoco for Windows 4.5作CCA分析。CCA典范對應(yīng)分析時,浮游植物物種經(jīng)篩選,至少在一個監(jiān)測點(diǎn)數(shù)量占該監(jiān)測點(diǎn)總量>12.5%進(jìn)入排序[14]。對應(yīng)分析前,先將各月主要浮游植物物種數(shù)據(jù)作去趨勢間接梯度分析(Detrended correspondence analysis,DCA),每個排序軸梯度長度值(Lengths of gradient)均<3。盡管最長梯度<3,但線性模型可視為單峰模型特例。因此,選擇單峰模型作CCA分析。分析前,浮游植物物種數(shù)據(jù)和環(huán)境因子數(shù)據(jù)(除pH外)均作lg(x+1)形式轉(zhuǎn)換,分析結(jié)果用物種—環(huán)境因子關(guān)系雙序圖表示[15-16]。

    2 結(jié)果與分析

    2.1 物種組成

    在烏梁素海冰下水體共檢出浮游植物7門32屬37種,其中,硅藻門10屬12種,占32.43%;綠藻門8屬9種,占24.32%;裸藻門4屬6種,占16.22%;甲藻門4屬4種,占10.81%;藍(lán)藻門3屬3種,占8.11%;隱藻門2屬2種,占5.41%;褐藻門1屬1種,占2.70%。冰層中共檢出浮游植物7門34屬60種,其中,硅藻門16屬34種,占56.67%;裸藻門3屬8種,占13.33%;綠藻門5屬6種,占10.00%;藍(lán)藻門3屬5種,占8.33%;甲藻門4屬4種,占6.67%;隱藻門2屬2種,占3.33%;金藻門1屬1種,占1.67%。與烏梁素海夏季物種組成相比[10],硅藻門、綠藻門種類數(shù)仍占主導(dǎo)地位,但綠藻門所占比例下降,硅藻門和裸藻門上升。在冰層物種組成中,硅藻門類為優(yōu)勢種,達(dá)56.67%。

    2.2 優(yōu)勢種

    烏梁素海浮游植物冰封期各采樣點(diǎn)優(yōu)勢種及優(yōu)勢度見表1。

    表1 烏梁素海冰封期浮游植物優(yōu)勢種及優(yōu)勢度Table 1 Dominantspecies and dominance ofphytoplankton inWuliangsuhaiLakeduring

    水體共發(fā)現(xiàn)優(yōu)勢種10種,分別由硅藻門大柄鏈藻(Podosira maxima)、暗色圓篩藻(Coscinodiscus obscurus)、極小冠盤藻(Stephanodiscusminutulus)和星形冠盤藻(Stephanodiscusastraea),綠藻門水溪綠球藻(Chlorococcuminfusionum),甲藻門利馬原甲藻(Prorocentrum lima)、薄甲藻(Glenodiniumpulviscu?lus)、不顯著多甲藻(Peridiniuminconspicuum),藍(lán)藻門細(xì)長聚球藻(Synechococcuselongatus),隱藻門卵形紅胞藻(Rhodomonasovalis)組成,不同采樣點(diǎn)優(yōu)勢種組成存在差距。冰體內(nèi)浮游植物優(yōu)勢種共14種,分別為暗色圓篩藻(Coscinodiscus obscurus),極小冠盤藻(Stephanodiscusminutulus),大柄鏈藻(Podosira maxima),薄甲藻(Glenodiniumpulviscu?lus),不顯著多甲藻(Peridiniuminconspicuum),卵形紅胞藻(Rhodomonasovalis),利馬原甲藻(Proro?centrum lima),肘狀針桿藻(Synedra ulna vardanica),平片針桿藻(Synedratabulata),線形舟形藻(Navic?ulagraciloides),星形柄鏈藻(Podosirastelliger),塔瑪亞歷山大藻(Alexandriumtamarense),萎軟幾內(nèi)亞藻(Guinardiaflaccida),粗壯細(xì)鞘絲藻(Lyngby?aattenuata Fritsch)。

    各監(jiān)測點(diǎn)優(yōu)勢種差異性顯著(見表1),J11點(diǎn)主要以薄甲藻、不顯著多甲藻、利馬原甲藻、卵形紅細(xì)胞組成,其中薄甲藻占相對優(yōu)勢,優(yōu)勢度為0.267;I12點(diǎn)主要由硅藻門的暗色圓篩藻、極小冠盤藻組成,暗色圓篩藻占相對優(yōu)勢,優(yōu)勢度為0.229;L15點(diǎn)優(yōu)勢種最多,為8種,其中大柄鏈藻占絕對優(yōu)勢,優(yōu)勢度為0.243;O10、Q10、L11、V3點(diǎn)均以硅藻門極小冠盤藻占絕對優(yōu)勢,優(yōu)勢度分別為0.314、0.162、0.195、0.295。N13、Q8 和二點(diǎn)均以甲藻門的利馬原甲藻占絕對優(yōu)勢,分別為0.338、0.376、0.354;二點(diǎn)和Q10點(diǎn)優(yōu)勢種最少,均為3種。但不顯著多甲藻,薄甲藻,大柄鏈藻,極小冠盤藻等在各監(jiān)測點(diǎn)出現(xiàn)頻率較大,不顯著多甲藻在冰體出現(xiàn)頻率大于90%,優(yōu)勢度大于0.1,最大優(yōu)勢度為0.443,均與硅藻、甲藻生長最適溫度低有關(guān)[17]。

    2.3 冰封期冰體中浮游植物密度垂直變化

    調(diào)查期間,各冰層中浮游植物群落結(jié)構(gòu)特征差異顯著。由于各采樣點(diǎn)冰厚不同,無法對各采樣點(diǎn)藻類密度定量分析,僅根據(jù)冰厚百分比描述藻類密度變化趨勢(假定浮游植物群落在冰層中均勻線性分布)。

    監(jiān)測點(diǎn)J11、Q10從上層冰到底層冰浮游植物藻類密度呈先增后減趨勢;I12和O10監(jiān)測點(diǎn)呈先緩減后增趨勢;N13和二點(diǎn)兩個監(jiān)測點(diǎn)呈先減后保持穩(wěn)定趨勢;在V3監(jiān)測點(diǎn),藻類密度則表現(xiàn)為先減后增趨勢;Q8、L11、L15等監(jiān)測點(diǎn)較復(fù)雜,呈先減后升再緩慢減小趨勢(見圖2)。烏梁素海冰封期冰體中浮游植物藻類在豎直梯度上呈復(fù)雜分布規(guī)律,與冰體環(huán)境因子復(fù)雜性有關(guān)。從浮游植物門類分析,各監(jiān)測點(diǎn)主要浮游植物門類均為甲藻和硅藻,監(jiān)測點(diǎn)浮游植物群落藻類密度組成中占主體地位。垂直梯度分布規(guī)律和各監(jiān)測點(diǎn)藻類細(xì)胞總密度分布規(guī)律一致。

    2.4 冰封期浮游植物多樣性分析

    浮游植物Shannon-Weaver多樣性指數(shù)H可描述浮游植物群落構(gòu)成情況及水體污染對生物影響。某一群落種類多樣性越高,結(jié)構(gòu)越復(fù)雜,穩(wěn)定性越好[18-19]。Margalef物種豐富度指數(shù)D可反映植物群落與環(huán)境間關(guān)系,理論上浮游植物種類越多,個體數(shù)量分布越均勻,豐富度指數(shù)越大,指示環(huán)境越穩(wěn)定。均勻度指數(shù)J為實(shí)際多樣性指數(shù)與理論最大多樣性指數(shù)比值,反映物種個體數(shù)目分配均勻程度。由圖3可知,總體上多樣性指數(shù)H、均勻度指數(shù)J各監(jiān)測點(diǎn)差異較小。其中位于N13和L11兩監(jiān)測點(diǎn)多樣性指數(shù)均高于其他點(diǎn)位,此處監(jiān)測點(diǎn)位于蘆葦密集區(qū)域,水生態(tài)系統(tǒng)物種組成相對復(fù)雜,物種更豐富,穩(wěn)定性、均勻度相對較高。位于烏梁素海出水口V3、Q8和進(jìn)水口J11多樣性指數(shù)明顯偏低,說明浮游植物種類組成多樣性較低,監(jiān)測點(diǎn)附近水域內(nèi)浮游植物群落結(jié)構(gòu)相對不穩(wěn)定。其他各點(diǎn)多樣性指數(shù)接近。Margalef物種豐富度指數(shù)D對物種數(shù)依賴性較強(qiáng),可充分反映浮游植物在水體中分布情況和群落特征[20]。J11、V3和Q8由于分別位于入水口、人工漁場及出水口位置,因污染物外源性輸入及養(yǎng)殖業(yè)廢水污染,水質(zhì)惡化,浮游植物物種數(shù)減少,豐富度指數(shù)D顯著低于其他監(jiān)測點(diǎn)。

    2.5 冰封期浮游植物群落與環(huán)境因子CCA分析

    CCA分析中冰下水體和冰體中浮游植物種類代碼見表2~3。

    烏梁素海冰封期水體篩選出20種浮游植物及9個環(huán)境因子作CCA分析,結(jié)果見圖4~5。

    圖2 各采樣點(diǎn)冰層浮游植物密度Fig.2 A algae density of phytoplankton at each sampling point

    圖3 烏梁素海冰封期浮游植物多樣性指數(shù)Fig.3 Phytoplankton diversity index in the Wuliangsuhai Lake

    表2 CCA分析冰下水體浮游植物種類代碼Table 2 Codes of phytoplankton species for CCA analysis

    表3 CCA分析冰體中浮游植物種類代碼Table 3 Codes of phytoplankton species for CCA analysis

    圖4 冰下水體采樣點(diǎn)與環(huán)境因子CCA排序Fig.4 CCA biplot of sample stations and environmental variables in the Wuliangsuhai Lake

    圖5 冰下水體浮游植物群落與環(huán)境因子CCA排序Fig.5 CCA biplot of phytoplankton species and environmental variables in the Wuliangsuhai Lake

    兩個物種排序軸和環(huán)境排序軸相關(guān)系數(shù)均為0,表明排序結(jié)果可靠[21],軸1和軸2反映48.9%物種信息量。由圖4~5可知,排序軸1呈正相關(guān)因子主要有ORP,相關(guān)系數(shù)為0.3653,負(fù)相關(guān)因子主要為透明度(WD),相關(guān)系數(shù)為-0.4716,排序軸2呈正相關(guān)因子主要是TP、DO、TN和TDP,相關(guān)系數(shù)分別為 0.7484、0.6912、0.5483 和 0.5017,TP、DO、TN、TDP對烏梁素海浮游植物群落結(jié)構(gòu)影響較大。烏梁素海浮游植物在CCA第1、第2排序軸上分化較好,綠藻種類主要集中在第3、4象限,與各環(huán)境因子距離較遠(yuǎn),說明環(huán)境因子最適值較低;裸藻種類主要集中在第1、2象限,TP、TDP、DO、TN和對其均有影響,主要在TP因子較大時有最適值;甲藻種類主要集中在第2象限,與TDP、TP、DO和WD正相關(guān),硅藻種類4個象限均有分布,說明硅藻類對烏梁素海冰封期水質(zhì)有較強(qiáng)適應(yīng)性,其中藍(lán)藻門細(xì)長聚球藻(S17)與各環(huán)境因子距離較遠(yuǎn),說明各環(huán)境因子最適值較低,而具緣微囊藻S19截然相反,與TDP、DO、TP、TN、pH等環(huán)境指標(biāo)正相關(guān),在TP因子較大時有最適值。各采樣點(diǎn)浮游植物群落組成差異較大,J11樣點(diǎn)主要由甲藻門的薄甲藻、隱藻門的卵形細(xì)胞組成,受TP、DO、TDP、TN影響不大;I12樣點(diǎn)主要由硅藻門的暗色圓篩藻、極小冠盤藻組成,受TP、TN和氨氮等影響較大,受WD影響較小;二點(diǎn)樣點(diǎn)以甲藻門的利馬原甲藻占絕對優(yōu)勢,該采點(diǎn)受WD影響較大;O10采樣點(diǎn)以硅藻門極小冠盤藻和甲藻門利馬原甲藻占絕對優(yōu)勢,L11采樣點(diǎn)以硅藻門極小冠盤藻、星形冠盤藻占優(yōu)勢,V3主要以硅藻門極小冠盤藻為主,3個采樣點(diǎn)受、、PH影響較大??梢?,、、PH含量變化影響L11、O10、V3采樣點(diǎn)物種組成和數(shù)量。

    對烏梁素海冰體中經(jīng)篩選23種浮游植物及8個環(huán)境因子作統(tǒng)計分析,CCA分析結(jié)果見圖6~7。結(jié)果顯示兩個物種排序軸和兩個環(huán)境排序軸相關(guān)系數(shù)均為0,表明排序結(jié)果有可靠性[22]。由圖6~7可知,、TN、、COD對烏梁素海冰封期冰層浮游植物群落結(jié)構(gòu)影響較大,其中、與軸1正相關(guān),相關(guān)系數(shù)為0.4351、0.7159,COD與軸1負(fù)相關(guān),相關(guān)系數(shù)為-0.4143。TP、TN與軸2正相關(guān),相關(guān)系數(shù)分別為0.7484、0.5483。

    烏梁素海冰封期冰層浮游植物在CCA第1、2排序軸上分化較好,綠藻種類較少,主要集中在第3象限,與各環(huán)境因子距離較遠(yuǎn),說明環(huán)境因子最適值較低;裸藻、隱藻種類較少,主要集中在第1~2象限,與各環(huán)境指標(biāo)均為負(fù)相關(guān)。甲藻種類分布較廣,其中塔瑪亞歷山大藻在環(huán)境變量TDP較大時有最適值,利馬原甲藻在環(huán)境變量較大時有最適值;硅藻種類在4個象限均有分布,說明硅藻門對冰封期烏梁素海冰體環(huán)境有較強(qiáng)適應(yīng)性,其中H30舟形藻屬(Navicula)、H37平片針桿藻(Synedra tabulata)、H55卵圓雙眉藻(Amphora ovalis)在TN值較大時有最適值。

    圖6 冰體采樣點(diǎn)與環(huán)境因子CCA排序Fig.6 CCA biplot of sample stations and environmental variables in the Wuliangsuhai Lake ice

    圖7 冰體浮游植物群落與環(huán)境因子CCA排序Fig.7 CCA biplot of phytoplankton species and environmental variables in the Wuliangsuhai Lake ice

    3 討論與結(jié)論

    3.1 冰封期冰體及冰下浮游植物群落結(jié)構(gòu)特點(diǎn)

    烏梁素海冰封期水體浮游植物以硅藻門種類最多,占總種類數(shù)32.43%;其次為綠藻和裸藻,分別占24.32%和16.22%。而在冰層中硅藻門占56.67%;裸藻門、綠藻門次之,分別占13.33%、10.00%。與烏梁素海夏季藻類物種組成相比,硅藻門、綠藻門種類數(shù)仍占主導(dǎo)地位,但綠藻門所占比例下降,下降24.68%。硅藻門和裸藻門物種比例有所上升,分別上升9.43%、8.22%。在冰層物種組成中變化幅度顯著,占主導(dǎo)地位的硅藻門變化幅度達(dá)56.67%。冬季低溫淘汰藍(lán)藻門和部分綠藻門種類,適宜低溫硅藻門種類迅速生長,冰下水體環(huán)境為硅藻提供低溫和低光照良好條件,硅藻在極低溫度條件下仍可進(jìn)行光合作用,繁殖速度超過其他藻類。烏梁素海冰下及冰體內(nèi)硅藻數(shù)量占絕對優(yōu)勢,與極地冰區(qū)[23]、河流冰下水體[24]硅藻數(shù)量占優(yōu)一致。

    從藻類豐度數(shù)據(jù)分析結(jié)果可知,硅藻和甲藻為烏梁素海冰封期優(yōu)勢種。其中甲藻門浮游植物物種種類不多,但數(shù)量多、分布廣,以利馬原甲藻和不顯著多甲藻為常見種類。甲藻具有由纖維素組成的細(xì)胞壁,分上甲和下甲,繁殖以細(xì)胞分裂為主。甲藻在烏梁素海成為優(yōu)勢種與海水越冬池冰下浮游植物優(yōu)勢種研究結(jié)果一致[25],但與南極冰區(qū)和河流冰下水體中藻類特征不同[23-24],可能與甲藻不適合高緯度地區(qū)氣候條件或河流冰下水體水動力條件有關(guān)。

    因此,烏梁素海冰下水體中的浮游植物密度中硅藻和甲藻普遍較大,烏梁素海冰封期水體和冰體內(nèi)浮游植物優(yōu)勢種組成呈硅藻——甲藻型,而烏梁素海春夏季浮游植物優(yōu)勢種組成呈硅藻——綠藻型[26]。表明非冰封期和冰封期水體浮游植物群落結(jié)構(gòu)特征不同。烏梁素海冰封期藻類數(shù)量和種類均少于非冰封期,主要因?yàn)槎舅疁?、光照等自然條件不利于藻類生長,符合藻類生物特性。但烏梁素海冰中藻類屬、種量高于冰下水體,冰體中硅藻種類數(shù)量顯著高于冰下水體,其他藻類則變化微小。因此,硅藻在低溫條件下新陳代謝強(qiáng)度明顯高于其他浮游植物,可通過減少細(xì)胞中水分和增加細(xì)胞中糖類、脂肪等物質(zhì)降低冰點(diǎn),增加抗寒能力,在低溫條件下具有較寬生態(tài)位,適宜低溫生存條件[27]。

    由于各采樣點(diǎn)所處位置及環(huán)境因子差異,冰體與水體中藻類多樣性指標(biāo)并無顯著性差異。各采樣點(diǎn)Margalef豐富度指數(shù)(D)、Shannon-Wiener多樣性指數(shù)(H)、Pielou均勻度指數(shù)(J)變化趨勢一致。采樣點(diǎn)多樣性指標(biāo)相關(guān)性分析表明,多樣性指數(shù)(H)受均勻度指數(shù)(J)和豐富度指數(shù)(D)變化影響,其中物種豐富度指數(shù)(D)(R=0.713,P=0.000)影響較顯著。從多樣性指數(shù)分析結(jié)果看,烏梁素海冰封期冰下水體總體浮游植物多樣性、穩(wěn)定性及均勻度較好[19]。

    3.2 冰封期影響冰體和冰下浮游植物環(huán)境因子

    在水生態(tài)系統(tǒng)中,環(huán)境對浮游植物特征變化有主導(dǎo)作用,而浮游植物分布及種類對環(huán)境也有指示作用,不同時期環(huán)境條件改變可直接或間接影響浮游植物特征。從浮游植物物種與環(huán)境因子關(guān)系CAA分析來看,TP為影響烏梁素海冰下水體浮游植物群落分布限制性環(huán)境因子。與明水期松花江哈爾濱段浮游植物分布與環(huán)境因子中總磷具有顯著相關(guān)性的研究結(jié)果一致[28],而與冰下松花江哈爾濱段與浮游植物密切相關(guān)的環(huán)境因子為溶解氧、硝酸鹽氮不同[24]。由于烏梁素海冬季受化學(xué)需氧量、總磷污染加重,湖泊內(nèi)部水體交換遲緩,與北方河流冬季特點(diǎn)不同,水動力條件對于增加水體中溶解氧含量,提高水體自凈能力十分重要,而溶解氧是決定水體環(huán)境重要因子。冰層中NH4+、TP、TN等環(huán)境因子變化對浮游植物群落在冰層中分布影響較大。原因是冰生長過程中層理結(jié)構(gòu)發(fā)生變化,對各類污染物有不同排斥效應(yīng)。水結(jié)冰過程中,污染物會從冰體中不同程度析出后進(jìn)入水體;通過檢測冰體和水體中浮游植物發(fā)現(xiàn),水結(jié)冰過程對浮游植物存在排斥效應(yīng),因水結(jié)冰過程中對浮游植物和不同環(huán)境因子排斥程度不同,影響冰體和水體浮游植物分布環(huán)境因子發(fā)生變化。烏梁素海冰封期冰體浮游植物群落與環(huán)境因子CCA分析中,各環(huán)境指標(biāo)與第二軸呈負(fù)相關(guān),且耐低溫甲藻門和硅藻門在4個象限均有分布,說明其對烏梁素海環(huán)境冰封期冰層環(huán)境條件有較好適應(yīng)性。冰層內(nèi)浮游植物與環(huán)境因子關(guān)系研究十分鮮見,有待后續(xù)研究。另外,硅藻門和甲藻門對烏梁素海冰封期冰層中主要污染因子NH4+、NO2--N、TP、TN變化響應(yīng)不顯著。烏梁素海冰封期冰下水體10個監(jiān)測點(diǎn)平均水溫0.438℃,采樣時烏梁素海湖面氣溫為-18.63℃。根據(jù)熱力學(xué)平衡理論,烏梁素海冰封期冰體溫度應(yīng)在0.438~-18.63℃(或更低)。溫度梯度是否對烏梁素海冰封期冰層中浮游植物群落分布造成影響,有待進(jìn)一步探究。本研究在冰封期條件下,對烏梁素海浮游植物群落特征與環(huán)境因子響應(yīng)關(guān)系作初步研究。冰封期與非冰封期不同,主要表現(xiàn):①冰蓋阻礙水體大氣復(fù)氧,溶解氧與水生態(tài)系統(tǒng)關(guān)系密切,含量直接影響浮游植物生消規(guī)律,溶解氧含量高有利于浮游植物生長;②冰蓋阻斷水體與外界風(fēng)力接觸,而風(fēng)力可驅(qū)動水流變化,水動力條件影響浮游植物群落結(jié)構(gòu)特征,流速可改變水體表層混合強(qiáng)度,不利于浮游植物滯留,低流速條件下,利于浮游植物繁殖。③冰表面積雪影響太陽凈輻射進(jìn)入冰下水體,光是影響單細(xì)胞浮游植物生長關(guān)鍵因子,適宜光照強(qiáng)度可延長浮游植物生長期,利于浮游植物新陳代謝和光合產(chǎn)物形成。

    [參考文獻(xiàn)]

    [1] 趙志剛,徐奇友,羅亮,等.添加碳源對松浦鏡鯉養(yǎng)殖池塘魚體生長及水質(zhì)影響[J].東北農(nóng)業(yè)大學(xué)學(xué)報,2013,44(9):105-112.

    [2] Liu X,Huang B,Huang Q,et al.Seasonal phytoplankton response to physical processes in the southern Yellow Sea[J].Journal of Sea Research,2015,95:45-55.

    [3] 謝斌,張碩,李莉,等.海州灣海洋牧場浮游植物群落結(jié)構(gòu)特征及其與水質(zhì)參數(shù)的關(guān)系[J].環(huán)境科學(xué)學(xué)報,2017,37(1):121-129.

    [4] Shi X,Zhao X,Zhang M,et al.The responses of phytoplankton communities to elevated CO2show seasonal variations in the high?ly eutrophic lake Taihu[J].Canadian Journal of Fisheries and Aquatic Sciences,2016,73(5):727-736.

    [5] 許海,朱廣偉,秦伯強(qiáng),等.氮磷比對水華藍(lán)藻優(yōu)勢形成的影響[J].中國環(huán)境科學(xué),2011,31(10):1676-1683.

    [6] Francesco P,Jokela J,Castiglioni S,et al.Water-borne pharma?ceuticals reduce phenotypic diversity and response capacity of nat?ural phytoplankton communities[J].Plos One,2017,12(3):1-18.

    [7] 胡俊,胡鑫,米瑋潔,等.多沙河流夏季浮游植物群落結(jié)構(gòu)變化及水環(huán)境因子影響分析[J].生態(tài)環(huán)境學(xué)報,2016,25(12):1974-1982.

    [8] 姜忠峰.烏梁素海綜合需水分析及生態(tài)系統(tǒng)健康評價[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué),2011.

    [9] 姜慧琴.烏梁素海營養(yǎng)鹽在冰體中的空間分布及其在凍融過程中釋放規(guī)律的試驗(yàn)研究[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué),2011.

    [10] 李興,張樹禮,李暢游,等.烏梁素海浮游植物群落特征分析[J].生態(tài)環(huán)境學(xué)報,2012,21(11):1865-1869.

    [11] 國家環(huán)境保護(hù)總局水和廢水監(jiān)測分析方法編委會.水和廢水監(jiān)測分析方法[M].北京:中國環(huán)境科學(xué)出版社,2002.

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

    [13] 金相燦,屠清瑛.湖泊富營養(yǎng)化調(diào)查規(guī)范[M].北京:中國環(huán)境科學(xué)出版社,1990.

    [14] Licursi M,Sierra M V,Gómez N.Diatom assemblages from a tur?bid coastal plain estuary:Río de la Plata(South America)[J].Jour?nal of Marine Systems,2006,62(1-2):35-45.

    [15] Lopes M R M,Bicudo C E,Ferragut M C.Short tern spatial and temporal variation of phytoplankton in a shallow tropical oligotro?phic reservoir,southeast Brazil[J].Hydrobiologia,2005,542(1):235-247.

    [16] Muylaert K,Sabbe K,Vyverman W.Spatial and temporal dynam?ics of phytoplankton communities in a freshwater tidal estuary(Schelde,Belgium)[J].Estuarine Coastal and Shelf Science,2000,50(5):673-687.

    [17] 譚香,夏小玲,程曉莉,等.丹江口水庫浮游植物群落時空動態(tài)及其多樣性指數(shù)[J].環(huán)境科學(xué),2011,32(10):2875-2882.

    [18] 楊虹,汪益嬪,薄芳芳,等.淀山湖浮游植物群落時空分布動態(tài)初步研究[J].華東師范大學(xué)學(xué)報:自然科學(xué)版,2010(6):54-63.[19]楊虹.淀山湖浮游植物群落時空分布及與環(huán)境因子關(guān)系的研究[D].上海:華東師范大學(xué),2010.

    [20] 吳波,陳德輝,吳瓊,等.黃浦江浮游植物群落結(jié)構(gòu)及其對水環(huán)境的指示作用[J].武漢植物學(xué)研究,2007,25(5):467-472.

    [21] 夏品華,馬建榮,李存雄,等.紅楓湖水庫冬春季浮游生物群落與環(huán)境因子的典范對應(yīng)分析[J].環(huán)境科學(xué)研究,2011,24(4):378-386.

    [22] 孟順龍,陳家長,范立民,等.2007年太湖五里湖浮游植物生態(tài)學(xué)特征[J].湖泊科學(xué),2009,21(6):845-854.

    [23] 陳興群,迪克曼.南極威德爾海陸緣固冰區(qū)葉綠素a及硅藻的分布[J].海洋學(xué)報,1989,11(4):501-509.

    [24] 李喆,霍堂斌,唐富江,等.松花江哈爾濱段冰下浮游生物群落結(jié)構(gòu)與環(huán)境因子的相關(guān)分析[J].水產(chǎn)學(xué)雜志,2014,27(6):44-50.

    [25] 祁峰,劉谞,趙艷紅,等.海水越冬池冰下浮游生物群落結(jié)構(gòu)特征的研究[J].水產(chǎn)科學(xué),2013,32(1):7-14.

    [26] 李建茹,李暢游,李興,等.烏梁素海浮游植物群落特征及其與環(huán)境因子的典范對應(yīng)分析[J].生態(tài)環(huán)境學(xué)報,2013,22(6):1032-1040.

    [27] 李興,李建茹,李暢游.內(nèi)蒙古烏梁素海浮游植物優(yōu)勢種的生態(tài)位分析[J].水生態(tài)學(xué)雜志,2017,38(6):40-47.

    [28] 李曉鈺,于洪賢,馬成學(xué).松花江哈爾濱段浮游植物群落典范對應(yīng)分析及多樣性分析[J].東北林業(yè)大學(xué)學(xué)報,2013,41(10):103-107.

    猜你喜歡
    冰體烏梁素海甲藻
    甲藻甾醇在海洋和湖泊古環(huán)境重建中的應(yīng)用
    彈體高速侵徹冰體研究
    高速彈體侵徹冰材料過程數(shù)值模擬研究
    冰體質(zhì)量和撞擊角度對船首結(jié)構(gòu)碰撞性能的影響
    養(yǎng)只恐龍當(dāng)寵物
    學(xué)與玩(2019年8期)2019-10-29 03:32:16
    基于船-水-冰耦合技術(shù)的撞擊參數(shù)對船冰碰撞性能的影響
    疑似甲藻水華監(jiān)測及應(yīng)急管理探究
    十字路口的烏梁素海
    中國三峽(2017年4期)2017-06-06 10:44:22
    烏梁素海地區(qū)近50年氣候變化特征分析
    淺談烏梁素海的環(huán)境保護(hù)與水土保持工作措施
    中國水利(2015年16期)2015-02-28 15:14:47
    国产一区二区三区综合在线观看 | 啦啦啦啦在线视频资源| 日韩电影二区| 国产成人精品久久久久久| 日韩制服骚丝袜av| 性色avwww在线观看| 嘟嘟电影网在线观看| 日韩一区二区视频免费看| 亚洲中文av在线| 亚洲精品久久久久久婷婷小说| 免费不卡的大黄色大毛片视频在线观看| 黄色怎么调成土黄色| 狂野欧美白嫩少妇大欣赏| 97在线人人人人妻| 免费av不卡在线播放| 麻豆国产97在线/欧美| 中国国产av一级| 美女cb高潮喷水在线观看| 久热久热在线精品观看| 内射极品少妇av片p| 成人毛片a级毛片在线播放| 久久99热这里只有精品18| 国产亚洲午夜精品一区二区久久| 色婷婷久久久亚洲欧美| 亚洲国产毛片av蜜桃av| 美女脱内裤让男人舔精品视频| 亚洲精品日韩在线中文字幕| 国产精品一区二区在线不卡| 一级av片app| 蜜桃在线观看..| 日韩欧美 国产精品| 99久久精品热视频| 国产精品久久久久久av不卡| av线在线观看网站| 免费观看在线日韩| 精品午夜福利在线看| 日本一二三区视频观看| 亚洲av成人精品一区久久| 大香蕉久久网| 精品一区二区免费观看| 免费观看的影片在线观看| 国产爽快片一区二区三区| 在线精品无人区一区二区三 | 国产精品久久久久久精品古装| 在线观看免费视频网站a站| 亚洲电影在线观看av| 日日撸夜夜添| 久久国内精品自在自线图片| 久久婷婷青草| 亚洲精品,欧美精品| 成人18禁高潮啪啪吃奶动态图 | 国产精品人妻久久久影院| 一本久久精品| 在线观看一区二区三区激情| 午夜福利网站1000一区二区三区| 精品久久久久久久久av| 国产视频内射| 91狼人影院| 免费高清在线观看视频在线观看| 中国国产av一级| 色视频www国产| 日韩一区二区视频免费看| 精品人妻偷拍中文字幕| 国产亚洲最大av| 国产中年淑女户外野战色| 成人一区二区视频在线观看| 欧美三级亚洲精品| 日日撸夜夜添| 国产 精品1| 人妻 亚洲 视频| 日韩av在线免费看完整版不卡| 午夜激情久久久久久久| 少妇人妻精品综合一区二区| 亚洲国产日韩一区二区| 中国三级夫妇交换| 久久久久国产精品人妻一区二区| 欧美一级a爱片免费观看看| 免费av不卡在线播放| 99热国产这里只有精品6| 久久韩国三级中文字幕| 丝瓜视频免费看黄片| 黄色日韩在线| 亚洲国产精品专区欧美| 免费观看的影片在线观看| 免费黄频网站在线观看国产| 国产乱人视频| 欧美国产精品一级二级三级 | 秋霞伦理黄片| 亚洲精品久久久久久婷婷小说| 日韩av免费高清视频| 久久久色成人| 精品一品国产午夜福利视频| 欧美日韩视频精品一区| 能在线免费看毛片的网站| 丝袜喷水一区| 男女啪啪激烈高潮av片| av国产精品久久久久影院| 高清毛片免费看| 成人亚洲欧美一区二区av| 亚洲精品乱久久久久久| 国产人妻一区二区三区在| 另类亚洲欧美激情| 亚洲欧美精品专区久久| videossex国产| 最近最新中文字幕大全电影3| 久久久久久久久久人人人人人人| 大码成人一级视频| 久久久久国产精品人妻一区二区| 看免费成人av毛片| 中国三级夫妇交换| 色视频www国产| 亚洲国产精品999| 日韩中字成人| 精品一区二区三区视频在线| 亚洲av中文字字幕乱码综合| 亚洲成人中文字幕在线播放| 亚洲婷婷狠狠爱综合网| 又黄又爽又刺激的免费视频.| 夜夜看夜夜爽夜夜摸| 中文精品一卡2卡3卡4更新| 精品人妻熟女av久视频| 亚洲无线观看免费| 亚洲国产精品专区欧美| 国产精品嫩草影院av在线观看| 一级爰片在线观看| 国产午夜精品久久久久久一区二区三区| 日韩成人av中文字幕在线观看| 晚上一个人看的免费电影| .国产精品久久| 日韩三级伦理在线观看| 亚洲精品第二区| 在线精品无人区一区二区三 | 国产日韩欧美在线精品| 成人毛片a级毛片在线播放| 国产精品精品国产色婷婷| 国产v大片淫在线免费观看| 五月开心婷婷网| 国产精品一及| 少妇人妻精品综合一区二区| 黑丝袜美女国产一区| 欧美变态另类bdsm刘玥| 91aial.com中文字幕在线观看| 亚洲精品久久午夜乱码| 男人添女人高潮全过程视频| 观看av在线不卡| 国产免费一区二区三区四区乱码| 午夜视频国产福利| 欧美bdsm另类| 亚洲av免费高清在线观看| 伊人久久国产一区二区| 亚洲精品久久久久久婷婷小说| 观看美女的网站| 欧美成人精品欧美一级黄| 久久国产精品大桥未久av | 高清欧美精品videossex| 午夜免费观看性视频| 欧美精品一区二区免费开放| 亚洲欧美清纯卡通| 国产精品99久久久久久久久| 亚洲国产精品专区欧美| 欧美丝袜亚洲另类| 九九在线视频观看精品| 一级毛片 在线播放| 婷婷色综合www| 麻豆乱淫一区二区| 精品久久久噜噜| 纵有疾风起免费观看全集完整版| 天美传媒精品一区二区| 国产亚洲av片在线观看秒播厂| 人妻夜夜爽99麻豆av| 免费久久久久久久精品成人欧美视频 | 国产精品久久久久久久久免| 天堂中文最新版在线下载| 我要看日韩黄色一级片| 97在线人人人人妻| 黄片无遮挡物在线观看| 在现免费观看毛片| 伊人久久精品亚洲午夜| 国产男人的电影天堂91| 亚洲av不卡在线观看| 在线天堂最新版资源| 亚洲一区二区三区欧美精品| 一级毛片 在线播放| 一区二区三区乱码不卡18| 性色av一级| 少妇的逼水好多| av.在线天堂| 爱豆传媒免费全集在线观看| 美女福利国产在线 | av视频免费观看在线观看| 人妻夜夜爽99麻豆av| 亚洲伊人久久精品综合| av卡一久久| 欧美日韩亚洲高清精品| 日韩国内少妇激情av| 免费观看a级毛片全部| 天美传媒精品一区二区| 99热6这里只有精品| 国产色婷婷99| 国产高清不卡午夜福利| 亚洲av.av天堂| 最近手机中文字幕大全| 久久久久久久亚洲中文字幕| 亚洲国产高清在线一区二区三| av免费观看日本| 在线观看美女被高潮喷水网站| 简卡轻食公司| 午夜免费观看性视频| 在线观看av片永久免费下载| 2018国产大陆天天弄谢| 亚洲国产最新在线播放| 99久久人妻综合| 免费大片黄手机在线观看| 成人黄色视频免费在线看| 视频中文字幕在线观看| 久热久热在线精品观看| 日本黄色日本黄色录像| 欧美日韩视频精品一区| 国精品久久久久久国模美| 如何舔出高潮| 免费大片黄手机在线观看| 男女下面进入的视频免费午夜| 精品少妇久久久久久888优播| 欧美3d第一页| 亚洲最大成人中文| av国产免费在线观看| 亚洲国产欧美人成| 色婷婷久久久亚洲欧美| 在线观看人妻少妇| 人人妻人人澡人人爽人人夜夜| 一区在线观看完整版| 啦啦啦视频在线资源免费观看| 99热这里只有是精品在线观看| 久久av网站| av免费观看日本| 少妇猛男粗大的猛烈进出视频| 精品99又大又爽又粗少妇毛片| 99国产精品免费福利视频| 最黄视频免费看| 美女福利国产在线 | 日韩,欧美,国产一区二区三区| 3wmmmm亚洲av在线观看| 91久久精品国产一区二区成人| 国产爽快片一区二区三区| 亚洲成色77777| 一本色道久久久久久精品综合| 男男h啪啪无遮挡| 国产精品一区www在线观看| 日韩电影二区| 国产精品无大码| 91精品国产国语对白视频| 高清毛片免费看| 三级经典国产精品| 久久99蜜桃精品久久| 91在线精品国自产拍蜜月| 欧美丝袜亚洲另类| 亚洲在久久综合| 少妇人妻 视频| 亚洲欧美日韩卡通动漫| 午夜激情福利司机影院| 午夜视频国产福利| 国产黄色免费在线视频| 人人妻人人添人人爽欧美一区卜 | 国产成人午夜福利电影在线观看| 男女啪啪激烈高潮av片| 久久综合国产亚洲精品| 建设人人有责人人尽责人人享有的 | 蜜臀久久99精品久久宅男| 99久久精品国产国产毛片| 国产在线免费精品| 少妇裸体淫交视频免费看高清| 啦啦啦视频在线资源免费观看| 中文在线观看免费www的网站| 亚洲精品日韩在线中文字幕| 国语对白做爰xxxⅹ性视频网站| 国产成人一区二区在线| 我要看日韩黄色一级片| 国产在线男女| 多毛熟女@视频| 国产高清三级在线| 精品一区二区免费观看| 丝瓜视频免费看黄片| 国产真实伦视频高清在线观看| 男人和女人高潮做爰伦理| 尾随美女入室| 多毛熟女@视频| 国产精品麻豆人妻色哟哟久久| 日本av免费视频播放| 精品人妻熟女av久视频| 亚洲电影在线观看av| 男男h啪啪无遮挡| 亚洲在久久综合| 国产视频首页在线观看| 人妻 亚洲 视频| 精品久久久久久电影网| 91精品国产九色| 国产极品天堂在线| 最近的中文字幕免费完整| 永久网站在线| 毛片女人毛片| 看十八女毛片水多多多| 女人十人毛片免费观看3o分钟| 男人添女人高潮全过程视频| xxx大片免费视频| 国产乱人偷精品视频| 亚洲精品日韩av片在线观看| 色视频www国产| 黄色配什么色好看| 3wmmmm亚洲av在线观看| 日日摸夜夜添夜夜爱| 国产欧美日韩一区二区三区在线 | 简卡轻食公司| 欧美激情国产日韩精品一区| 日日撸夜夜添| 国产中年淑女户外野战色| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 亚洲国产精品成人久久小说| 欧美人与善性xxx| 亚洲综合色惰| 人妻少妇偷人精品九色| 一个人看视频在线观看www免费| 最新中文字幕久久久久| 一本一本综合久久| 国产成人一区二区在线| 日日撸夜夜添| 色视频www国产| 国产在线免费精品| 大片电影免费在线观看免费| 久久精品久久久久久久性| 欧美激情极品国产一区二区三区 | 舔av片在线| 亚洲人成网站在线播| 3wmmmm亚洲av在线观看| 亚洲色图综合在线观看| 少妇精品久久久久久久| 国产真实伦视频高清在线观看| 国产精品久久久久久av不卡| 久久久精品免费免费高清| 在线观看一区二区三区| 国产精品久久久久久精品电影小说 | 日韩av不卡免费在线播放| 一区二区三区四区激情视频| 一个人免费看片子| 男女边摸边吃奶| 国产精品一区二区在线不卡| 亚洲人成网站在线播| 肉色欧美久久久久久久蜜桃| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 插逼视频在线观看| 久久ye,这里只有精品| 直男gayav资源| 亚洲激情五月婷婷啪啪| 卡戴珊不雅视频在线播放| 亚洲三级黄色毛片| 免费观看a级毛片全部| 亚洲精品日本国产第一区| 免费大片黄手机在线观看| 亚洲内射少妇av| 中文字幕制服av| 天天躁夜夜躁狠狠久久av| 卡戴珊不雅视频在线播放| 欧美日本视频| 国产精品一及| 人妻系列 视频| 夜夜看夜夜爽夜夜摸| 深爱激情五月婷婷| 老师上课跳d突然被开到最大视频| 三级国产精品欧美在线观看| 成人特级av手机在线观看| 久热久热在线精品观看| 人体艺术视频欧美日本| 一级毛片久久久久久久久女| 国产色婷婷99| 亚洲美女黄色视频免费看| 精品一区二区三区视频在线| 精品酒店卫生间| 日韩av免费高清视频| 精品99又大又爽又粗少妇毛片| av.在线天堂| av在线播放精品| av专区在线播放| 一区在线观看完整版| 欧美xxxx性猛交bbbb| 欧美区成人在线视频| 久久久久人妻精品一区果冻| 美女中出高潮动态图| 亚洲怡红院男人天堂| 久久99热6这里只有精品| 国产黄色免费在线视频| 性色av一级| 亚洲天堂av无毛| 欧美极品一区二区三区四区| 国产一区二区在线观看日韩| 日日啪夜夜爽| 午夜激情久久久久久久| 亚洲精华国产精华液的使用体验| 久久ye,这里只有精品| 久久av网站| 最近中文字幕高清免费大全6| 亚洲无线观看免费| 校园人妻丝袜中文字幕| 三级经典国产精品| av不卡在线播放| 中文天堂在线官网| 成年免费大片在线观看| 亚洲欧美日韩无卡精品| 亚洲精品乱码久久久v下载方式| 最新中文字幕久久久久| 亚洲图色成人| av又黄又爽大尺度在线免费看| 亚洲国产毛片av蜜桃av| 欧美少妇被猛烈插入视频| kizo精华| 国产高清有码在线观看视频| 岛国毛片在线播放| 成人无遮挡网站| av天堂中文字幕网| 国产综合精华液| 夜夜爽夜夜爽视频| 99久久精品热视频| 久久99热6这里只有精品| 身体一侧抽搐| 欧美成人a在线观看| 国产又色又爽无遮挡免| 97超视频在线观看视频| 啦啦啦啦在线视频资源| 免费看av在线观看网站| 国产一区二区三区综合在线观看 | 国产一区亚洲一区在线观看| 精品99又大又爽又粗少妇毛片| 中文字幕av成人在线电影| 在线观看免费高清a一片| 18+在线观看网站| 日韩制服骚丝袜av| 如何舔出高潮| 久久99蜜桃精品久久| 亚洲欧美一区二区三区国产| 欧美三级亚洲精品| 亚洲精品日本国产第一区| 亚洲一区二区三区欧美精品| 亚洲av不卡在线观看| 精品国产三级普通话版| 乱码一卡2卡4卡精品| 免费久久久久久久精品成人欧美视频 | 国产av码专区亚洲av| 男人爽女人下面视频在线观看| 99精国产麻豆久久婷婷| 毛片一级片免费看久久久久| 免费看光身美女| 成人影院久久| 岛国毛片在线播放| 99久国产av精品国产电影| 日韩av在线免费看完整版不卡| 日韩一区二区三区影片| 国产精品av视频在线免费观看| 极品少妇高潮喷水抽搐| 国国产精品蜜臀av免费| 一个人看视频在线观看www免费| 啦啦啦视频在线资源免费观看| 国产白丝娇喘喷水9色精品| 国产免费一区二区三区四区乱码| 寂寞人妻少妇视频99o| 99热全是精品| 国产精品av视频在线免费观看| 国产大屁股一区二区在线视频| 精品午夜福利在线看| 欧美高清成人免费视频www| 少妇裸体淫交视频免费看高清| 看免费成人av毛片| av不卡在线播放| 久久久午夜欧美精品| 晚上一个人看的免费电影| 日韩大片免费观看网站| 亚洲精华国产精华液的使用体验| 91在线精品国自产拍蜜月| 亚洲精品自拍成人| 国产精品一二三区在线看| 看非洲黑人一级黄片| tube8黄色片| 黑人猛操日本美女一级片| 久久精品国产亚洲av涩爱| 最后的刺客免费高清国语| av在线蜜桃| 少妇 在线观看| 99热国产这里只有精品6| 国产在线一区二区三区精| 啦啦啦在线观看免费高清www| 国产午夜精品一二区理论片| 久久精品久久精品一区二区三区| 欧美日韩综合久久久久久| 欧美精品亚洲一区二区| 少妇被粗大猛烈的视频| 亚洲第一av免费看| 色视频在线一区二区三区| 最近中文字幕2019免费版| 51国产日韩欧美| 国产午夜精品久久久久久一区二区三区| 久久久久久久国产电影| 国产亚洲欧美精品永久| 国产一区二区三区综合在线观看 | 国产高潮美女av| 国产精品国产三级国产专区5o| 精品人妻视频免费看| 在线观看人妻少妇| 丰满迷人的少妇在线观看| 久久国产精品大桥未久av | 夜夜爽夜夜爽视频| 亚洲精品久久午夜乱码| 身体一侧抽搐| 极品教师在线视频| 国产成人精品婷婷| 黄色配什么色好看| 下体分泌物呈黄色| 少妇被粗大猛烈的视频| h视频一区二区三区| 啦啦啦中文免费视频观看日本| 欧美日韩综合久久久久久| 免费黄色在线免费观看| 日韩不卡一区二区三区视频在线| 精品一区二区三卡| 女人十人毛片免费观看3o分钟| 日本免费在线观看一区| 热re99久久精品国产66热6| 丰满人妻一区二区三区视频av| 久热这里只有精品99| 在线观看免费高清a一片| 久久久久久九九精品二区国产| 男人狂女人下面高潮的视频| 亚洲av福利一区| 日韩国内少妇激情av| 热99国产精品久久久久久7| av网站免费在线观看视频| 麻豆乱淫一区二区| 亚洲熟女精品中文字幕| 亚洲精品国产av蜜桃| 亚洲欧美成人综合另类久久久| 久久久久久久国产电影| 99热国产这里只有精品6| 国产在线一区二区三区精| 丰满人妻一区二区三区视频av| 亚洲av在线观看美女高潮| 亚洲国产精品专区欧美| videos熟女内射| 成人国产av品久久久| 黄色一级大片看看| 在线 av 中文字幕| 99热国产这里只有精品6| h视频一区二区三区| 国产精品女同一区二区软件| 大香蕉久久网| 一本久久精品| 欧美日韩在线观看h| 日本欧美视频一区| 日韩强制内射视频| 亚洲国产日韩一区二区| 一区二区三区乱码不卡18| 91aial.com中文字幕在线观看| 乱码一卡2卡4卡精品| 少妇裸体淫交视频免费看高清| 熟女电影av网| 免费在线观看成人毛片| 搡女人真爽免费视频火全软件| 秋霞伦理黄片| 91久久精品国产一区二区成人| 久久久久久久久久久免费av| 国产在视频线精品| 精品人妻偷拍中文字幕| a 毛片基地| 亚洲欧美日韩无卡精品| 午夜精品国产一区二区电影| 亚洲综合色惰| 亚洲三级黄色毛片| 舔av片在线| 久久久久久久久久成人| 少妇猛男粗大的猛烈进出视频| h日本视频在线播放| 91精品伊人久久大香线蕉| 在线天堂最新版资源| 热99国产精品久久久久久7| 婷婷色麻豆天堂久久| 亚洲精品456在线播放app| 黄片wwwwww| av在线app专区| 又粗又硬又长又爽又黄的视频| 在线观看免费视频网站a站| 99九九线精品视频在线观看视频| 欧美+日韩+精品| 久久6这里有精品| h日本视频在线播放| 18禁裸乳无遮挡免费网站照片| 精品国产露脸久久av麻豆| 在线免费十八禁| 免费人妻精品一区二区三区视频| 超碰97精品在线观看| 亚洲一级一片aⅴ在线观看| 亚洲欧美日韩东京热| av女优亚洲男人天堂| 久久精品久久久久久噜噜老黄| 久久久精品免费免费高清| 亚洲精品国产av成人精品| 精华霜和精华液先用哪个| 亚洲丝袜综合中文字幕| 国产精品精品国产色婷婷| 日本黄大片高清| 特大巨黑吊av在线直播| 麻豆国产97在线/欧美| 久久女婷五月综合色啪小说| 妹子高潮喷水视频| 日本猛色少妇xxxxx猛交久久| 国产黄片美女视频| 久久久久性生活片| 各种免费的搞黄视频| 三级经典国产精品|