趙 茜，王興華，王 蒙，李格樂(lè)，王 輝，殷安齊

（1. 海南熱帶海洋學(xué)院 熱帶海洋生物資源利用與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室，海南 三亞 572022; 2. 海南熱帶海洋學(xué)院 海南省熱帶海洋漁業(yè)資源保護(hù)與利用重點(diǎn)實(shí)驗(yàn)室，海南 三亞 572022; 3. 海南熱帶海洋學(xué)院 海南省近岸海洋生態(tài)環(huán)境過(guò)程與碳匯重點(diǎn)實(shí)驗(yàn)室，海南 三亞 572022; 4. 海南熱帶海洋學(xué)院 海南省現(xiàn)代化海洋牧場(chǎng)工程研究中心，海南 三亞 572022; 5. 海南熱帶海洋學(xué)院 生態(tài)環(huán)境學(xué)院，海南 三亞 572022; 6. 湛江港（集團(tuán)）股份有限公司，廣東 湛江 524068; 7. 三亞市農(nóng)業(yè)技術(shù)推廣服務(wù)中心，海南 三亞 572019）

海洋浮游植物是海洋生態(tài)系統(tǒng)的重要組成部分，是物質(zhì)循環(huán)和能量流動(dòng)的關(guān)鍵過(guò)程[1-2]。浮游植物群落結(jié)構(gòu)的變化可以及時(shí)反映水體環(huán)境的變化以及整個(gè)水生生態(tài)系統(tǒng)的穩(wěn)定性[3-4]，水環(huán)境的變化也會(huì)改變浮游植物群落結(jié)構(gòu)的物種組成和豐度。浮游植物與其生活環(huán)境密切相關(guān)，營(yíng)養(yǎng)鹽水平、溫度和光照條件對(duì)浮游植物群落結(jié)構(gòu)影響很大，其影響水平能夠通過(guò)浮游植物的生長(zhǎng)和繁殖率表現(xiàn)出來(lái)[5-6]。因此，根據(jù)浮游植物群落結(jié)構(gòu)多樣性及其時(shí)空變化，能夠?qū)λw的營(yíng)養(yǎng)水平和污染狀態(tài)作出判斷，有助于深入了解海洋生態(tài)系統(tǒng)的健康狀況。

三亞灣（18°13′～18°17′N；109°22′～109°29′E）是緊挨三亞市區(qū)的六大海灣之一，是三亞市重點(diǎn)監(jiān)測(cè)港灣。隨著海南建設(shè)自由貿(mào)易港、打造國(guó)際旅游島等一系列國(guó)家戰(zhàn)略決策的實(shí)施，三亞市旅游業(yè)、加工業(yè)及港口航運(yùn)等發(fā)展迅速[7-9]。近年來(lái)，由于生活污水和船舶廢水的排放以及海岸工程的建設(shè)，三亞灣近岸海域污染日趨嚴(yán)重，部分區(qū)域出現(xiàn)海水渾濁和海灘變黑的現(xiàn)象[10]。此外，海灣的大量建筑物削弱了該海域的水動(dòng)力，使得該海域?qū)I(yíng)養(yǎng)鹽等污染物質(zhì)的凈化能力降低。目前，三亞灣有關(guān)浮游植物的調(diào)查已開(kāi)展多次：2004 年，楊志浩等[11]在三亞灣共鑒定出169 種網(wǎng)采浮游植物；2011 年，車(chē)志偉等[12]在春、夏、秋、冬4個(gè)季節(jié)對(duì)三亞灣海域網(wǎng)采浮游植物種類進(jìn)行研究和分析；2019 年，程賢松等[13]在春、秋兩季對(duì)三亞灣的浮游植物及水質(zhì)進(jìn)行網(wǎng)采取樣調(diào)查；2020 年4 月和2021 年4 月，洪鵬等[14]對(duì)三亞灣近岸海域浮游植物進(jìn)行2 次網(wǎng)采取樣調(diào)查。近年來(lái)，對(duì)三亞灣浮游植物的研究主要是用淺水Ⅲ型浮游生物網(wǎng)（77 μm）的方式進(jìn)行樣品采集，水采的方式較少。網(wǎng)采樣品中小尺寸浮游植物的遺漏導(dǎo)致該區(qū)域浮游植物豐度的結(jié)果偏低[15-16]；而水采樣品中小尺寸浮游植物豐度較高，誤差較小，水采對(duì)不同尺寸浮游植物統(tǒng)計(jì)更全面[16]，更適用于三亞灣近岸海域浮游植物的定量研究。本研究基于2021 年夏冬兩季三亞灣近岸海域水采浮游植物與水環(huán)境因子調(diào)查，分析浮游植物群落結(jié)構(gòu)特征及其與環(huán)境因子之間的關(guān)系，以期為三亞灣近岸海域的生態(tài)系統(tǒng)保護(hù)提供基礎(chǔ)數(shù)據(jù)和支持。

