涠洲岛造礁石珊瑚群落变化特征及其环境影响因子

杨振雄,张敬怀,吕向立,郑琰晶,姜广甲,3,*,时小军,陶 伟,李秀芹,吕意华

1 国家海洋局南海环境监测中心,广州 510300 2 自然资源部南沙珊瑚礁生态系统野外科学观测研究站,广州 510300 3 广西南海珊瑚礁研究重点实验室(广西大学),南宁 530004

涠洲岛不仅被认为是南海珊瑚为适应高温迁移的潜在避难所,还是研究人类活动影响珊瑚礁生态的理想实验室[1-2]。近年来梁文等[3-5]、黄晖等[6]、王欣等[7]、陈天然等[8]等学者从不同角度对涠洲岛进行了珊瑚礁生态调查,并研究了石珊瑚群落生态变化的影响因子。例如,黄晖等[6]2005年6个站位的调查发现石珊瑚死亡率高,推断涠洲岛珊瑚遭到了长期的、特别是近两年的持续破坏,并认为原因是人为因素导致的海水水质变差、沉积物增多,以及海水温度异常和人为直接破坏或破坏性的捕鱼作业。梁文等[3]2007、2008年6个站位的调查发现涠洲岛受1998年的“厄尔尼诺”引起的全球性极高温事件、2008年1—2月极低温事件及不合理人类活动胁迫影响较大。陈天然等[8]通过测定滨珊瑚骨骼生长率变化验证了余克服等[9]提出的夏季异常高SST的持续上升可能会使珊瑚生长处于非常敏感边缘的推论,证实了涠洲岛珊瑚会对气候变暖和极端事件进行响应,并指出涠洲岛的人类活动胁迫降低了珊瑚耐受气候变化的阈值,降低了其适应气候变化和自我进化的能力。Yu等[1]假设涠洲岛是一个理想的珊瑚避难所(珊瑚礁具有复杂形态结构、珊瑚覆盖率高、种类丰富和复原能力强的特点),通过2015年的6个站位的现状调查进行评估,结果发现涠洲岛珊瑚礁在逐渐加强的人类活动影响下快速退化,其作为潜在珊瑚避难所的功能性在减弱。全球性变暖、近岸海域富营养化和草食性动物的减少常被认为是珊瑚礁退化的主要原因[10]。虽然与全球气候相比,人类活动的影响范围较小,但对于地理单元小且相对独立的涠洲岛来说,其环境承载能力有限,人类活动干扰超过了珊瑚礁生态系统的容忍限度,珊瑚礁生态系统会发生退化甚至崩溃[8]。因此,研究涠洲岛珊瑚礁现状及其面临的威胁对于深入研究区域珊瑚礁生态特征及其影响因素具有重要的意义。

尽管涠洲岛珊瑚礁历史调查资料较多,但对于全面评价涠洲岛珊瑚群落特征和演替规律还需要更多的现场数据支撑。现有的生态调查多以分析珊瑚礁区群落特征为主,对石珊瑚的群落结构特征、优势种群及其与环境因子关系的定量分析相对较少,而人类活动影响多以定性描述为主,缺乏系统的调查和定量的评价[2]。可用于定量化趋势评价的数据较少,历史资料还不能完全反映当前涠洲岛珊瑚礁生态资源及其面临的威胁。同时,很多学者认为涠洲岛珊瑚退化是过度捕捞、污染、物理损伤(踏踩、抛锚、水下工程等)和大尺度环境变化(如气候异常)等原因共同造成的[2],但其不同时期的主导影响因素需要进一步研究。

为摸清我国珊瑚礁资源空间分布格局和生态状况,自然资源部2019年启动了包括涠洲岛在内的全国珊瑚礁基线调查,在涠洲岛海域布设了15个站位45条断面,对涠洲岛珊瑚礁资源进行生态精细化调查。基于涠洲岛珊瑚礁基线调查数据,分析石珊瑚种群生态特征及其与环境因子的相关性,并结合历史数据,探求10年内(2009—2019年)涠洲岛珊瑚退化的主导因素,为涠洲岛珊瑚礁资源开发、保护和管理提供基础数据及技术支撑。

