1. 引言
湿地是地球上重要的生存环境和生态系统,被誉为“地球之肾”、“生命的摇篮”、“文明的发源地”和“物种的基因库”,湿地与森林、海洋一起并列为全球三大生态系统,它在蓄洪防旱、调节气候、控制土壤侵蚀、促淤造陆、降解环境污染等方面都发挥着极其重要的作用 [1] [2] [3] [4] [5]。天津滨海湿地类型丰富,近年来随着经济的快速发展,由于一系列自然及人为如盲目开垦、污染物排放、水资源短缺和气候变化等因素,导致滨海湿地资源面临着污染和退化的威胁 [6]。相关研究表明,导致天津滨海湿地生态系统健康恶化的主要原因为赤潮灾害频发、城市化进程加速带来的人口增长和大型围填海工程的建设等 [7]。
天津大港滨海湿地生态系统位于天津市滨海新区大港以东,面积为106.37 km2,岸线长度为9.69 km,是天津典型的湿地生态系统。2014年7月,天津市海洋局(现天津市规划和自然资源局)将其选划为海洋生态红线区。天津大港滨海湿地地貌类型为粉沙淤泥质平原海岸,沿岸地势低平,潮滩和水下浅滩宽缓,主要植被为米草。近年来,随着大量海洋工程的开发建设,该海域的海洋环境发生了变化,其原有生态功能也受到了一定影响。大型底栖动物是海洋生态系统的重要组成部分,它们通过摄食、掘穴和建管等活动与周围环境发生着相互影响,其群落结构是认识海洋环境特点、预测环境质量的重要指标,研究它们的结构和功能在理论和实践上都具有重要的意义 [8] [9]。因此,对天津大港滨海湿地生态系统附近海域海洋生态环境现状进行调查研究,旨为天津大港滨海湿地生态系统的环境保护、修复及资源的合理利用提供相关科学依据。
2. 材料与方法
2.1. 调查站位
本研究于2021年8月对天津大港滨海湿地生态系统附近海域进行现场调查,共布设6个调查站位(图1)。
2.2. 调查项目
海水水质调查项目为:pH、溶解氧(DO)、化学需氧量(COD)、无机氮(亚硝酸盐–氮、硝酸盐–氮、氨–氮)、活性磷酸盐、总氮、总磷、叶绿素a;海洋生物调查项目为:大型底栖生物。
2.3. 样品采集与处理
2021年8月对大港滨海湿地生态系统附近海域进行水质和海洋生物调查,样品的采集、保存、运输和分析均按照《海洋监测规范》 [10] 和《海洋调查规范》 [11] 执行。
2.4. 评价方法
2.4.1. 单因子质量指数法
1) 采用单因子指数法进行评价,其标准指数计算公式为 [12]:
(1)
式中:Pi为第i项污染物的污染指数;Ci为第i项污染物的实测浓度;Coi为第i项污染物的评价标准。单因子污染指数大于1.0,则表明水质超出规定的评价标准。
2) pH值评价指数按下式:
(2)
其中:
;
;
式中:SpH为pH的污染指数;pH为本次调查实测值;pHsu为海水pH标准的上限值;pHsd为海水pH标准的下限值。
3) 根据溶解氧(DO)的特点,溶解氧评价按下式:
(3)
(4)
式中:
;DO为溶解氧的实测浓度;DOf为饱和溶解氧的浓度;DOs为溶解氧的评价标准
值;T为水温(℃)。
2.4.2. 多项水质参数综合指数法
多项水质参数综合评价方法可较全面地反映水质整体状况。采用单因子质量指数的算术平均法计算多项水质参数综合指数,计算公式为 [13]:
(5)
式中:Qi为多项水质参数综合指数;Sij为单因子质量指数;n为单因子个数。
根据计算结果综合评价并分级:Qi < 0.75为清洁,0.75 < Qi < 1.00为轻度污染,1.00 < Qi < 1.25为中度污染,Qi > 1.25为重度污染。
2.4.3. 有机污染指数评价法
有机污染指数的计算公式为 [14]:
(6)
式中:A为有机污染指数;CODi、DINi、DIPi和DOi分别为化学需氧量、无机氮、无机磷和溶解氧的实测值;COD0、DIN0、DIP0和DO0分别为化学需氧量、无机氮、无机磷和溶解氧的海水水质标准值。
根据计算结果综合评价并分级:A < 0为良好,0 ≤ A < 1为较好,1 ≤ A < 2为开始受到污染,2 ≤ A < 3为轻度污染,3 ≤ A < 4为中度污染,A ≥ 4为重度污染。
2.4.4. 富营养化指数法
富营养化指数能较好地反映水质富营养化程度,计算公式为 [15] [16]:
(7)
式中:E为富营养化指数;COD、DIN、DIP分别为化学需氧量、无机氮、无机磷和溶解氧的实测值;
根据计算结果综合评价并分级:0 £ E < 0.5为贫营养,0.5 £ E < 1.0为中营养,1.0 £ E < 3.0为富营养,E ³ 3为高富营养。
2.4.5. 大型底栖生物评价方法
大型底栖生物群落特征指数计算公式如下 [17]:
1) 香农–韦弗(Shannon-Weaver)多样性指数:
(8)
式中:H'——种类多样性指数;n——样品中的种类总数;Pi——第i种的个体数(ni)与总个体数(N)的比值。
2) 均匀度(Pielou指数):
(9)
式中:J——均匀度;H'——种类多样性指数值;Hmax——为
多样性指数的最大值,S为样品中总种类数。
3) 物种优势度:
(10)
式中:Y——物种优势度;ni——所有测站第i种个体数;N——样品中的总个体数;fi——第i种的测站出现频率。如果Y ≥ 0.02,该种为优势种。
4)丰度(Margalef计算公式):
(11)
式中:d——丰度;S——样品中的种类总数;N——样品中的生物个体数。
2.5. 评价标准
调查海域内布设的6个站位,均位于大港滨海湿地海洋特别保护区内,要求水质符合第二类海水水质标准。海水质量状况评价参照《中华人民共和国海水水质标准》(GB3097-1997) [18],见表1。
![](Images/Table_Tmp.jpg)
Table 1. Seawater quality evaluation standard
表1. 海水水质评价标准
注:单位:mg/L,pH无量纲。
3. 评价结果
3.1. 海水质量评价
3.1.1. 单因子指数法
单因子指数法是将每个评价指标的实测浓度与评价标准进行比较来确定各个评价指标的水质类别,并以最差的水质类别作为断面水质类别 [19],能够较为直观地反映水质的超标指标,是目前最常用的水质评价方法之一。
调查海域pH变化范围在7.80~8.07之间,均值为7.91;pH最大值出现在4号站,最小值出现在1号站。监测海域所有站位海水pH均符合第一(二)类海水水质标准;按照各站位所处海洋功能区划所要求的海水水质标准进行评价,海水pH评价指数变化范围为0.23~1.00,无超标站位,见表2。
![](Images/Table_Tmp.jpg)
Table 2. Marine water quality investigation and evaluation results
表2. 海水水质调查与评价结果
注:盐度、pH无量纲;“—”表示未评价,其中盐度、总磷、总氮无评价标准,悬浮物无本底值,pH、DO不是污染因子,pH称为pH评价指数,DO称为DO评价指数。
溶解氧浓度变化范围在(6.67~7.29) mg∙L−1之间,均值为7.04 mg∙L−1;监测海域溶解氧浓度最大值出现在5号站,最小值出现在2号站。监测海域所有站位海水的溶解氧浓度均符合第一类海水水质标准。
化学需氧量变化范围在(1.90~2.89) mg∙L−1之间,均值为2.46 mg∙L−1;监测海域海水化学需氧量最大值出现在4号站,最小值出现在6号站。化学需氧量符合第一和第二类海水水质标准的站位比例分别为16.7%和83.3%;按照各站位所处的海洋功能区划所要求的海水水质标准进行评价,化学需氧量污染指数变化范围为0.63~0.96,无超标站位。
无机氮浓度变化范围在(0.354~0.439) mg∙L−1之间,均值为0.406 mg∙L−1;监测海域无机氮浓度最大值出现在2号站,最小值出现在4号站。无机氮符合第三、四类海水水质标准的站位比例分别为33.3%和66.7%;按照各站位所处的海洋功能区划所要求的海水水质标准进行评价,无机氮污染指数变化范围为1.18~1.46,超标站位比例为100%。
活性磷酸盐浓度变化范围在(0.00252~0.00393) mg∙L−1之间,均值为0.00328 mg∙L−1;监测海域活性磷酸盐最大值出现在1号站,最小值出现在5号站。监测海域海水活性磷酸盐全部符合第一类海水水质标准。按照各站位所处的海洋功能区划所要求的海水水质标准进行评价,活性磷酸盐污染指数变化范围为0.08~0.13,无超标站位。
总氮浓度变化范围在(0.00103~0.00126) mg∙L−1之间,均值为0.00111 mg∙L−1;监测海域总氮最大值出现在3号站,最小值出现在6号站。
总磷浓度变化范围在(0.0391~0.0752) mg∙L−1出现在3号站,最小值出现在2号站。
叶绿素a浓度变化范围在(5.97~17.0) µg∙L−1之间,均值为11.5 µg∙L−1;监测海域叶绿素a最大值出现在4号站,最小值出现在2号站。
2021年8月调查海域海水水质超标因子主要为无机氮,这与王娟娟 [20]、徐冠球 [21] 等在该区域的研究结果相一致;所有站位无机氮均超出第二类海水水质标准要求,符合第三、四类海水水质标准的站位比例分别为:33.3%和66.7%。除此之外,其他各项监测因子均能满足其相应的海洋功能区划要求。
3.1.2. 多项水质参数综合指数法
多项水质参数综合指数法是以单因子污染指数法为基础,将每个水质指标的实测浓度与评价标准之比作为单个污染指数,再将各单个污染指数相加计算算术平均值,进而求出综合污染指数,从而判断水体的污染程度 [19]。由表3可知,天津大港滨海湿地生态系统附近海域多项水质参数综合指数的范围为0.54~0.76,平均值为0.65。除1号站位为轻度污染(0.75 < Qi < 1.00)外,其他站位评价结果均为清洁(Qi < 0.75)。所有站位水质参数综合指数均符合第二类海水水质标准要求,总体评价结果为清洁。
![](Images/Table_Tmp.jpg)
Table 3. Evaluation results of multiple water quality parameters by comprehensive index method
表3. 多项水质参数综合指数法评价结果
3.1.3. 水质有机污染指数法
水质有机污染指数能够反映水体的有机污染状况,该方法选取COD、DIN、DIP和DO四个指标作为评估因子,侧重于对水质的有机污染状况进行评估。由表4可知,天津大港滨海湿地生态系统附近海域有机污染指数范围为0.64~1.09,平均值为0.87。除1号站位评价结果为开始污染(1 ≤ A < 2)外,其他各站位评价结果均为较好(0 ≤ A < 1)。调查海域水质有机污染指数法总体评价结果为较好。
![](Images/Table_Tmp.jpg)
Table 4. Evaluation results of organic pollution index of water quality
表4. 水质有机污染指数评价结果
3.1.4. 富营养化指数法
富营养化指数由冈市友利于1972年提出 [15],是对水体富营养化程度进行评价的指标,富营养化会造成浮游生物的过度繁殖,与海水中赤潮、绿潮的发生密切相关 [22]。由表5可知,天津大港滨海湿地生态系统附近海域有机污染指数范围为0.53~1.00,平均值为0.73。其中1号站评价结果为富营养(1.0 < E < 3.0),其他各站位评价结果均为中营养(0.5 £ E < 1.0)。调查海域富营养化指数法总体评价结果为中营养状态。
![](Images/Table_Tmp.jpg)
Table 5. Evaluation results of eutrophication index method
表5. 富营养化指数法评价结果
3.2. 大型底栖生物评价
3.2.1. 大型底栖生物种类组成
2021年8月调查海域共鉴定出大型底栖生物6大类38种(含未定名种),其中环节动物15种,占种类组成的39.5%;软体动物14种,占种类组成的36.8%;节肢动物4种,占种类组成的10.5%;螠形动物和棘皮动物各2种,均占种类组成的5.3%,纽形动物1种,占种类组成的2.6% (图2,表6)。由此可知环节动物、软体动物和节肢动物是大型底栖生物种类组成的主要部分。
![](//html.hanspub.org/file/6-2830176x23_hanspub.png?20220701091547712)
Figure 2. Percentage diagram of species composition of macrobenthos in the sea area under investigation
图2. 调查海域大型底栖生物种类组成百分比示意图
![](Images/Table_Tmp.jpg)
Table 6. List of macrobenthic species in the sea area under investigation
表6. 调查海域大型底栖生物种名录
3.2.2. 优势种
2021年8月调查海域大型底栖生物优势种及其优势度见表7,调查海域大型底栖生物优势种有6种(优势度 ≥ 0.02),其中软体动物2种,环节动物、节肢动物、棘皮动物和螠形动物各1种,具体为短吻铲荚螠、小荚蛏、长岛角螺、绒毛细足蟹、棘刺锚参和不倒翁虫。其中短吻铲荚螠优势度最高,在所有测站中出现的频率为16.7%,是本次调查海域大型底栖生物群落的第一优势种,平均密度为155.0 ind∙m−2,占总密度的13.3%;小荚蛏为第二优势种,在所有测站中出现频率为50.0%,平均密度为20.0 ind∙m−2,占总密度的5.2%;长岛角螺为第三优势种,出现频率为83.3%,平均密度为58.0 ind∙m−2,占总密度的24.9%。
![](Images/Table_Tmp.jpg)
Table 7. Dominant species and dominance of macrobenthos in the sea area under investigation
表7. 调查海域大型底栖生物优势种及其优势度
3.2.3. 种群密度和生物量分布
2021年8月调查海域大型底栖生物平均密度为194 ind∙m−2,各站位数量波动范围在(100~310) ind∙m−2之间,站位间的种群密度分布差异较小,密度最高值出现在3号站位,最低值出现在6号站位(图3);大型底栖生物平均生物量为16.2076 g∙m−2,各站位数量波动范围在(7.8060~46.5795) g∙m−2之间,站位间生物量差异较明显,调查海域大型底栖生物生物量在4号站位最高,在1号站位生物量最低(图3)。
![](//html.hanspub.org/file/6-2830176x24_hanspub.png?20220701091547712)
Figure 3. Schematic diagram of density and biomass distribution of benthos in the sea area under investigation
图3. 调查海域底栖生物密度、生物量分布示意图
3.2.4. 群落特征参数
2021年8月份底栖生物多样性指数平均为2.66,各站位波动范围在2.33~3.40之间;均匀度指数平均值为0.75,各站位波动范围在0.65~0.89之间;站位优势度均值0.61,各站位波动范围在0.44~0.70之间;丰度均值1.43,各站位波动范围在1.14~1.87之间(表8)。根据《滨海湿地生态监测技术规程》 [23] 可知,2021年8月监测海域大型底栖生物多样性指数均处于中等水平。
![](Images/Table_Tmp.jpg)
Table 8. Characteristic parameters of macrobenthos community in the investigated sea area
表8. 调查海域大型底栖生物群落特征参数
4. 结论
水质评价结果表明,2021年8月调查海域海水水质超标因子主要为无机氮,与2015年该区域的相关研究结果 [20] [21] 相一致;除此之外,其他各项监测因子均能满足其相应的海洋功能区划要求。
大型底栖生物评价结果表明,2021年8月调查海域大型底栖生物多样性指数均处于中等水平。与2015年 [24] 相比,调查海域大型底栖生物种类有所增加;优势种由单一优势种凸壳肌蛤增加为短吻铲荚螠、小荚蛏、长岛角螺等6种;大型底栖生物平均密度和平均生物量有所降低,多样性指数有所增加。
基金项目
北海局海洋科技项目(202115)。