河口生态系统处于陆地、河流和海洋相互作用的界面上，具有一定的开放性、敏感性和脆弱性^[1]。河口生态系统在人为活动和气候变化的影响下，容易受到富营养化、有害藻华、缺氧和pH值变化等方面的威胁^[2]。我国长江口^{[3, 4]}、珠江口^[5]等大型河口近岸海域仍然面临水质污染严重、海洋生态功能受损、滨海湿地生态退化等问题。相比于大型河流，中小型河流河口在人为活动的影响下，更容易出现海洋生态问题，河口生态环境的保护十分重要。广西北部湾沿岸多为中小河流，河口常有滨海湿地、红树林等生态系统，具有重要的生态价值。广西北海冯家江是北海唯一的内陆潮汐河流，冯家江及江口东面海域是北海滨海国家湿地公园，拥有200 hm²原生红树林、丰富的鸟类和鱼类资源。冯家江口西面海域为北海银滩，具有较高的旅游观光价值^{[6, 7]}。冯家江口及其周边海域旅游开发活动活跃，人与海洋环境关系密切，海洋生态环境质量对该区域的发展有着重要意义。冯家江流域水环境未进行治理之前，大量未经处理的养殖废水、雨污水、农业废水等被直接排入冯家江中，使得冯家江及邻近海域生态系统承受巨大的环境压力，面临着水质恶化、生态退化和生产力下降等生态环境问题^[8]。2018年广西北海市开展冯家江流域水环境治理工作，治理效果和治理后研究区海洋生态环境质量是人们关注的重点，但在这方面一直缺乏科学数据作为支撑，因此有必要对冯家江口及其邻近海域生态环境进行综合性评价。

河口及其邻近海域同时受河流和海洋影响，生态系统组成复杂，其性质由该区域水文物理、地质、化学环境以及生物所决定^[9]。该区域的生态环境容易受人类活动影响，对其进行海洋生态环境综合评价具有一定复杂性。海洋生态环境评价是目前海洋研究的重要内容之一，随着对海洋环境评价方法的不断探索和研究，海洋环境评价方法也从单一指标评价发展到海洋生态环境综合评价。现在国外^[10-12]和国内^[13-15]采用的海洋生态环境综合评价模型均属于多参数评价体系，能够较全面地评估海域的生态环境质量。赖俊翔等^{[16, 17]}曾采用《近岸海洋生态健康评价指南》中的方法对广西钦州湾和北仑河口近岸海域进行生态环境健康评价，结果表明广西钦州湾和北仑河口近岸海域均处于亚健康状态。目前关于北海冯家江及其周边海域生态环境质量的研究较少，在前人的调查研究中，通常只是采用单一指标对海区中的水质^[18]、沉积物^[19]或生物生态^[20]进行评价，未对海区生态环境进行综合评价。

本文结合前人对研究区海洋生态环境的调查研究和评价结果，采用指标体系法^{[21, 22]}，结合水环境、沉积物环境和生物生态3个方面构建三级评价模型，对研究区海洋生态环境进行综合评价, 拟为冯家江河口生态保护、滨海湿地管理和区域旅游开发提供科学数据支撑。

1 材料与方法 1.1 研究区概况

研究区位于广西北海银海区冯家江口及其邻近海域，地理坐标为109°6′-109°14′ E，21°22′-21°26′ N。研究区包含冯家江河口生态系统、湿地生态系统和红树林生态系统3种典型的近岸海洋生态系统。该区域属海洋性季风气候，具有典型的亚热带特色。年平均气温22.6℃，年平均降水量1 832.6 mm，降雨主要集中在每年的5-9月，冯家江上游有鲤鱼地水库，在水库调节下，河流径流量季节变化较小。区域风向季节变化显著，冬季盛吹北风，夏季盛吹偏南风。研究区潮汐类型为不正规的全日潮。

1.2 样品及分析