1 材料與方法

1.1 樣品采集與分析 參考《海洋調(diào)查規(guī)范》[17]，結(jié)合三亞灣近岸海域?qū)嶋H情況，在水深小于10 m處設(shè)置6 個(gè)站位（圖1）。S1位于三亞灣與天涯灣分界處，處于圍填海狀態(tài)；S2位于三亞灣派出所附近，人為干擾較少，可作為背景值；S3位于三亞景利萊文旅商業(yè)中心附近，人為干擾較多；S4位于海月廣場(chǎng)附近的淺橋下，人為干擾較為嚴(yán)重；S5位于三亞港碼頭處，受航運(yùn)情況、船舶廢水的影響較多；S6位于半山半島帆船港處，處于圍填海狀態(tài)。三亞地處熱帶區(qū)域，四季劃分不明顯，于2021年6月（夏季）和12 月（冬季）進(jìn)行采樣。

1.1.1 浮游植物樣品 參考《近岸海域環(huán)境監(jiān)測(cè)技術(shù)規(guī)范》[18]，使用玻璃采水器在表層水深0.5 m處采集水樣，并將水樣立即轉(zhuǎn)入1 L的采樣瓶，用中性福爾馬林溶液固定樣品至終濃度為1%～2%，帶回實(shí)驗(yàn)室靜置至少48 h后采用虹吸法吸去上清液，濃縮至25 mL，獲得浮游植物定量樣品，在蔡司顯微鏡Axiolab 5下進(jìn)行種類鑒定和細(xì)胞計(jì)數(shù)。

1.1.2 水環(huán)境樣品 用1.5 L采水器采集混合水樣，使用多參數(shù)水質(zhì)測(cè)定儀現(xiàn)場(chǎng)測(cè)定溫度（T）、鹽度（S）、酸堿度（pH）、溶解氧（DO）；將采集后的水樣裝入溶解氧瓶和采樣瓶中，置于裝有冰袋的泡沫箱中低溫保存。參考《海洋監(jiān)測(cè)規(guī)范 第4部分:海水分析》（GB 17 378.4-2007）[19]，將帶回實(shí)驗(yàn)室的水樣進(jìn)行化學(xué)需氧量（COD）、五日生化需氧量（BOD5）、氨氮（NH3-N）、亞硝酸鹽氮（NO2-N）、硝酸鹽氮（NO3-N）和無(wú)機(jī)磷（DIP）的測(cè)定，各指標(biāo)的分析方法見(jiàn)表1。

表1 水質(zhì)分析方法

1.2 數(shù)據(jù)處理

1.2.1 優(yōu)勢(shì)種的確定 根據(jù)每個(gè)種類的優(yōu)勢(shì)度（Y）[20]確定優(yōu)勢(shì)種。

式中，Ni表示所測(cè)樣品中第i種物種的個(gè)數(shù)，N表示樣品中所有種類的總個(gè)數(shù)，fi表示第i種在各站點(diǎn)出現(xiàn)的頻率。Y＞0.02的種類為優(yōu)勢(shì)種。

1.2.2 多樣性的分析 采用生物多樣性指數(shù)（Shannon-Wiener[21]；Pielou[22]；Margalef[23]）進(jìn)行多樣性分析。