1 研究方法

1.1 研究区域和调查方法

涠洲岛(109°00′—109°15′E,21°00′—21°10′N)位于广西壮族自治区南部,北部湾海域中部(图1)。调查日期为2019年5月26日—6月1日。采用截线样带调查法(Line intercept transect),以500—1000 m/条间距,布设垂直于岸线的共计15个站位45条样带(每个站位布设3条样带,每条样带50 m),水深范围为2—12 m。拍摄水下样线录像和水下特写照片,提取样线底质类型、石珊瑚覆盖率、石珊瑚属种组成、硬珊瑚补充量、大型藻类覆盖率、泥沙覆盖率等数据。石珊瑚鉴定通过分析水下特写照片和样品细微结构。同步开展水温、盐度、溶解氧、悬浮物、无机氮、活性磷酸盐等环境因子调查,样品采集和分析参照《海洋监测规范》(GB17378—2007)、《海洋调查规范》(GB/T12763—2007)。

图1 调查站位示意图Fig.1 Skecth map of sampling stations

1.2 数据分析

造礁石珊瑚的优势属种选取采用重要值法[1-2]。同时,为全面评价石珊瑚与环境因子间的关系,本文还采用以下的标准选取参与环境因子分析的石珊瑚作为优势种,即:该物种在各样点出现的频度≥12.5%,且在至少一个样点的平均覆盖率≥1%(阈值确定的标准为各物种频度和平均覆盖率的排序中值)。采用Shannon-Wienner指数(H′)、Pielou均匀度指数(J)和Simpson多样性指数(d)分析珊瑚群落多样性特征[11]。使用Raunkiaer频度法分析珊瑚群落稳定性[12],分级标准为:A为频度1%—20%的比例,B为频度21%—40%的比例,C为频度41%—60%的比例,D为频度61%—80%的比例,E为频度81%—100%的比例。通常样方内各频度等级种群占总物种数比例关系为A>B>C≥或≤D级

采用去趋势对应分析方法(DCA)分析珊瑚群落在不同区域的相互关系；采用聚类分析(Cluster)初步评估不同环境因子中各珊瑚群落分布的分类关系；采用主成分分析法(PCA)识别影响站位分组的主要环境因子；鉴于珊瑚生物群落在水下鉴定可能有一定误差,本文对石珊瑚优势属、优势种与环境因子的关系进行分析:使用DCA来确定生物数据适用于线性模型还是单峰型模型(若DCA中所有排序轴的值均小于3,则选择线性模型),若为线性模型则使用冗余分析RDA获得珊瑚群落与环境变量的关系(利用CANOCO软件自动向前选择程序对环境因子进行逐一筛选,并利用Monte Carlo检测,其结果可判断第二轴的解释量是否影响第一轴解释量的显著性),挑选出对物种分布显著的环境因子(P<0.05))；使用多元线性回归(Linear Regression)拟合群落生态特征值与环境因子的回归方程。物种矩阵经过lg(x+1)转换,环境数据除pH外进行lg(x+1)转换。数据分析使用CANOCO 4.5、SPSS 16、PRIMER 5.0等软件。

2 结果与分析

2.1 石珊瑚种群结构特征

研究表明,涠洲岛造礁石珊瑚有9科22属42种,以裸肋珊瑚科种类最多(60.5%)。总体来说,石珊瑚覆盖率为3.8%—19.8%,平均为10.5%,西南、西北和东南部覆盖率较高(图2)。石珊瑚主要分布在5—10 m的水深范围内,珊瑚覆盖率在4.8%—42.8%之间(平均值:15.3%)；<5 m和>10 m的水深区域珊瑚覆盖率相差不大,范围分别为在2.0%—25.8%(平均值:8.4%)和0.6%—23.0%(平均值:8.7%)。石珊瑚主要优势属为角蜂巢珊瑚属(Favitesp.)、滨珊瑚属(Poritesp.)、角孔珊瑚属(Gonioporasp.)和陀螺珊瑚属(Turbinariasp.)等,形态以团块状为主。主要优势种为澄黄滨珊瑚(Poriteslutea)、秘密角蜂巢珊瑚(Favitesabdita)、斯氏角孔珊瑚(Gonioporastutchburyi)、海孔角蜂巢珊瑚(Faviteshalicora)(表1),其中澄黄滨珊瑚、秘密角蜂巢珊瑚的重要值最高,其属级重要值百分比分别为9.7%和8.5%。