样品来自国家海洋局北海海洋环境监测中心站2021年9月在广西北海冯家江口及其邻近海域开展的海洋环境质量调查，对研究区20个水质站位、10个沉积物站位、12个生态站位和3个潮间带断面(图 1)进行样品采集和分析。水质要素包括海水温度、盐度、pH值、溶解氧、化学需氧量、硝酸盐、亚硝酸盐、铵盐、活性磷酸盐、石油类、汞、镉、铅、铜和锌等15项。沉积物要素包括石油类、有机碳、铜、铅、镉、锌和汞7项。生物生态调查内容包括浮游植物、大型浮游动物、大型底栖生物和潮间带生物。样品采集和分析均按《海洋监测规范》(GB 17378-2007)^[23]中的要求执行，其中浮游植物现场调查采用浅水Ⅲ型浮游生物网(网口面积0.1 m²，网口直径37 cm，网长140 cm)由海底至海面作垂直拖网1次；浮游动物现场调查采用浅水Ⅰ型浮游生物网(网口面积0.2 m²，网口直径50 cm，网长145 cm)由海底至海面垂直拖网1次；底栖生物现场调查采用开口面积为0.05 m²的抓斗式采泥器采集，每站采样2次。

1.3 数据分析 1.3.1 水质要素评价

采用单因子标准指数法^[24]对海水中pH值、溶解氧、化学需氧量、无机氮、活性磷酸盐、石油类、汞、镉、铅、铜、锌等水质要素进行评价，根据研究区所在海洋功能区划管理要求，评价标准采用《海水水质标准》(GB 3097-1997)^[25]中的二类海水水质标准。重金属是影响水质的重要因素之一，利用重金属污染指数法^[22]对海水中的汞、镉、铅、铜、锌等5种重金属污染水平进行评价。为了解研究区海水水质要素的空间分布规律，采用Ocean Data View软件绘制海水温度、盐度、pH值、溶解氧、化学需氧量、无机氮、活性磷酸盐、石油类和重金属污染指数的空间分布图^[26]。

1.3.2 沉积物要素评价

采用单因子标准指数法对沉积物中的石油类、有机碳、汞、镉、铅、铜和锌进行评价，评价标准采用《海洋沉积物质量》(GB 18668-2002)^[27]中的一类海洋沉积物标准。另外，采用瑞典学者Hankanson提出的沉积物重金属潜在生态风险指数法^[28]，对沉积物中5种重金属的潜在生态风险进行评价。潜在生态危害指数能反映沉积物中某一种重金属的生态危害污染程度，根据潜在生态危害系数法，先计算海区沉积物中的汞、镉、铅、铜和锌等5种重金属的潜在生态危害指数，然后再结合毒性系数计算出各调查站位沉积物中多种重金属的潜在生态风险指数(RIs)。重金属RIs反映的是研究区或调查站位多种重金属整体潜在生态危害程度。重金属RIs值越高，潜在生态风险越大。综合考虑前人研究中沉积物的采样范围与采样时间^{[19, 29, 30]}，本文采用1988年在廉州湾和铁山港湾调查中的沉积物重金属平均值作为沉积物重金属背景参考值^[30]。沉积物重金属含量背景参考值和毒性系数见表 1。

1.3.3 生物生态评价

研究区的浮游植物、浮游动物、底栖生物和潮间带生物评价采用生物多样性指数(H′)进行分析，生物多样性指数(H′)采用Shannon-Wiener多样性指数^[31]公式进行计算：

$ H^{\prime}=-\sum\limits_{i=1}^s P_i \log _2 P_i, $

式中，H′为生物多样性指数，s为生物总种类数，P_i为第i种生物的个体数与总个体数的比值。

生物多样性指数(H′)适用于海区环境污染程度评价，当H′值为3-4时表示环境清洁，H′值为2-3时表示环境受轻度污染，H′值为1-2时表示环境受中度污染，H′值< 1时表示环境受严重污染^[32]。

1.3.4 生态环境综合评价指标体系和模型

为了能够更加系统全面地评价研究区海洋生态环境质量现状，本文采用指标体系法对研究区进行综合评价。范海梅等^[22]曾采用指标体系法对长江口及其邻近海域生态环境进行综合评价，本文在其构建的三级评价标准体系的基础上，将溶解氧、化学需氧量和有机碳纳入评价标准体系，对研究区海洋生态环境进行综合评价。构建的三级评价标准体系将评价指标分为目标层、准则层和指标层3层。目标层为水环境、沉积物环境和生物生态，准则层为选择的典型指标，指标层为调查中样品具体分析要素，各项指标见表 2。构建指标体系模型的难处在于对指标进行分级和赋值评价，本文参考前人的相关研究成果^[22]，确定准则层各项评价指标的标准值范围和对应的赋值(表 2)。