式中：H′表示Shannon-Wiener多樣性指數(shù)，J表示Pielou均勻度指數(shù)，d表示Margalef豐富度指數(shù)，S表示樣品的總種類數(shù)，Pi表示第i種的個(gè)數(shù)（Ni）與總個(gè)數(shù)（N）的比值（Ni/N）。

1.2.3 統(tǒng)計(jì)分析 運(yùn)用Excel 2019和Origin 2021進(jìn)行統(tǒng)計(jì)分析和繪圖。采用SPSS 26.0對(duì)各項(xiàng)數(shù)據(jù)進(jìn)行正態(tài)性檢驗(yàn)，符合正態(tài)分布且方差齊性時(shí)，采用單因素ANOVA方差分析組間差異性，反之則用Wilcoxon秩和檢驗(yàn)（兩組之間）和Kruskal-Wallis檢驗(yàn)（多組之間）分析組間差異性。采用Canoco 5.0對(duì)優(yōu)勢(shì)浮游植物種類作去趨勢(shì)對(duì)應(yīng)分析（DCA），DCA大于4時(shí)，選擇典范對(duì)應(yīng)分析（CCA）來(lái)解釋優(yōu)勢(shì)浮游植物和環(huán)境因子之間的關(guān)系；DCA小于3時(shí)，應(yīng)選擇冗余分析（RDA）；DCA介于3～4之間，兩種分析模型均可[24]。

2 結(jié)果與分析

2.1 浮游植物種類組成 本研究共鑒定出浮游植物158 種（含變種、變型，不含未定種），隸屬于6門(mén)64 屬。硅藻種類最多，42 屬116 種，占總種類數(shù)的73.42%；甲藻17 屬36 種，占總種類數(shù)的22.78%。其中，角毛藻屬（Chaetoceros）種類最多，共19 種，占總種類數(shù)的12.03%；其次是根管藻屬（Rhizosolenia）和原多甲藻屬（Protoperidinium），均為10 種，各占總種類數(shù)的6.33%；其他61 屬的種類數(shù)均小于10 種。夏季共發(fā)現(xiàn)浮游植物130 種，其中硅藻99 種、甲藻25 種；冬季共發(fā)現(xiàn)浮游植物111 種，其中硅藻80 種、甲藻26 種。硅藻與甲藻種類數(shù)占比在95%左右，可見(jiàn)該海域浮游植物群落主要由硅藻和甲藻構(gòu)成。夏季S2站位的浮游植物種類最多，為88 種，S6最少，為36 種；冬季S4站位的浮游植物種類最多，為80 種，S6最少，為21 種（圖2）。

圖2 浮游植物種類組成

2.2 浮游植物細(xì)胞豐度 如圖3所示，在夏季和冬季，S4站位浮游植物豐度均最高，分別為39.04×105個(gè)·L-1和20.46×105個(gè)·L-1；S6站位浮游植物豐度最低，分別為4.33×105個(gè)·L-1和0.60×105個(gè)·L-1。在季節(jié)方面，夏季浮游植物平均豐度（15.07×105個(gè)·L-1）顯著高于冬季（7.06×105個(gè)·L-1）。

圖3 浮游植物細(xì)胞豐度

2.3 浮游植物優(yōu)勢(shì)種 調(diào)查海域的浮游植物優(yōu)勢(shì)種季節(jié)變化明顯：夏季（表2）優(yōu)勢(shì)種13 種，鏈狀曲殼藻（Achnanthes catenate）為第一優(yōu)勢(shì)種（Y=0.123），細(xì)胞豐度為13.28×105個(gè)·L-1，占總豐度的14.71%，其豐度最大值出現(xiàn)在S4站位（8.63×105個(gè)·L-1），僅在S5站位未發(fā)現(xiàn)該藻種。冬季（表3）優(yōu)勢(shì)種10 種，翼根管藻（Rhizosolenia alata）為第一優(yōu)勢(shì)種（Y=0.109），細(xì)胞豐度為4.62×105個(gè)·L-1，占總豐度的10.89%，其豐度最大值出現(xiàn)在S3站位（1.78×105個(gè)·L-1），最小值出現(xiàn)在S6站位（0.003×105個(gè)·L-1）。整體來(lái)說(shuō)，三亞灣近岸海域硅藻優(yōu)勢(shì)最為明顯。此外，冬季一些物種的豐度占比雖高于0.02，但因其出現(xiàn)頻率低，未成為優(yōu)勢(shì)種，如鏈狀曲殼藻（A. catenate）和中肋骨條藻（S. costatum）等。

