2019年4月和9月，参考《海洋调查规范》（GB12763.9—2007），在三亚湾海域设置12个调查站位（图1），地理坐标见表1，依据GB12763.6—2007《海洋监测规范》规定，对三亚湾浮游植物进行网采、贮存、运输和预处理。浅水区域使用浅水Ⅲ型浮游生物网（网长140 cm，网口内径37 cm，网口面积0.1 m²）进行网采，距底层2 m的位置开始往上垂直拖网3次采集样品。水深超过50 m的调查站位，由50 m至表层垂直拖网3次采集样品，现场测定气温、水温、水深、盐度、透明度、溶解氧、pH值和电导率。取1 L网采样品用终浓度为2%的甲醛固定后带回实验室备用，剩余样品随后静置沉淀浓缩至50 mL，取0.1 mL浓缩样品在LEICADM6000B生物显微镜下进行种类鉴定和计数。

图  1  三亚湾浮游植物采样站位示意图

Figure 1.  The schematic diagram of phytoplankton survey stations in Sanya Bay

取0.1 mL浓缩样品在生物显微镜下计数、鉴定。Shannon-Wiener多样性指数(H′)、Pielou均匀度指数(J′)、Margalef丰富度指数(D)和优势种确定（优势度Y）的计算公式^[8-11]如下：

式中：H′为多样性指数；S为种类数；P_i=ni/N（ni是第i个物种的个体数，N是全部物种的个体数）；fi为第i种在各个站位出现的频率，P_i为第i种占各个站位细胞总量的比例。

调查期间三亚湾春季平均水温26.5 ℃，秋季平均水温28.87 ℃。调查水域春秋两季pH值均为8.28。春季盐度平均为31，秋季盐度高达35.31。春季电导率为20.53～20.87 μs·cm⁻¹，平均值为20.75 μs·cm⁻¹。秋季电导率有所增高，平均值达到24.31 μs·cm⁻¹。春季溶解氧为8.97～10.5 mg·L⁻¹，平均值为9.88 mg·L⁻¹，明显高于秋季，秋季的平均值5.65 mg·L⁻¹（表2）。

三亚湾调查海域春季鉴定出浮游植物6门82种，秋季鉴定出4门85种。其中，春、秋季硅藻门分别占浮游植物总种类数的67.07%和78.82%，占据优势类群。春季鉴定出甲藻门22种，秋季16种；春秋两季均鉴定出蓝藻门1种和金藻门1种。除此之外，春季还鉴定出定鞭藻门2种。秋季各调查站位鉴定的浮游植物种类均比春季多（图2）。

图  2  2019年三亚湾各站位浮游植物物种数

Figure 2.  Number of phytoplankton species at each station in Sanya Bay in 2019

三亚湾春季浮游植物细胞平均密度为5.13×10³个·m⁻³；秋季浮游植物细胞平均密度为99.8×10³个·m⁻³。秋季平均细胞密度是春季的20倍左右（表3）。除春季站位4-3以定鞭藻门为优势类群，其余站位均以硅藻门为优势种类。春秋两季硅藻细胞密度均明显高于其他藻类，春季占细胞总密度的78.31%，明显低于秋季（96.08%），但春季甲藻细胞密度明显高于秋季（图3）。

图  3  各门浮游植物空间分布

Figure 3.  Spatial distribution of phytoplankton

优势度决定优势种，将优势度≥0.02的浮游植物种类作为该海域的优势种类。春季调查期间三亚湾浮游植物优势种共4种，分别为角毛藻属（Chaetoceros sp.）、印度角毛藻（Chaetoceros indicus）、三角新角藻（Neoceratium tripos）和薄壁几内亚藻(Guinardia flaccida)。角毛藻属（Chaetoceros sp.）为第一优势种，优势度为0.163，平均密度为0.84×10³个·m⁻³（表4），在SYW3-2号站位密度最大；印度角毛藻为第二优势种，优势度为0.091，平均密度为0.62×10³个·m⁻³，在SYW4-1号站位密度最大（图4）。秋季调查期间该海域浮游植物优势种共9种，分别为菱形海线藻(Thalassionema nitzschioides)、热带骨条藻(Skeletonema tropicum)、霍氏半管藻(Hemiaulus hauckii)、日本星杆藻(Asterionella japonica)、丹麦细柱藻(Leptocylindrus danicus)、劳氏角毛藻(Chaetoceros lorenzianus)、柔弱拟菱形藻(Pseudo-nitzschia delicatissima)、窄隙角毛藻(Chaetoceros affinis)和异角毛藻(Chaetoceros diversus)。菱形海线藻为第一优势种，优势度为0.184，平均密度为18.3×10³个·m⁻³（表4），在2-1号站位密度最大；热带骨条藻为第二优势种，优势度为0.074，平均密度为80.0×10³个·m⁻³，在2-1号站位密度最大（表4，图4），优势种显微镜下放大400倍的生物学特征见图5。