对准则层指标进行赋值后，可以根据以下模型计算环境综合评价指数，确定海洋生态环境综合状况。评价模型如下：

$ E=\frac{W+S+B}{3}, $

式中，E为环境综合评价指数，W为水环境指数，取溶解氧含量赋值(W₁)、化学需氧量赋值(W₂)、无机氮含量赋值(W₃)、活性磷酸盐含量赋值(W₄)、重金属污染指数赋值(W₅)和石油类含量赋值(W₆)的平均值；S为沉积物环境指数，取沉积物的重金属生态风险指数赋值(S₁)、有机碳赋值(S₂)和石油类含量赋值(S₃)的平均值；B为生物生态指数，取浮游植物生物多样性指数赋值(B₁)、浮游动物生物多样性指数赋值(B₂)、底栖生物生物多样性指数赋值(B₃)和潮间带生物生物多样性指数赋值(B₄)的平均值。根据得分进行环境评价，参考前人研究成果^[22]，水环境指数、沉积物环境指数、生物生态指数和综合评价指数分级及环境评价见表 3。评价指数的数值为1-4，数值越大环境质量越好。

2 结果与分析 2.1 水质评价

研究区海水中的温度、盐度、pH值、溶解氧、化学需氧量、无机氮、活性磷酸盐、石油类和重金属调查结果见表 4，各水质要素的空间分布见图 2。由图 2可知，研究区海水中的盐度、pH值和溶解氧由江口向外呈递增趋势，溶解氧在江口附近区域出现高值区。海水中化学需氧量、无机氮、活性磷酸盐和石油类含量由江口向外呈递减趋势，其中无机氮与活性磷酸盐在出江口处迅速下降。海水中重金属污染在冯家江口邻近海域较低，在东南面海区污染水平相对较高。由1, 2, 3, 4, 5, 6, 9, 10号调查站位组成的断面能反映各水质要素由江口向外的变化趋势(图 3)，以5号站为界线，江口与邻近海域海水中各水质要素分布存在明显差异。

采用单因子标准指数法进行评价，结果见表 4。研究区海水中pH值、化学需氧量、石油类、汞、镉、铅和铜均符合二类海水水质标准。溶解氧在冯家江口的1, 2, 3号站超出二类海水水质标准，其中2, 3号站满足三类标准，1号站满足四类标准。无机氮和活性磷酸盐在1-4号站都超出二类海水水质标准，冯家江口4个站中海水无机氮超出四类标准，活性磷酸盐在1-3号站超出四类标准，在4号站满足四类标准。锌在17号站超出二类海水水质标准，满足三类标准。综合来看，冯家江口水体中的溶解氧含量较低，无机氮和活性磷酸盐含量较高，冯家江口海水受地表水的影响较大，如果按照地表水的标准进行评价，则1-4号站均能满足地表水水质标准中的Ⅳ类标准。

采用重金属污染指数法进行评价，结果表明各站位海水重金属污染指数(0.07-0.32)较低，平均值为0.15，研究区海水重金属污染水平低。

2.2 沉积物环境质量评价

研究区沉积物中石油类、有机碳和重金属的含量不高，石油类和有机碳的平均含量分别为121.3×10^-6和0.45%，汞、镉、铅、铜和锌等重金属平均含量分别为0.037×10^-6，0.14×10^-6，5.40×10^-6，2.29×10^-6和17.8×10^-6。单项标准指数法评价结果显示，沉积物中有机碳、石油类、铜、铅、锌、镉、汞等评价指标均符合一类海洋沉积物标准(表 5)。沉积物中的重金属潜在生态危害污染程度分析结果表明，铜、铅、锌、镉、汞等5种重金属的潜在生态危害指数均远小于40，生态潜在危害指数Hg > Cd > Pb > Cu > Zn。海区各调查站位的重金属RIs值为14.34-56.53，平均值为29.61，研究区沉积物重金属潜在生态风险低。

2.3 生物生态评价