表2 夏季浮游植物優(yōu)勢(shì)種及優(yōu)勢(shì)度

表3 冬季浮游植物優(yōu)勢(shì)種及優(yōu)勢(shì)度

2.4 浮游植物群落結(jié)構(gòu)指數(shù) 三亞灣近岸海域浮游植物群落結(jié)構(gòu)指數(shù)如圖4所示，Shannon-Wiener多樣性指數(shù)的變化范圍為2.08～5.16，均值為4.14，最高值出現(xiàn)在夏季的S2站位，最低值出現(xiàn)在冬季的S6站位，水質(zhì)多處于無(wú)污染-輕污染狀態(tài)。Pielou均勻度指數(shù)在0.47～0.82之間變化，均值為0.70。均勻度指數(shù)最高值出現(xiàn)在夏季的S2站位，最低值出現(xiàn)在冬季的S4站位，水質(zhì)多處于輕污染-中污染狀態(tài)。Margalef豐富度指數(shù)變化范圍為1.26～4.45，均值為3.13，最高值為夏季的S2站位，最低值出現(xiàn)在冬季的S6站位，水質(zhì)多處于中污染狀態(tài)。經(jīng)Wilcoxon秩和檢驗(yàn)和Kruskal-Wallis檢驗(yàn)得出，各項(xiàng)浮游植物群落結(jié)構(gòu)指數(shù)時(shí)空差異性不明顯。

2.5 浮游植物與環(huán)境因子相關(guān)性分析

2.5.1 環(huán)境因子的分布特征 三亞灣近岸海域水環(huán)境因子狀況見(jiàn)表4，經(jīng)Wilcoxon秩和檢驗(yàn)，除S、COD、BOD5、NO3-N和DIP（P＞0.05）以外，其他環(huán)境因子均存在顯著性差異。其中，T（P＜0.01）、pH（P＜0.05）和NH3-N（P＜0.01）的夏季數(shù)值顯著高于冬季，DO（P＜0.01）和NO2-N（P＜0.05）的冬季數(shù)值顯著高于夏季。水質(zhì)符合二類海水水質(zhì)標(biāo)準(zhǔn)。

表4 三亞灣近岸海域水環(huán)境因子狀況

2.5.2 冗余分析（RDA） 選取夏季、冬季各站位出現(xiàn)頻率≥20%，且至少取一個(gè)站位的相對(duì)豐度高于1%的10 種浮游植物進(jìn)行去趨勢(shì)對(duì)應(yīng)分析[25-26]，梯度長(zhǎng)度分別為1.7和2.3，均小于3，采用冗余分析對(duì)優(yōu)勢(shì)浮游植物豐度與環(huán)境因子進(jìn)行RDA排序（圖5），藻類代碼見(jiàn)表5。

表5 進(jìn)行RDA排序的浮游植物種類代碼

圖5 夏季（A）和冬季（B）優(yōu)勢(shì)浮游植物豐度與環(huán)境因子的RDA排序

利用向前引入法逐步篩選環(huán)境因子，經(jīng)Monte Carlo置換檢驗(yàn)，結(jié)果表明，夏季COD（F=15.2，P=0.002）、DO（F=12.0，P=0.002）、NO3-N（F=9.3，P=0.002）、S（F=6.5，P=0.012）和NO2-N（F=3.5，P=0.028）這5 個(gè)環(huán)境因子達(dá)到顯著水平，是影響夏季優(yōu)勢(shì)浮游植物豐度的主要環(huán)境因子；冬季NH3-N（F=8.9，P=0.002）、COD（F=3.5，P=0.002）、S（F=5.1，P=0.002）、NO3-N（F=4.0，P=0.002）和pH（F=2.9，P=0.012）這5個(gè)環(huán)境因子達(dá)到顯著水平，是影響冬季優(yōu)勢(shì)浮游植物豐度的主要環(huán)境因子。