图  4  2019年三亚湾优势种空间分布

Figure 4.  Spatial distribution of dominant species in Sanya Bay in 2019

图  5  三亚湾浮游植物优势种显微图片

Figure 5.  Photomicrograph of dominant species of phytoplankton in Sanya Bay

浮游植物多样性指数是反映其种类的多寡和各个种类数量差异的函数关系，均匀度则反映其种类数量的分布情况，可以作为生态监测的参数。计算结果表明，春季三亚湾浮游植物的多样性指数和均匀度平均值分别为3.33和0.80，秋季多样性指数和均匀度分别为4.05和0.79。秋季调查时，所有站位的多样性指数均大于3，表明三亚湾海域生境环境优良。春季丰富度指数在0.98～4.18之间，平均为2.26。秋季丰富度指数在2.13～5.41之间，平均为3.715（表5）。

三亚湾海域12个站点水体的总氮、硝态氮、亚硝态氮、氨氮、总磷、COD测定结果（图6）显示，三亚湾春季总氮含量平均值为4.35 mg·L⁻¹，低于秋季的6.45 mg·L⁻¹。春秋两季硝酸盐氮含量平均值分别为0.83 mg·L⁻¹和0.89 mg·L⁻¹，差异不大。春秋两季亚硝酸盐氮含量平均值均为0.008 mg·L⁻¹。秋季氨氮含量平均值为1.91 mg·L⁻¹，高于春季的平均值（0.53 mg·L⁻¹）。春秋两季总磷含量平均值分别为2.09 mg·L⁻¹和0.80 mg·L⁻¹，秋季含量明显下降。春秋两季COD平均值分别为2.81 mg·L⁻¹和5.02 mg·L⁻¹，秋季的含量高于春季的。

图  6  2019年三亚湾各调查站位水质相关指标测定结果

Figure 6.  Quality-related indicators of seawater at various survey stations in Sanya Bay in 2019

根据采样调查数据，利用Canoco 5软件对浮游植物群落和环境因子关系进行对应分析。选择出浮游植物丰度前10的物种，首先对数据矩阵进行降趋势对应分析(DCA)，根据分析结果中排序轴数值大小特征，判断下一步分析方法。分析结果表明，春秋两季数据采用冗余分析（RDA）较合适。通过蒙特卡洛方法对春秋两季浮游植物群落与环境因子RDA对应关系分析结果。基于三亚湾浮游植物丰度前10的藻株进行相关分析（图7），浮游植物物种主要分布在第一和第四象限，春季总磷、氨氮、水温和pH对浮游植物群落结构影响较大，其中棕囊藻属与总磷和氨氮呈正相关关系，优势藻种中的印度角毛藻、具槽直链藻、丹麦角毛藻和薄壁几内亚藻与水温和pH 2个环境因子呈正相关关系。大角管藻、角毛藻属、三角新角藻、螺端根管藻和膜状舟形藻与总磷、氨氮、硝态氮呈负相关关系。秋季影响浮游植物最大的因素COD、电导率、溶解氧和盐度，水温次之；在水质方面，总氮、总磷和硝态氮影响最大。柔弱拟菱形藻与电导率、COD和总氮呈正相关，与pH、水温、硝态氮和总磷呈负相关；异角毛澡、窄隙角毛澡、菱形海链藻和柔弱根管藻与pH、水温和硝态氮正相关，与COD、电导率和总氮呈负相关；霍氏半根藻、劳氏角毛藻、日本星杆藻、丹麦细柱藻和热带骨条藻与总磷和盐度呈正相关，与溶解氧和亚硝态氮呈负相关。

图  7  三亚湾春季（a）和秋季（b）浮游植物与环境因子RDA分析

Figure 7.  RDA analysis of phytoplankton and environmental factors in Sanya Bay in spring (a) and autumn (b)

