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

    流沙灣浮游植物群落特征季節(jié)變化及其與養(yǎng)殖活動(dòng)的關(guān)系*

    2018-12-19 08:13:38蘇家齊朱長(zhǎng)波李俊偉陳素文頡曉勇
    漁業(yè)科學(xué)進(jìn)展 2018年6期
    關(guān)鍵詞:沙灣養(yǎng)殖區(qū)浮游

    蘇家齊 朱長(zhǎng)波 李俊偉 李 婷 陳素文 頡曉勇 張 博

    ?

    流沙灣浮游植物群落特征季節(jié)變化及其與養(yǎng)殖活動(dòng)的關(guān)系*

    蘇家齊 朱長(zhǎng)波①李俊偉 李 婷 陳素文 頡曉勇 張 博

    (農(nóng)業(yè)農(nóng)村部南海漁業(yè)資源開(kāi)發(fā)利用重點(diǎn)實(shí)驗(yàn)室 中國(guó)水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所 廣東省漁業(yè)生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室 廣州 510300)

    流沙灣是中國(guó)海水珍珠“南珠”的主產(chǎn)區(qū)和廣東省重要的貝類(lèi)養(yǎng)殖區(qū)。為評(píng)估灣內(nèi)養(yǎng)殖活動(dòng)的環(huán)境效應(yīng),于2015~2016年對(duì)流沙灣海區(qū)進(jìn)行了夏(8月)、秋(11月)、冬(2月)、春(5月) 4個(gè)季節(jié)的浮游植物和海水理化因子調(diào)查。共檢出浮游植物171種,包括硅藻門(mén)43屬122種、甲藻門(mén) 10屬44種、藍(lán)藻門(mén)2屬2種、金藻門(mén)2屬2種和裸藻門(mén)1屬1種。流沙灣內(nèi)灣浮游植物細(xì)胞豐度為(0.05~79.04)×104個(gè)/L,夏季>春季>秋季>冬季,且夏季豐度遠(yuǎn)大于其他三季,內(nèi)灣和外灣差異不顯著。春季須狀角毛藻()、紅海束毛藻()和明壁圓篩藻()為主要優(yōu)勢(shì)種,夏季優(yōu)勢(shì)種主要為中肋骨條藻(),秋季優(yōu)勢(shì)種主要為擬彎角毛藻()、奇異棍形藻()、洛氏角毛藻(),冬季主要優(yōu)勢(shì)種為威氏圓篩藻()、柔弱根管藻()、島脆桿藻()。其中,奇異棍形藻為春、秋、冬季的優(yōu)勢(shì)種,紅海束毛藻為春、夏、秋季的優(yōu)勢(shì)種。冗余分析表明,流沙灣浮游植物優(yōu)勢(shì)種在春季與水溫和氨氮密切相關(guān),夏季受透明度顯著影響,而在秋、冬季受亞硝酸氮含量影響顯著。流沙灣不同養(yǎng)殖區(qū)的浮游植物多態(tài)性和豐度有明顯季節(jié)差異。與2012年相比,流沙灣外灣的養(yǎng)殖覆蓋率提高了近50%,內(nèi)灣的珍珠貝養(yǎng)殖減少了近90%,養(yǎng)殖品種、規(guī)模和分布格局都發(fā)生了明顯變化,目前流沙灣的浮游植物群落特征正是對(duì)其變化的一種響應(yīng)。夏季魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖提高了水域營(yíng)養(yǎng)鹽水平并降低了浮游植物多樣性。大規(guī)模的扇貝養(yǎng)殖則導(dǎo)致了扇貝養(yǎng)殖區(qū)浮游植物豐度的降低。

    流沙灣;浮游植物;聚類(lèi)分析;養(yǎng)殖活動(dòng)

    水產(chǎn)養(yǎng)殖業(yè)屬環(huán)境依賴型產(chǎn)業(yè),水環(huán)境狀況與養(yǎng)殖品種的生長(zhǎng)繁育息息相關(guān)。養(yǎng)殖活動(dòng)也會(huì)影響 周邊水域環(huán)境和生態(tài)系統(tǒng),從而引發(fā)新的生態(tài)效應(yīng)。養(yǎng)殖活動(dòng)與環(huán)境之間的相互作用具有多樣性和復(fù)雜性,不同的養(yǎng)殖模式和養(yǎng)殖種類(lèi)的環(huán)境效應(yīng)也存在較大差異。高密度養(yǎng)殖設(shè)施會(huì)影響海水流動(dòng),造成水文情況的改變。同時(shí),養(yǎng)殖過(guò)程中,投藥、施肥等養(yǎng)殖管理活動(dòng)也會(huì)造成一定的環(huán)境壓力。在魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖投飼過(guò)程中產(chǎn)生的大量殘餌和糞便會(huì)給環(huán)境帶來(lái)大量的營(yíng)養(yǎng)物質(zhì),過(guò)高的養(yǎng)殖密度易造成水體的富營(yíng)養(yǎng)化和污染,進(jìn)而改變周?chē)h(huán)境的生物多樣性,甚至造成養(yǎng)殖魚(yú)類(lèi)的死亡(Grigorakis, 2011; Edwards, 2015)。濾食性貝類(lèi)以有機(jī)碎屑和浮游生物為食,規(guī)?;B(yǎng)殖可引起養(yǎng)殖區(qū)懸浮顆粒物的衰減,影響浮游生物動(dòng)力學(xué),并加速顆粒有機(jī)物的沉積 (張繼紅等, 2009)。相對(duì)于魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖,貝類(lèi)養(yǎng)殖對(duì)浮游生物具有調(diào)控作用,一定程度上可以修復(fù)或延緩水域的富營(yíng)養(yǎng)化(張繼紅等, 2013; 唐啟升等, 2016; 蔣增杰等, 2012)。

    流沙灣(22.36°~20.50°N, 109.80°~110.02°E)位于廣東省湛江市雷州半島西南部,為西北向呈葫蘆形的半封閉海灣,面積約為69 km2,是中國(guó)最大的海水珍珠育苗與養(yǎng)殖基地,海水珍珠產(chǎn)量占全國(guó)的60%以上。自2005年以來(lái),流沙灣周邊養(yǎng)殖戶們先后引入了墨西哥灣扇貝()、卵形鯧鲹()、華貴櫛孔扇貝()等養(yǎng)殖品種,至此,流沙灣形成了多品種與模式共存的養(yǎng)殖格局。隨著扇貝等新興品種養(yǎng)殖生產(chǎn)的擴(kuò)大,流沙灣的養(yǎng)殖區(qū)域和養(yǎng)殖布局也隨之發(fā)生改變。2015年流沙灣的養(yǎng)殖面積已達(dá)到16 km2,占其總面積的近1/4。與2012年相比,扇貝養(yǎng)殖面積由1132 hm2增加到1400 hm2,珍珠貝養(yǎng)殖面積由272 hm2下降至約100 hm2,網(wǎng)箱養(yǎng)殖面積約為95 hm2,其中,新增深水網(wǎng)箱約27 hm2。

    流沙灣是粵西典型的養(yǎng)殖海灣,自2010年,灣內(nèi)養(yǎng)殖品種和布局一直在不斷調(diào)整中,養(yǎng)殖格局的變化必然引起水域環(huán)境的響應(yīng)。多元的養(yǎng)殖品種和養(yǎng)殖方式共同作用,及其與環(huán)境因素相互作用,可能對(duì)流沙灣的浮游植物的群落結(jié)構(gòu)特征產(chǎn)生潛在影響。本研究針對(duì)流沙灣水域開(kāi)展浮游植物的調(diào)查,并對(duì)珍珠貝、扇貝和魚(yú)類(lèi)網(wǎng)箱等不同養(yǎng)殖類(lèi)型海域的浮游植物群落結(jié)構(gòu)和環(huán)境因子進(jìn)行對(duì)比分析,探討?zhàn)B殖活動(dòng)對(duì)流沙灣浮游植物群落的影響,進(jìn)而為流沙灣養(yǎng)殖結(jié)構(gòu)優(yōu)化和生態(tài)環(huán)境評(píng)價(jià)提供一定的參考。

    1 材料與方法

    流沙灣海域養(yǎng)殖品種呈多元化,不同品種的養(yǎng)殖方式與養(yǎng)殖周期差異較大,主要養(yǎng)殖品種及養(yǎng)殖周期見(jiàn)表1。

    表1 流沙灣水域養(yǎng)殖情況

    1.1 采樣站位點(diǎn)的設(shè)定

    分別于2015年8月(夏季)、11月(秋季)和2016年2月(冬季)、5月(春季)在流沙灣內(nèi)港、外港布設(shè) 12個(gè)站位,進(jìn)行4個(gè)航次的采樣調(diào)查。采樣站位如圖1所示。其中,扇貝養(yǎng)殖區(qū)有1、2、3和6號(hào)站位,珠母貝養(yǎng)殖區(qū)有7號(hào)站。8號(hào)站位位于灣口魚(yú)類(lèi)網(wǎng)箱 養(yǎng)殖區(qū),4和9號(hào)站位分別位于內(nèi)、外港主航道, 5號(hào)站位為外海邊界,10和11號(hào)站位為內(nèi)灣采樣點(diǎn),12號(hào)站位為外灣近岸采樣點(diǎn)。

    1.2 浮游植物的采集與分析

    浮游植物樣品的采集與定量參照《海洋調(diào)查規(guī)范》(GB17378.7-2007)。定性樣品采用淺水Ⅲ型浮游生物網(wǎng),自距底2 m深度向海面垂直拖曳,定量樣品采用有機(jī)玻璃采水器收集表層及底層的水樣(各1 L)進(jìn)行混合,現(xiàn)場(chǎng)用魯戈氏液固定后帶回實(shí)驗(yàn)室,用浮游生物計(jì)數(shù)框鏡檢、計(jì)數(shù)。

    圖1 流沙灣養(yǎng)殖布局及采樣站位圖

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

    海水溶解氧(DO)、水溫(T)、鹽度(S)和pH采用便攜式水質(zhì)分析儀(YSI-556型)現(xiàn)場(chǎng)測(cè)定,透明度(SD)采用薩式盤(pán)法測(cè)定。水樣中,亞硝酸鹽(NO2-N)采用萘乙二胺分光光度法,氨氮(NH4-N)采用靛酚藍(lán)法,硝酸鹽(NO3-N)采用鋅(Zn)-鎘(Cd)還原法,活性磷酸鹽(PO4-P)采用磷鉬藍(lán)分光光度法測(cè)定,總氮(TN)和總磷(TP)采用過(guò)硫酸鉀氧化法測(cè)定。葉綠素a(Chl-)顆粒有機(jī)物(POM)的測(cè)定方法分別為丙酮萃取法和灼燒稱重法。以上方法均按照《海洋監(jiān)測(cè)規(guī)范》(GB17378.4-2007)指定方法進(jìn)行。

    1.4 多樣性分析方法

    浮游植物的多樣性、均勻度分別采用Shannon- Weaver多樣性指數(shù)(′)和Pielou均勻度指數(shù)(),計(jì)算公式如下:

    式中,為物種數(shù),P為第物種在全部采樣中的比例,n為第種的個(gè)體數(shù),為所有種的個(gè)體總數(shù),為樣品種類(lèi)總數(shù),f為第種在各站點(diǎn)的監(jiān)測(cè)頻率?!?.02視為優(yōu)勢(shì)種。

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

    采用PRIMER 6.0軟件對(duì)浮游植物豐度進(jìn)行4次方根轉(zhuǎn)換,建立Bray-Curtis相似性聚類(lèi)分析,并進(jìn)行非度量多維尺度排序(nMDS)。用Canoco4.5軟件分析環(huán)境因子對(duì)浮游植物物種分布的影響,相關(guān)數(shù)據(jù)除pH外,經(jīng)過(guò)lg(+1)轉(zhuǎn)換。本研究將各季出現(xiàn)的優(yōu)勢(shì)種作為排序物種,對(duì)排序物種數(shù)據(jù)進(jìn)行除趨勢(shì)對(duì)應(yīng)分析(DCA),發(fā)現(xiàn)4個(gè)軸中梯度最大值<3,故選擇RDA進(jìn)行物種?環(huán)境相關(guān)性分析。采樣站位平面圖用ArcMap 10.2軟件繪制。

    2 結(jié)果

    2.1 浮游植物種類(lèi)組成及優(yōu)勢(shì)種的時(shí)空變化特征

    2015~2016年四季流沙灣浮游植物的種類(lèi)組成見(jiàn)表2。從表2可以看出,夏、秋、冬、春4個(gè)航次共鑒定出浮游植物5門(mén)58屬,共171種(含變種)。其中,硅藻42屬122種,占總種數(shù)的71.35%;甲藻10屬44種,占總種數(shù)的25.73%;金藻2屬2種;藍(lán)藻 2屬2種;裸藻1屬1種。在冬季,硅藻種數(shù)最少,而甲藻種數(shù)最多。流沙灣浮游植物優(yōu)勢(shì)種見(jiàn)表3。如表3所示,調(diào)查期間發(fā)現(xiàn)浮游植物優(yōu)勢(shì)種共16種,秋季最多,夏季最少。其中,奇異棍形藻為春、秋、冬季優(yōu)勢(shì)種,紅海束毛藻為春、夏、秋季優(yōu)勢(shì)種。在夏季,中肋骨條藻在各采樣站位占據(jù)絕對(duì)優(yōu)勢(shì),其豐度占夏季浮游植物總豐度的69%。

    2.2 浮游植物豐度及群落結(jié)構(gòu)分析

    在調(diào)查期間,流沙灣海域各采樣點(diǎn)的浮游植物豐度相差較大。春季,浮游植物豐度的變化范圍為0.05′104~13.15′104/L。夏季各采樣點(diǎn)游植物豐度較高,其變化范圍為2.02′104~79.04′104個(gè)/L。秋季浮游植物豐度變化范圍為0.11′104~7.45′104個(gè)/L,冬季浮游植物豐度較低,范圍為0.15′104~4.53′104個(gè)/L。流沙灣四季平均豐度分別為0.10′104、0.12′104、0.16′104、2.52′104個(gè)/L,季節(jié)變化趨勢(shì)為夏>春>秋>冬,且夏季內(nèi)、外灣豐度相差不大。甲藻的豐度在冬季最高,占冬季總豐度的6.3%,在其余季節(jié)不超過(guò)1%。藍(lán)藻主要為紅海束毛藻,其豐度在春季最高,占總豐度的35.8%。硅藻豐度在春、冬、秋和夏各占總豐度的63.6%、81.2%、82.3%和92.7%。

    表2 2015~2016年四季流沙灣浮游植物的種類(lèi)組成

    Tab.2 Seasonal phytoplankton composition in Liusha Bay during 2015~2016

    表3 流沙灣浮游植物優(yōu)勢(shì)種

    Tab.3 Dominant species of phytoplankton in Liusha Bay

    流沙灣全年浮游植物nMDS排序結(jié)果如圖2所示,不同季節(jié)間浮游植物群落組成差異顯著。因此,可將浮游植物群落大致按照季節(jié)劃分為4個(gè)大類(lèi)群,而各季節(jié)不同站位的群落組成相似度有一定差異。分別對(duì)春、夏、秋和冬4個(gè)航次的浮游群落進(jìn)行Bray- Curtis聚類(lèi)分析,結(jié)果如圖3所示,在50%的相似度水平上,春季各站位被分為5個(gè)類(lèi)群,分別為處于航道及附近的1、2、8和9號(hào)站位,位于內(nèi)灣的10和11號(hào)站位,位于扇貝養(yǎng)殖區(qū)的3號(hào)站位,位于珠母貝養(yǎng)殖區(qū)的7號(hào)站位和位于外灣的其他站位。夏季各站位可以分為4個(gè)類(lèi)群:位于魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖區(qū)及其附近的8和9號(hào)站位、位于內(nèi)灣的10號(hào)站位、位于內(nèi)灣近岸的11號(hào)站位和位于外灣的其他站位。在秋季,各站位也被分為4個(gè)類(lèi)群:位于魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖區(qū)及其附近的8號(hào)站位、位于流沙灣分界點(diǎn)的5號(hào)站位、位于內(nèi)灣的11和10號(hào)站位和位于外灣的其他站位。冬季各站位可以分為3個(gè)類(lèi)群:位于內(nèi)灣的10和11號(hào)站位和位于航道以北的2、3和12號(hào)站位及其他站位。

    2.3 多樣性指數(shù)與均勻度的時(shí)空分布

    流沙灣四季多樣性指數(shù)(¢)和均勻度指數(shù)()冬季最高、夏季最低(圖4)。春季,各站位的浮游植物多樣性指數(shù)范圍為4.21~2.22,其中,2號(hào)站位最高,為4.21,3、9、1和7號(hào)站位較低,分別為3.10、3.09、2.73和2.22。夏季,各站位的浮游植物多樣性指數(shù)范圍為2.66~0.23,其中,2號(hào)站位最高,為2.66,3、6、9、11、8和10號(hào)站位¢均不超過(guò)1。秋季,各站位的浮游植物多樣性指數(shù)范圍為3.60~2.02,其中,8和9號(hào)站位最高,為3.60和3.59,5和3號(hào)站位較低,分別為2.12和2.02。冬季各站位¢相差不大,范圍為3.15~2.50,其中,9號(hào)站位最高、4號(hào)站位最低。而分別為春季0.469~0.879、夏季0.522~0.048、秋季0.440~0.749、冬季0.520~0.706。

    2.4 不同季節(jié)浮游植物群落結(jié)構(gòu)的影響因素分析

    調(diào)查期間,流沙灣4個(gè)季節(jié)環(huán)境因子的平均值見(jiàn)表4。對(duì)優(yōu)勢(shì)種與環(huán)境因子間進(jìn)行冗余分析(Redundancy analysis, RDA)及Monte Carlo置換檢驗(yàn)。如圖5所示,檢驗(yàn)結(jié)果顯示,春季溫度(=4.05,=0.022)和氨氮(=3.03,=0.05)、夏季透明度(=58.34,=0.002)、秋季亞硝酸氮(=4.55,=0.006)、冬季亞硝酸氮(=3.91,=0.022)與流沙灣浮游植物優(yōu)勢(shì)種密切相關(guān)并影響顯著。透明度是影響夏季流沙灣浮游植物優(yōu)勢(shì)種的關(guān)鍵因子,對(duì)浮游植物物種變化的解釋度高達(dá)85%;亞硝酸氮對(duì)秋、冬季浮游植物物種變化的解釋程度分別為34%和28%。春季溫度和氨氮對(duì)浮游植物物種變化的解釋度分別為25%和15%。

    圖2 浮游植物群落的多維尺度分析

    圖3 流沙灣12個(gè)采樣點(diǎn)4個(gè)季節(jié)的聚類(lèi)

    圖4 流沙灣浮游植物多樣性指數(shù)(a)和均勻度指數(shù)( b)

    2.5 流沙灣不同養(yǎng)殖區(qū)浮游植物豐度、多樣性指數(shù)及可溶性營(yíng)養(yǎng)鹽含量

    比較了外灣對(duì)照區(qū)(4和5號(hào)站位)、扇貝主養(yǎng)區(qū)(1、3和6號(hào)站位)、珍珠貝養(yǎng)殖區(qū)(7號(hào)站位)、網(wǎng)箱魚(yú)類(lèi)養(yǎng)殖區(qū)(8號(hào)站位)和內(nèi)灣區(qū)(9和10號(hào)站位)的浮游植物豐度和多樣性指數(shù)(圖6)。結(jié)果顯示,在浮游植物豐度方面,魚(yú)類(lèi)小網(wǎng)箱養(yǎng)殖區(qū)浮游植物豐度在秋季和冬季高于其他區(qū)域。扇貝養(yǎng)殖區(qū)的浮游植物豐度較低,在夏、秋、冬三季低于外灣對(duì)照區(qū),除秋季外,低于珍珠貝養(yǎng)殖區(qū)。除夏季外,內(nèi)灣區(qū)的浮游植物豐度高于外灣對(duì)照區(qū)。在多樣性方面,春、夏兩季貝類(lèi)養(yǎng)殖區(qū)浮游植物多樣性小于外海對(duì)照區(qū)。在夏季,魚(yú)類(lèi)小網(wǎng)箱養(yǎng)殖區(qū)多樣性最低。秋季魚(yú)類(lèi)養(yǎng)殖區(qū)與內(nèi)灣區(qū)多樣性高于外灣。冬季各區(qū)多樣性指數(shù)均接近。

    外灣對(duì)照區(qū)、扇貝主養(yǎng)區(qū)、珍珠貝養(yǎng)殖區(qū)、網(wǎng)箱魚(yú)類(lèi)養(yǎng)殖區(qū)和內(nèi)灣區(qū)的可溶性無(wú)機(jī)氮(DIN)和磷酸鹽(PO4-P)含量見(jiàn)圖7。從圖7可以看出,流沙灣內(nèi)灣區(qū)和魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖區(qū)的DIN含量在春季低于其他區(qū)域,而在夏、秋季則高于其他區(qū)域。秋季各區(qū)的DIN含量相差不大。冬季魚(yú)類(lèi)網(wǎng)箱區(qū)的DIN為各區(qū)域之間最低。同時(shí),流沙灣外海對(duì)照區(qū)全年的氨氮濃度普遍低于其他海區(qū)。在PO4-P方面,除冬季外,魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖區(qū)和內(nèi)灣區(qū)的磷含量高于其他三區(qū)。與此同時(shí),外海對(duì)照區(qū)的PO4-P濃度則普遍低于其他海區(qū)。

    表4 流沙灣主要水質(zhì)參數(shù)

    Tab.4 Main parameters of water quality in Liusha Bay

    圖5 流沙灣浮游植物與環(huán)境因子冗余分析排序

    圖6 流沙灣不同區(qū)域浮游植物浮游植物豐度及多樣性指數(shù)

    圖7 流沙灣不同養(yǎng)殖區(qū)可溶性無(wú)機(jī)氮和磷酸鹽含量

    1:外海對(duì)照區(qū);2:扇貝養(yǎng)殖區(qū);3:珍珠貝養(yǎng)殖區(qū);4:魚(yú)類(lèi)小網(wǎng)箱養(yǎng)殖區(qū);5:內(nèi)灣區(qū);C:春季;X:夏季;Q:秋季;D:冬季

    1: No-culturing area in outer bay; 2: Scallop culturing area; 3: Pearl oyster culturing area; 4: Fish culturing area; 5: Inner bay; C: Spring; X: Summer; Q: Autumn; D: Winter

    3 討論

    3.1 流沙灣浮游植物群落結(jié)構(gòu)的季節(jié)變化特征及其影響因素分析

    浮游植物在地球生態(tài)系統(tǒng)中起著重要的功能和作用。浮游植物作為海洋食物網(wǎng)的基礎(chǔ),提供約占全球總產(chǎn)出一半的初級(jí)生產(chǎn)量,在O2的制造、CO2的固定和儲(chǔ)存、海洋營(yíng)養(yǎng)物質(zhì)循環(huán)等生態(tài)過(guò)程中有著不可或缺的貢獻(xiàn)(Not, 2012)。浮游植物群落結(jié)構(gòu)對(duì)生態(tài)系統(tǒng)中環(huán)境因子的改變有直接的響應(yīng),其時(shí)空變化特征與水文、溫度、營(yíng)養(yǎng)元素等環(huán)境因子關(guān)系密切(王符菁等, 2015)。Bode等(2017)對(duì)拉科努尼亞沿海地區(qū)的浮游植物多樣性進(jìn)行調(diào)查,發(fā)現(xiàn)河水徑流和海灣上升流共同影響浮游植物的群落結(jié)構(gòu)組成,同時(shí),浮游植物的群落結(jié)構(gòu)也受到季節(jié)的影響。Wu等(2014)對(duì)大亞灣海區(qū)浮游植物群落與環(huán)境因子的關(guān)系進(jìn)行了分析,結(jié)果顯示,環(huán)境溫度和營(yíng)養(yǎng)等因素對(duì)浮游植物的季節(jié)演替率有很大的影響。

    在本次調(diào)查中,流沙灣海域浮游植物群落的細(xì)胞豐度夏季最高,平均豐度為2.52×104個(gè)/L,遠(yuǎn)高于其他3個(gè)季節(jié)。春季、秋季和冬季從高到低依次為0.16×104、0.12×104、0.1×104個(gè)/L。各季節(jié)浮游植物群落的主要優(yōu)勢(shì)種出現(xiàn)了明顯的更替,夏季主要優(yōu)勢(shì)種為中肋骨條藻,占總豐度的69.5%,秋季優(yōu)勢(shì)種較多,但占總豐度的比例不高,其中,主要優(yōu)勢(shì)種擬彎角毛藻、奇異棍形藻、洛氏角毛藻的豐度之和僅占總豐度的43%。冬季主要優(yōu)勢(shì)種為威氏圓篩藻、柔弱根管藻、島脆桿藻,三者占總豐度的60%。春季以須狀角毛藻、紅海束毛藻和明壁圓篩藻為主要優(yōu)勢(shì)種,3種微藻的豐度之和共占流沙灣浮游植物總豐度的67%。多維尺度排序的脅強(qiáng)值為0.15,表明本次排序結(jié)果可以較好地解釋物種的群落結(jié)構(gòu)分組情況(Manly, 2005)。結(jié)果顯示,調(diào)查站位在相同季節(jié)浮游植物群落組成的相似度較高,即浮游植物群落結(jié)構(gòu)存在明顯的季節(jié)差異。

    在流沙灣海域,奇異棍形藻為秋、冬、春季優(yōu)勢(shì)種,紅海束毛藻為夏、秋、春季優(yōu)勢(shì)種。其中,奇異棍形藻為廣溫沿岸種(王符菁等, 2015),其在夏季之后在流沙灣海域廣泛存在。而紅海束毛藻對(duì)溫度的要求較高(丁昌玲等, 2016),在水溫較低的冬季只存在于內(nèi)灣近岸的11號(hào)站位,未成為冬季優(yōu)勢(shì)種。流沙灣冬季平均水溫為18.2℃,期間出現(xiàn)的主要優(yōu)勢(shì)種如威氏圓篩藻、柔弱根管藻等也屬于不耐高溫的溫帶近岸或廣布種。在春、秋兩季,出現(xiàn)較多的主要優(yōu)勢(shì)種如須狀角毛藻、擬彎角毛藻等則屬于偏暖性近岸種。在水溫最高的夏季,耐高溫的中肋骨條藻和喜高溫的紅海束毛藻成為流沙灣海區(qū)的主要優(yōu)勢(shì)種。表明環(huán)境溫度對(duì)流沙灣浮游植物的季節(jié)演替率有很大的影響。喬倩等(2016)發(fā)現(xiàn),中肋骨條藻可以有效利用多種有機(jī)氮源,對(duì)無(wú)機(jī)氮具有較低的親和能力。夏季流沙灣養(yǎng)殖活動(dòng)頻繁,為海區(qū)帶來(lái)較多的有機(jī)氮源,同時(shí),DIN含量在4個(gè)季節(jié)中最低,也為中肋骨條藻的大量繁殖提供了有利條件。

    3.2 規(guī)?;KB(yǎng)殖對(duì)浮游植物的影響

    多樣性指數(shù)和均勻度是反映水生態(tài)穩(wěn)定性的重要指標(biāo),也是評(píng)價(jià)群落是否成熟和穩(wěn)定的特征之一(陳露等, 2015)。在夏季各采樣點(diǎn)中,8號(hào)站的多樣性指數(shù)為0.29,在春、秋、冬季均大于3,分別為3.71、3.59、3.01,顯示出良好的多態(tài)性,表明魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖在夏季對(duì)浮游植物多態(tài)性影響顯著。夏季,養(yǎng)殖魚(yú)類(lèi)生長(zhǎng)迅速,因而也是網(wǎng)箱養(yǎng)殖活動(dòng)的高峰期,傳統(tǒng)養(yǎng)殖網(wǎng)箱密集分布通常會(huì)造成水體富營(yíng)養(yǎng)化。周志誠(chéng)等(2014)發(fā)現(xiàn),流沙灣網(wǎng)箱養(yǎng)殖區(qū)水體中氮、磷濃度自5~9月逐漸增加,10月開(kāi)始下降。認(rèn)為這是網(wǎng)箱養(yǎng)魚(yú)投餌、養(yǎng)殖生物體集中排泄以及浮游生物活動(dòng)共同作用的結(jié)果。

    高密度的濾食性貝類(lèi)養(yǎng)殖會(huì)降低浮游植物的多樣性。張玲等(2015)對(duì)大鵬澳牡蠣養(yǎng)殖區(qū)及其鄰近海域進(jìn)行了采樣調(diào)查,研究表明,牡蠣養(yǎng)殖對(duì)浮游植物群落結(jié)構(gòu)造成了一定影響,能夠顯著降低浮游植物豐度、種類(lèi)多樣性和均勻度。董婧等(2003)在對(duì)小窯灣的牡蠣、扇貝養(yǎng)殖區(qū)調(diào)查中發(fā)現(xiàn),在高密度貝類(lèi)筏式養(yǎng)殖區(qū),隨著貝類(lèi)攝食強(qiáng)度的增強(qiáng),浮游植物數(shù)量、種類(lèi)都大量減少,浮游植物群落多樣性指數(shù)偏低,豐度指數(shù)偏低,均勻度指數(shù)因季節(jié)不同而有較大差異。本次調(diào)查中,貝類(lèi)養(yǎng)殖區(qū)平均多樣性指數(shù)在春、夏兩季低于外灣對(duì)照區(qū),在秋季則低于內(nèi)灣區(qū)。表明,流沙灣貝類(lèi)養(yǎng)殖對(duì)浮游植物多態(tài)性產(chǎn)生了一定的影響。本次調(diào)查還發(fā)現(xiàn),位于扇貝養(yǎng)殖區(qū)的3和6號(hào)站位在夏、秋季,多樣性指數(shù)均明顯低于其他站位。同時(shí),1和2號(hào)站位則表現(xiàn)出很高的多態(tài)性,表明在水流通暢的灣口,扇貝養(yǎng)殖活動(dòng)對(duì)浮游植物多樣性影響不大。

    浮游植物的生長(zhǎng)同時(shí)受到海水營(yíng)養(yǎng)鹽的上行效應(yīng)以及浮游動(dòng)物濾食性雙殼類(lèi)及魚(yú)類(lèi)下行效應(yīng)的影響。魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖產(chǎn)生的殘餌和糞便會(huì)給海水環(huán)境帶來(lái)大量的有機(jī)營(yíng)養(yǎng)物,而貝類(lèi)主要以NH4-N和無(wú)機(jī)磷酸鹽的方式向水中排泄N或P。魚(yú)、貝養(yǎng)殖可以增加海水中營(yíng)養(yǎng)鹽的含量,促進(jìn)藻類(lèi)的生長(zhǎng)。研究表明,海灣扇貝()和合浦珠母貝的養(yǎng)殖活動(dòng)會(huì)影響海區(qū)DO、COD和P、TN的變化(徐東會(huì)等, 2017; 王增煥等, 2011)。流沙灣外海對(duì)照區(qū)無(wú)養(yǎng)殖活動(dòng),全年的氨氮濃度普遍低于養(yǎng)殖區(qū),同時(shí),磷酸鹽含量往往也低于養(yǎng)殖區(qū)。僅在貝類(lèi)代謝水平降低且無(wú)魚(yú)類(lèi)小網(wǎng)箱養(yǎng)殖活動(dòng)的冬季,外海對(duì)照區(qū)磷酸鹽含量高于其他區(qū)域。同時(shí),在魚(yú)類(lèi)養(yǎng)殖活動(dòng)旺盛的夏季,小網(wǎng)箱區(qū)的DIN和磷酸鹽含量明顯高于其他海區(qū),而在無(wú)養(yǎng)殖活動(dòng)的冬季則低于其他海區(qū)。表明流沙灣海區(qū)的多種養(yǎng)殖活動(dòng)在養(yǎng)殖季節(jié)對(duì)水質(zhì)影響較大。

    濾食性的扇貝和珍珠貝等貝類(lèi)能夠給浮游植物帶來(lái)攝食壓力。流沙灣扇貝養(yǎng)殖區(qū)浮游植物豐度,除春季高于外海對(duì)照區(qū)外,全年都低于其他海區(qū)。而在無(wú)貝類(lèi)攝食壓力的魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖區(qū)和內(nèi)灣區(qū),浮游植物豐度則高于外海對(duì)照區(qū)。而在珍珠貝區(qū),浮游植物豐度高于扇貝區(qū),可能是因?yàn)檎渲樨悈^(qū)養(yǎng)殖規(guī)模小,對(duì)浮游植物帶來(lái)的攝食壓力較低的緣故。此外,不同貝類(lèi)對(duì)浮游植物的選擇性攝食也可能會(huì)引起浮游植物群落結(jié)構(gòu)的改變。櫛孔扇貝()主動(dòng)選擇攝食甲藻和定鞭金藻(),對(duì)骨條藻的選擇性不高。在中肋骨條藻為優(yōu)勢(shì)種的夏季,櫛孔扇貝胃含物最高的品種為海鏈藻()(喬芮, 2015)。而中肋骨條藻可以很好的作為合浦珠母貝幼貝和珠母貝()的餌料(Kuwatani, 1964; Martínez-Fernández, 2007)。夏季,流沙灣合浦珠母貝對(duì)中肋骨條藻的攝食作用可以在一定程度上抑制中肋骨條藻的優(yōu)勢(shì),導(dǎo)致夏季合浦珠母貝養(yǎng)殖區(qū)多樣性指數(shù)高于扇貝養(yǎng)殖區(qū)。

    3.3 養(yǎng)殖布局對(duì)浮游植物的影響

    隨著水產(chǎn)養(yǎng)殖活動(dòng)規(guī)模的擴(kuò)大,其對(duì)近岸環(huán)境的影響和污染問(wèn)題也越來(lái)越突出(鮑旭騰等, 2012)。傳統(tǒng)網(wǎng)箱養(yǎng)殖廣泛投喂人工飼料和天然雜魚(yú)等餌料,大量的殘餌是海域的污染物質(zhì)來(lái)源之一,易造成水體營(yíng)養(yǎng)指數(shù)偏高、底部沉積物堆積、養(yǎng)殖水體底層貧氧現(xiàn)象嚴(yán)重等老化待征。而貝類(lèi)筏式養(yǎng)殖可以阻擋海水流動(dòng),降低養(yǎng)殖海區(qū)水體的流速、改變水流方向,使得有機(jī)物質(zhì)和營(yíng)養(yǎng)鹽被截留(苗衛(wèi)衛(wèi)等, 2007)。同時(shí),即使在營(yíng)養(yǎng)鹽、水溫等富營(yíng)養(yǎng)化條件滿足的情況下,流量、流速和水位等水動(dòng)力條件也可能成為水體富營(yíng)養(yǎng)化的制約因子(劉洪波等, 2013)。許棟等(2017)對(duì)南渡江下游進(jìn)行二維水動(dòng)力及水生態(tài)數(shù)值模擬,結(jié)果顯示,枯水期水質(zhì)富營(yíng)養(yǎng)化風(fēng)險(xiǎn)遠(yuǎn)高于豐水期。

    與往年調(diào)查相比,流沙灣的浮游植物種數(shù)和硅藻種數(shù)有所恢復(fù),與2008年持平(表5)。然而,豐度和多樣性指數(shù)(¢)在夏季卻出現(xiàn)了較大的波動(dòng)。吳雪等(2017)認(rèn)為,一般¢<1時(shí),表明浮游植物群落可能受到其他環(huán)境因素的干擾;均勻度>0.3,則表明海區(qū)內(nèi)浮游植物多樣性較好。即以¢>1且>0.3作為浮游植物多樣性較好的評(píng)價(jià)標(biāo)準(zhǔn)。在本次調(diào)查中,夏季內(nèi)灣8~11號(hào)站位、外灣扇貝養(yǎng)殖區(qū)3、6號(hào)站位的生物多樣性指數(shù)和豐富度均低于標(biāo)準(zhǔn),表明該位點(diǎn)浮游植物多樣性差、結(jié)構(gòu)不穩(wěn)定。郭永堅(jiān)等(2015)在2012年對(duì)流沙灣的調(diào)查中發(fā)現(xiàn),魚(yú)貝養(yǎng)殖區(qū)的多樣性指數(shù)、均勻度指數(shù)較高,且網(wǎng)箱魚(yú)類(lèi)養(yǎng)殖區(qū)的多樣性指數(shù)一直為各區(qū)最高,與本次調(diào)查得出的結(jié)論完全相反。上述現(xiàn)象可能是因?yàn)榱魃碁仇B(yǎng)殖設(shè)施增多,對(duì)海水流動(dòng)性的阻礙增大,引起局部水體污染造成的。

    張靜等(2013)研究表明,流沙灣潮流運(yùn)動(dòng)為帶有旋轉(zhuǎn)流的往復(fù)流運(yùn)動(dòng),落潮流速大于漲潮流速,且余速較強(qiáng),海水自凈能力很強(qiáng)。羅昭林等(2014)的水動(dòng)力建模表明,流沙灣內(nèi)、外灣口水體流速差異顯著,退潮時(shí),內(nèi)灣口的水流速度趨于靜止?fàn)顟B(tài)(圖8)。2015年,外灣扇貝筏式養(yǎng)殖面積進(jìn)一步擴(kuò)大,勢(shì)必對(duì)流沙灣的水交換產(chǎn)生更大的影響。朱宜平等(2010)通過(guò)圍隔裝置研究3種流速對(duì)浮游動(dòng)植物生長(zhǎng)的影響,發(fā)現(xiàn)浮游植物生物量隨著流速增長(zhǎng)出現(xiàn)先升高后降低的趨勢(shì),并在總量上小于靜水環(huán)境。魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖活動(dòng)易造成水體的富營(yíng)養(yǎng)化,水動(dòng)力不足易造成內(nèi)灣透明度降低,營(yíng)養(yǎng)鹽水平升高,其與夏季高溫的共同作用導(dǎo)致耐高溫和鹽度適應(yīng)能力較強(qiáng)的中肋骨條藻在內(nèi)灣大面積增殖。此外,本次調(diào)查發(fā)現(xiàn),紅海束毛藻首次在流沙灣海區(qū)發(fā)現(xiàn),并在夏季普遍出現(xiàn)于灣外各采樣點(diǎn),秋季則遍布于內(nèi)外灣,冬季僅存在于水溫最高的11號(hào)站位,春季則出現(xiàn)在內(nèi)灣的10、11號(hào)站位和外灣的3、4號(hào)站位。紅海束毛藻體積小,生存溫度在18℃~32℃之間,在北部灣海域大量分布(Boatman, 2017; 丁昌玲等, 2016),幾乎不受貝類(lèi)攝食壓力的影響,其在流沙灣全年的時(shí)空分布出現(xiàn)了明顯的溫度相關(guān),也間接地表明,流沙灣水動(dòng)力受到養(yǎng)殖設(shè)施較大的干擾。

    與2012年相比,流沙灣無(wú)機(jī)磷和總顆粒有機(jī)物的含量有所上升,其中,顆粒有機(jī)物范圍由原先的2.64~6.18(郭永堅(jiān)等, 2015)上升到本次的5.01~10.22,表明隨著水動(dòng)力不足,部分海域可能出現(xiàn)了富營(yíng)養(yǎng)化趨勢(shì)。此外,本次調(diào)查中,溶解氧水平低于2012年的調(diào)查結(jié)果,表明流沙灣可能已出現(xiàn)養(yǎng)殖負(fù)荷超過(guò)系統(tǒng)承載力的狀況。水動(dòng)力不足的情況下,在養(yǎng)殖活動(dòng)強(qiáng)度大的5~11月,應(yīng)注意赤潮和水中氧含量過(guò)低情況的發(fā)生。

    表5 近年來(lái)流沙灣浮游植物群落特征比較

    Tab.5 Comparison of characteristics of phytoplankton community in Liusha Bay in recent years

    圖8 流沙灣漲/落潮流場(chǎng)(羅昭林等, 2014)

    Bao XT, Xu H, Zhang JH,. Best management practices for controlling aquaculture non-point pollution. South China Fisheries Science, 2012, 8(3): 79–86 [鮑旭騰, 徐皓, 張建華, 等. 水產(chǎn)養(yǎng)殖面源污染控制的最佳管理實(shí)踐. 南方水產(chǎn)科學(xué), 2012, 8(3): 79–86]

    Boatman TG, Lawson T, Geider RJ. A key marine diazotroph in a changing Ocean: The interacting effects of temperature, CO2and light on the growth ofIMS101. PLoS One, 2017, 12(1): e0168796

    Bode A, Varela M, Prego R,. The relative effects of upwelling and river flow on the phytoplankton diversity patterns in the ria of A Coru?a (NW Spain). Marine Biology, 2017, 164(4): 93

    Chen L, Li CH, Wei XL,. Influence of nitrogen and phosphorus enrichment on phytoplankton community structure in Nansha area in summer. South China Fisheries Science, 2015, 11(5): 56–66 [陳露, 李純厚, 魏小嵐, 等. 南沙海域夏季氮磷添加模擬實(shí)驗(yàn)中浮游植物群落的變化. 南方水產(chǎn)科學(xué), 2015, 11(5): 56–66]

    Ding CL, Sun J, Xue B. Winter and summer diazotrophic cyanobacteria in the western South China Sea. Haiyang Xuebao, 2016, 38(4): 84–94 [丁昌玲, 孫軍, 薛冰. 冬夏季南海西部固氮藍(lán)藻. 海洋學(xué)報(bào), 2016, 38(4): 84–94]

    Dong J, Bi YP, Wang WB,. Effect of high density bivalve rafts cultivation on phytoplankton community. Marine Fisheries Research, 2003, 24(3): 50–54 [董婧, 畢遠(yuǎn)溥, 王文波, 等. 小窯灣高密度貝類(lèi)筏式養(yǎng)殖對(duì)浮游植物群落的影響. 海洋水產(chǎn)研究, 2003, 24(3): 50–54]

    Edwards P. Aquaculture environment interactions: Past, present and likely future trends. Aquaculture, 2015, 447: 2–14

    Grigorakis K, Rigos G. Aquaculture effects on environmental and public welfare-The case of Mediterranean mariculture. Chemosphere, 2011, 85(6): 899–919

    Guo YJ, Luo ZL, Zhu CB,. Influence of aquaculture on characteristics of phytoplankton community in Liusha Bay. South China Fisheries Science, 2015, 11(2): 57–65 [郭永堅(jiān), 羅昭林, 朱長(zhǎng)波, 等. 水產(chǎn)養(yǎng)殖對(duì)流沙灣浮游植物群落特征的影響. 南方水產(chǎn)科學(xué), 2015, 11(2): 57–65]

    Jiang ZJ, Fang JG, Mao YZ,. Environmental effect of marine fish cage aquaculture and integrated multi-trophic aquaculture. Environmental Science and Management, 2012, 37(1): 120–124 [蔣增杰, 方建光, 毛玉澤, 等. 海水魚(yú)類(lèi)網(wǎng)箱養(yǎng)殖的環(huán)境效應(yīng)及多營(yíng)養(yǎng)層次的綜合養(yǎng)殖. 環(huán)境科學(xué)與管理, 2012, 37(1): 120–124]

    Kuwatani Y. Rearing of the pearl oyster,(Dunker), on the diet of. Nippon Suisan Gakkaishi, 1964, 30(2): 104–113

    Liu HB, Zhu ML, Wang JZ,Influence of reservoir hydrodynamic on eutrophication of water body. Journal of Water Resources and Water Engineering, 2013, 24(2): 19–21 [劉洪波, 朱夢(mèng)羚, 王精志, 等. 水庫(kù)水動(dòng)力對(duì)水體富營(yíng)養(yǎng)化影響. 水資源與水工程學(xué)報(bào), 2013, 24(2): 19–21]

    Luo ZL, Zhu CB, Guo YJ,. Distribution characteristics of C, N and P in Liusha Bay surface sediment and their pollution assessment. South China Fisheries Science, 2014, 10(3): 1–8 [羅昭林, 朱長(zhǎng)波, 郭永堅(jiān), 等. 流沙灣表層沉積物中碳、氮、磷的分布特征和污染評(píng)價(jià). 南方水產(chǎn)科學(xué), 2014, 10(3): 1–8]

    Manly BF. Multivariate statistical methods: A primer. Chapman and Hall, 2005: 1–224

    Martínez-Fernández E, Southgate PC. Use of tropical microalgae as food for larvae of the black-lip pearl oyster. Aquaculture, 2007, 263(1–4): 220–226

    Miao WW, Jiang M. Environmental impacts and sustainable development of aquaculture in China. Journal of Agro- Environment Science, 2007, 26(s1): 319–323 [苗衛(wèi)衛(wèi), 江敏. 我國(guó)水產(chǎn)養(yǎng)殖對(duì)環(huán)境的影響及其可持續(xù)發(fā)展. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2007, 26(s1): 319–323]

    Not F, Siano R, Kooistra WHCF,. Chapter one-diversity and ecology of eukaryotic marine phytoplankton. Advances in Botanical Research, 2012, 64: 1–53

    Qiao Q, Wang CH, Guo X. Effects of nitrogen sources on the growth of. Chinese Journal of Ecology, 2016, 35(8): 2110–2116 [喬倩, 王朝暉, 郭鑫. 不同氮源對(duì)中肋骨條藻()生長(zhǎng)的影響. 生態(tài)學(xué)雜志, 2016, 35(8): 2110–2116]

    Qiao R. Microscope and phytoplankton pigments analysis of the shellfish feeding habits. Master¢s Thesis of Shanghai Ocean University, 2015 [喬芮. 基于鏡檢和浮游植物色素分析的貝類(lèi)食性研究. 上海海洋大學(xué)碩士研究生學(xué)位論文, 2015]

    Tang QS, Han D, Mao YZ,. Species composition, non-fed rate and trophic level of Chinese aquaculture. Journal of Fishery Sciences of China, 2016, 23(4): 729–758 [唐啟升, 韓冬, 毛玉澤, 等. 中國(guó)水產(chǎn)養(yǎng)殖種類(lèi)組成、不投餌率和營(yíng)養(yǎng)級(jí). 中國(guó)水產(chǎn)科學(xué), 2016, 23(4): 729–758]

    Wang FJ, Lin YS, Cao WQ,. The relationship between nutrients and phytoplankton community structure in northern Beibu Gulf. Journal of Tropical Oceanography, 2015, 34(6): 73–85 [王符菁, 林元燒, 曹文清, 等. 北部灣北部浮游植物群落結(jié)構(gòu)及其與營(yíng)養(yǎng)鹽的關(guān)系. 熱帶海洋學(xué)報(bào), 2015, 34(6): 73–85]

    Wang ZH, Ke CL, Wang XN,.Marine environmental quality assessment of shellfish culture in Liusha Bay. South China Fisheries Science, 2011, 7(3): 24–30 [王增煥, 柯常亮, 王許諾, 等. 流沙灣貝類(lèi)養(yǎng)殖海域環(huán)境質(zhì)量評(píng)價(jià). 南方水產(chǎn)科學(xué), 2011, 7(3): 24–30]

    Wu ML, Wang YS, Wang YT,. Influence of environmental changes on phytoplankton pattern in Daya Bay, South China Sea. Revista De Biologia Marina Y Oceanografia, 2014, 49(2): 323–337

    Wu X, Li XL, Yang JL,. Ecological characteristics of phytoplankton in scallop culture area in Laizhou Bay. Fisheries Science, 2017, 36(3): 347–352 [吳雪, 李希磊, 楊俊麗, 等. 萊州灣扇貝養(yǎng)殖區(qū)浮游藻類(lèi)的生態(tài)特征. 水產(chǎn)科學(xué), 2017, 36(3): 347–352]

    Xu D, Wang D, Ji CN,. Numerical simulation of two-dimensional hydrodynamics and aquatic ecology of the lower reach of Nandujiang River in Hainan Island, China. Research of Environmental Sciences, 2017, 30(2): 214–223 [許棟, 王迪, 及春寧, 等. 海南南渡江下游二維水動(dòng)力及水生態(tài)數(shù)值模擬. 環(huán)境科學(xué)研究, 2017, 30(2): 214–223]

    Xu DH, Chen BJ, Jiang T,. Temporal and spatial variation characteristics of water nutrients and eutrophication evaluation in aquaculture water areas of Sanggou Bay. Progress in Fishery Sciences, 2017, 38(3): 1–11 [徐東會(huì), 陳碧鵑, 江濤, 等. 桑溝灣養(yǎng)殖海域營(yíng)養(yǎng)鹽時(shí)空分布特征及富營(yíng)養(yǎng)化評(píng)價(jià). 漁業(yè)科學(xué)進(jìn)展, 2017, 38(3): 1–11]

    Zhang CX, Chen HY, Sun SL,. Temporal and spatial distribution of phytoplankton in Liusha BayActa Ecologica Sinica, 2012, 32(5): 1527–1537 [張才學(xué), 陳慧妍, 孫省利, 等. 流沙灣海草床海域浮游植物的時(shí)空分布及其影響因素. 生態(tài)學(xué)報(bào), 2012, 32(5): 1527–1537]

    Zhang J, Liu Y, Ma ZH,. Three-dimensional numerical simulation of tide and tidal currents in Liusha Bay and the study on its sea environmental capacity. Marine Science Bulletin, 2013, 32(1): 34–44 [張靜, 劉玉, 馬志華, 等. 流沙灣海域潮汐潮流的三維數(shù)值模擬和海域環(huán)境容量研究. 海洋通報(bào), 2013, 32(1): 34–44]

    Zhang JH, Fang JG, Wang W. Progress in studies on ecological carrying capacity of mariculture for filter-feeding shellfish. Journal of Fishery Sciences of China, 2009, 16(4): 626–632 [張繼紅, 方建光, 王巍. 淺海養(yǎng)殖濾食性貝類(lèi)生態(tài)容量的研究進(jìn)展. 中國(guó)水產(chǎn)科學(xué), 2009, 16(4): 626–632]

    Zhang L, Li ZJ, Chen YF,. Impact of oyster aquaculture on phytoplankton population structure in Dapeng Cove. Oceanologia et Limnologia Sinica , 2015, 46(3): 549–555 [張玲, 李政菊, 陳飛羽, 等. 大鵬澳牡蠣養(yǎng)殖對(duì)浮游植物種群結(jié)構(gòu)的影響研究. 海洋與湖沼, 2015, 46(3): 549–555]

    Zhang L, Wu T, Xu D,Feeding on feces, diets and cage farming sediments of Pby Japanese scallop.Progress in Fishery Sciences, 2013, 34(2): 96–102 [張繼紅, 吳桃, 徐東, 等. 蝦夷扇貝對(duì)牙鲆的糞便、殘餌及網(wǎng)箱養(yǎng)殖區(qū)沉積物的攝食. 漁業(yè)科學(xué)進(jìn)展, 2013, 34(2): 96–102]

    Zhou ZC, Shen YC, Xue HB. Temporal and spatial variation of nitrogen and phosphorus and their load estimation in cage-culture waters in Liusha Bay, Guangdong. Journal of Guangdong Ocean University, 2014, 34(3): 65–69 [周志誠(chéng), 申玉春, 薛海波. 廣東流沙灣網(wǎng)箱養(yǎng)魚(yú)區(qū)水體氮磷的時(shí)空變化及負(fù)荷估算. 廣東海洋大學(xué)學(xué)報(bào), 2014, 34(3): 65–69]

    Zhu YP, Zhang HP, Li FP,. Enclosure experiments about the hydrodynamics effects on the plankton.Environmental Science, 2010, 31(1): 69–75 [朱宜平, 張海平, 李飛鵬, 等. 水動(dòng)力對(duì)浮游生物影響的圍隔研究. 環(huán)境科學(xué), 2010, 31(1): 69–75]

    Phytoplankton Community Characteristics in Different Seasons and Their Relationship with Aquaculture in Liusha Bay

    SU Jiaqi, ZHU Changbo①, LI Junwei, LI Ting, CHEN Suwen, XIE Xiaoyong, ZHANG Bo

    (Key Laboratory of South China Sea Fishery Resources Exploitation & Utilization, Ministry of Agriculture and Rural Affairs; South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangdong Provincial Key Laboratory of Fishery Ecology and Environment, Guangzhou 510300)

    Liusha Bay is one of the main aquaculture areas of marine mollusks in Guangdong Province, characterized by the main producing areas of South Pearl. To assess the environmental effects of aquaculture activities, the phytoplankton composition, abundance, community structure, diversity, and water physicochemical factors were investigated in each season from August 2015 to May 2016. In total, 171 phytoplankton species were detected, with 122 species of diatom belonging to 43 genera, followed by dinoflagellates represented by 44 taxa (10 genera), blue-green algae by two taxa (two genera), chrysophyceae by two taxa (two genera), and euglena by one taxa (one genus). The range of phytoplankton density was 0.05×104~79.04×104cells/L, which ranked as: summer > spring > autumn > winter. In general, the phytoplankton density in the outer bay was similar to that in the inner bay. The abundance of phytoplankton in summer was distinctly higher than that in the other three seasons. The most common dominant species were,,andin spring, Sin summer,,, andin autumn, andandin winter.was the dominant species in spring, autumn, and winter.was the dominant species in all seasons except for autumn. Redundancy analysis suggested that the main variables affecting the dominant species were water temperature and nitrate (NO3-N) in spring, transparency in summer, and nitrite (NO2-N) in autumn and winter. The phytoplankton diversity and abundance of different cultured zones presented seasonal differences. Compared to 2012, the total coverage of the aquaculture zone in 2015 increased by approximately 50% in the outer bay. Furthermore, coverage of the pearl oyster farming zone decreased by about 90% in the inner bay. Both the distribution of the aquaculture zone and the scale of the mariculture species were significantly changed, and the characteristics of the phytoplankton community responded to these changes. In summer, cage fish farming could efficiently increase nutrient content in seawater and decrease phytoplankton diversity. Furthermore, phytoplankton abundance was decreased by large-scale scallop culture.

    Liusha Bay; Phytoplankton; Cluster analysis; Aquaculture activity

    ZHU Changbo, E-mail: changbo@scsfri.ac.cn

    蘇家齊, 朱長(zhǎng)波, 李俊偉, 李婷, 陳素文, 頡曉勇, 張博. 流沙灣浮游植物群落特征季節(jié)變化及其與養(yǎng)殖活動(dòng)的關(guān)系. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(6): 11–23

    Su JQ, Zhu CG, Li JW, Li T, Chen SW, Xie XY, Zhang B. Phytoplankton community characteristics in different seasons and their relationship with aquaculture in Liusha Bay. Progress in Fishery Sciences, 2018, 39(6): 11–23

    * 廣東省海洋漁業(yè)科技與產(chǎn)業(yè)發(fā)展專(zhuān)項(xiàng)項(xiàng)目(Z2015013)和中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)南海水產(chǎn)研究所專(zhuān)項(xiàng)資金項(xiàng)目(2013ZD01)共同資助[This work was supported by Special Fund for Marine Fisheries Research and Extension of Guangdong Province (Z2015013), and Scientific Research Funds for Central Non-Profit Institutes, South China Sea Fisheries Research Institute (2013ZD01)]. 蘇家齊,E-mail: sujiaqi2005@126.com

    朱長(zhǎng)波,副研究員,E-mail: changbo@scsfri.ac.cn

    2017-10-24,

    2017-12-26

    10.19663/j.issn2095-9869.20171024002

    S932.7

    A

    2095-9869(2018)06-0011-13

    (編輯 馬璀艷)

    猜你喜歡
    沙灣養(yǎng)殖區(qū)浮游
    廣西北部灣近岸養(yǎng)殖區(qū)表層沉積物和養(yǎng)殖生物體重金屬污染狀況及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)*
    江蘇近海紫菜養(yǎng)殖區(qū)的空間擴(kuò)張模式研究
    浞河浮游藻類(lèi)的調(diào)查研究與水質(zhì)評(píng)價(jià)
    《浮游》
    流行色(2019年10期)2019-12-06 08:13:26
    番禺沙灣古鎮(zhèn)保護(hù)更新的分析與研究
    江西建材(2018年1期)2018-04-04 05:25:54
    浮游藻類(lèi)在水體PAHs富集與降解中的研究
    養(yǎng)殖活動(dòng)對(duì)超微型浮游生物分布影響的研究*
    海洋與湖沼(2014年6期)2014-03-19 12:17:34
    圍墾對(duì)椒江口夏季浮游植物群落結(jié)構(gòu)和多樣性的影響
    多源多相復(fù)雜儲(chǔ)集層預(yù)測(cè)——以春風(fēng)油田新近系沙灣組一段1砂組為例
    斷塊油氣田(2014年5期)2014-03-11 15:33:44
    小沙灣黃河特大橋針對(duì)高低塔的減震設(shè)計(jì)
    老司机影院成人| 亚洲伊人久久精品综合| 美女xxoo啪啪120秒动态图| 美女国产高潮福利片在线看| 亚洲av综合色区一区| 秋霞伦理黄片| 成人免费观看视频高清| 久久精品国产自在天天线| 黑人欧美特级aaaaaa片| 久久99热6这里只有精品| 久久精品熟女亚洲av麻豆精品| 亚洲中文av在线| 亚洲美女黄色视频免费看| 欧美另类一区| 99热全是精品| 久久久久久久久久成人| 18在线观看网站| 夫妻午夜视频| 亚洲欧美一区二区三区国产| 大又大粗又爽又黄少妇毛片口| 高清午夜精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品久久午夜乱码| 国产爽快片一区二区三区| 亚洲,欧美,日韩| 日韩av在线免费看完整版不卡| 日韩av免费高清视频| 我要看黄色一级片免费的| 十八禁高潮呻吟视频| 伦理电影大哥的女人| 日本午夜av视频| 综合色丁香网| 晚上一个人看的免费电影| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线观看三级黄色| 亚洲av免费高清在线观看| 精品亚洲成国产av| 国产免费福利视频在线观看| 激情五月婷婷亚洲| 国产精品久久久久久精品电影小说| 一级爰片在线观看| 22中文网久久字幕| 在线观看国产h片| 亚洲av免费高清在线观看| 狠狠精品人妻久久久久久综合| 99国产精品免费福利视频| 日本与韩国留学比较| 国产视频首页在线观看| 亚洲精品第二区| 国产淫语在线视频| 国产精品国产三级专区第一集| 91久久精品国产一区二区三区| 欧美精品亚洲一区二区| 国产毛片在线视频| 成人国产麻豆网| 亚洲,一卡二卡三卡| 欧美丝袜亚洲另类| 性高湖久久久久久久久免费观看| 最近最新中文字幕免费大全7| 日韩欧美精品免费久久| 亚洲三级黄色毛片| 欧美亚洲 丝袜 人妻 在线| 9色porny在线观看| 女的被弄到高潮叫床怎么办| 大香蕉久久网| 97超视频在线观看视频| 国产黄色视频一区二区在线观看| 亚洲精品一二三| 汤姆久久久久久久影院中文字幕| 日本av手机在线免费观看| 亚洲成人一二三区av| 欧美人与善性xxx| 日本av免费视频播放| 人人妻人人添人人爽欧美一区卜| 亚洲精品乱码久久久v下载方式| 女人久久www免费人成看片| 人人妻人人澡人人看| 国产成人免费无遮挡视频| 中文字幕制服av| www.av在线官网国产| 久久精品久久久久久噜噜老黄| 91久久精品国产一区二区成人| 国产爽快片一区二区三区| 久久女婷五月综合色啪小说| 一级片'在线观看视频| 亚洲欧美中文字幕日韩二区| 亚洲无线观看免费| 人人妻人人澡人人看| 欧美日韩国产mv在线观看视频| 精品亚洲成a人片在线观看| 一级毛片黄色毛片免费观看视频| 亚洲欧洲日产国产| 三级国产精品片| 伦精品一区二区三区| 91aial.com中文字幕在线观看| 日韩成人av中文字幕在线观看| 国产一级毛片在线| 99精国产麻豆久久婷婷| 男女边吃奶边做爰视频| 欧美 日韩 精品 国产| 国产成人精品福利久久| 在线观看三级黄色| 在线观看www视频免费| 亚洲在久久综合| 国产熟女午夜一区二区三区 | 欧美三级亚洲精品| 多毛熟女@视频| 99热这里只有是精品在线观看| 欧美激情国产日韩精品一区| 亚洲av.av天堂| 99久国产av精品国产电影| 男女无遮挡免费网站观看| 国产永久视频网站| 精品少妇久久久久久888优播| 亚洲国产成人一精品久久久| 成年美女黄网站色视频大全免费 | 成年人免费黄色播放视频| 99热国产这里只有精品6| 丰满少妇做爰视频| 少妇被粗大的猛进出69影院 | 极品人妻少妇av视频| 久久国内精品自在自线图片| 熟女av电影| 少妇被粗大猛烈的视频| 日韩欧美精品免费久久| 免费人成在线观看视频色| 18在线观看网站| 国产黄片视频在线免费观看| 又黄又爽又刺激的免费视频.| 美女大奶头黄色视频| 在线观看国产h片| 特大巨黑吊av在线直播| 成人毛片60女人毛片免费| 激情五月婷婷亚洲| 国产乱人偷精品视频| 国产成人精品久久久久久| 男女边吃奶边做爰视频| 高清毛片免费看| 亚洲欧美精品自产自拍| 观看美女的网站| 亚洲美女搞黄在线观看| 亚洲激情五月婷婷啪啪| 嫩草影院入口| 考比视频在线观看| 日韩不卡一区二区三区视频在线| 国产精品99久久99久久久不卡 | 婷婷色综合www| 男女无遮挡免费网站观看| 久久99精品国语久久久| 免费观看无遮挡的男女| 夜夜爽夜夜爽视频| 欧美日韩精品成人综合77777| 日韩欧美精品免费久久| 中国美白少妇内射xxxbb| 美女主播在线视频| 成年人午夜在线观看视频| 久久av网站| 视频区图区小说| 国产精品久久久久成人av| 日韩一区二区视频免费看| 街头女战士在线观看网站| 亚洲欧美日韩卡通动漫| 免费少妇av软件| 大香蕉久久成人网| 大片免费播放器 马上看| 国产免费又黄又爽又色| 欧美人与性动交α欧美精品济南到 | 一级黄片播放器| 九草在线视频观看| 日本wwww免费看| 色婷婷久久久亚洲欧美| av有码第一页| 搡老乐熟女国产| 国产爽快片一区二区三区| 人妻制服诱惑在线中文字幕| √禁漫天堂资源中文www| 在线观看www视频免费| 成人亚洲精品一区在线观看| 亚洲综合色网址| 一边摸一边做爽爽视频免费| av.在线天堂| 国产日韩一区二区三区精品不卡 | 街头女战士在线观看网站| 精品少妇久久久久久888优播| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区| 日日摸夜夜添夜夜爱| 日本欧美视频一区| 一区二区av电影网| 国产日韩欧美视频二区| 亚洲成人av在线免费| 男女边吃奶边做爰视频| 精品亚洲成a人片在线观看| 亚洲精品日韩在线中文字幕| 少妇猛男粗大的猛烈进出视频| 欧美97在线视频| 亚洲色图 男人天堂 中文字幕 | 麻豆乱淫一区二区| 满18在线观看网站| 嫩草影院入口| 亚洲国产精品999| 亚洲av.av天堂| 国产精品国产三级国产专区5o| 国产成人精品在线电影| 99九九线精品视频在线观看视频| 久久国产精品男人的天堂亚洲 | 久久 成人 亚洲| 男人添女人高潮全过程视频| 国产伦精品一区二区三区视频9| 国产精品久久久久久精品电影小说| 国产欧美亚洲国产| 国产片内射在线| 亚洲精品乱码久久久v下载方式| 日韩人妻高清精品专区| 亚洲综合精品二区| www.av在线官网国产| 国产男女内射视频| 亚洲不卡免费看| 国产精品 国内视频| 免费观看无遮挡的男女| 亚洲无线观看免费| 色哟哟·www| 免费黄网站久久成人精品| 欧美精品人与动牲交sv欧美| 高清不卡的av网站| 亚洲精华国产精华液的使用体验| 大码成人一级视频| 欧美 日韩 精品 国产| 欧美精品一区二区大全| 日韩欧美精品免费久久| 黄色一级大片看看| 激情五月婷婷亚洲| 青春草视频在线免费观看| 欧美人与性动交α欧美精品济南到 | 在线观看免费日韩欧美大片 | 久久99精品国语久久久| 国产亚洲精品久久久com| 亚洲人成网站在线观看播放| 777米奇影视久久| 久久久久久久国产电影| 91aial.com中文字幕在线观看| 日本欧美国产在线视频| 美女脱内裤让男人舔精品视频| 欧美性感艳星| 人妻 亚洲 视频| 精品国产一区二区久久| 久久久久久久大尺度免费视频| 成人毛片60女人毛片免费| 久久国产亚洲av麻豆专区| 80岁老熟妇乱子伦牲交| 国产男人的电影天堂91| 午夜福利视频在线观看免费| 搡女人真爽免费视频火全软件| 亚洲精品第二区| 精品熟女少妇av免费看| 18禁在线无遮挡免费观看视频| 亚洲精品一二三| av又黄又爽大尺度在线免费看| 国产一级毛片在线| 91国产中文字幕| 男人爽女人下面视频在线观看| 亚洲av二区三区四区| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美成人综合另类久久久| 日韩一区二区三区影片| 狠狠婷婷综合久久久久久88av| 久久久久久久精品精品| 在线播放无遮挡| 亚洲综合色网址| 久久精品国产鲁丝片午夜精品| 欧美精品一区二区免费开放| 秋霞在线观看毛片| 母亲3免费完整高清在线观看 | 一本—道久久a久久精品蜜桃钙片| 欧美亚洲日本最大视频资源| 久久99精品国语久久久| a 毛片基地| 最近中文字幕2019免费版| 久久精品久久精品一区二区三区| 欧美人与善性xxx| √禁漫天堂资源中文www| 丁香六月天网| 精品久久久久久久久亚洲| 狂野欧美白嫩少妇大欣赏| 永久网站在线| 我要看黄色一级片免费的| 女人久久www免费人成看片| 成人18禁高潮啪啪吃奶动态图 | 色吧在线观看| 国产亚洲精品第一综合不卡 | 大片免费播放器 马上看| 欧美激情国产日韩精品一区| 9色porny在线观看| 尾随美女入室| 亚洲欧洲精品一区二区精品久久久 | 日本-黄色视频高清免费观看| 在线观看免费高清a一片| 国产伦理片在线播放av一区| 国产视频首页在线观看| 午夜免费男女啪啪视频观看| 亚洲国产欧美在线一区| 欧美精品高潮呻吟av久久| 免费播放大片免费观看视频在线观看| 999精品在线视频| 亚洲一区二区三区欧美精品| 亚洲欧美一区二区三区国产| 日韩人妻高清精品专区| 最近2019中文字幕mv第一页| 亚洲精品aⅴ在线观看| 一区二区三区四区激情视频| 亚洲精品一区蜜桃| 99re6热这里在线精品视频| 久久99热6这里只有精品| 亚洲性久久影院| 亚洲av成人精品一二三区| 久久午夜福利片| 午夜视频国产福利| 欧美日本中文国产一区发布| 看非洲黑人一级黄片| 在线观看免费高清a一片| 极品少妇高潮喷水抽搐| 人体艺术视频欧美日本| 一区二区三区四区激情视频| 国产日韩欧美亚洲二区| 天天操日日干夜夜撸| 午夜免费鲁丝| 丝袜喷水一区| 国产在线一区二区三区精| 亚洲国产日韩一区二区| 国产精品三级大全| 99久久人妻综合| 精品人妻偷拍中文字幕| 亚洲情色 制服丝袜| 丝袜在线中文字幕| 人妻制服诱惑在线中文字幕| 亚州av有码| 国产日韩欧美视频二区| 久久久国产精品麻豆| 中国三级夫妇交换| 美女视频免费永久观看网站| 色网站视频免费| 十分钟在线观看高清视频www| 麻豆精品久久久久久蜜桃| 亚洲人成77777在线视频| 亚洲国产av影院在线观看| 国产成人aa在线观看| 日本免费在线观看一区| 欧美性感艳星| 亚洲精品乱码久久久v下载方式| 日韩av在线免费看完整版不卡| 王馨瑶露胸无遮挡在线观看| 亚洲精品久久成人aⅴ小说 | 街头女战士在线观看网站| 一级毛片我不卡| 午夜视频国产福利| 极品人妻少妇av视频| 免费看av在线观看网站| 中文字幕人妻丝袜制服| 丰满迷人的少妇在线观看| 青春草视频在线免费观看| 欧美日韩成人在线一区二区| 亚洲精品久久午夜乱码| 亚洲不卡免费看| 国产精品无大码| 2021少妇久久久久久久久久久| 亚洲,欧美,日韩| 伦理电影免费视频| 免费黄频网站在线观看国产| 亚洲精品日韩在线中文字幕| 春色校园在线视频观看| 性高湖久久久久久久久免费观看| 国产日韩一区二区三区精品不卡 | 久久国产精品大桥未久av| 制服人妻中文乱码| 国产深夜福利视频在线观看| 纵有疾风起免费观看全集完整版| 久久青草综合色| 男人添女人高潮全过程视频| 亚洲av男天堂| 少妇人妻精品综合一区二区| www.av在线官网国产| √禁漫天堂资源中文www| 女的被弄到高潮叫床怎么办| 日韩av不卡免费在线播放| av在线播放精品| 国产免费又黄又爽又色| 一本色道久久久久久精品综合| 国产乱来视频区| 大陆偷拍与自拍| 精品少妇久久久久久888优播| 热99久久久久精品小说推荐| 高清av免费在线| av在线老鸭窝| 最黄视频免费看| 一区二区三区乱码不卡18| 日韩熟女老妇一区二区性免费视频| 免费久久久久久久精品成人欧美视频 | 久久久国产一区二区| 亚洲美女搞黄在线观看| 亚洲精品日韩在线中文字幕| 精品一区二区三区视频在线| 免费观看av网站的网址| 大片免费播放器 马上看| 少妇高潮的动态图| 亚洲av电影在线观看一区二区三区| 特大巨黑吊av在线直播| 午夜福利网站1000一区二区三区| 啦啦啦中文免费视频观看日本| 日韩电影二区| 少妇人妻久久综合中文| 妹子高潮喷水视频| 制服丝袜香蕉在线| 91精品伊人久久大香线蕉| 日韩在线高清观看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 另类精品久久| 精品人妻偷拍中文字幕| 亚洲欧美清纯卡通| 亚洲中文av在线| 久久影院123| 日韩精品免费视频一区二区三区 | 男女边摸边吃奶| 亚洲欧美清纯卡通| 女人久久www免费人成看片| 男男h啪啪无遮挡| 免费高清在线观看日韩| 亚洲精品日本国产第一区| 母亲3免费完整高清在线观看 | 日本爱情动作片www.在线观看| av有码第一页| 大香蕉97超碰在线| 日韩一本色道免费dvd| 国产精品国产三级国产专区5o| 亚洲熟女精品中文字幕| 极品少妇高潮喷水抽搐| 色哟哟·www| 亚洲性久久影院| 啦啦啦视频在线资源免费观看| 亚洲人成77777在线视频| 国内精品宾馆在线| 亚洲av中文av极速乱| 免费看av在线观看网站| 少妇人妻久久综合中文| 国产精品成人在线| 免费日韩欧美在线观看| 亚洲精品亚洲一区二区| 亚洲av不卡在线观看| 亚洲怡红院男人天堂| 国产亚洲av片在线观看秒播厂| 五月玫瑰六月丁香| 欧美老熟妇乱子伦牲交| 精品少妇内射三级| 精品国产露脸久久av麻豆| 下体分泌物呈黄色| 国产 一区精品| 久久久久久人妻| 色吧在线观看| 18+在线观看网站| 日本爱情动作片www.在线观看| 欧美最新免费一区二区三区| 亚洲欧美中文字幕日韩二区| 亚洲熟女精品中文字幕| 老女人水多毛片| 亚洲国产av影院在线观看| 欧美一级a爱片免费观看看| 久久ye,这里只有精品| 在线免费观看不下载黄p国产| 亚洲精品aⅴ在线观看| 国产精品国产三级国产专区5o| 欧美bdsm另类| 国产成人精品一,二区| 精品一区在线观看国产| 亚洲中文av在线| 99久国产av精品国产电影| freevideosex欧美| 成年av动漫网址| 日韩中文字幕视频在线看片| 尾随美女入室| 欧美日韩国产mv在线观看视频| 中文天堂在线官网| 特大巨黑吊av在线直播| 国产男人的电影天堂91| 人人妻人人爽人人添夜夜欢视频| 一区二区三区免费毛片| 男女啪啪激烈高潮av片| 纯流量卡能插随身wifi吗| 各种免费的搞黄视频| 少妇的逼水好多| 免费观看性生交大片5| 中文字幕人妻丝袜制服| 久久99热6这里只有精品| 日韩制服骚丝袜av| 欧美变态另类bdsm刘玥| 美女福利国产在线| 交换朋友夫妻互换小说| 女的被弄到高潮叫床怎么办| 中文精品一卡2卡3卡4更新| 国产伦理片在线播放av一区| 日韩欧美精品免费久久| 2022亚洲国产成人精品| 亚洲精品一区蜜桃| av一本久久久久| 国产白丝娇喘喷水9色精品| 成年人午夜在线观看视频| 如日韩欧美国产精品一区二区三区 | 黄色配什么色好看| 国产精品熟女久久久久浪| 中文欧美无线码| 老女人水多毛片| 亚洲国产精品一区三区| 精品国产露脸久久av麻豆| 午夜激情福利司机影院| 菩萨蛮人人尽说江南好唐韦庄| 久热这里只有精品99| 啦啦啦视频在线资源免费观看| 我的老师免费观看完整版| 国产一级毛片在线| 国产男人的电影天堂91| 欧美日本中文国产一区发布| 在线 av 中文字幕| 一区在线观看完整版| 日本欧美视频一区| 啦啦啦中文免费视频观看日本| 永久网站在线| 香蕉精品网在线| 国产一区亚洲一区在线观看| 国产成人一区二区在线| 我要看黄色一级片免费的| av播播在线观看一区| 丝袜喷水一区| 黄色怎么调成土黄色| 国产免费福利视频在线观看| 飞空精品影院首页| 啦啦啦在线观看免费高清www| 国产探花极品一区二区| 亚洲精品日本国产第一区| 五月伊人婷婷丁香| 久久婷婷青草| 亚洲成人手机| 亚洲精品乱码久久久久久按摩| 久久精品国产亚洲av天美| 激情五月婷婷亚洲| 99精国产麻豆久久婷婷| 日本黄色片子视频| 黑人高潮一二区| 不卡视频在线观看欧美| 国产精品嫩草影院av在线观看| 国产精品无大码| 精品久久久噜噜| 亚洲欧美日韩卡通动漫| 国产永久视频网站| 亚洲精品一区蜜桃| 一级二级三级毛片免费看| 免费观看av网站的网址| 最近的中文字幕免费完整| 只有这里有精品99| 国产精品一区www在线观看| 日韩一区二区三区影片| 一级,二级,三级黄色视频| 日韩亚洲欧美综合| 国产精品一区二区三区四区免费观看| 国产免费视频播放在线视频| 国产男女超爽视频在线观看| 91精品伊人久久大香线蕉| 人体艺术视频欧美日本| 国产精品久久久久成人av| 欧美精品人与动牲交sv欧美| 在线观看www视频免费| 看非洲黑人一级黄片| 亚洲国产精品国产精品| av又黄又爽大尺度在线免费看| 午夜视频国产福利| 黄色怎么调成土黄色| 成人亚洲精品一区在线观看| 亚洲综合精品二区| 伊人亚洲综合成人网| 性色av一级| 精品国产一区二区久久| 国产成人精品一,二区| 男女高潮啪啪啪动态图| 91精品三级在线观看| 午夜福利视频精品| 寂寞人妻少妇视频99o| 最近的中文字幕免费完整| 日本色播在线视频| 久久精品国产鲁丝片午夜精品| 亚洲综合色网址| 日韩在线高清观看一区二区三区| 久久97久久精品| 插逼视频在线观看| 久久精品夜色国产| 亚洲天堂av无毛| 免费人成在线观看视频色| av卡一久久| 中文欧美无线码| 亚洲精品久久成人aⅴ小说 | 精品久久久噜噜| 丝袜美足系列| 日本av免费视频播放| 久久这里有精品视频免费| 午夜老司机福利剧场| 欧美日韩视频精品一区| 亚洲丝袜综合中文字幕| 国语对白做爰xxxⅹ性视频网站| 国产 一区精品| 亚洲人成网站在线观看播放| 免费观看的影片在线观看| 国产国语露脸激情在线看| a级片在线免费高清观看视频| 国产精品熟女久久久久浪| 国产精品秋霞免费鲁丝片| 99国产精品免费福利视频| 欧美xxxx性猛交bbbb| 国产欧美另类精品又又久久亚洲欧美| 日韩熟女老妇一区二区性免费视频|