图2 造礁石珊瑚覆盖度、主要珊瑚种群及泥沙覆盖分布Fig.2 Distribution of scleractinian coral coverage, major coral community, and sediment coverage

表1 涠洲岛不同水深范围内的主要珊瑚礁优势种Table 1 The dominant species of coral reefs in different water depths of Weizhou Island

多样性指数H′范围为2.64—4.00,均值为3.30,有2/3的海域H′大于3,在涠洲岛西南、东南和东北海域较高,而在南湾(人类活跃区,WZ01)和西部(风浪作用强烈及海岸与海底侵蚀活跃区,WZ09—WZ10)附近海域较低。主要石珊瑚优势种Spearman秩相关分析结果(图3)显示,正相关种对数略多于负相关,正相关对数占总对数的50.5%,负相关对数占48.4%,正负关联比为1.05,极显著或显著正相关、负相关的占15.8%,说明总体上石珊瑚群落种间联结性较松散,仅部分物种之间相关性较强,如秘密角蜂巢珊瑚与海孔角蜂巢珊瑚显著正相关,澄黄滨珊瑚与斯氏角孔珊瑚极显著正相关等。

图3 珊瑚主要优势种对的Spearman秩相关分析Fig.3 Spearman rank correlation analysis of main species in coral communitiesPor.lut :澄黄滨珊瑚 Porites lutea:Fav.abd :秘密角蜂巢珊瑚 Favites abdita；Gon.stu :斯氏角孔珊瑚 Goniopora stutchburyi；Fav.hal:海孔角蜂巢珊瑚 Favites halicora；Tur.pel :盾形陀螺珊瑚 Turbinaria peltata；Fav.pen:五边角蜂巢珊瑚 Favites pentagona；Mer.amp:阔裸肋珊瑚 Merulina ampliata；Ple.ver:多孔同星珊瑚 Plesiastrea versipora；Por.lob:团块滨珊瑚 Porites lobata；Mon.sp.:蔷薇珊瑚 sp. Montipora sp.；Ech.asp:粗糙刺叶珊瑚 Echinophyllia aspera；Gon.asp :粗糙菊花珊瑚 Goniastrea aspera；Gon.sp.:角孔珊瑚 sp. Goniopora sp.；Fav.fav:黄廯蜂巢珊瑚 Favia favus；Por.sp.:滨珊瑚sp. Porites sp.

2.2 石珊瑚群落稳定性、群落分布及其环境分组特征

Raunkiaer频度分析结果(图4)表明,A>B>C≥D>E,A级所占比例最高,48.2%；B级次之,26.8%；C级为10.7%,D级为8.9%,E级所占的比例最低,为5.4%。群落的Raunkiaer标准频度直方图不呈反J型,且A级和E级的比值为8.9,说明珊瑚群落受到严重干扰而处于不稳定或不成熟状态。

图4 珊瑚分布频度和珊瑚群落物种Raunkiaer标准频度分级Fig.4 Frequency distribution & Satandard frequency histogram of coral communitiesA为频度1%—20%的比例,B为频度21%—40%的比例,C为频度41%—60%的比例,D为频度61%—80%的比例,E为频度81%—100%的比例

对优势石珊瑚群落(统计到属)进行DCA排序(图5),发现DCA的第一排序轴值长度为1.733,且所有排序轴的值均不大于1.733。滨珊瑚属(第一优势类群)和扁脑珊瑚属卡方距离较近,表明具有较相似的生长环境条件,二者具有潜在的竞争关系,滨珊瑚属和扁脑珊瑚主要分布在区域Ⅱ(WZ01、WZ02、WZ03和WZ10),即在涠洲岛的东西部,而角蜂巢珊瑚属(第二优势类群)主要分布在区域Ⅰ(WZ04、WZ05),即涠洲岛北部。