圖5表明，夏季優(yōu)勢(shì)浮游植物主要集中在第一、四象限，與環(huán)境因子有明顯關(guān)系，可以分為3類：第一類包括圓海鏈藻（T. rotula）、柔弱擬菱形藻（P. delicatissima）、丹麥細(xì)柱藻（L. danicus）、勞氏角毛藻（C. lorenzianus）、中肋骨條藻（S. costatum），與DO、pH、BOD5呈正相關(guān)，與NH3-N、DIP、NO2-N、S和NO3-N呈負(fù)相關(guān)；第二類是小球藻（C.sp.），與T和COD呈負(fù)相關(guān)；第三類包括骨條藻（S.sp.）、角毛藻（C.spp.）、鏈狀曲殼藻（A.catenate）和日本星桿藻（A. japonica），與T、COD呈正相關(guān)。冬季優(yōu)勢(shì)浮游植物主要集中在第二、三象限，也分為3類：第一類包括角毛藻（C.spp.）、柔弱擬菱形藻（P. delicatissima）、勞氏角毛藻（C.lorenzianus）、丹麥細(xì)柱藻（L. danicus）、尖刺擬菱形藻（P. pungens），與T、DO呈正相關(guān)，與NO2-N、NO3-N和DIP呈負(fù)相關(guān)；第二類為小球藻（C.sp.），與COD、BOD5、NH3-N呈正相關(guān)，與pH、S呈負(fù)相關(guān)；第三類包括翼根管藻（R. alata）、暹羅角毛藻（C. siamense）、扁面角毛藻（C. compressus）、并基角毛藻（C. decipiens），與pH、S正相關(guān)，與COD、BOD5、NH3-N呈負(fù)相關(guān)。

結(jié)果表明，硅藻門(mén)在三亞灣近岸海域浮游植物群落結(jié)構(gòu)中起主要貢獻(xiàn)作用，RDA分析結(jié)果表明，S、pH、DO、COD、NO2-N、NO3-N和NH3-N是夏冬兩季浮游植物群落結(jié)構(gòu)的關(guān)鍵驅(qū)動(dòng)因子。

3 討 論

3.1 浮游植物群落結(jié)構(gòu)組成 本研究結(jié)果表明，三亞灣近岸海域夏冬兩季浮游植物共158 種，兩季皆以硅藻為主，其次為甲藻，兩者共占總種類數(shù)的96.20%，其他門(mén)類僅占極少數(shù)，該海域浮游植物群落結(jié)構(gòu)的生態(tài)類型為硅藻-甲藻型，與其他海灣的浮游植物群落結(jié)構(gòu)研究結(jié)果一致[27-29]。關(guān)于三亞灣海域浮游植物長(zhǎng)期變化的研究報(bào)道均是使用淺水Ⅲ型浮游生物網(wǎng)進(jìn)行采集樣品[11-14]，而基于水采方法的浮游植物群落結(jié)構(gòu)比傳統(tǒng)的網(wǎng)采方法更能反映浮游植物的實(shí)際情況，因?yàn)榫W(wǎng)采樣品不能將水體中全部粒徑的浮游植物采獲，算是1種半定量的方法。本研究中水采的浮游植物豐度的數(shù)量級(jí)為105個(gè)·L-1，而程賢松[13]、車(chē)志偉[12]、楊志浩[11]、洪鵬[14]等網(wǎng)采樣品的數(shù)量級(jí)在105個(gè)·m-3左右。在鏡檢樣品時(shí)發(fā)現(xiàn)，粒徑小的浮游植物遠(yuǎn)遠(yuǎn)多于粒徑大的，不僅是部分硅藻，絕大多數(shù)甲藻的尺寸都小于77 μm。網(wǎng)采樣品中浮游植物種類數(shù)、豐度與甲藻種類數(shù)均呈現(xiàn)出水采高于網(wǎng)采的規(guī)律[11,30-32]，此次水采樣品的總豐度及甲藻種類數(shù)均高于之前的網(wǎng)采樣品，所以采樣方式可能是導(dǎo)致該海域浮游植物豐度數(shù)量級(jí)與前人研究不同的原因。