研究区生物调查统计结果见表 6，浮游植物、浮游动物、底栖生物和潮间带生物的种类数、密度和生物量在各调查站位的变化范围比较大。海区浮游植物密度较大，平均密度为297.34×10⁶ cells/m³；浮游动物平均密度为205 ind./m³，平均生物量为188.5 mg/m³；底栖生物平均密度为513.6 ind./m²，平均生物量为198.86 g/m²；潮间带生物平均密度为860.4 ind./m²，平均生物量为352.50 g/m²。海区浮游植物和浮游动物的多样性指数平均值分别为2.31和2.28，指示海洋环境为轻度污染，海区底栖生物和潮间带生物的生物多样性指数平均值分别为1.77和1.97，指示海洋环境为中度污染。

从各个调查站位的生物分布来看，冯家江口及其邻近海域存在一定的差异，具体见图 4。位于冯家江口2号和3号站的浮游植物种类数分别为6种和9种，密度分别为48.32×10⁶ cells/m³和13.4×10⁶ cells/m³，远低于研究区的平均值18种和297.34×10⁶ cells/m³。2号和3号站的浮游动物种类数分别为1种和2种，密度分别为10 ind./m³和20 ind./m³，也远低于研究区的平均值9种和205 ind./m³。浮游动物的种类数、密度和生物量分布由江口向外呈递增趋势(图 4)。底栖生物的分布与浮游植物和浮游动物不同，江口和邻近海域底栖生物的种类没有明显差异，而密度和生物量江口明显高于邻近海域。潮间带生物的种类数、密度和生物量分布与底栖生物相似，具有江口高于邻近海域的特征。冯家江口2号和3号站浮游动物的生物多样性指数较低，分别仅有0和1，与邻近海域差异较大。浮游植物、底栖生物和潮间带生物的生物多样性指数在江口和邻近海域差别不大(表 6，图 4)。

2.4 生态环境综合评价

对研究区各调查站位的准则层指标进行赋值和综合评价，结果见表 7。由赋值结果可知，在水环境方面，冯家江口的2号和3号站位中的溶解氧、化学需氧量、无机氮和活性磷酸盐的赋值结果明显低于其他站位；在沉积物环境方面，2号和3号站沉积物中重金属的潜在生态风险指数赋值低于其他站位；在生物生态方面，2号和3号的浮游动物多样性指数赋值低于其他站位。其他指标赋值结果在江口和邻近海域没有明显差异。

研究区的水环境指数为3.5，水环境良，其中溶解氧、化学需氧量、活性磷酸盐赋值结果为4，无机氮、重金属污染指数和油类赋值结果为3。研究区的沉积物环境指数为3.7，沉积物环境优，其中重金属潜在生态风险指数和有机碳的赋值结果为4，石油类的赋值结果为3。研究区的生物生态指数为2.5，生物生态环境中等，其中浮游植物和浮游动物的多样性指数赋值为3，底栖生物和潮间带生物多样性指数赋值仅为2，相对较低。研究区的生态环境综合评价指数为3.2，海区生态环境质量为优。

综合评价结果显示，位于冯家江口的2号和3号站的生态环境综合评价指数较低，分别仅有2.4和2.8，生态环境质量分别为中等和良好。邻近海域的生态环境综合指数均大于3，生态环境质量为优。各调查站位的水环境指数(Indices of MwENV)、沉积物环境指数(Indices of SdENV)、生物生态指数(Indices of BeCHA)和综合评价指数(Integrated indices)见图 5。其中2号站和3号站的水环境指数都较低，2号站的生物生态指数较低，2号、3号和6号站的综合评价指数有明显的递增趋势，表明由冯家江口向外生态环境质量逐渐变好。

3 讨论 3.1 水环境和沉积物环境影响因素分析

广西北海冯家江口及其邻近海域水环境总体评价为良好，冯家江口水环境污染相对严重，超标的水质要素主要有溶解氧、无机氮和活性磷酸盐。与2016年海区环境质量调查结果^[18]相比，2021年研究区海水中无机氮的超标率由30%降到20%，石油类超标率由60%降到0，化学需氧量超标率由25%降为0，研究海区经过治理后水环境质量得到改善。

