• <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
    亚洲熟妇中文字幕五十中出| 99久久综合精品五月天人人| 男男h啪啪无遮挡| 久热爱精品视频在线9| 国产av在哪里看| 18禁国产床啪视频网站| 这个男人来自地球电影免费观看| 桃色一区二区三区在线观看| 精品国产国语对白av| 丝袜在线中文字幕| 免费av毛片视频| 亚洲人成77777在线视频| 青草久久国产| 啦啦啦韩国在线观看视频| 精品不卡国产一区二区三区| 国产主播在线观看一区二区| av超薄肉色丝袜交足视频| 老司机午夜福利在线观看视频| 中文字幕av电影在线播放| 中文亚洲av片在线观看爽| 91麻豆精品激情在线观看国产| 欧美黑人欧美精品刺激| 成人手机av| 伦理电影免费视频| 日本黄色视频三级网站网址| 免费av毛片视频| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 岛国在线观看网站| 日本 av在线| 午夜亚洲福利在线播放| 色av中文字幕| 亚洲自偷自拍图片 自拍| 人人妻,人人澡人人爽秒播| 99精品欧美一区二区三区四区| 老鸭窝网址在线观看| 中出人妻视频一区二区| 欧美日韩精品网址| 少妇被粗大的猛进出69影院| 亚洲五月天丁香| 老司机福利观看| 搡老熟女国产l中国老女人| 女人被躁到高潮嗷嗷叫费观| 欧美日本中文国产一区发布| 国产成人一区二区三区免费视频网站| 精品国产乱子伦一区二区三区| 亚洲中文字幕日韩| 一二三四社区在线视频社区8| 国产亚洲欧美在线一区二区| 精品国产一区二区久久| 国产伦一二天堂av在线观看| av中文乱码字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 久久久国产精品麻豆| 两性夫妻黄色片| 欧美大码av| 少妇 在线观看| 一区二区日韩欧美中文字幕| 日韩高清综合在线| 久久午夜综合久久蜜桃| 老司机在亚洲福利影院| 久久中文看片网| 人成视频在线观看免费观看| 国产欧美日韩综合在线一区二区| 人妻久久中文字幕网| 欧美一级a爱片免费观看看 | 曰老女人黄片| 亚洲激情在线av| 国产精品久久视频播放| 91老司机精品| 黑人巨大精品欧美一区二区蜜桃| 午夜福利免费观看在线| 97超级碰碰碰精品色视频在线观看| 久久精品国产综合久久久| 久久欧美精品欧美久久欧美| 91av网站免费观看| 欧洲精品卡2卡3卡4卡5卡区| 动漫黄色视频在线观看| 久久久久久久午夜电影| 欧美中文日本在线观看视频| 日韩精品中文字幕看吧| 中文字幕人成人乱码亚洲影| 男人的好看免费观看在线视频 | 精品久久久久久成人av| 侵犯人妻中文字幕一二三四区| 国产野战对白在线观看| 99国产精品免费福利视频| 精品免费久久久久久久清纯| 精品国产一区二区久久| 中文字幕最新亚洲高清| 欧美激情久久久久久爽电影 | 国产成人av教育| 精品久久蜜臀av无| 国产在线精品亚洲第一网站| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成人久久性| 欧美黑人精品巨大| or卡值多少钱| 一区二区三区精品91| 在线观看免费午夜福利视频| 欧美另类亚洲清纯唯美| 91精品国产国语对白视频| 亚洲成人久久性| 久久久久久大精品| 免费久久久久久久精品成人欧美视频| 在线观看www视频免费| 国产精品二区激情视频| 老司机午夜十八禁免费视频| 天天一区二区日本电影三级 | 国产成人av激情在线播放| 国产一级毛片七仙女欲春2 | 久久久国产成人精品二区| 亚洲一区二区三区不卡视频| 在线观看免费视频网站a站| 一级片免费观看大全| 国产片内射在线| 亚洲中文字幕日韩| 露出奶头的视频| 久久久久久亚洲精品国产蜜桃av| 九色亚洲精品在线播放| 在线天堂中文资源库| 国产成人av教育| 一本大道久久a久久精品| 男人舔女人下体高潮全视频| 制服诱惑二区| 99久久99久久久精品蜜桃| 久久国产乱子伦精品免费另类| а√天堂www在线а√下载| 免费久久久久久久精品成人欧美视频| 激情视频va一区二区三区| 国产一区二区三区综合在线观看| 午夜福利在线观看吧| 国产av又大| 又黄又粗又硬又大视频| 色av中文字幕| 亚洲天堂国产精品一区在线| 色播在线永久视频| 亚洲精品在线美女| 午夜福利一区二区在线看| 国产av精品麻豆| 日韩有码中文字幕| 日日摸夜夜添夜夜添小说| 亚洲精华国产精华精| 国产熟女午夜一区二区三区| 99在线视频只有这里精品首页| 女人被狂操c到高潮| 最近最新中文字幕大全电影3 | 久久香蕉国产精品| 国产成人影院久久av| 午夜久久久在线观看| 欧美精品亚洲一区二区| 国产不卡一卡二| 母亲3免费完整高清在线观看| 欧美丝袜亚洲另类 | 精品国产超薄肉色丝袜足j| 黄片小视频在线播放| 免费女性裸体啪啪无遮挡网站| 久久精品国产99精品国产亚洲性色 | 亚洲精品国产色婷婷电影| 亚洲熟女毛片儿| 亚洲男人天堂网一区| 一边摸一边抽搐一进一小说| 一边摸一边抽搐一进一出视频| 亚洲少妇的诱惑av| 90打野战视频偷拍视频| 欧美乱妇无乱码| 亚洲国产精品久久男人天堂| 亚洲欧美精品综合一区二区三区| 久久精品国产清高在天天线| 禁无遮挡网站| 免费一级毛片在线播放高清视频 | 别揉我奶头~嗯~啊~动态视频| 欧美大码av| 成人特级黄色片久久久久久久| 亚洲熟妇中文字幕五十中出| 国产国语露脸激情在线看| 波多野结衣巨乳人妻| 人人妻,人人澡人人爽秒播| 久久久精品欧美日韩精品| 久久九九热精品免费| 久久久久久久久免费视频了| 久久精品91蜜桃| 国产男靠女视频免费网站| 亚洲国产精品sss在线观看| 91成人精品电影| 一级a爱片免费观看的视频| 国产一区在线观看成人免费| 美女免费视频网站| 男人舔女人下体高潮全视频| 日本 av在线| 国产私拍福利视频在线观看| 午夜福利,免费看| АⅤ资源中文在线天堂| 一二三四社区在线视频社区8| 黄色视频不卡| 国产一区二区在线av高清观看| 高清黄色对白视频在线免费看| 可以在线观看的亚洲视频| 久久久久久久久免费视频了| 精品一区二区三区视频在线观看免费| 免费在线观看完整版高清| 大型av网站在线播放| 国产精品美女特级片免费视频播放器 | 国产亚洲欧美精品永久| 亚洲一区二区三区不卡视频| 亚洲精品一区av在线观看| 黄色视频,在线免费观看| 久久精品国产99精品国产亚洲性色 | 亚洲av五月六月丁香网| 久久国产精品男人的天堂亚洲| 人人妻人人澡欧美一区二区 | 女人被狂操c到高潮| 在线播放国产精品三级| 丝袜在线中文字幕| 香蕉丝袜av| aaaaa片日本免费| 国产亚洲精品综合一区在线观看 | 国产高清videossex| 91麻豆av在线| 欧美日韩精品网址| 亚洲国产精品sss在线观看| 嫩草影院精品99| 中文字幕人妻熟女乱码| 国产1区2区3区精品| 在线观看免费视频日本深夜| 国产欧美日韩一区二区三| 一级,二级,三级黄色视频| 午夜久久久在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久国产精品男人的天堂亚洲| 久久中文字幕一级| 十八禁人妻一区二区| 中文字幕色久视频| 看黄色毛片网站| 夜夜看夜夜爽夜夜摸| 精品乱码久久久久久99久播| 麻豆av在线久日| 欧美中文综合在线视频| 国产亚洲精品一区二区www| 不卡av一区二区三区| 免费看a级黄色片| 国产精品爽爽va在线观看网站 | 一个人免费在线观看的高清视频| 欧美日韩福利视频一区二区| 精品一区二区三区av网在线观看| 成人永久免费在线观看视频| 日韩 欧美 亚洲 中文字幕| 老鸭窝网址在线观看| 啦啦啦观看免费观看视频高清 | 91精品三级在线观看| 久久久久久久久久久久大奶| 午夜成年电影在线免费观看| 一区二区三区激情视频| 欧美黑人欧美精品刺激| 97超级碰碰碰精品色视频在线观看| 色婷婷久久久亚洲欧美| 99re在线观看精品视频| 在线观看午夜福利视频| 国产在线精品亚洲第一网站| 一进一出抽搐gif免费好疼| 国产精品av久久久久免费| 午夜福利成人在线免费观看| 一级黄色大片毛片| 国产精品免费一区二区三区在线| 精品熟女少妇八av免费久了| 久久人妻av系列| 亚洲精品在线美女| 久久久久久国产a免费观看| 中文字幕另类日韩欧美亚洲嫩草| 操美女的视频在线观看| 99久久99久久久精品蜜桃| 欧美日本亚洲视频在线播放| 国产精品久久久久久人妻精品电影| 亚洲欧美精品综合久久99| 老司机在亚洲福利影院| 亚洲欧美日韩无卡精品| 电影成人av| 国产精品 欧美亚洲| 中文亚洲av片在线观看爽| 国产成人精品在线电影| 操美女的视频在线观看| 免费无遮挡裸体视频| 欧美日本中文国产一区发布| 老司机午夜福利在线观看视频| 日本黄色视频三级网站网址| 淫秽高清视频在线观看| 久久久久精品国产欧美久久久| 一级黄色大片毛片| 99精品欧美一区二区三区四区| 后天国语完整版免费观看| 成人精品一区二区免费| 波多野结衣av一区二区av| 一个人观看的视频www高清免费观看 | 午夜亚洲福利在线播放| 久久伊人香网站| 午夜激情av网站| 美女高潮喷水抽搐中文字幕| 久久性视频一级片| 啦啦啦观看免费观看视频高清 | 桃色一区二区三区在线观看| 亚洲最大成人中文| 国产欧美日韩一区二区三| 欧美人与性动交α欧美精品济南到| 日日爽夜夜爽网站| 久久人人97超碰香蕉20202| 精品国产乱码久久久久久男人| 人人妻人人澡人人看| 国产伦一二天堂av在线观看| 日本 av在线| 禁无遮挡网站| 日日夜夜操网爽| 国产男靠女视频免费网站| 1024视频免费在线观看| 男男h啪啪无遮挡| 中文字幕人妻丝袜一区二区| 成人18禁高潮啪啪吃奶动态图| 搡老岳熟女国产| 两人在一起打扑克的视频| 亚洲 国产 在线| 制服诱惑二区| 悠悠久久av| 国产色视频综合| 一二三四在线观看免费中文在| 搡老妇女老女人老熟妇| 国产亚洲精品一区二区www| 久久精品影院6| 久久人人爽av亚洲精品天堂| 久热这里只有精品99| 人妻久久中文字幕网| 亚洲欧美日韩无卡精品| 国产97色在线日韩免费| 亚洲精品一区av在线观看| 精品国产一区二区久久| 亚洲精品在线观看二区| 精品一区二区三区四区五区乱码| 欧美+亚洲+日韩+国产| 在线观看日韩欧美| 女人被躁到高潮嗷嗷叫费观| 欧美日韩一级在线毛片| 国产高清视频在线播放一区| 欧美人与性动交α欧美精品济南到| 亚洲精品在线美女| 久久性视频一级片| 老司机靠b影院| tocl精华| 免费在线观看影片大全网站| 欧美色视频一区免费| 亚洲中文字幕日韩| 成年人黄色毛片网站| 精品少妇一区二区三区视频日本电影| 日韩高清综合在线| 精品一区二区三区视频在线观看免费| 91国产中文字幕| 在线观看免费视频日本深夜| 亚洲aⅴ乱码一区二区在线播放 | 少妇熟女aⅴ在线视频| 黄频高清免费视频| 成人国产综合亚洲| 午夜福利在线观看吧| 亚洲欧美精品综合一区二区三区| 又大又爽又粗| 亚洲成人免费电影在线观看| 久久午夜亚洲精品久久| 在线观看www视频免费| 午夜福利视频1000在线观看 | 精品无人区乱码1区二区| www.自偷自拍.com| 欧美丝袜亚洲另类 | 欧美一区二区精品小视频在线| 免费搜索国产男女视频| 亚洲 欧美一区二区三区| 男人舔女人下体高潮全视频| 国产麻豆69| 精品不卡国产一区二区三区| 一个人免费在线观看的高清视频| 级片在线观看| 搡老岳熟女国产| 久久香蕉激情| 又紧又爽又黄一区二区| 成人免费观看视频高清| 九色国产91popny在线| 亚洲国产看品久久| 精品一区二区三区四区五区乱码| 精品国产超薄肉色丝袜足j| 欧美在线一区亚洲| 国产激情欧美一区二区| 在线永久观看黄色视频| 成年女人毛片免费观看观看9| 18禁观看日本| 日韩高清综合在线| 老熟妇乱子伦视频在线观看| 99国产综合亚洲精品| 一级毛片女人18水好多| 国产精品香港三级国产av潘金莲| 99久久国产精品久久久| 激情视频va一区二区三区| 国产av精品麻豆| 一级毛片高清免费大全| 国产一卡二卡三卡精品| 日韩欧美国产一区二区入口| 老鸭窝网址在线观看| 一本大道久久a久久精品| 欧美激情极品国产一区二区三区| √禁漫天堂资源中文www| 国产三级在线视频| 中国美女看黄片| 欧美乱色亚洲激情| 亚洲久久久国产精品| 少妇的丰满在线观看| ponron亚洲| 桃红色精品国产亚洲av| 99国产精品一区二区三区| 国产黄a三级三级三级人| 国产精品 欧美亚洲| 国产精品久久久av美女十八| 久久香蕉激情| 黄网站色视频无遮挡免费观看| 19禁男女啪啪无遮挡网站| 黄色丝袜av网址大全| 日韩成人在线观看一区二区三区| tocl精华| 亚洲五月色婷婷综合| 亚洲av日韩精品久久久久久密| 日韩大尺度精品在线看网址 | 国产精品秋霞免费鲁丝片| 久久久久久免费高清国产稀缺| 69av精品久久久久久| 国产成人欧美| 免费高清视频大片| 中文字幕色久视频| 少妇粗大呻吟视频| 母亲3免费完整高清在线观看| 夜夜躁狠狠躁天天躁| 精品久久久久久,| 成人手机av| 色婷婷久久久亚洲欧美| 琪琪午夜伦伦电影理论片6080| 久久伊人香网站| 99国产极品粉嫩在线观看| 久久亚洲真实| 中文字幕高清在线视频| 免费av毛片视频| 精品欧美一区二区三区在线| 91字幕亚洲| 黄片大片在线免费观看| 精品国产超薄肉色丝袜足j| 久久午夜亚洲精品久久| 黑人巨大精品欧美一区二区蜜桃| 欧美成人免费av一区二区三区| 国产精品久久久久久亚洲av鲁大| 亚洲第一欧美日韩一区二区三区| 久久精品国产亚洲av香蕉五月| 久久久久久久久免费视频了| 国产一区二区三区在线臀色熟女| 欧美乱码精品一区二区三区| 日本 欧美在线| 亚洲中文字幕日韩| 精品一区二区三区四区五区乱码| 美女免费视频网站| 精品第一国产精品| 国产一区二区三区在线臀色熟女| 9热在线视频观看99| 国产成人影院久久av| 午夜视频精品福利| 十分钟在线观看高清视频www| 最好的美女福利视频网| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 欧美黄色淫秽网站| 两个人免费观看高清视频| 亚洲欧美日韩高清在线视频| 国产一区二区三区在线臀色熟女| 亚洲色图av天堂| tocl精华| 亚洲欧美日韩无卡精品| 国产aⅴ精品一区二区三区波| 欧美日本中文国产一区发布| 中文字幕精品免费在线观看视频| 国产成人影院久久av| 大型黄色视频在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 在线免费观看的www视频| 中文字幕高清在线视频| 色播亚洲综合网| 亚洲一区二区三区色噜噜| 免费不卡黄色视频| 日韩欧美一区视频在线观看| 中文字幕av电影在线播放| 搡老妇女老女人老熟妇| 久久精品91蜜桃| 老司机午夜十八禁免费视频| 午夜a级毛片| 怎么达到女性高潮| 一夜夜www| 国产亚洲av高清不卡| 国产麻豆成人av免费视频| 99re在线观看精品视频| 免费观看精品视频网站| 亚洲精品中文字幕一二三四区| 亚洲欧美日韩另类电影网站| 在线观看免费午夜福利视频| 一a级毛片在线观看| 9色porny在线观看| 中文字幕另类日韩欧美亚洲嫩草| 成人特级黄色片久久久久久久| x7x7x7水蜜桃| 欧美黑人精品巨大| 麻豆一二三区av精品| 亚洲一码二码三码区别大吗| 亚洲中文字幕一区二区三区有码在线看 | 国产精品国产高清国产av| 18禁裸乳无遮挡免费网站照片 | 午夜福利欧美成人| 亚洲在线自拍视频| 久久精品影院6| 中文字幕高清在线视频| 免费人成视频x8x8入口观看| 在线观看www视频免费| 国产国语露脸激情在线看| 国产xxxxx性猛交| 精品国产乱子伦一区二区三区| 丝袜美腿诱惑在线| 亚洲 欧美 日韩 在线 免费| 在线av久久热| 欧美激情 高清一区二区三区| 中文字幕高清在线视频| 亚洲国产精品成人综合色| 身体一侧抽搐| 成人三级做爰电影| 久久香蕉国产精品| 97人妻天天添夜夜摸| 一级作爱视频免费观看| avwww免费| 脱女人内裤的视频| 亚洲aⅴ乱码一区二区在线播放 | 日韩中文字幕欧美一区二区| 欧美丝袜亚洲另类 | 麻豆一二三区av精品| 一个人观看的视频www高清免费观看 | 亚洲片人在线观看| 窝窝影院91人妻| 51午夜福利影视在线观看| 色哟哟哟哟哟哟| 首页视频小说图片口味搜索| 中文字幕人妻丝袜一区二区| www日本在线高清视频| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一品国产午夜福利视频| 中文字幕高清在线视频| 日韩高清综合在线| 亚洲自拍偷在线| 我的亚洲天堂| 亚洲av第一区精品v没综合| 亚洲国产欧美日韩在线播放| 免费久久久久久久精品成人欧美视频| 色综合站精品国产| 精品少妇一区二区三区视频日本电影| 夜夜夜夜夜久久久久| 美女高潮喷水抽搐中文字幕| 色婷婷久久久亚洲欧美| 国产精品美女特级片免费视频播放器 | 久久中文看片网| 黄片播放在线免费| 狠狠狠狠99中文字幕| 欧美一区二区精品小视频在线| 国产成人啪精品午夜网站| 国产一区二区激情短视频| 色老头精品视频在线观看| 久久久国产成人免费| 日韩精品免费视频一区二区三区| 丁香六月欧美| 这个男人来自地球电影免费观看| 波多野结衣高清无吗| 亚洲国产精品久久男人天堂| 日韩视频一区二区在线观看| 妹子高潮喷水视频| 大码成人一级视频| 涩涩av久久男人的天堂| 亚洲欧美日韩高清在线视频| 色播亚洲综合网| 久久久久亚洲av毛片大全| 精品第一国产精品| 国产成人精品久久二区二区91| 老汉色∧v一级毛片| 美女扒开内裤让男人捅视频| 亚洲精品国产精品久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 免费不卡黄色视频| 91精品三级在线观看| 亚洲色图av天堂| 99久久精品国产亚洲精品| 中文字幕人妻丝袜一区二区| 欧美日韩精品网址| 精品一区二区三区四区五区乱码| 亚洲男人天堂网一区| 成人国产综合亚洲| 国内精品久久久久久久电影| 人人妻人人爽人人添夜夜欢视频| 国产精品国产高清国产av| 国内精品久久久久久久电影| 国产精品二区激情视频| 亚洲色图av天堂| 国产区一区二久久| 啦啦啦 在线观看视频| 侵犯人妻中文字幕一二三四区| 久久精品国产99精品国产亚洲性色 | 一区二区三区精品91| 精品国产一区二区三区四区第35| 看片在线看免费视频| 在线观看午夜福利视频| 97人妻精品一区二区三区麻豆 | 亚洲国产精品合色在线| 非洲黑人性xxxx精品又粗又长| 深夜精品福利| 国产精品1区2区在线观看.|