图5 主要珊瑚种群DCA分布图Fig.5 The DCA distribution of major coral community 三角形表珊瑚优势属；Favi.sp:角蜂巢珊瑚属 Favites sp.；Por.sp:滨珊瑚属 Porites sp.；Gon.sp:角孔珊瑚属 Goniopora sp.；Tur.sp:陀螺珊瑚属Turbinaria.sp.；Fav.sp:蜂巢(盘星)珊瑚属 Favia sp.；Pla.sp:扁脑珊瑚属 Platygyra sp.

Cluster分析结果(图6)表明,涠洲岛石珊瑚礁区生物群落生境可主要分为两类:一类位于南湾湾内的WZ01站,即受人类活动影响最大的生境；另一类则在南湾湾内之外(除WZ01站外)的涠洲岛其他海域生境。利用PCA对各个站位进行分组,以分析站位中起主要作用的环境因子。结果表明(图6,第一第二主成分共解释了总变异的99.8%)环境因子主要可以分为3组:组1包括站位WZ03—WZ06、WZ08、WZ10、WZ11、WZ14,涵盖了53.3%的区域。组2包括站位WZ09、WZ12和WZ13。组3为站位WZ01。各组的结果与聚类分析Cluster相似,其中,组1与营养盐、大型海藻覆盖率、透明度等相关性最高。组3位于南湾,与其他站位显著差异,且与营养盐负相关最高。

图6 珊瑚群落调查站位分布聚类和环境因子主成分分析Fig.6 Distribution clustering & PCA of the environmental factors in coral communites箭头表环境因子,圆圈为站位号;Dept:水深 Depth;Tran:透明度 Transparency;Tem:温度 Temperature;Sal:盐度 Salinity;COD:化学需氧量 Chemical oxygen demand;SS:悬浮物含量 Suspended solids content;Oil:石油类 Oils;DIN:无机氮Inorganic nitrogen;Chla:叶绿素a Chlorophyll-a;SC:泥沙覆盖率 Sediment coverage；AC:大型海藻覆盖率 Coverage of macroalgae

2.3 珊瑚种群与环境因子的关系

对珊瑚群落生态特征值与主要环境因子(表2)进行相关性和多元线性回归分析,结果表明(表3):珊瑚覆盖率与水体悬浮物含量(-0.53,P<0.05)显著负相关,与水深(0.53,P<0.05)显著正相关；珊瑚种类数与铵盐(0.60,P<0.05)和透明度(0.58,P<0.05)显著正相关,与泥沙覆盖率(-0.62,P<0.05)显著负相关；珊瑚多样性指数与珊瑚种类数(0.87,P<0.01)极显著正相关,也与铵盐(0.59,P<0.05)显著正相关,与泥沙覆盖率(-0.52,P<0.05)显著负相关；硬珊瑚补充量与石油类(-0.56,P<0.05)显著负相关。

表2 涠洲岛海域主要海洋环境因子Table 2 The main marine environmental factors in waters of Weizhou Island

表3 珊瑚群落特征与环境因子形成的回归方程Table 3 Regression equations between coral community characteristics and the environmental factors

利用RDA分析物种与环境因子的关系(DCA各排序轴均小于3),发现对石珊瑚优势属经Monte Carlo检测(Number of permutations排序次数499),筛选出7个主要环境因子(解释量比例为73.8%),为水深、悬浮物含量、无机氮、叶绿素、温度、pH、盐度、化学需氧量,其中为水深、悬浮物含量是显著影响因子(P<0.05)；对珊瑚优势种经Monte Carlo检测(Number of permutations排序次数499),筛选出8个主要环境因子(解释量比例为75.8%),为水深、悬浮物含量、pH、无机氮、磷酸盐、叶绿素、泥沙覆盖率、大型藻类覆盖率,其中为水深、悬浮物含量是显著影响因子(P<0.05)。