冯家江属于潮汐河流，所在流域过去排污情况严重。得益于2018年冯家江水环境治理工程，大量的排污口和养殖排水被整治清理，冯家江口海水中的石油类和化学需氧量含量明显降低。受城镇生活污染物、人畜粪便及农业生产中化肥施用等影响，冯家江流域地表径流的氮磷营养盐含量较高^[34]，地表径流将较多的氮磷营养盐带入冯家江中，导致冯家江口的无机氮和活性磷酸盐超标。研究区1-4号站位位于冯家江口，该区域受河流和潮汐的共同影响，由江口向外水体交换能力逐渐增强。地表径流带来的无机氮和活性磷酸盐被不断稀释，因此由江口向外呈现逐渐降低趋势(图 3)。

冯家江口海水中溶解氧、无机氮和活性磷酸盐在部分站位超出四类海水水质标准，但这种超标现象影响范围较小，溶解氧、无机氮和活性磷酸盐含量在出江口位置迅速恢复到正常水平，邻近海域的各项水质要素均能满足二类海水水质标准，因此研究区整体水环境仍为良好。

沉积物是持久性有机污染物、石油类和重金属主要的汇，沉积物若受到污染将会对海区的浮游生物、底栖生物乃至红树林生态系统产生不利影响。本文的分析结果表明研究区的沉积物环境处于优良状态，调查结果显示沉积物中有机碳、石油类、铜、铅、锌、镉、汞含量均较低，这与前人的研究结果^[19]相同，冯家江口及其邻近海域沉积物仍保持清洁状态，重金属潜在生态风险低。

3.2 生物多样性影响因素分析

由于淡水和陆源物质的注入，河口及其邻近水域形成了独特的河口类型海洋生态系统。一般情况下河口区生物的种类组成较为复杂，多样性指数较高，不过冯家江口及其邻近海域的生物多样性指数并不高(表 6)。调查结果显示研究区的浮游植物和浮游动物的生物多样性整体处于中等水平，而底栖生物和潮间带生物的生物多样性整体处于较低水平。

生物多样性与海洋环境中的水质、沉积物质量关系密切，环境污染会对海区的生物多样性造成不利影响。环境对生物多样性的影响主要是对敏感种的影响，当环境受到污染，不耐污种会消失，耐污种大量繁殖，导致种类单一，个体数量较多，生物多样性下降^[35]。研究海区的大型底栖生物优势种为珠带拟蟹守螺(Cerithidea cingulata)和斜肋齿蜷(Sermyla riqueti)，潮间带生物的优势种为斜肋齿蜷、珠带拟蟹守螺、奥莱彩螺(Clithon oualaniensis)和毛掌活额寄居蟹(Diogenes penicillatus)，耐污种的优势明显，这与该海区2017年的调查结果^[20]相似。冯家江流域经过治理后，水环境和沉积物环境质量得到改善，但是大型底栖生物和潮间带生物种群结构仍以耐污种为主，生物多样性仍未得到恢复。这可能与底栖生物对环境的敏感性较低有关，因为对环境变化不敏感，底栖生物种群结构变化可能滞后。底栖生物多样性指数较低的情况在北部湾近岸海域较为常见，庞碧剑等^[36]研究发现北部湾近岸海域底栖生物多样性指数平均值仅为1.32，如果采用底栖生物的多样性指数进行污染评价，污染程度评价结果会比用浮游植物多样性指数作出的评价结果污染程度更严重，因此认为底栖生物多样性指数不适于用来评价像北部湾这样较为清洁的海域。

海区浮游植物和浮游动物的生物多样性整体处于中等水平，与邻近海域相比，冯家江口浮游植物和浮游动物都具有种类少、密度低的特征(图 4)。浮游动物在江口内外的分布差异特征更加明显，江口浮游动物生物多样性指数也明显低于邻近海域。研究区浮游植物和浮游动物没有明显的优势种，冯家江口也未出现耐污种(指示环境受到污染)。冯家江口海水中的盐度、pH值、营养盐、溶解氧等要素变化较大(图 3)。浮游动物对这种环境变化反应较为敏感，表现为种类数、密度和生物多样性指数由河口向外呈递增的趋势，这种分布特征与北部湾浮游动物的空间分布规律一致^[37]。

3.3 综合评价结果与主要影响因素分析

广西北海冯家江口及其邻近海域生态环境综合评价结果整体为优，其中冯家江口海洋生态环境质量相对较差，由江口向外生态环境质量呈逐渐变好趋势。冯家江口海水溶解氧含量较低，无机氮和活性磷酸盐含量高，浮游动物生物多样性较低，导致江口海洋生态环境综合评价结果只达到中等或良好水平。