三亞灣近岸海域夏冬兩季優(yōu)勢(shì)種差異較大，交替明顯，18 種優(yōu)勢(shì)種中只有5 種是共有優(yōu)勢(shì)種，均為廣溫廣布性種，優(yōu)勢(shì)種中近岸廣布性種占主要地位，并在夏冬兩季均有出現(xiàn)。夏季第一優(yōu)勢(shì)種為鏈狀曲殼藻（A. catenate），鏡檢結(jié)果表明，其主要以6～20 個(gè)的群體形態(tài)出現(xiàn)，易志勇[33]研究發(fā)現(xiàn)，鏈狀曲殼藻（A. catenate）在低氮濃度下更易于形成鏈狀群體，當(dāng)?shù)獫舛葹?.5 mg·L-1時(shí)，鏈狀曲殼藻（A. catenate）生長(zhǎng)情況最佳。三亞灣近岸海域的氮元素水平均處于較低的狀態(tài)，鏈狀曲殼藻（A.catenate）以群體形態(tài)出現(xiàn)與易志勇的研究是一致的。冬季的第一優(yōu)勢(shì)種是翼根管藻（R. alata），夏季數(shù)量不多，卻在冬季成為第一優(yōu)勢(shì)種，羅先香等[31]對(duì)黃河口及鄰近海域的調(diào)查發(fā)現(xiàn)，翼根管藻（R. alata）對(duì)水環(huán)境變化的敏感程度相對(duì)較高，與S呈顯著正相關(guān)關(guān)系，與T、COD和pH呈顯著負(fù)相關(guān)關(guān)系，與本研究中的翼根管藻（R. alata）的豐度分布情況大致相同。

3.2 浮游植物群落的多樣性 Shannon-Wiener多樣性指數(shù)用于判斷群落結(jié)構(gòu)的復(fù)雜程度，浮游植物的種類越多或者分布越均勻，值越大，群落結(jié)構(gòu)越復(fù)雜[34]。均勻度指數(shù)用于判斷群落是否成熟穩(wěn)定[35]，豐富度指數(shù)的值越大，群落結(jié)構(gòu)越穩(wěn)定，生態(tài)系統(tǒng)的抗干擾能力越強(qiáng)[36]。本研究結(jié)果表明，夏季各站位浮游植物群落的H′、J和D的平均值分別為4.44、0.73和3.33，冬季分別為3.84、0.67和2.93，說(shuō)明該海域浮游植物多樣性較高，物種較為豐富，群落結(jié)構(gòu)較復(fù)雜、穩(wěn)定性較好，生態(tài)系統(tǒng)抗干擾能力較強(qiáng)。

綜合對(duì)比圖2和圖4可以發(fā)現(xiàn)，夏、冬兩季物種數(shù)與豐富度指數(shù)變化大致相同。兩者在S2站位均達(dá)到峰值，不同之處出現(xiàn)在夏季的S4站位，雖然此站位的物種數(shù)處于較高水平，但由于硅藻（豐度占比99.12%）大量繁殖導(dǎo)致豐度遠(yuǎn)遠(yuǎn)高于其他站位，使得該站位豐富度指數(shù)較低，出現(xiàn)偏差。單一藻類優(yōu)勢(shì)過(guò)高，會(huì)影響群落結(jié)構(gòu)指數(shù)的變化趨勢(shì)[37]，夏季鏈狀曲殼藻（A. catenate）的豐度占比為14.71%，冬季翼根管藻（R. alata）的豐度占比為10.89%，是引起該海域多樣性指數(shù)與豐富度指數(shù)變化趨勢(shì)不同的原因。三亞旅游業(yè)淡、旺季特征明顯，4 月至10 月是旅游淡季，而1 月至3 月和11 月至12 月為旺季[38]。S2站位相對(duì)于其他站位人流量最少，受人為影響最小，夏季該站位的浮游植物種類數(shù)、多樣性指數(shù)和豐富度指數(shù)均處于高峰期，表明該站位物種種類多、多樣性好，分布均勻。

