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

    嶗山近岸浮游植物群落結(jié)構(gòu)季節(jié)變化及其環(huán)境影響因素*

    2020-02-13 07:00:42陳碧鵑孫雪梅王曉曉陳聚法崔正國(guó)曲克明
    漁業(yè)科學(xué)進(jìn)展 2020年1期
    關(guān)鍵詞:嶗山浮游海域

    柴 然 馮 娟 陳碧鵑 夏 斌 孫雪梅 王曉曉 陳聚法 崔正國(guó) 曲克明

    嶗山近岸浮游植物群落結(jié)構(gòu)季節(jié)變化及其環(huán)境影響因素*

    柴 然1,2馮 娟1陳碧鵑2,3夏 斌2,3孫雪梅2,3①王曉曉4陳聚法2崔正國(guó)2曲克明2

    (1. 青島大學(xué)環(huán)境科學(xué)與工程學(xué)院 青島 266071;2. 中國(guó)水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所 農(nóng)業(yè)農(nóng)村部海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室 青島 266071;3. 青島海洋科學(xué)與技術(shù)試點(diǎn)國(guó)家實(shí)驗(yàn)室海洋生態(tài)與環(huán)境科學(xué)功能實(shí)驗(yàn)室 青島 266071;4. 青島市嶗山區(qū)自然資源局 青島 266106)

    為了解山東省青島市嶗山近岸海域浮游植物群落結(jié)構(gòu)變化特征及其與環(huán)境因子的關(guān)系,于2016~2017年的春、夏、秋3個(gè)季節(jié)在該海域設(shè)置15個(gè)監(jiān)測(cè)站位進(jìn)行調(diào)查,并同步監(jiān)測(cè)海區(qū)環(huán)境因子。結(jié)果顯示,共檢出浮游植物2門(mén)31屬69種。其中,硅藻門(mén)(Bacillariophyta)為絕對(duì)優(yōu)勢(shì)門(mén)類(lèi),共25屬60種,甲藻門(mén)(Pyrrophyta)有6屬9種。調(diào)查海區(qū)浮游植物細(xì)胞豐度變化范圍為(732.10~ 1142.19)×104cell/m3,平均為937.15×104cell/m3。浮游植物細(xì)胞豐度季節(jié)變化明顯,2016年最高峰出現(xiàn)在秋季;2017年最高峰出現(xiàn)在夏季;空間分布上,春、夏2個(gè)季節(jié)呈現(xiàn)由北向南遞減的趨勢(shì)。浮游植物群落結(jié)構(gòu)特征指數(shù)分析結(jié)果顯示,調(diào)查海域浮游植物的多樣性指數(shù)和豐富度指數(shù)均較好,表明該海域浮游植物群落結(jié)構(gòu)比較穩(wěn)定。聚類(lèi)分析結(jié)果顯示,春季浮游植物細(xì)胞豐度季節(jié)變化不大,2016和2017年相似性高達(dá)70%左右,但夏、秋季的浮游植物細(xì)胞豐度季節(jié)變化較大,相似性?xún)H為40%左右,這主要與環(huán)境因子變化有關(guān)。浮游植物群落結(jié)構(gòu)與環(huán)境因子相關(guān)性分析表明,浮游植物細(xì)胞豐度平均值與化學(xué)需氧量(COD)的相關(guān)系數(shù)最高,最高值為0.536,其次與之相關(guān)的雙因子組為pH與COD、COD與無(wú)機(jī)氮(DIN)、溶解氧(DO)與COD。因此,嶗山近岸海區(qū)浮游植物細(xì)胞豐度的關(guān)鍵限制因子為COD、pH、DIN和DO。

    浮游植物;嶗山近岸;群落結(jié)構(gòu);環(huán)境因子

    在海洋生態(tài)系統(tǒng)中,浮游植物是一種含有色素或色素體的自養(yǎng)單細(xì)胞微生物,其特點(diǎn)是生命力強(qiáng)、繁殖周期短、對(duì)環(huán)境變化反應(yīng)敏感,其作為食物鏈中的初級(jí)生產(chǎn)者,承擔(dān)著海洋生態(tài)系統(tǒng)中能量轉(zhuǎn)化與流動(dòng)和物質(zhì)循環(huán)的重要環(huán)節(jié)。海洋環(huán)境的變化可以影響浮游植物群落的種類(lèi)組成、結(jié)構(gòu)和細(xì)胞豐度等(王文海等, 2011)。因此,浮游植物可作為環(huán)境監(jiān)測(cè)評(píng)估中的指示生物,其群落結(jié)構(gòu)特征是評(píng)價(jià)海域營(yíng)養(yǎng)水平、污染狀況、資源現(xiàn)狀以及生產(chǎn)潛力和海洋生態(tài)系統(tǒng)穩(wěn)定性的重要指標(biāo)。

    嶗山近岸為山東青島嶗山區(qū)周邊開(kāi)放海域,北接即墨區(qū)近岸海域,西接青島市區(qū)近海,王哥莊灣、仰口灣分布于其北部海域,南部分布有太清宮灣、流清河灣、沙子口灣以及多條河流入???。海灣總面積約為188 km2,最大水深為14 m,生物資源豐富,灣內(nèi)與外海水體交換情況良好,為山東省魁蚶()等貝類(lèi)資源放流的重點(diǎn)海域(張波等, 2015)。近岸海域是海陸之間相互作用最強(qiáng)烈的區(qū)域之一,是許多海洋生物重要的索餌場(chǎng)和洄游通道,物質(zhì)循環(huán)和能量傳遞較為活躍(宋秀凱等, 2016)。近年來(lái),人類(lèi)對(duì)近岸海域的開(kāi)發(fā)利用不斷深入,使得水質(zhì)污染日趨嚴(yán)重,生物多樣性日益降低(喬芮等, 2015)。在近海環(huán)境不斷惡化的脅迫下,浮游植物群落結(jié)構(gòu)常處于不穩(wěn)定狀態(tài)(郭術(shù)津等, 2014),因此,需要針對(duì)近岸海域生態(tài)系統(tǒng)開(kāi)展深入研究。

    近年來(lái),山東半島近海海域浮游植物群落結(jié)構(gòu)方面的研究多集中在萊州灣和膠州灣(陳碧鵑等, 2001)。研究表明,萊州灣浮游植物數(shù)量變動(dòng)與溫度呈顯著的正相關(guān),與鹽度呈顯著的負(fù)相關(guān),但與單一營(yíng)養(yǎng)鹽之間的相關(guān)關(guān)系并不顯著(寧璇璇等, 2011; 于瀟等, 2018)。膠州灣浮游植物生物量(細(xì)胞豐度、葉綠素質(zhì)量濃度)在不同月份的變化與海水溫度、營(yíng)養(yǎng)鹽結(jié)構(gòu)和營(yíng)養(yǎng)鹽濃度呈顯著正相關(guān)(王妍等, 2014)。為客觀反映近岸海域生態(tài)環(huán)境變化對(duì)浮游植物群落演變的影響,本研究對(duì)2016~2017年春、夏和秋3個(gè)季節(jié)嶗山近岸浮游植物的群落結(jié)構(gòu)特征進(jìn)行了分析,并探討其與關(guān)鍵環(huán)境因子的相關(guān)性,以期為研究人類(lèi)活動(dòng)對(duì)近岸海域生態(tài)系統(tǒng)演化的影響提供基礎(chǔ)數(shù)據(jù),同時(shí),為嶗山近岸海域開(kāi)展貝類(lèi)資源增殖放流和漁業(yè)資源保護(hù)提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 監(jiān)測(cè)區(qū)域

    本次調(diào)查的采樣時(shí)間為2016年春季(5月)、夏季(8月)、秋季(11月)和2017年的春季(5月)、夏季(9月)、秋季(11月),采樣地點(diǎn)為山東青島市嶗山灣的沿岸近海海域(120.507°~120.707°E, 36.083°~36.142°N),共布設(shè)15個(gè)監(jiān)測(cè)站位(圖1)。

    圖1 2016~2017年嶗山近岸浮游植物采樣站位

    1.2 樣品采集

    浮游植物的采集參照《海洋調(diào)查規(guī)范》(GB/T12763.6-2007)的方法。具體操作方法:在每個(gè)站位,通過(guò)淺水Ⅲ型浮游生物網(wǎng)在采樣點(diǎn)進(jìn)行垂直拖網(wǎng)采樣,所采集生物樣品采用含量為5%的甲醛溶液進(jìn)行固定,用顯微鏡對(duì)浮游植物進(jìn)行分類(lèi)鑒定與數(shù)量統(tǒng)計(jì)?,F(xiàn)場(chǎng)采樣時(shí),使用多參數(shù)水質(zhì)監(jiān)測(cè)儀(YSI 556型,美國(guó))同步測(cè)定海水表層溫度()、表層鹽度()、化學(xué)需氧量(COD)、溶解氧(DO)、pH等指標(biāo),同步采集水樣,帶回實(shí)驗(yàn)室測(cè)定無(wú)機(jī)氮(DIN)、活性磷酸鹽(DIP)和葉綠素(Chl-)等。海水樣品的采集、儲(chǔ)存和環(huán)境參數(shù)測(cè)定嚴(yán)格按照《海洋監(jiān)測(cè)規(guī)范》(GB/T17378-2007)中的方法進(jìn)行。

    1.3 數(shù)據(jù)分析

    1.3.1 優(yōu)勢(shì)種的確定 采用種類(lèi)優(yōu)勢(shì)度()公式:

    1.3.2 群落結(jié)構(gòu)特征指數(shù)的計(jì)算 Margalef豐富度指數(shù)公式:

    Pielou均勻度指數(shù)()公式:

    1.4 群落結(jié)構(gòu)分析

    本研究通過(guò)生物統(tǒng)計(jì)學(xué)分析軟件PRIMER 7.0,對(duì)2016~2017年嶗山近岸浮游植物群落結(jié)構(gòu)進(jìn)行分析,以浮游植物種類(lèi)、細(xì)胞豐度及監(jiān)測(cè)站位組成原始數(shù)據(jù)矩陣,所監(jiān)測(cè)的站點(diǎn)為矩陣的樣本,各站點(diǎn)浮游植物種類(lèi)和豐度為矩陣的變量,將2016~2017年嶗山近岸浮游植物原始豐度數(shù)據(jù)進(jìn)行平方根轉(zhuǎn)化,以Bray-Curtis相似性系數(shù)為基礎(chǔ),構(gòu)建嶗山近岸各站站位的等級(jí)相似矩陣。同時(shí),采用BIOENV分析2016~ 2017年嶗山近岸浮游植物群落的細(xì)胞豐度與環(huán)境因子之間的相關(guān)性(周紅等, 2003)。嶗山近岸海域浮游植物豐度平面分布圖采用Surfer 8.0軟件繪制。

    2 結(jié)果與討論

    2.1 浮游植物種類(lèi)組成

    2016~2017年,嶗山近岸調(diào)查海域共檢出浮游植物2門(mén)31屬69種。其中,硅藻門(mén)(Bacillariophyta)共有25屬60種,占浮游植物的86.95%;以圓篩藻屬()的種類(lèi)數(shù)最多,達(dá)11種;其次是角毛藻屬(),達(dá)9種。甲藻門(mén)(Pyrrophyta)有6屬9種,占浮游植物的13.04%;其中,角藻屬()有2種。由此可見(jiàn),2016~2017年嶗山近岸海域浮游植物群落構(gòu)成以硅藻門(mén)和甲藻門(mén)為主,硅藻門(mén)在物種數(shù)量和細(xì)胞豐度上占絕對(duì)優(yōu)勢(shì)。浮游植物的生態(tài)類(lèi)型主要是廣溫、廣鹽的廣布性種類(lèi),如中肋骨條藻()、丹麥細(xì)柱藻()、柔弱角毛藻()、尖刺偽菱形藻()、扁面角毛藻()、星臍圓篩藻()等;近岸廣溫性種類(lèi),如派格棍形藻()、旋鏈角毛藻()、勞氏角毛藻()和冰河擬星桿藻()等;溫帶沿岸性種類(lèi),如柔弱偽菱形藻()、丹麥細(xì)柱藻()、泰晤士旋鞘藻()、柔弱幾內(nèi)亞藻()和旋鏈角毛藻()等;外洋性種類(lèi),如虹彩圓篩藻()、偏心圓篩藻()等;浮游性種類(lèi),如中華齒狀藻()和翼根管藻印度變形(f.)等;暖水性種類(lèi),如薄壁幾內(nèi)亞藻()等。

    表1為2016~2017年嶗山近岸海域浮游植物群落優(yōu)勢(shì)種組成具有明顯的季節(jié)演替現(xiàn)象。春季主要優(yōu)勢(shì)種為虹彩圓篩藻(0.048~0.174)、梭角角藻(,=0.029~0.089)、星臍圓篩藻(=0.081~0.149)等。圓篩藻屬在各站位出現(xiàn)頻率為100%,且豐度較高。角毛藻屬以北溫帶廣布性種為主,如柔弱角毛藻等。夏季主要優(yōu)勢(shì)種有虹彩圓篩藻(=0.024~0.131)、勞氏角毛藻(=0.039~0.117)、柔弱偽菱形藻(= 0.02~0.133)、星臍圓篩藻(=0.02~0.102)、中肋骨條藻(=0.02~0.099)等。夏季出現(xiàn)的優(yōu)勢(shì)種以廣布性類(lèi)群與近岸廣溫性類(lèi)群為主。秋季主要優(yōu)勢(shì)種有虹彩圓篩藻(=0.147~0.241)、卡氏角毛藻(,=0.036~0.053)、星臍圓篩藻(=0.137~ 0.221)、翼根管藻印度變形(=0.032~0.203)等。浮游性類(lèi)群在秋季浮游植物群落中占有一定的地位,但優(yōu)勢(shì)種類(lèi)出現(xiàn)不多。溫帶沿岸性物種數(shù)量有一定程度的升高,但優(yōu)勢(shì)地位卻有所下降,對(duì)秋季嶗山近岸海域浮游植物群落構(gòu)成的貢獻(xiàn)有一定程度的下降。與夏季相比,近岸廣溫性類(lèi)群中的一些物種如勞氏角毛藻依然作為優(yōu)勢(shì)種出現(xiàn)。

    表1 2016~2017年調(diào)查海域浮游植物優(yōu)勢(shì)種及優(yōu)勢(shì)度

    Tab.1 Dominant species and predominance of phytoplankton in the survey area in 2016 and 2017

    2.2 浮游植物豐度的時(shí)空分布

    2.2.1 平面分布 從圖2可以看出,調(diào)查期間,同一季度不同年份中,嶗山近岸浮游植物豐度的平面分布均有不同。2016、2017年春季,浮游植物豐度的平面分布均呈現(xiàn)從北部向南部逐漸減少的趨勢(shì)。高值區(qū)出現(xiàn)在王哥莊灣一帶,其數(shù)值范圍為(44.19~ 70.16)×104cell/m3。2017年5月,南部沿岸浮游植物細(xì)胞豐度整體大于2016年同期。夏季浮游植物豐度均遠(yuǎn)高于5月,在兩年夏季中,最大峰值出現(xiàn)的海域不同。2016年夏季,浮游植物東部近岸數(shù)量較大且分布較為均勻,南岸除在南窯河及石老人附近海域豐度較高外,其余整體豐度較低。2017年夏季,浮游植物豐度分布仍呈現(xiàn)由北向南遞減的趨勢(shì),但東部近岸水域豐度整體低于上年同期,南部浮游植物的數(shù)量分布相對(duì)較均勻,高值區(qū)出現(xiàn)在王哥莊、太清宮口水域。在2016和2017年秋季,浮游植物的分布格局與春、夏季有明顯的差異,南部近岸整體浮游植物豐度明顯增高,與北部、東部近岸持平,高值區(qū)在南窯河與石老人海水浴場(chǎng)一帶。綜上所述,嶗山近岸浮游植物豐度平面分布較密集的海區(qū)主要出現(xiàn)在王哥莊灣、南窯河與石老人海水浴場(chǎng)一帶,春、夏季浮游植物豐度分布特征為北部沿岸大于南部沿岸;秋季南、北部沿岸分布較均勻。

    2.2.2 季節(jié)變化 浮游植物細(xì)胞豐度季節(jié)變化明顯,2016年秋季>夏季>春季,而2017年夏季>春季>秋季。據(jù)中國(guó)海灣志記載(中國(guó)海灣志編纂委員會(huì), 1993),嶗山灣浮游植物群落上、下半年分別出現(xiàn)1個(gè)數(shù)量高峰,宋秀凱等(2016)的調(diào)查結(jié)果與這一變化趨勢(shì)相一致。由于本次只調(diào)查了3個(gè)季節(jié),無(wú)法體現(xiàn)出雙峰變化,但仍然可以看出,上半年有1個(gè)高峰期,且夏季浮游植物豐度大于春季。

    圖2 嶗山近岸浮游植物豐度平面分布(×103 cell/m3)

    從圖3a可以看出,調(diào)查海域春季的細(xì)胞豐度平均值為496×104cell/m3,略高于同處于黃海北部的膠州灣(175×104cell/m3) (王妍等, 2014)。此時(shí),甲藻數(shù)量占比達(dá)最高,與渤海灣(周然等, 2013)、膠州灣(孫曉霞等, 2011)春季的甲藻數(shù)量增多的結(jié)論相一致。夏季細(xì)胞豐度平均值為937×104cell/m3,略高于2012年的調(diào)查結(jié)果[(200~500)×104cell/m3](孫雪梅等, 2013)。該季節(jié)的浮游植物以硅藻為主,豐度占總量的95.6%~97.0%,甲藻占總量的3.0%~4.4%(圖3b)。在采集到的浮游植物的總物種數(shù)中,仍以硅藻為主,約占總物種數(shù)的97.0%~98.8%,甲藻占總物種數(shù)的1.2%~3.0%。在該季度監(jiān)測(cè)期間,浮游植物結(jié)構(gòu)變化相對(duì)穩(wěn)定。秋季細(xì)胞豐度平均值為621×104cell/m3,大于宋秀凱等(2016)的2012~2014年對(duì)嶗山近岸浮游植物細(xì)胞豐度調(diào)查值(26.6×104cell/m3),說(shuō)明近幾年浮游植物數(shù)量略有提升,但仍遠(yuǎn)遠(yuǎn)低于膠州灣(2900×104cell/m3)(孫曉霞等, 2011)與萊州灣(6623.22×104cell/m3)(寧璇璇等, 2011)浮游植物細(xì)胞豐度值。

    2.3 豐富度指數(shù)、多樣性指數(shù)、均勻度指數(shù)的季節(jié)變化

    Shannon-Winner多樣性指數(shù)用于判斷群落結(jié)構(gòu)的復(fù)雜程度(趙其彪等, 2015),均勻度指數(shù)用于判斷群落是否成熟穩(wěn)定(楊曉改等, 2014)。多樣性指數(shù)>1、均勻度指數(shù)>0.3,則認(rèn)為該群落多樣性較好(李廣樓等, 2006)。2016~2017年浮游植物群落結(jié)構(gòu)參數(shù)變化情況見(jiàn)圖4。從圖4可以看出,各季節(jié)的多樣性指數(shù)平均值為1.15~2.50,年平均值為2.16,各月均大于1。各季節(jié)均勻度指數(shù)的平均值為0.68~0.90,年平均值為0.82,各月均大于0.3。可以認(rèn)為,2016~2017年嶗山近岸浮游植物群落結(jié)構(gòu)較為穩(wěn)定,組成物種數(shù)較多,生物多樣性較好,物種類(lèi)生物量分布比較均勻。從全年來(lái)看,多樣性指數(shù)與均勻度指數(shù)呈現(xiàn)較為明顯的季節(jié)變化。2016年春季多樣性指數(shù)為2.46,夏季略微升高,達(dá)到2.50,最低值出現(xiàn)在秋季,為2.24。而2017年春季多樣性指數(shù)為2.34,夏季較春季有所下降,為2.26,在秋季降至最低,為1.15。兩年中均在秋季呈驟降趨勢(shì)。從圖4可以看出,均勻度指數(shù)的季節(jié)變化趨勢(shì)與多樣性指數(shù)相同,最低值都出現(xiàn)在秋季,但季節(jié)間的差異不顯著。

    圖3 2016~2017年嶗山近岸浮游植物細(xì)胞豐度變化(a)與各類(lèi)群數(shù)量百分比(b)

    圖4 2016~2017年嶗山近岸豐富度指數(shù)、多樣性指數(shù)、均勻度指數(shù)季節(jié)變化

    從圖4可以看出,豐富度指數(shù)呈較明顯的變化趨勢(shì),反映出不同季節(jié)浮游植物種類(lèi)數(shù)量差異較大。2016年為秋季>夏季>春季,而2017年豐富度指數(shù)整體下降,最高值出現(xiàn)在夏季,秋季最低。綜合對(duì)比圖3a與圖4可以看出,浮游植物細(xì)胞豐度與多樣性指數(shù)的變化趨勢(shì)一致。即在2016年3個(gè)季節(jié)都穩(wěn)步上升,浮游植物細(xì)胞豐度峰值出現(xiàn)在11月,同時(shí),本月的豐富度指數(shù)也是全年最高值,這與硅藻的大量繁殖有關(guān)。春季豐富度指數(shù)與多樣性指數(shù)變化趨勢(shì)不同,這通常與單一藻類(lèi)優(yōu)勢(shì)過(guò)高有關(guān)。2017年9月,浮游植物細(xì)胞豐度、多樣性指數(shù)和豐富度指數(shù)均處于高峰期,即物種種類(lèi)多、各物種多樣性好、分布均勻,這可能是因?yàn)橄募臼歉∮沃参锓N類(lèi)交替的季節(jié),從而使浮游植物呈現(xiàn)細(xì)胞豐度高且多樣性較好的特點(diǎn)。

    2.4 浮游植物群落結(jié)構(gòu)季節(jié)變化聚類(lèi)分析

    不同季節(jié)浮游植物群落的聚類(lèi)分析結(jié)果見(jiàn)圖5。從圖5可以看出,浮游植物群落相似性在40%水平,主要分為2個(gè)組群。第1組群為2017年的9、11月,其余月份為第2組群,在第2組群中,2016年5月與2017年5月相似度高達(dá)70%。第1組群2017年9、11月的特征種為虹彩圓篩藻和勞氏角毛藻,該組群種類(lèi)組成以近岸廣溫性種為主,外洋性種類(lèi)的細(xì)胞數(shù)量在這一組群中占有一定優(yōu)勢(shì)。2017年春季與夏季的相似度較低。從圖4還可以看出,2017年豐富度指數(shù)9月明顯高于5月。這可能是由不同藻類(lèi)生長(zhǎng)依賴(lài)于不同的影響因子所致(許海等, 2012)。如新月菱形藻、舟形藻的主要影響因子為磷酸鹽,中肋骨條藻的主要影響因子為COD(王欣等, 2014)。由于嶗山近岸海域5~9月關(guān)鍵環(huán)境因子有所改變,部分浮游植物因環(huán)境改變而減少或死亡,一些適合新環(huán)境的其他浮游植物迅速生長(zhǎng)繁殖,從而使5、9月嶗山近岸浮游植物群落構(gòu)成相似度較低。第2組群中,2016年 5月與2017年5月浮游植物群落組成相似度較高,該組群的特征種為高齒狀藻、虹彩圓篩藻、派洛棍形藻、梭角角藻、星臍圓篩藻等。

    圖5 2016~2017年嶗山近岸浮游植物群落的CLUSTER聚類(lèi)分析

    2.5 浮游植物群落結(jié)構(gòu)與環(huán)境因子關(guān)系

    嶗山近岸調(diào)查環(huán)境因子的平均值見(jiàn)表2。從表2可以看出,溫度的變化范圍為13.35℃~26.76℃,鹽度的變化范圍為31.74~32.11,DO變化范圍為5.75~ 9.30 mg/L,pH的變化范圍為7.97~8.31,COD的變化范圍為0.36~1.27 mg/L,DIP的變化范圍為5.75~ 17.53 μg/L,DIN的變化范圍為30.28~182.44 μg/L,Chl-的變化范圍為0.40~1.58 μg/L。

    將環(huán)境因子和浮游植物豐度進(jìn)行回歸統(tǒng)計(jì)學(xué)分析結(jié)果見(jiàn)表3。從表3可以看出,嶗山近岸浮游植物細(xì)胞豐度的變化與COD的相關(guān)性最高,相關(guān)系數(shù)為0.536。產(chǎn)生主要影響的雙因子組為pH與COD,相關(guān)系數(shù)為0.411;COD與DIN,相關(guān)系數(shù)為0.400;DO與COD,相關(guān)系數(shù)為0.321。

    表2 2016~2017年嶗山近岸水域環(huán)境參數(shù)

    Tab.2 Environmental parameters in waters in offshore area of Laoshan in 2016 and 2017

    表3 2016~2017年嶗山近岸浮游植物群落與環(huán)境因子的相關(guān)關(guān)系

    Tab.3 The Pearson correlation between phytoplankton community and environmental factors in offshore area of Laoshan in 2016 and 2017

    COD是評(píng)價(jià)水質(zhì)污染常用的指標(biāo),大致反映了水體中有機(jī)質(zhì)的含量(王欣等, 2014)。調(diào)查期間,COD呈明顯的季節(jié)變化。孫偉等(2016)研究表明,嶗山近岸COD濃度主要受海流擴(kuò)散及浮游植物消耗影響。2017年秋季,COD高值原因可能是該季度風(fēng)力較強(qiáng),海流擴(kuò)散作用高于其他季度所致。此外,王哥莊灣、南窯河、石老人海水浴場(chǎng)一帶由于徑流入海,COD陸源排放量較大,藻類(lèi)更適應(yīng)在富含有機(jī)顆粒營(yíng)養(yǎng)物質(zhì)的水體中生長(zhǎng),因此,細(xì)胞豐度較高。

    在本次調(diào)查中,海水環(huán)境是主要測(cè)定參數(shù)之一,pH的變化范圍較大,一年中最低值出現(xiàn)在夏季,最高值出現(xiàn)在秋季,而海洋環(huán)境是弱酸、弱堿及鹽形成的體系,海水中pH的高低很大程度上影響生物的代謝過(guò)程。近岸的營(yíng)養(yǎng)鹽主要來(lái)源于降水及生活廢水排海、沉積物再懸浮釋放和周邊海域的補(bǔ)充(過(guò)鋒等, 2015)。夏季風(fēng)力較弱,海底沉積物懸浮釋放量相對(duì)較少,且洋流擴(kuò)散作用弱,使外海對(duì)嶗山近岸海域營(yíng)養(yǎng)鹽補(bǔ)給量小(孫偉等, 2016),因此,夏季N、P營(yíng)養(yǎng)鹽含量低于其他季節(jié)。浮游植物是水域中DO的供應(yīng)者,一般浮游植物大量繁殖時(shí),水中DO含量增高,但同時(shí)DO與溫度呈負(fù)相關(guān),水溫升高使海洋動(dòng)物生長(zhǎng)繁殖加快,消耗大量的DO(劉長(zhǎng)東等, 2016),因此,夏季雖有浮游植物補(bǔ)充DO含量,但總體消耗量較大,導(dǎo)致海水DO含量較低。

    結(jié)合聚類(lèi)分析可以看出,環(huán)境因子的改變使浮游植物的群落結(jié)構(gòu)受到擾動(dòng),進(jìn)而造成群落相似性的改變。硅藻等能適應(yīng)較高pH、COD環(huán)境,與DO表現(xiàn)出負(fù)相關(guān)(王華等, 2016),如中肋骨條藻、卡氏角毛藻、薄壁幾內(nèi)亞藻等(朱旭宇等, 2017)。角毛藻、柔弱角毛藻、旋鏈角毛藻與pH和DO呈顯著正相關(guān),與COD呈顯著負(fù)相關(guān)(羅先香等, 2018),由此可見(jiàn),環(huán)境因子與浮游植物群落結(jié)構(gòu)關(guān)系密切,在實(shí)際應(yīng)用中可通過(guò)監(jiān)控環(huán)境因子的變化預(yù)測(cè)浮游植物種群的分布情況,為貝類(lèi)資源增殖放流確定合理時(shí)間與地點(diǎn)。

    3 結(jié)論

    2016~2017春、夏、秋3個(gè)季節(jié)共檢出浮游植物2門(mén)31屬69種,硅藻門(mén)25屬60種,甲藻門(mén)6屬9種。浮游植物優(yōu)勢(shì)種隨季節(jié)變化明顯,春季的主要優(yōu)勢(shì)種為虹彩圓篩藻、梭角角藻、星臍圓篩藻等;夏季主要優(yōu)勢(shì)種有虹彩圓篩藻、勞氏角毛藻、柔弱偽菱形藻等;秋季主要優(yōu)勢(shì)種有虹彩圓篩藻、卡氏角毛藻、星臍圓篩藻等,3個(gè)季節(jié)優(yōu)勢(shì)種交替明顯。

    浮游植物細(xì)胞豐度年內(nèi)變化明顯,2016年最高峰出現(xiàn)在秋季,次高峰為夏季,2017年最高峰出現(xiàn)在夏季,次高峰為秋季。在調(diào)查中,嶗山近岸海域的浮游植物多樣性指數(shù)和豐富度指數(shù)均較好,群落結(jié)構(gòu)比較穩(wěn)定。

    浮游植物群落結(jié)構(gòu)聚類(lèi)分析結(jié)果顯示,相似度在40%水平聚類(lèi)為2大類(lèi)群,其中,相似度為70%水平的為2016年5月和2017年5月,這可能與2016年5月和2017年5月的關(guān)鍵環(huán)境因子COD變化不大相關(guān)。

    對(duì)環(huán)境因子和浮游植物豐度進(jìn)行相關(guān)性分析,結(jié)果顯示,浮游植物豐度最主要的環(huán)境影響因子為COD,相關(guān)系數(shù)為0.536,相關(guān)的雙因子為pH與COD組、COD與DIN組和DO與COD組。

    Chen BJ, Chen JF, Cui Y,. The ecological characteristics of phytoplankton in the Eastern marineculture area of Laizhou Bay. Marine Fisheries Research, 2001, 22(3): 64–70 [陳碧鵑, 陳聚法, 崔毅, 等. 萊州灣東部養(yǎng)殖區(qū)浮游植物的生態(tài)特征. 漁業(yè)科學(xué)進(jìn)展, 2001, 22(3): 64–70]

    China Gulf Records Compilation Committee. China gulf records- volume 4-southern Shandong Peninsula and gulf of Jiangsu Province. Beijing: Ocean Press, 1993 [中國(guó)海灣志編纂委員會(huì). 中國(guó)海灣志–第四分冊(cè)–山東半島南部和江蘇省海灣, 北京: 海洋出版社1993]

    Guo F, Chen JF, Qu KM,. Analysis and evaluation on nutritional status and organic pollution condition in the southern offshore area of Huangdao District of Qingdao City. Progress in Fishery Sciences, 2015, 36(4): 44–50 [過(guò)鋒, 陳聚法, 張艷, 等. 黃島區(qū)南部海域營(yíng)養(yǎng)水平和有機(jī)污染狀況分析與評(píng)價(jià). 漁業(yè)科學(xué)進(jìn)展, 2015, 36(4): 44–50]

    Guo SJ, Li YQ, Zhang CX,. Phytoplankton community in the Bohai Sea and its relationship with environmental factors. Marine Science Bulletin, 2014, 33(1): 95–105 [郭術(shù)津, 李彥翹, 張翠霞, 等. 渤海浮游植物群落結(jié)構(gòu)及與環(huán)境因子的相關(guān)性分析. 海洋通報(bào), 2014, 33(1): 95–105]

    Li GL, Chen BJ, Cui Y,. Ecological characteristics of phytoplankton in the Laizhou Bay. Journal of Fishery Sciences of China, 2006, 13(2): 292–299 [李廣樓, 陳碧鵑, 崔毅,等. 萊州灣浮游植物的生態(tài)特征. 中國(guó)水產(chǎn)科學(xué), 2006, 13(2): 292–299]

    Liu CD, Yi J, Guo XF,Phytoplankton community structure in artificial reef are around Lidao, Rongcheng, and its relationship with environmental factors. Periodical of Ocean University of China (Natural Science), 2016, 46(3): 50–59 [劉長(zhǎng)東, 易堅(jiān), 郭曉峰, 等. 榮成俚島人工魚(yú)礁區(qū)浮游植物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 46(3): 50–59]

    Luo XX, Shan Y, Yang JQ. Distribution of phytoplankton community and its relationship with water environment in the Yellow River Estuary. Periodical of Ocean University of China(Natural Science), 2018, 48(4): 16–23 [羅先香, 單宇, 楊建強(qiáng). 黃河口及鄰近海域浮游植物群落分布特征及與水環(huán)境的關(guān)系. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 48(4): 16–23]

    Margalef DR. Perspectives in eclogical theory. Chicago: University of Chicago Press, 1968, 1–111

    Ning XX, Ji L, Wang G,. Phytoplankton community in the nearshore waters of Laizhou Bay in 2009. Transactions of Oceanology and Limnology, 2011(3): 97–104 [寧璇璇, 紀(jì)靈, 王剛, 等. 2009年萊州灣近岸海域浮游植物群落的結(jié)構(gòu)特征. 海洋湖沼通報(bào), 2011(3): 97–104]

    Pielou EC. An introduction to mathematical ecology. New York: Wiley-Inter-Science, 1969

    Qiao R, Luan QS, Wang J,. Distribution and controlling environmental factors of phytoplankton communities in the Laoshan Bay during spring in 2013. Marine Sciences, 2015, 39(7): 13–21 [喬芮, 欒青杉, 王俊, 等. 2013年春季嶗山灣浮游植物群落及其環(huán)境控制. 海洋科學(xué), 2015, 39(7): 13–21]

    Shannon PT, Stroupe E, Tory EM. Batch and continuous thickening. Basic theory. Solids flux for rigid spheres. Industrial and Engineering Chemistry Fundamentals, 1963, 3(2): 203–211

    Sun W, Tang XC, Xu YD,. Characteristics of nutrient sand eutrophication assessment of the Laoshan Bay, Qingdao. Transactions of Oceanology and Limnology, 2016(6): 47–54 [孫偉, 湯憲春, 徐艷東, 等. 青島嶗山灣營(yíng)養(yǎng)鹽分布特征及富營(yíng)養(yǎng)化評(píng)價(jià). 海洋湖沼通報(bào), 2016(6): 47–54]

    Song XK, Tang XC, Cheng L,. Spatio-temporal distribution of net-collected phytoplankton and its relationship with environmental factors in Laoshan Bay. Oceanologia et Limnologia Sinica, 2016, 47(6): 1205–1213 [宋秀凱, 湯憲春, 程玲, 等. 嶗山灣網(wǎng)采浮游植物時(shí)空分布特征及其與環(huán)境因子關(guān)系. 海洋與湖沼, 2016, 47(6): 1205–1213]

    Sun XM, Xia B, Guo F,. Structure characteristics of phytoplankton community and the relationship with environmental factors in offshore area of Laoshan, Qingdao. Porgress in Fishery Sciences, 2013, 34(3): 46–53 [孫雪梅, 夏斌, 過(guò)鋒, 等. 青島嶗山近岸海域浮游植物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系. 漁業(yè)科學(xué)進(jìn)展, 2013, 34(3): 46–53]

    Sun XX, Sun S, Wu YY,. Long-term changes of phytoplankton community structure in the Jiaozhou Bay. Oceanologia et Limnologia Sinica, 2011, 42(5): 639–646 [孫曉霞, 孫松, 吳玉霖, 等. 膠州灣網(wǎng)采浮游植物群落結(jié)構(gòu)的長(zhǎng)期變化. 海洋與湖沼, 2011, 42(5): 639–646]

    Wang H, Yang SP, Fang CZ,. Canonical correspondence analysis of relationship between characteristics of phytoplankton community and environmental factors in Dianchi Lake. China Environmental Science, 2016, 36(2): 544–552 [王華, 楊樹(shù)平, 房晟忠, 等. 滇池浮游植物群落特征及與環(huán)境因子的典范對(duì)應(yīng)分析.中國(guó)環(huán)境科學(xué), 2016, 36(2): 544–552]

    Wang WH, Wu SY, Feng AP. Discussion on healthy bay and index for bay health. Coastal Engineering, 2011, 30(3): 90– 97 [王文海, 吳桑云, 豐愛(ài)平. 試論健康海灣與海灣健康指標(biāo). 海岸工程, 2011, 30(3): 90–97]

    Wang X, Sheng HX, Tang YL,. Relationships between the structure of phytoplankton community and environmental factors in three artificial reef areas of Laoshan Bay. Prorgress in Fishery Sciences, 2014, 35(4): 7–12 [王欣, 盛化香, 唐衍力, 等. 嶗山灣人工魚(yú)礁區(qū)浮游植物群落結(jié)構(gòu)與環(huán)境因子的關(guān)系. 漁業(yè)科學(xué)進(jìn)展, 2014, 35(4): 7–12]

    Wang Y, Zhang Y, Wang YY,. Spatial and temporal characteristics of phytoplankton and its relation with the environmental factors in Jiaozhou Bay. Journal of Safety and Environment, 2013, 13(1): 163–170 [王妍, 張永, 王玉玨, 等. 膠州灣浮游植物的時(shí)空變化特征及其與環(huán)境因子的關(guān)系. 安全與環(huán)境學(xué)報(bào), 2013, 13(1): 163–170]

    Xu H, Qin BQ, Zhu GW. Nutrient limitation of cyanobacterial growth in different regions of Lake Taihu in summer.China Environmental Science, 2012, 32(12): 2230–2236 [許海, 秦伯強(qiáng), 朱廣偉. 太湖不同湖區(qū)夏季藍(lán)藻生長(zhǎng)的營(yíng)養(yǎng)鹽限制研究. 中國(guó)環(huán)境科學(xué), 2012, 32(12): 2230–2236]

    Yang XG, Xue Y, Zan XX,. Community structure of phytoplankton in Haizhou Bay and adjacent waters and its relationships with environmental factors. Chinese Journal of Applied Ecology, 2014, 25(7): 2123–2131 [楊曉改, 薛瑩, 昝肖肖, 等. 海州灣及其鄰近海域浮游植物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系. 應(yīng)用生態(tài)學(xué)報(bào), 2014, 25(7): 2123–2131]

    Yu X, Li XL, Yang JL,. Study on the community changes of phytoplankton in scallop culture area of Laizhou Bay. Journal of Anhui Agricultural Sciences, 2018, 46(18): 65–70 [于瀟, 李希磊, 楊俊麗, 等. 萊州灣扇貝養(yǎng)殖區(qū)浮游植物的群落變化研究. 安徽農(nóng)業(yè)科學(xué), 2018, 46(18): 65–70]

    Zhou H, Zhang ZN. Rationale of the multivariate statistical software PRIMER and its application in benthic community ecology. Journal of Ocean University of Qingdao(Natural Science), 2003, 33(1): 58–64 [周紅, 張志南. 大型多元統(tǒng)計(jì)軟件PRIMER的方法原理及其在底棲群落生態(tài)學(xué)中的應(yīng)用. 青島海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2003, 33(1): 58–64]

    Zhao QB, Sun J, Li D,. Seasonal changes of the phytoplankton along hypoxia area and adjacent waters in the East China Sea. Acta Ecologica Sinica, 2015, 35(7): 2366–2379 [趙其彪, 孫軍, 李丹, 等. 東海低氧區(qū)及鄰近水域浮游植物的季節(jié)變化. 生態(tài)學(xué)報(bào), 2015, 35(7): 2366–2379]

    Zhou R, Peng ST, Qin XB,. Phytoplankton assemblages and their relation to environmental factors by multivariate statistic analysis in Bohai Bay. Environmental Science, 2013, 13(1): 864–873 [周然, 彭士濤, 覃雪波, 等. 渤海灣浮游植物與環(huán)境因子關(guān)系的多元分析. 環(huán)境科學(xué), 2013,13(1): 864–873]

    Zhu XY, Xu HH, Xu X,. Community composition of net- collected phytoplankton and its relation to environmental factors near Lianyungang sea area. Journal of Applied Oceanography, 2017, 36(3): 387–394 [朱旭宇,許海華,許嫻,等.連云港鄰近海域網(wǎng)采浮游植物分布特征及其影響因素. 應(yīng)用海洋學(xué)學(xué)報(bào), 2017, 36(3): 387–394]

    Seasonal Variations in the Phytoplankton Community Structure and Their Environmental Impact Factors in the Offshore Area of Laoshan, Qingdao

    CHAI Ran1,2, FENG Juan1, CHEN Bijuan2,3, XIA Bin2,3, SUN Xuemei2,3①, WANG Xiaoxiao4, CHEN Jufa2, CUI Zhengguo2, QU Keming2

    (1. School of Environmental Science and Engineering, Qingdao University, Qingdao 266071; 2. Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Qingdao 266071; 3. Marine Ecology and Environmental Science Laboratory, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266071; 4. Qingdao Laoshan Marine and Fisheries Bureau, Qingdao 266106)

    To explore the seasonal variations in the phytoplankton community structure and their relationship with environmental factors in the offshore area of Laoshan, we designed 15 monitoring stations and made six investigative voyages in the spring, summer, and autumn during 2016~2017 while monitoring other environmental factors simultaneously. We detected 69 species belonging to 31 genera and two phyla of phytoplankton, most of which belong to temperate coastal species. Diatoma was the dominant phylum comprising 25 genera and 60 species, and dinoflagellates comprised 6 genera and 9 species. The seasonal variations in the dominant species were obvious. The most dominant species were,,andThe total number of phytoplankton ranged from 732.10×104to 1142.19×104cell/m3, with an average of 937.15×104cell/m3. The number of phytoplankton cells changed significantly during the year. The highest peak in 2016 appeared in autumn and the second highest peak in summer; the highest peak in 2017 appeared in summer and the second highest peak in autumn. The spatial distribution showed a decreasing trend from north to south, and the high-value areas appeared in Wanggezhuang Bay, Nanyao River, and Shilaoren Beach. The diversity and abundance of phytoplankton in the surveyed area are at a high level, and the phytoplankton community structure in this area is relatively stable, indicated by the analysis of the phytoplankton community structure index. The result of cluster analysis showed that the interannual variations in phytoplankton cell abundance are unobvious in spring, and the similarity is approximately 70% both in 2016 and 2017. However, the interannual variations in phytoplankton cell abundance are relatively obvious in summer and autumn, and the similarity is only approximately 40%, which is mainly related to environmental factors. The correlation analysis between phytoplankton community structure and environmental factors showed that the correlation coefficient between cell abundance average and COD is the highest (0.536) and the two factors that were most related to environmental factors were pH and COD, COD and DIN, and DO and COD. Therefore, the key limiting factors for phytoplankton cell abundance in the offshore area of Laoshan are COD, pH, DIN, and DO.

    Phytoplankton; Offshore area of Laoshan; Community structure; Environmental factors

    S931.2

    A

    2095-9869(2020)01-0021-10

    10.19663/j.issn2095-9869.20190125003

    * 國(guó)家重點(diǎn)研發(fā)計(jì)劃–生態(tài)環(huán)境損害鑒定評(píng)估業(yè)務(wù)化技術(shù)研究(2016YFC0503600)和國(guó)家基金委–山東省聯(lián)合基金項(xiàng)目“海洋生態(tài)與環(huán)境科學(xué)”(U1606404)共同資助 [This work was supported by the National Key Research and Development Program-Business Technology Research on Eco-Environmental Damage Assessment (2016YFC0503600), and the Joint Project of State Fund Committee and Shandong Province "Marine Ecology and Environmental Science" (U1606404)]. 柴 然,E-mail: 1215113055@qq.com

    孫雪梅,E-mail: sunxm@ysfri.ac.cn

    2019-01-25,

    2019-03-21

    http://www.yykxjz.cn/

    柴然, 馮娟, 陳碧鵑, 夏斌, 孫雪梅, 王曉曉, 陳聚法, 崔正國(guó), 曲克明. 嶗山近岸浮游植物群落結(jié)構(gòu)季節(jié)變化及其環(huán)境影響因素. 漁業(yè)科學(xué)進(jìn)展, 2020, 41(1): 21–30

    Chai R, Feng J, Chen BJ, Xia B, Sun XM, Wang XX, Chen JF, Cui ZG, Qu KM. Seasonal variations in the phytoplankton community structure and their environmental impact factors in the offshore area of Laoshan, Qingdao. Progress in Fishery Sciences, 2020, 41(1): 21–30

    SUN Xuemei, E-mail: sunxm@ysfri.ac.cn

    (編輯 陳嚴(yán))

    猜你喜歡
    嶗山浮游海域
    遺落海域
    浞河浮游藻類(lèi)的調(diào)查研究與水質(zhì)評(píng)價(jià)
    埕島海域海上獨(dú)立樁拆除方案探討
    《浮游》
    流行色(2019年10期)2019-12-06 08:13:26
    嶗山茶事
    商周刊(2018年12期)2018-07-11 01:27:22
    登嶗山
    編制嶗山
    商周刊(2017年12期)2017-06-22 12:01:57
    廣東省海域使用統(tǒng)計(jì)分析
    嶗山又日新
    商周刊(2017年24期)2017-02-02 01:42:54
    浮游藻類(lèi)在水體PAHs富集與降解中的研究
    少妇人妻精品综合一区二区| 亚洲色图综合在线观看| 人妻夜夜爽99麻豆av| 成人高潮视频无遮挡免费网站| 91精品一卡2卡3卡4卡| 国产成人精品一,二区| 免费看av在线观看网站| av在线播放精品| 黄色怎么调成土黄色| 男女国产视频网站| 精品久久久久久久久av| 国产精品国产三级国产av玫瑰| 日韩av免费高清视频| 亚洲欧洲日产国产| 啦啦啦在线观看免费高清www| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩精品成人综合77777| 哪个播放器可以免费观看大片| av线在线观看网站| 一级毛片久久久久久久久女| 欧美日韩视频高清一区二区三区二| 久久精品综合一区二区三区| 成人国产av品久久久| 久久久精品免费免费高清| 男女无遮挡免费网站观看| 夜夜爽夜夜爽视频| 精品久久久噜噜| 欧美变态另类bdsm刘玥| 人人妻人人看人人澡| 欧美老熟妇乱子伦牲交| 免费看不卡的av| 国产精品人妻久久久久久| kizo精华| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲5aaaaa淫片| 国产亚洲5aaaaa淫片| 亚洲天堂av无毛| 亚洲欧洲国产日韩| 国产伦精品一区二区三区视频9| 深夜a级毛片| 亚洲真实伦在线观看| 亚洲自拍偷在线| 哪个播放器可以免费观看大片| 国产成人91sexporn| 国产男女内射视频| 日韩av免费高清视频| 精品久久国产蜜桃| 狂野欧美激情性xxxx在线观看| 偷拍熟女少妇极品色| 国精品久久久久久国模美| 91精品伊人久久大香线蕉| 毛片女人毛片| 国产精品一区www在线观看| 日韩一区二区三区影片| 美女高潮的动态| 一级黄片播放器| 少妇人妻精品综合一区二区| 女的被弄到高潮叫床怎么办| 伊人久久国产一区二区| 亚洲人成网站高清观看| 成年女人在线观看亚洲视频 | av在线亚洲专区| 成人综合一区亚洲| 国产成人午夜福利电影在线观看| 国产精品99久久久久久久久| 看免费成人av毛片| 欧美区成人在线视频| 熟女电影av网| 日韩一本色道免费dvd| 久久精品久久久久久噜噜老黄| 七月丁香在线播放| 啦啦啦在线观看免费高清www| 中国国产av一级| 综合色av麻豆| 卡戴珊不雅视频在线播放| 最后的刺客免费高清国语| 久久久久网色| 久久99热6这里只有精品| 国产欧美日韩一区二区三区在线 | 少妇的逼水好多| 久久久久久久国产电影| 99热这里只有精品一区| 乱码一卡2卡4卡精品| 亚洲欧美一区二区三区国产| 男插女下体视频免费在线播放| 青青草视频在线视频观看| 精品一区在线观看国产| 我要看日韩黄色一级片| 天天躁日日操中文字幕| 久久久色成人| 色视频www国产| 别揉我奶头 嗯啊视频| 免费少妇av软件| 免费看不卡的av| 久久女婷五月综合色啪小说 | 久久久久久久午夜电影| 69人妻影院| 国产淫片久久久久久久久| 女的被弄到高潮叫床怎么办| 中文字幕久久专区| 伦理电影大哥的女人| 美女脱内裤让男人舔精品视频| 久久久久久九九精品二区国产| 亚洲av男天堂| 亚洲精品国产色婷婷电影| 亚洲国产精品成人综合色| 国产人妻一区二区三区在| 国产v大片淫在线免费观看| 极品少妇高潮喷水抽搐| 欧美日本视频| 免费大片黄手机在线观看| 国产午夜福利久久久久久| 欧美日韩视频精品一区| 大码成人一级视频| 一区二区三区四区激情视频| 99久国产av精品国产电影| 婷婷色av中文字幕| 国产成人a区在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产精品女同一区二区软件| 一级毛片我不卡| 最近2019中文字幕mv第一页| 久热久热在线精品观看| 亚洲图色成人| 色5月婷婷丁香| 成人免费观看视频高清| 国产高潮美女av| 精华霜和精华液先用哪个| 久久久精品94久久精品| 最近中文字幕高清免费大全6| 老司机影院成人| 超碰97精品在线观看| 亚洲国产精品国产精品| 国产伦精品一区二区三区四那| 欧美激情久久久久久爽电影| 成年女人看的毛片在线观看| 亚洲精品久久久久久婷婷小说| 日韩不卡一区二区三区视频在线| 大片电影免费在线观看免费| 丝袜美腿在线中文| 亚洲四区av| 超碰97精品在线观看| 国产色婷婷99| 亚洲经典国产精华液单| 日韩一本色道免费dvd| 国产精品一及| 偷拍熟女少妇极品色| 色婷婷久久久亚洲欧美| 中文字幕免费在线视频6| 免费人成在线观看视频色| 成年版毛片免费区| 男女啪啪激烈高潮av片| 日韩伦理黄色片| 中国国产av一级| 五月天丁香电影| 亚洲精品乱码久久久v下载方式| 亚洲av中文字字幕乱码综合| 国产黄片美女视频| 亚洲精品久久久久久婷婷小说| 亚洲熟女精品中文字幕| av又黄又爽大尺度在线免费看| 亚洲国产精品成人综合色| 精品人妻熟女av久视频| 在线亚洲精品国产二区图片欧美 | 一级片'在线观看视频| 国产成人福利小说| av国产久精品久网站免费入址| 亚洲人与动物交配视频| 别揉我奶头 嗯啊视频| 成人午夜精彩视频在线观看| 欧美日韩综合久久久久久| av在线观看视频网站免费| 婷婷色综合大香蕉| 日韩av免费高清视频| 久久韩国三级中文字幕| 99热全是精品| av在线老鸭窝| 不卡视频在线观看欧美| 欧美最新免费一区二区三区| 久久久午夜欧美精品| 日韩大片免费观看网站| 午夜免费鲁丝| 中文乱码字字幕精品一区二区三区| 91精品国产九色| 日韩一区二区三区影片| 一级毛片电影观看| 免费人成在线观看视频色| 在线观看一区二区三区| 国产视频首页在线观看| 亚洲人与动物交配视频| 91久久精品国产一区二区成人| 能在线免费看毛片的网站| 国产一级毛片在线| 亚洲天堂国产精品一区在线| 久久精品国产亚洲av涩爱| av网站免费在线观看视频| 亚洲欧洲国产日韩| 视频中文字幕在线观看| 午夜激情福利司机影院| 97在线人人人人妻| 伦理电影大哥的女人| 亚洲精品,欧美精品| 黑人高潮一二区| 国产成人免费无遮挡视频| 水蜜桃什么品种好| 久久精品国产鲁丝片午夜精品| 亚洲在久久综合| 亚洲真实伦在线观看| 精品午夜福利在线看| 亚洲aⅴ乱码一区二区在线播放| 99久国产av精品国产电影| 国产爽快片一区二区三区| 91狼人影院| 国产爱豆传媒在线观看| 国产成年人精品一区二区| 啦啦啦中文免费视频观看日本| 国产片特级美女逼逼视频| 亚洲精品国产av蜜桃| 欧美高清成人免费视频www| 最近中文字幕2019免费版| 成人黄色视频免费在线看| 国产亚洲最大av| 精品一区二区三区视频在线| 亚洲精品乱码久久久v下载方式| 亚洲欧美成人精品一区二区| 成人毛片60女人毛片免费| 最新中文字幕久久久久| 在线观看美女被高潮喷水网站| 青春草国产在线视频| 中文字幕久久专区| 亚洲色图av天堂| 日韩欧美精品v在线| 丰满乱子伦码专区| 丝袜脚勾引网站| 又黄又爽又刺激的免费视频.| 一个人看的www免费观看视频| 少妇被粗大猛烈的视频| 日本av手机在线免费观看| 久久久久性生活片| 大又大粗又爽又黄少妇毛片口| 亚洲欧美一区二区三区黑人 | 一二三四中文在线观看免费高清| 一级黄片播放器| 欧美潮喷喷水| 国产黄片美女视频| 免费观看a级毛片全部| 日本猛色少妇xxxxx猛交久久| 插阴视频在线观看视频| 国产精品久久久久久久久免| 另类亚洲欧美激情| 午夜福利网站1000一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 建设人人有责人人尽责人人享有的 | 亚洲在线观看片| 国产黄片美女视频| 成人国产麻豆网| 久久99热这里只频精品6学生| 韩国av在线不卡| 成年女人在线观看亚洲视频 | 美女内射精品一级片tv| 成年av动漫网址| 亚洲色图av天堂| 99热国产这里只有精品6| 校园人妻丝袜中文字幕| 久久久久久久大尺度免费视频| 国产精品.久久久| 亚洲av男天堂| 国产高清有码在线观看视频| 各种免费的搞黄视频| 国精品久久久久久国模美| 亚洲精品国产av蜜桃| 国产黄片美女视频| 国产亚洲91精品色在线| 国产黄色免费在线视频| 男女那种视频在线观看| 日韩一区二区视频免费看| 高清午夜精品一区二区三区| 一区二区av电影网| 嫩草影院入口| 亚洲国产成人一精品久久久| 国产精品99久久久久久久久| 亚洲真实伦在线观看| 噜噜噜噜噜久久久久久91| 国产美女午夜福利| 美女被艹到高潮喷水动态| 亚洲精品亚洲一区二区| 狂野欧美激情性xxxx在线观看| 麻豆精品久久久久久蜜桃| 国产又色又爽无遮挡免| av在线观看视频网站免费| 一级毛片久久久久久久久女| 婷婷色av中文字幕| 亚洲精品一二三| 国产亚洲5aaaaa淫片| 日本av手机在线免费观看| 国产91av在线免费观看| 人体艺术视频欧美日本| xxx大片免费视频| 男插女下体视频免费在线播放| 国产 一区 欧美 日韩| 欧美xxⅹ黑人| 91久久精品国产一区二区三区| 亚洲精品成人av观看孕妇| 久久精品熟女亚洲av麻豆精品| 亚洲欧美一区二区三区黑人 | 综合色av麻豆| 天堂俺去俺来也www色官网| 人妻制服诱惑在线中文字幕| 日韩制服骚丝袜av| 91在线精品国自产拍蜜月| 丝瓜视频免费看黄片| 最近中文字幕2019免费版| 中文精品一卡2卡3卡4更新| 久久精品夜色国产| 18+在线观看网站| 国产伦精品一区二区三区视频9| 免费观看的影片在线观看| 欧美97在线视频| 国产乱人视频| 亚洲最大成人av| 久久久久久久大尺度免费视频| 黄色日韩在线| xxx大片免费视频| 97精品久久久久久久久久精品| 伦精品一区二区三区| 国内揄拍国产精品人妻在线| 91精品伊人久久大香线蕉| 最近最新中文字幕大全电影3| 国产午夜福利久久久久久| 亚洲内射少妇av| 色视频www国产| 黄色一级大片看看| 久久久久久久久大av| 亚洲国产精品成人久久小说| 赤兔流量卡办理| 亚洲自拍偷在线| 69人妻影院| 国产在线一区二区三区精| 肉色欧美久久久久久久蜜桃 | 免费黄频网站在线观看国产| 天天躁夜夜躁狠狠久久av| 国产综合精华液| 国产乱人偷精品视频| 国产白丝娇喘喷水9色精品| 97超视频在线观看视频| 国产成人aa在线观看| 18+在线观看网站| 欧美日韩视频精品一区| 六月丁香七月| 久热这里只有精品99| 亚洲自拍偷在线| 91久久精品国产一区二区成人| 国产综合懂色| 免费看a级黄色片| 久久精品久久精品一区二区三区| 亚洲国产精品国产精品| 日韩欧美 国产精品| 午夜亚洲福利在线播放| 99热这里只有精品一区| 在线天堂最新版资源| 91久久精品国产一区二区三区| 嘟嘟电影网在线观看| 精品国产一区二区三区久久久樱花 | av线在线观看网站| 精品熟女少妇av免费看| 国产高清国产精品国产三级 | 国产久久久一区二区三区| 亚洲第一区二区三区不卡| 精品久久国产蜜桃| 亚洲精品中文字幕在线视频 | 九色成人免费人妻av| 欧美 日韩 精品 国产| 老司机影院毛片| 亚洲国产欧美人成| 欧美成人精品欧美一级黄| 大陆偷拍与自拍| 视频中文字幕在线观看| av在线老鸭窝| 亚洲精品久久久久久婷婷小说| 最后的刺客免费高清国语| 18+在线观看网站| 亚洲av.av天堂| 黄色视频在线播放观看不卡| 日韩国内少妇激情av| 亚洲精品乱久久久久久| 国产精品一二三区在线看| 看十八女毛片水多多多| 亚洲内射少妇av| 欧美日韩亚洲高清精品| 亚洲电影在线观看av| 欧美性感艳星| 欧美一级a爱片免费观看看| 欧美xxⅹ黑人| 亚洲精品国产色婷婷电影| 久久久久久久大尺度免费视频| 日韩国内少妇激情av| av免费在线看不卡| 一级黄片播放器| 日本免费在线观看一区| 听说在线观看完整版免费高清| 国产精品久久久久久久久免| 1000部很黄的大片| 校园人妻丝袜中文字幕| 国产精品麻豆人妻色哟哟久久| 亚洲精品视频女| 色播亚洲综合网| 一级爰片在线观看| 深夜a级毛片| 一级毛片我不卡| 不卡视频在线观看欧美| 日本猛色少妇xxxxx猛交久久| 国产欧美日韩一区二区三区在线 | 亚洲天堂av无毛| 在线免费十八禁| 久久精品国产亚洲av涩爱| 日本猛色少妇xxxxx猛交久久| 亚洲美女搞黄在线观看| 交换朋友夫妻互换小说| 丝袜美腿在线中文| 在线播放无遮挡| 99久久中文字幕三级久久日本| 青春草国产在线视频| 99精国产麻豆久久婷婷| 亚洲av一区综合| 日韩视频在线欧美| 久久精品国产自在天天线| 午夜激情福利司机影院| 麻豆成人av视频| 国产一区亚洲一区在线观看| 一级av片app| 国内精品美女久久久久久| 欧美zozozo另类| 黄色欧美视频在线观看| 免费观看av网站的网址| 精品一区在线观看国产| 欧美 日韩 精品 国产| 男人舔奶头视频| 一区二区三区免费毛片| 日本三级黄在线观看| 亚洲成人av在线免费| 久久久精品欧美日韩精品| 亚洲成人久久爱视频| 人妻制服诱惑在线中文字幕| av又黄又爽大尺度在线免费看| 免费看日本二区| 亚洲av欧美aⅴ国产| 国产成人a∨麻豆精品| 免费看不卡的av| 九色成人免费人妻av| 国产乱人偷精品视频| 中文欧美无线码| 成年人午夜在线观看视频| 亚洲国产日韩一区二区| 男女国产视频网站| 99久久人妻综合| 日本午夜av视频| 免费在线观看成人毛片| 色网站视频免费| 啦啦啦中文免费视频观看日本| 80岁老熟妇乱子伦牲交| 七月丁香在线播放| 久热久热在线精品观看| 国产乱人视频| 国产成人a∨麻豆精品| 亚洲性久久影院| 国产欧美另类精品又又久久亚洲欧美| 18禁动态无遮挡网站| 一级爰片在线观看| 成人亚洲欧美一区二区av| 激情 狠狠 欧美| 国产综合懂色| 精品久久久久久久久亚洲| 亚洲欧美日韩另类电影网站 | 久久精品国产亚洲av天美| 一区二区三区精品91| av女优亚洲男人天堂| 在线a可以看的网站| 69人妻影院| 日韩欧美一区视频在线观看 | 国产欧美日韩精品一区二区| www.色视频.com| 免费看日本二区| 2022亚洲国产成人精品| 亚洲最大成人手机在线| 久久女婷五月综合色啪小说 | 国产精品国产三级国产专区5o| 永久网站在线| 国产成人免费无遮挡视频| 国产成人午夜福利电影在线观看| 丰满人妻一区二区三区视频av| 内地一区二区视频在线| 亚洲精品,欧美精品| 男人爽女人下面视频在线观看| .国产精品久久| 最近中文字幕高清免费大全6| 欧美高清成人免费视频www| 又粗又硬又长又爽又黄的视频| 亚洲久久久久久中文字幕| 国产一区二区三区av在线| 国产中年淑女户外野战色| 永久免费av网站大全| 插阴视频在线观看视频| 网址你懂的国产日韩在线| 亚洲av一区综合| 国产精品国产av在线观看| 99热网站在线观看| 国产v大片淫在线免费观看| av在线观看视频网站免费| 日韩免费高清中文字幕av| 看十八女毛片水多多多| 黄色怎么调成土黄色| 欧美老熟妇乱子伦牲交| 国国产精品蜜臀av免费| 成人一区二区视频在线观看| 一本色道久久久久久精品综合| 久久精品国产a三级三级三级| 亚洲精品日韩av片在线观看| 久久精品国产亚洲网站| 男人舔奶头视频| 国产人妻一区二区三区在| 在线看a的网站| 色婷婷久久久亚洲欧美| 人体艺术视频欧美日本| 国产精品久久久久久av不卡| 久久久久国产网址| 青青草视频在线视频观看| 人人妻人人看人人澡| 国国产精品蜜臀av免费| 国产色婷婷99| 国产成年人精品一区二区| 国产高清三级在线| 一级毛片 在线播放| 国产大屁股一区二区在线视频| 99热国产这里只有精品6| 女的被弄到高潮叫床怎么办| 黄色日韩在线| 性插视频无遮挡在线免费观看| 亚洲在久久综合| 嫩草影院精品99| 少妇的逼水好多| 欧美另类一区| 成人欧美大片| 九色成人免费人妻av| 美女被艹到高潮喷水动态| 国产国拍精品亚洲av在线观看| 日韩视频在线欧美| 啦啦啦在线观看免费高清www| 最近最新中文字幕大全电影3| 久久久久久久久大av| 又爽又黄a免费视频| 一级毛片久久久久久久久女| 午夜免费观看性视频| av一本久久久久| 国产精品不卡视频一区二区| 97在线视频观看| 久久久欧美国产精品| 夫妻午夜视频| av专区在线播放| 久久精品综合一区二区三区| av专区在线播放| 国产精品三级大全| 人妻 亚洲 视频| 亚洲国产欧美人成| 国产精品人妻久久久久久| 亚洲av电影在线观看一区二区三区 | 国产欧美亚洲国产| 观看美女的网站| 久久久久国产网址| 在线观看免费高清a一片| 女人久久www免费人成看片| 午夜免费鲁丝| 伊人久久精品亚洲午夜| 丝袜美腿在线中文| 亚洲经典国产精华液单| 欧美日韩国产mv在线观看视频 | 久久女婷五月综合色啪小说 | av黄色大香蕉| 色吧在线观看| tube8黄色片| 777米奇影视久久| 成人一区二区视频在线观看| 国产成人福利小说| 久久久欧美国产精品| 亚洲av电影在线观看一区二区三区 | 三级国产精品欧美在线观看| 久久久久精品性色| 国产乱人视频| 三级国产精品片| 久久国产乱子免费精品| 色视频www国产| 26uuu在线亚洲综合色| 亚洲精品一区蜜桃| 亚洲av男天堂| 大又大粗又爽又黄少妇毛片口| 精品久久国产蜜桃| 日日摸夜夜添夜夜添av毛片| 少妇裸体淫交视频免费看高清| 亚洲在久久综合| 欧美zozozo另类| 一级毛片黄色毛片免费观看视频| 五月天丁香电影| 国产乱人视频| 欧美丝袜亚洲另类| 三级国产精品欧美在线观看| 最近2019中文字幕mv第一页| 国产精品成人在线| 狠狠精品人妻久久久久久综合| 亚洲国产av新网站| 天堂网av新在线| 久久久精品94久久精品| 欧美老熟妇乱子伦牲交| 啦啦啦中文免费视频观看日本| 黄色日韩在线| 麻豆成人午夜福利视频| 日产精品乱码卡一卡2卡三|