为了研究影响海区综合生态环境质量的主要因素，本文对海水温度、pH值、盐度、溶解氧、化学需氧量、无机氮、活性磷酸盐、重金属污染指数、石油类，沉积物中的重金属风险指数、有机碳、石油类，生物生态中的浮游植物生物多样性指数、浮游动物生物多样性指数、底栖生物生物多样性指数之间的相关性进行分析，结果见表 8。

相关性分析结果显示各评价指标之间具有较强的相关性。由表 8可知海水中化学需氧量、无机氮和活性磷酸盐三者相互之间具有显著的正相关性，化学需氧量、无机氮和活性磷酸盐分别与pH值、盐度和溶解氧具有显著的负相关性。无机氮和活性磷酸盐由地表径流带入海区，化学需氧量能反映水体中还原性物质的污染程度，三者具有显著正相关性，说明地表径流也将部分还原性物质带到海区。无机氮和活性磷酸盐由河口向邻近海域逐渐扩散，浓度逐渐降低，而在河口位置海水由于地表水的注入，pH值和盐度由邻近海域向河口也呈现逐渐降低趋势，因此海水中的无机氮、活性磷酸盐分别与pH值、盐度具有显著的负相关性。海水中还原性物质含量越高，则需要消耗越多的溶解氧，使得水体中溶解氧降低，因此海水中化学需氧量与溶解氧具有一定的负相关性。

相关性分析结果显示浮游动物的生物多样性指数与海水中的化学需氧量、无机氮、活性磷酸盐和石油类具有显著负相关性，与海水的盐度和pH值具有显著正相关性，与溶解氧和浮游植物多样性指数具有正相关性，这表明浮游动物的多样性可能受多种因素影响。营养盐可能是影响浮游动物生物多样性的因素之一，入海河流通常会携带大量的营养盐入海，使得河口区富营养化程度较高，对近岸海域生态系统产生一定影响^[38]。冯家江口海水营养盐含量较高，改变了河口水生生物的生境，一定程度导致浮游动物种类、密度和生物多样性指数降低^{[39, 40]}。河口浮游动物群落结构会受到温度、盐度、食物和动物摄食等因素的影响, 盐度是决定浮游动物分布的关键性因素之一^[41]。本文中冯家江口浮游动物生物多样性指数较低，江口海水中的盐度、pH值、溶解氧等水环境因子也可能是浮游动物生物多样性降低的影响因素。与浮游动物不同，浮游植物和底栖生物多样性指数与各评价指标均不存在相关性，这也说明了浮游植物和底栖生物对河流引起的生境变化反应不敏感。

相关性分析结果还显示海水中的石油类、沉积物重金属生态风险指数与海水中的化学需氧量、无机氮和活性磷酸盐具有正相关性，与海水中的pH值、盐度和溶解氧具有负相关性。这说明研究区海水中的石油类和沉积物中的重金属也可能受地表径流的影响。

综合来看，冯家江口及其邻近海域生态环境质量整体为优，但在河口位置海水仍存在水体富营养化和缺氧的风险，河口水生生物多样性也处于中低水平。冯家江口生态保护需要重视河口营养盐的输入控制和海洋生物多样性的恢复。

4 结论

(1) 广西北海冯家江口及其邻近海域整体水环境良好。冯家江口海水中溶解氧含量较低，无机氮和活性磷酸盐含量较高，个别站位超出四类海水水质标准，邻近海域海水水质均能满足二类海水水质标准。沉积物环境为优，沉积物中有机碳、石油类和重金属指标都符合一类海洋沉积物标准，其中重金属含量低，潜在生态风险较低。

(2) 海区浮游植物和浮游动物多样性处于中等水平，底栖生物和潮间带生物的生物多样性处于较低水平。与其他类型生物相比，浮游动物对水环境变化反应较为敏感，冯家江口的浮游动物生物多样性较差。

(3) 综合评价结果表明，海区生态环境质量整体为优，其中冯家江口的生态环境质量相对较差，由江口向外呈逐渐向好趋势。

(4) 冯家江口水体中溶解氧含量较低，无机氮、活性磷酸盐含量较高，浮游动物多样性较差，导致江口海洋生态环境综合评价结果相对较差，只达到中等或良好水平。