夏冬兩季的浮游植物群落差異性明顯。夏季浮游植物豐度和種類數(shù)均顯著高于冬季。S2站位受人類活動(dòng)干擾較少，浮游植物種類較豐富，在夏冬兩季均為最高值，且物種分布均勻，群落穩(wěn)定性大，水質(zhì)狀況整體優(yōu)良。S4站位的兩季硅藻種類數(shù)相差不大，但冬季甲藻種類數(shù)高于夏季，導(dǎo)致該站位種類數(shù)季節(jié)差異與其他站位相反，即夏季種類數(shù)低于冬季。經(jīng)差異性檢驗(yàn)后發(fā)現(xiàn)，S1站位與S6站位無(wú)顯著性差異（P=0.805）。S1站位處以大壩作為三亞灣與天涯灣的分界，S6站位位于半山半島帆船港港內(nèi)。二者均處于圍填海狀態(tài)，僅有圍填海的程度不同。圍海可以興建養(yǎng)殖池塘、圍海制鹽等，填?？梢杂糜诟劭诖a頭建設(shè)、城鎮(zhèn)及工業(yè)用地發(fā)展等[39]，而圍填海這一用海方式能夠引發(fā)圍填海附近海域生物多樣性降低、優(yōu)勢(shì)種演替和群落結(jié)構(gòu)變化[40]，所以S1、S6站位的物種多樣性均低于其他站位可能與圍填海工程有關(guān)，且處于港灣內(nèi)部的S6站位受到的影響更為嚴(yán)重，浮游植物種類和豐度均處于最低水平，該站位浮游植物群落復(fù)雜程度低，穩(wěn)定性差，物種豐富度低，生態(tài)系統(tǒng)抗干擾能力弱。S4站位與S1、S6站位均存在極顯著性差異（P＜0.01）。S4站位位于海上淺橋的橋下，人流量較大，有機(jī)污染程度較高，COD和浮游植物豐度在該站位均達(dá)到最高。COD是反映水體有機(jī)和無(wú)機(jī)可氧化物質(zhì)污染的常用指標(biāo)，較高的有機(jī)物含量，有利于浮游植物的生長(zhǎng)[37]，與該站位浮游植物豐度達(dá)到最高的結(jié)果相一致。

3.3 浮游植物群落與環(huán)境因子的關(guān)系 浮游植物的生長(zhǎng)受溫度、鹽度、營(yíng)養(yǎng)鹽等諸多環(huán)境因子的綜合影響，各環(huán)境因子及其間的交互作用共同影響浮游植物群落結(jié)構(gòu)[41]，RDA分析結(jié)果顯示，夏季浮游植物主要受S、DO、COD、NO2-N和NO3-N的影響，冬季主要受S、pH、COD、NO3-N和NH3-N的影響。

S是影響浮游植物營(yíng)養(yǎng)鹽吸收和繁殖的關(guān)鍵因素之一[42]，本研究中，S的季節(jié)變化不明顯，兩季的S在S5站位均為最低值，可能與船舶廢水、生活污水等的排放有關(guān)。林庚明等[43]通過(guò)主成分分析發(fā)現(xiàn)S對(duì)楚科奇海浮游植物分布有重要影響；卜世勛等[44]發(fā)現(xiàn)S是影響2020 年撫寧海灣扇貝養(yǎng)殖區(qū)的浮游植物群落結(jié)構(gòu)的主要環(huán)境因子之一，與本研究結(jié)果相符。

三亞灣近岸海域水質(zhì)呈弱堿性，而弱堿性的水體環(huán)境有助于浮游植物吸收CO2，從而提高光合作用的效率[45]。本研究中發(fā)現(xiàn)，夏季的pH值是顯著高于冬季的，甲藻的種類數(shù)與豐度占比在冬季均有所升高，據(jù)Chust等[46]認(rèn)為，CO2的增加對(duì)硅藻應(yīng)該是略微有益的，對(duì)甲藻非常有益，因此，推測(cè)pH對(duì)冬季的浮游植物群落結(jié)構(gòu)組成影響較夏季明顯。