石珊瑚优势属及其优势种对海洋环境的需求基本相似(图7),均受水深、悬浮物含量、无机氮浓度的影响较大,特别是悬浮物含量与大部分石珊瑚显著负相关。主要优势属角蜂巢属及优势种海孔角蜂巢珊瑚、秘密角蜂巢珊瑚与水深显著正相关,与叶绿素a、无机氮正相关,但与泥沙覆盖率、pH显著负相关,说明角蜂巢珊瑚属对泥沙覆盖较差的耐受力；滨珊瑚属和角孔珊瑚属及其优势种澄黄滨珊瑚、斯氏角孔珊瑚与泥沙覆盖率相关性较低,表明它们对泥沙覆盖有一定的耐受力,但与大型海藻覆盖率(主要为匍扇藻(Lobophoravariegate))具有显著的负相关关系。涠洲岛大型海藻覆盖率范围为0.8%—39.3%,均值为6.8%,仅在WZ05(16%>>11%)、WZ10站(39.3%>>5.2%)显著高于珊瑚覆盖率,反映了局部区域大型海藻与珊瑚可能存在较强的竞争关系。

图7 珊瑚类群、优势珊瑚种与环境因子RDA分析Fig.7 The RDA analysis of coral community, dominant species and environmental factors黑色细箭头表示珊瑚物种,红色拉丁简写表前几位优势种,蓝色粗箭头表环境因子,黑色圆圈为站位号

3 讨论

3.1 涠洲岛石珊瑚群落变化趋势

历史数据分析表明[1-7],近30年涠洲岛平均石珊瑚覆盖率呈现下降趋势(图8),由20世纪80—90年代的60%以上下降到目前的10%左右。原珊瑚种群主要分布区(西南部、东南部和北部)退化严重,尤以西南部(涠洲岛主要人类活动区,特别是南湾)变化最为明显,由1991年的70%下降至2019年的9.6%。石珊瑚优势类群发生了较大的变化,块状珊瑚逐渐取代分枝状珊瑚成为优势类群,角蜂巢珊瑚、滨珊瑚、角孔珊瑚等是近年来的主要优势群落[1-2]。实际上,珊瑚形态在一定程度上反应生态适应性[15-16],生长迅速、空间竞争能力强、抵抗干扰能力弱的分枝形态的造礁石珊瑚群更容易遭受严重破坏,而生长缓慢、生态幅宽广、对栖息环境要求不高的团块状珊瑚群受到的破坏较小,更适应生长。尽管活珊瑚覆盖率递减和优势群落更替趋势显著,Yu等[1]在2015年发现在涠洲岛珊瑚礁区仍有大量的硬珊瑚补充量,本研究也发现在涠洲岛西南(W3站附近,补充量为6.51个/m2)、北部(WZ05站附近,补充量为5.80个/m2)海域仍有较多的硬珊瑚补充量,表明涠洲岛珊瑚礁在局部区域仍有较大的恢复潜力。

图8 近30年涠洲岛全岛及各区域造礁石珊瑚平均覆盖率变化趋势Fig.8 The trend of scleractinian coral coverage in all regions of Weizhou Island for nearly 30 years

研究发现,造礁石珊瑚多样性指数(H′)较梁文等[17](H′均值为2.40)高,但多样性分布趋势相似,可能由于本研究以种统计多样性指数,而梁文等[17]使用的是属级统计多样性指数。尽管分析结果表明涠洲岛的多样性指数较高,但其群落稳定性较差,原因可能是近年来人类活动干扰(如营养盐输入等)扰乱种间生态位互补机制,使得物种多样性对群落稳定性的正效应削弱[13],如澄黄滨珊瑚与斯氏角孔珊瑚种间联结极显著正相关,显示二者生态位可能互补可共同生长,但耐受能力更强[18]的澄黄滨珊瑚覆盖度和分布范围明显多于斯氏角孔珊瑚,说明对环境扰动具有较强耐受力的种类在一定程度限制其他种群的生长,使得其他种群竞争力减弱,进而降低整个群落的稳定性。

3.2 主要环境因子对涠洲岛石珊瑚群落的影响

3.2.1悬浮泥沙与泥沙沉降物