三亚湾春秋两个航次调查结果显示，春季共鉴定出浮游植物6门82种，秋季共鉴定出浮游植物4门85种，共发现167种浮游植物，其中春秋两季不同种共145种，共有22种相同浮游植物种类。春秋两季浮游植物种类均以硅藻为主，甲藻次之，与历年三亚湾浮游植物调查研究的结果一致^[12]。同时，杨志浩等的研究发现，三亚湾浮游植物硅藻不仅在种类数上占优势，其细胞密度上也同样占据优势，各季节的浮游植物均以硅藻为优势种，本研究的调查结果显示，春秋两季均以硅藻为优势种，硅藻细胞密度占优势，与之相一致^[13]。定鞭藻门仅在春季出现，秋季各站位浮游植物种类高多于春季。春季藻细胞平均密度为5.13×10³个·m⁻³；秋季藻细胞平均密度为99.8×10³个·m⁻³。秋季藻细胞密度是春季的20倍左右。

春季调查期间三亚湾浮游植物优势种共4种，分别为角毛藻属、印度角毛藻、三角新角藻和薄壁几内亚藻。其中角毛藻属为第一优势种。春季三亚湾受外海水温度影响较大^[14-15]，因此,一些外海广温性种则占据优势。角毛藻属于广温性种，是近岸海域常见的赤潮藻种^[16]。温度盐度是影响角毛藻分布的主要因素^[17]，角毛藻的大量生长是水体富营养化的指征之一^[18]。秋季优势种共9种，分别为菱形海线藻、热带骨条藻、霍氏半管藻、日本星杆藻、丹麦细柱藻、劳氏角毛藻、柔弱拟菱形藻、窄隙角毛藻和异角毛藻。菱形海线藻为第一优势种。春季三亚湾浮游植物的多样性指数、均匀度和丰富度指数平均值分别为3.33、0.80和2.26，秋季多样性指数、均匀度和丰富度分别为4.05、0.79和3.43。综合多样性指数和均匀度指数来看，三亚湾春秋两季水体环境均为优良等级，浮游植物种类多样性高，且结构稳定^[19]。

历年来研究显示，三亚湾主要的影响环境因子包括温度、盐度、气象条件和营养盐含量等^[20]。本研究水质调查结果显示，三亚湾总氮、氨氮和COD含量平均值秋季高于春季，总氮、氨氮和COD含量的增高是水体富营养化的特征指标^[21]，说明秋季浮游植物生境较春季差。海水中COD含量的增高说明了该海域受到有机污染物的污染，COD含量越高污染越严重，COD含量的增多是导致赤潮藻类大量繁殖的因素^[22]。说明可能是由于旅游业的开发，而对三亚湾的环境保护措施不足，导致三亚湾的水质受到污染。秋季总磷含量明显下降，硝酸盐氮和亚硝酸盐氮含量两季差异不明显。

RDA分析结果显示，春季影响浮游植物最主要的环境因子是总磷、氨氮、水温和pH，其中，棕囊藻属与总磷和氨氮呈正相关关系，氮磷元素是影响浮游藻类生长的主要元素^[23-24]。优势藻种中的印度角毛藻、具槽直链藻、丹麦角毛藻和薄壁几内亚藻与水温和pH2个环境因子呈正相关关系。水温变化会进一步导致水体的物理、化学性质和生物的活动发生改变，以此对浮游植物的组成及分布产生影响^[25]。研究报道，在河口浮游植物的调查中发现，pH是影响其分布组成的重要环境因子^[26-27]。秋季影响浮游植物最大的因素COD、电导率、溶解氧和盐度，水温次之。盐度的变化主要是影响浮游植物的渗透压，进而对浮游植物的群落组成产生影响^[28]。在水质方面，总氮、总磷和硝态氮影响最大。柔弱拟菱形藻与电导率、COD和总氮呈正相关，与pH、水温、硝态氮和总磷呈负相关。异角毛藻、窄隙角毛藻、菱形海链藻和柔弱根管藻与pH、水温和硝态氮正相关，与COD、电导率和总氮呈负相关；霍氏半根藻、劳氏角毛藻、日本星杆藻、丹麦细柱藻和热带骨条藻与总磷和盐度呈正相关，与溶解氧和亚硝态氮呈负相关。

总体来说，三亚湾海水中的氮磷营养盐、pH、盐度、COD和海水的温度是影响该海域浮游藻类的主要因素。三亚湾浮游藻类多样性调查发现，该海域浮游植物生存环境优良，但存在个别赤潮藻种。三亚湾属于旅游开放区，应降低旅游业对三亚湾水体的影响以及减少当地生活污水的排入，保护海域水体环境。