近岸海域是敏感區(qū)域，在多種作用的影響下，會(huì)導(dǎo)致DO呈現(xiàn)出不同的特點(diǎn)[42]，如水質(zhì)中營(yíng)養(yǎng)鹽含量較為豐富時(shí)，浮游植物大量繁殖，光合作用增強(qiáng)，導(dǎo)致溶解氧含量上升[47]。本研究中，夏季的DO含量與優(yōu)勢(shì)浮游植物呈顯著正相關(guān)關(guān)系，與O′boyle等[48]對(duì)愛(ài)爾蘭河口和近岸沿海水域浮游植物生物量積累影響因素的研究結(jié)果相一致。同時(shí)，DO也與溫度呈負(fù)相關(guān)關(guān)系，夏季雖有大量的浮游植物進(jìn)行光合作用，產(chǎn)生DO，但水溫的升高有利于海洋動(dòng)物的生長(zhǎng)繁殖，消耗水中的DO[30]，從而導(dǎo)致夏季的DO含量低于冬季。

COD是表示水體還原性物質(zhì)含量的指標(biāo)，也是衡量水體有機(jī)物含量的指標(biāo)[37]。本研究中，COD無(wú)明顯的季節(jié)性變化，但空間變化較為顯著。S4站位位于海月廣場(chǎng)附近，人流量較其他站位大，還處于海岸工程的橋下，COD含量均顯著高于其他站位。人類活動(dòng)輸入大量的營(yíng)養(yǎng)物質(zhì)[49]，導(dǎo)致有機(jī)物質(zhì)含量較高，為浮游植物的大量繁殖提供了條件，特別是甲藻，其豐度和種類數(shù)在此站位均達(dá)到最高值。甲藻與硅藻的生存條件相似，這與柴然等[37]對(duì)嶗山近岸浮游植物群落的研究結(jié)果相一致。

氮、磷是浮游植物生長(zhǎng)中常見(jiàn)的限制營(yíng)養(yǎng)元素，可以調(diào)節(jié)浮游植物的生長(zhǎng)，影響浮游植物的種類和豐度[50]。鄒景忠等[51]提出營(yíng)養(yǎng)鹽閾值標(biāo)準(zhǔn)為：水質(zhì)中DIN含量為0.2～0.3 mg·L-1，DIP含量為0.04 mg·L-1。本研究中，三亞灣近岸海域兩個(gè)季節(jié)的DIN和DIP含量均低于閾值，證明了該海域?yàn)樨殸I(yíng)養(yǎng)海域。結(jié)果表明，優(yōu)勢(shì)浮游植物豐度與無(wú)機(jī)氮呈顯著負(fù)相關(guān)關(guān)系，與無(wú)機(jī)磷的相關(guān)性不顯著。有研究結(jié)果表明，硅藻和甲藻對(duì)氮磷的需求是不同的，其中，硅藻對(duì)氮元素的需求較大，甲藻對(duì)磷元素的需求較大[52]。夏冬兩季硅藻豐度占比分別為97.81%和96.88%，對(duì)氮元素的需求更大，磷元素與甲藻豐度均處于較低水平，相關(guān)性不明顯。

綜上所述，三亞灣近岸海域浮游植物以硅藻門(mén)為主要種類，該海域浮游植物豐度存在明顯的季節(jié)性差異。該海域浮游植物群落多樣性較高，穩(wěn)定性較好，生態(tài)系統(tǒng)具有較強(qiáng)的抗干擾能力。兩季水環(huán)境因子表明該海域水質(zhì)處于二類海水水質(zhì)狀態(tài)，相關(guān)部門(mén)應(yīng)該采取相關(guān)措施，降低人為活動(dòng)及海岸工程建設(shè)對(duì)三亞灣近岸海域水質(zhì)的影響，保護(hù)海域生態(tài)環(huán)境。