涠洲岛海域水体悬浮物均值为5.1 mg/L,但在南湾(WZ01站)的悬浮物含量较高(>10 mg/L),与近10年涠洲岛悬浮物值变化范围(4.5—21.4 mg/L)一致[19-20]。虽然整体上悬浮物含量不高,但发现研究区域礁石上泥沙覆盖情况较严重(图2、图9)。全岛53.3%的区域处于严重覆盖,主要是分布在西南、西北和东南近岸海域,又以西南南湾最为显著。

图9 涠洲岛珊瑚礁区泥沙覆盖Fig.9 The sediment cover in coral reef of Weizhou Island

大部分造礁石珊瑚对悬浮沉积物的耐受能力差,悬浮沉积物会在一定程度上抑制珊瑚的生长和发育[21-22]。例如,统计结果表明泥沙覆盖率与透明度呈极显著负相关(-0.735),与珊瑚多样性(-0.515)、珊瑚种类数(-0.615)显著负相关。同时,珊瑚褪色率、死亡率和覆盖率与水体悬浮物有显著的相关性[6]。尽管如此,珊瑚对悬浮环境也有一定的适应能力。珊瑚对悬浮沉积物清除能力与珊瑚的种类、形态、生活习性甚至生长方向有很大关系,也受悬浮沉积物本身的数量以及种类影响[23],如滨珊瑚就能较强耐受悬浮沉积物[24]。一些珊瑚种类为了适应低辐射、低波长光谱等浑浊环境,共生藻会通过增加叶绿素含量及类囊体膜面积的方式来增加对有限光能的吸收[25],如盾形陀螺珊瑚生长更偏向于低波长光谱质量[26],也可通过改变生长形态或骨骼形态适应高悬浮物环境,如柳珊瑚[27]和圆菊珊瑚Montastreaannularis[28],说明具有较强耐受能力的团块状石珊瑚在高悬浮物环境下仍可以普遍生长。

涠洲岛较高的泥沙覆盖率可能与涠洲海岸风化和侵蚀及北部湾风升流周期性对涠洲岛沿岸影响有关[29],特别是春夏季时底层环流搅动所产生的向上覆水补充加上风浪作用造成的紊流扩散,不仅是泥沙沉积物扩散的动力,还是营养盐P的重要补充来源[30]。此外,近年来涠洲岛开展的客货码头、石化栈桥等多项涉海大型工程建设项目,施工期会产生砂石碎屑、粉尘等入海工程废物。同时,潜水观光、海水浴场等旅游项目也可能对海水和底质有较大的扰动。

3.2.2营养盐

珊瑚长期暴露在富营养化海水中,会影响其新陈代谢、捕食和共生藻密度。一般珊瑚礁区的营养盐较贫瘠,硝酸盐为0.1—0.5 μmol/L,铵盐为0.2—0.5 μmol/L,无机磷小于0.3 μmol/L[31]。涠洲岛海域水体活性磷酸盐均值为0.09 μmol/L,无机氮均值为11.61 μmol/L,其中硝酸盐均值为9.26 μmol/L,占比79.7%。铵盐均值为1.73 μmol/L(占比14.9%)。铵盐主要来源于生物体排泄物和有机物分解,微生物活动会导致氨氮在沉积物中积累并以无机氮的形式存于底质中[32]。研究区域内硝酸盐含量远高于铵盐含量,可能与两者时空分布有关。涠洲岛海域受季风影响较大,春夏季节转换时海水温度不高,生物作用较低,冬季贮存硝酸盐尚未被浮游植物大量消耗,水体硝酸盐现存量较高,但细菌分解有机质的速率较低,导致水体中的铵盐较少。

结合涠洲历史调查数据[33-34]可知(图10),近10年来,涠洲岛无机氮(硝酸盐、铵盐等)含量均远高于一般热带珊瑚礁区含量,也大于限制浮游植物生长的营养盐浓度阈值(DIN=1 μmol/L)[35],近年来有增加的趋势。例如,在WZ11、WZ05站,无机氮超二类海水标准(>13 μmol/L),石珊瑚覆盖率和种类数均较低。WZ11站位邻近南湾养殖海域,而WZ05站则发现了大型藻类—匍扇藻(Lobophoravariegate)过度生长,该大型藻易获取营养物质快速生长,能抑制浮浪幼虫的附着或使浮浪幼虫产生回避行为[36]。朱葆华等[37]发现0.001 mol/L的氨氮或硝氮就能使3种珊瑚共生藻的释放明显增多。雷新明等[38]发现高浓度(>15 μmol/L)的硝酸盐、铵盐、磷酸盐及其营养盐组合浓度(>10 μmol/L)均对共生藻光合效率有不同程度的抑制作用,且随着营养盐浓度的升高抑制作用加强。本研究发现硝酸盐浓度可能对造礁石珊瑚生长不利,统计学分析也表明无机氮与大部分珊瑚群落成负相关,而铵盐与珊瑚种类数、珊瑚多样性指数均成显著正相关,说明适量的铵盐有利于珊瑚生长。同时发现,10年来涠洲岛溶解氧的变化范围为4.90—8.12 mg/L[19-20],本研究结果与该范围保持一致,但有55.6%的海域发现超第一类海水水质标准情况。其轻度低氧的现象可能与溶解氧季节性变化或水体富营养化有关[39]。细菌等微生物在低氧环境下固氮作用增强[40],产生氨氮,导致铵盐与珊瑚群落具有较高的相关性。

图10 近30年涠洲岛年均无机氮和无机磷含量及N/P变化趋势Fig.10 The trend of annual average of inoraganic nitrogen, inoraganic phosphorus,and N/P ratio in Weizhou Island for nearly 30 years

涠洲岛磷的存在形态以无机磷为主(约占83.3%),其分布与北部湾季风过程和生物活动密切相关[30]。尽管N/P比自2011年起就高于30,显示涠洲岛海域营养盐磷限制显著,但磷酸盐浓度仍可以满足该海域浮游动植物的需要[33]。研究发现涠洲岛磷酸盐浓度较低,统计分析表明其与溶解氧、叶绿素极显著负相关(-0.65和-0.94,P<0.01),这与海水中的磷循环过程一致,大部分颗粒态的磷最终沉积到水底,但在缺氧或无氧的状态下,微生物分解后的磷很快再矿化为溶解的磷,磷重新悬浮的比例增加[41]。

无机氮、溶解氧超标站位主要分布于涠洲岛西部和西南海域,该海域有南湾及蕉坑口养殖区、鳄鱼山景区、涠洲岛客运码头等人类活动区。统计分析表明南湾及其附近人类活动密集区与涠洲岛其他区域珊瑚环境完全不同,是叶绿素a的高值区(含量5.91 mg/m3,处于中营养水平)和悬浮物高值区(含量在10 mg/L以上),南湾是涠洲岛主要养殖用海,也是南珠的主要养殖区,养殖品种包括扇贝、石斑鱼、鲍鱼和珍珠[42]。养殖活动造成水体富营养化而导致沉积物缺氧或无氧状态,微生物的活动加速了无机盐从底质向上覆水的释放,加快了水体营养盐的循环过程[41]。

4 结论

据2019年5—6月涠洲岛海域珊瑚礁生态调查数据,分析了近10年来涠洲岛造礁石珊瑚的生态特征及其主导影响因子。涠洲岛造礁石珊瑚覆盖率显著降低,珊瑚形态组成趋于团块状化,尽管珊瑚群落多样性较高,但群落结构较不稳定,群落分布较松散,部分优势种种间竞争较激烈,主要优势种群是耐受干扰种类。悬浮物含量是影响石珊瑚群落最显著的环境因子,石珊瑚优势种群在不同水深中分布差异显著,泥沙覆盖率、营养盐对不同石珊瑚种群影响差异较大,反映不同珊瑚种群有不同的适应环境能力,而大型海藻覆盖率在局部区域也对优势珊瑚形成较强的竞争关系。短期内,涠洲岛珊瑚礁退化的主导因素是人类海洋开发活动,例如,营养盐和泥沙沉积物的增加与涠洲岛近海养殖业、生活排污、海岸工程及海岸侵蚀等密切相关。