Study on the ecosystem model of Daya Bay I. A preliminary approach on energy flow model
-
摘要:
主要根据1984~1986年和1986~1987年在大亚湾进行的环境、资源和生态调查资料,应用Ecopath with Ecosim(EwE)软件,构建大亚湾海域生态系统初步能量流动模型。文中根据大亚湾游泳动物的食物组成特点,把该海域生态系划分15个功能组,分别是海洋哺乳动物、肉食性鱼类、底栖捕食鱼类、滤食性鱼类、草食性鱼类、蟹类、虾类、头足类、底栖动物、水母、浮游动物、珊瑚、沉水植物、浮游植物和有机碎屑,功能组的划分基本能覆盖大亚湾海域生态系统的能量流动过程。经EwE软件模拟,结果表明:大亚湾海域生态系统的营养级范围为1~3.88级;各营养级的能量转换效率分别为7.2%,11.2%,8.7%,2.9%,可用构建金字塔形状来描述营养流动的转换效率;大亚湾生态系统的总能量传递效率为8.9%,略低于林德曼转换效率(10%左右),可能是由于在该海域大量的沉水植物(马尾藻)未能被充分利用而腐烂所造成;在能量流动过程中,直接来源于碎屑的比例占总流量的48%,而直接来源于初级生产者的比例为52%。
-
关键词:
- 营养模型 /
- 能量流动 /
- Ecopath with Ecosim软件 /
- 生态系统 /
- 大亚湾
Abstract:Ecopath with Ecosim (EwE) is designed for straightforward construction, parameterization and analysis of mass-balance trophic models of aquatic ecosystems. Base on the data from environments, resources and ecology survey in Daya Bay during 1984~1986 and 1986~1987, the energy flow model of ecosystem in Daya Bay is constructed using the EwE package. According to the characteristic of the food composition of swimming animal in Daya Bay, the ecosystem comprised 15 function groups, which are marine mammals, carnivorous fish, benthic-feeding fish, zooplanktivorous, herbivorous fish, crabs, shrimps, cephalopods, zoobenthos, jellyfish, zooplankton, corals, submerged plant, phytoplankton and detritus, respectively. And the function groups can cover the main trophic flow of the ecosystem in Daya Bay. The result shows that the trophic level of the function groups varied from 1.0 to 3.88. The transfer efficiencies for each trophic level are 7.2%, 11.2%, 8.7%, 2.9%, respectively. The transfer efficiencies can be used for constructing a figure presenting the trophic flows in form of a pyramid. The total transfer efficiency of Daya Bay′s ecosystem is 8.9%, less than that of Lindeman (about 10%), perhaps a lot of submerged plant (Sargassum) in Daya Bay can not be utilized fully, and decayed and broke down. In course of the energy flow, the proportion of total flow originating from detritus is 48%, and from primary producer is 52%.
-
Keywords:
- trophic model /
- energy flow /
- Ecopath with Ecosim package /
- ecosystem /
- Daya Bay
-
硫酸盐还原菌(sulfate-reducing bacteria,SRB)在生态系统中是属于土著微生物类群,是兼性厌氧菌,广泛存在于各种环境中,是自然界硫循环中SO42-异化还原的专门承担者,将SO42-还原同时降解有机物获取能量,产物为对周围环境有害的有毒气体H2S。目前,国内外对海洋沉积物中SRB研究多集中在海洋工程、石油勘探、管道铺设等金属防腐方面[1-2]。近年来因病害使国内外海水养殖产量出现大幅度的滑坡,其主要原因是养殖环境恶化,尤其是底质的污染严重,大部分养殖海区投喂天然饵料,饵料系数较高,养殖年限较长,网箱底部堆积大量富含有机物的沉积物,为SRB大量繁殖提供了有利条件,SRB繁殖会产生大量有毒气体H2S,造成水体污染并毒害养殖生物,使养殖水体环境状况进入恶性循环状态。目前,国内外对网箱养殖环境中SRB的相关研究报道不多[3],因此,开展网箱养殖海区沉积物中SRB生态特性和相关防治是很有必要的。本文通过研究大鹏澳网箱养殖海域的SRB,了解此环境中SRB的生态特性、SRB数量变化,以及SRB与其他环境因子之间的关系,为以后修复养殖水体环境,提高养殖生物产量,实现海水网箱养殖业的可持续发展提供理论依据。
大鹏澳是大亚湾西南部的一个小内湾,水面面积约1 400 km2。20世纪90年代以来,海水网箱养殖得到了迅速发展。目前,网箱养殖区水面面积约有20 km2,养殖网箱约4 200箱,养殖区平均水深4.5 m,污染相对比较严重[4-5]。
1. 材料与方法
1.1 样品采集
2006年11月,在大亚湾大鹏澳网箱养殖海域利用GPS定位5个站点,其中1、2、3号站位于鱼类网箱养殖区,4号站位于浮筏贝类养殖区,5号站作为对照站,采样站点如图 1所示。在每个站用柱状采泥器采集柱状沉积物样品2管,将沉积物置于内径为5 cm的PVC管内并用橡胶塞密封好,其中1管用于测定硫化物、氧化还原电位(oxidation-reduction potential,ORP)和pH,另1管用于测定SRB含量,并清晰标记样品垂直方向,其中柱状沉积物SRB计数按照1 cm层、5 cm层、10 cm层深度进行取样。
1.2 SRB培养基配方和计数方法
培养基配方为K2HPO4 0.5 g · L-1,NH4Cl 1.0 g · L-1,MgSO4 2.0 g · L-1,Na2SO4 0.5 g · L-1,CaCI2 0.1 g · L-1,酵母膏1.0 g · L-1,乳酸钠4 mL ·L-1依次加入净化陈海水中溶解。调节pH值约为7.0~7.5[6]。密封后用蒸气压力灭菌器121~125℃灭菌15~20 min后冷却至室温。称取硫酸亚铁铵1.2 g,抗坏血酸0.4 g,在无菌箱(室)内均匀地摊在离紫外线灯30 cm处灭菌30 min。在无菌操作下,把硫酸亚铁铵和维生素C溶解于事先准备好的40 mL无菌水中。按每100 mL培养基各加入1.0 mL硫酸亚铁铵溶液和1.0 mL维生素C溶液。
硫酸盐还原菌的计数(本标准适用于工业循环冷却水中硫酸盐还原菌的测定,也适用于原水、生活用水及粘泥中硫酸盐还原菌的测定)采用我国国家标准《工业循环冷却水中硫酸盐还原菌的测定》(GB/T14643.5-1993)规定的MPN法[7]。
1.3 沉积物硫化物、ORP、pH的测定
1.3.1 硫化物的测定
沉积物中硫化物的测定采用《海洋监测规范》[8]中的硫离子选择电极法。取5.00 g混匀的湿样置于50 mL烧杯中,加入20 mL抗氧化络合剂使用液,充分搅拌5 min,静置,待沉积物完全沉淀后,用倾斜法将上清液倒入50 mL量瓶中;再重复2次用15 mL抗氧化络合剂使用液按上述步骤浸取残留的沉积物,合并3次上清液,再用抗氧化络合剂使用液定容至刻度。测定时将上清液移入50 mL烧杯中,插入参比电极和活化后的硫离子选择电极,一边搅拌至读数稳定后记录结果。每次测定后均需用蒸馏水洗涤电极并用滤纸吸干,以备下一样品测定。
1.3.2 ORP的测定
沉积物ORP的测定参照《海洋监测规范》[8]中的电位计法。
1.3.3 pH的测定
沉积物pH的测定按照《海洋监测规范》[8]规定采用复合电极方法。
上述分析仪器为上海雷磁PHSJ-4A实验室pH/ORP测定仪,配有硫离子选择电极、参比电极、氧化还原电极、感温电极和pH复合电极。
2. 数据处理
利用Microsoft Excel对SRB的检出率、含量、变化幅度等按照不同站位、不同层次进行分析,利用数理统计软件SPSS对SRB与硫化物、ORP、pH做相关分析。
3. 结果与讨论
3.1 沉积物中的SRB的数量
2006年11月大亚湾大鹏澳网箱养殖海域沉积物中SRB的数量如表 1所示,沉积物中SRB数量波动范围在900~110 000 ind · g-1之间。
表 1 大鹏澳网箱养殖海域沉积物中SRB的数量Table 1. Content of SRB in cage culture area sediment in Dapeng′ao Coveind · g-1 垂直深度/cm
vertical depth1号站
station 12号站
station 23号站
station 34号站
station 4对照站
control station1 45 000 45 000 110 000 25 000 30 000 5 30 000 4 500 45 000 15 000 15 000 10 11 500 2 500 4 500 900 4 500 图 2显示了各个站位表层沉积物SRB的分布情况,SRB平面分布特征是鱼类网箱区>对照区>浮筏贝类区,网箱区的数量明显高于对照区和浮筏贝类区,这可能是由于网箱区内的网箱设置比较密集,饵料主要是冰鲜小杂鱼,残饵和粪便形成的生物沉积,为SRB的生长提供了有利的环境;浮筏贝类区SRB数量最少,可能是由于此区域的养殖生物为太平洋牡蛎,无需投饵,故底部富含有机质的沉积物较少,SRB含量也较少。
各个站位沉积物中SRB的垂直分布如图 3所示,可以明显看出SRB数量的分布特征为1 cm层>5 cm层>10 cm层,呈现由上到下逐渐减少的趋势。验证了海洋沉积物中硫酸盐还原菌数量是由上向下逐渐降低的结论[9],这也与高爱国等[10]在北极楚科奇海调查得到的SRB分布规律相似。
3.2 沉积物中的硫化物
2006年11月大亚湾大鹏澳网箱养殖海域5个站位沉积物中硫化物的含量如表 2所示,其硫化物含量波动范围为372.48~703.98 μg · g-1。
表 2 表层沉积物中SRB与其他环境因子Table 2. The amount of SRB and other factors of surface sediment站位
stationSRB数量
content of SRB硫化物含量/μg·g-1
sulfide content氧化还原电位
ORPpH 1号站 station 1 45 000 529.44 -226.5 8.12 2号站 station 2 45 000 527.76 -270.3 8.18 3号站 station 3 110 000 703.98 -384.5 8.16 4号站 station 4 25 000 372.48 -171.5 8.21 对照站 control station 30 000 445.18 -199.5 8.17 由图 4可以明显看出,鱼类网箱养殖区(1、2、3号站)沉积物中硫化物含量要高于贝类区和对照区,其主要原因是鱼类网箱区需要大量投饵,底部富含有机质,这一环境利于SRB大量繁殖,异化还原硫酸盐形成硫化物;浮筏贝类养殖区的硫化物含量最低,这是因为贝类养殖过程中不需投饵,底部沉积物有机质含量较少;另外,该海域贝类区养殖年限较短,水体底部的溶解氧和pH值均较高,不利于SRB的生长,因此,硫化物含量较低。而对照区的5号站,其硫化物含量相对较高于贝类养殖区,这可能是因为5号站位于大鹏澳湾口,在海流作用下,受网箱区沉积物往湾口外迁移的影响要大于贝类区有关。
3.3 沉积物中的ORP
各采样站位沉积物的氧化还原电位如表 2所示,其波动范围在-171.5~-384.5 mV之间。图 5显示了2006年11月大亚湾大鹏澳网箱养殖海域5个站位沉积物的ORP平面分布特征。可以看出,鱼类网箱区的ORP要低于贝类区和对照区,其主要原因是由于网箱区残饵和排泄粪便长期积累,网箱区底部富含有机质,微生物分解大量有机质,产生还原性的物质(如:S2-、NH4+、Fe2+等),从而导致氧化还原电位较低。
3.4 沉积物中的pH
图 6显示了2006年11月大亚湾大鹏澳网箱养殖海域各个站位沉积物的pH值,其特征为网箱区低于对照区和贝类区,影响沉积物pH变化的主要因素是呼吸作用[11],因为网箱区沉积物的有机质含量高于对照区和贝类区,所以网箱区沉积物中的细菌呼吸作用和有机质在细菌参与下的分解作用相对活跃,导致pH下降。
3.5 SRB与其他因子的相关关系
2006年11月大亚湾大鹏澳网箱养殖海域表层沉积物中SRB含量、硫化物含量、ORP、pH值如表 2所示。利用SPSS对SRB与其他环境因子所做的相关分析结果显示,SRB与硫化物的含量为显著性正相关(R=0.96, P < 0.05,表 3),这与李培英等[12]研究的浙江至闽北陆架积物硫酸盐还原菌与硫化物的相关关系的结果相一致。由于SRB在沉积物环境中主要电子受体是硫酸盐,当SRB大量存在时就会有更多的还原态的硫化物生成。SRB含量与Eh的关系显示极显著性负相关关系(R=-0.97,P < 0.01,表 3),这与陈皓文等[13]对北部湾东侧沉积物SRB研究结果相一致。ORP的高低主要取决于沉积物中的DO含量,当沉积物处于严重缺氧时,有机物最终分解为还原态物质,Eh为负值,而厌氧的低氧化还原状态环境适合SRB的大量繁殖,所以这一分析结果与理论相一致。而SRB与pH的关系系数为R=-2.71(P>0.05,表 3),两者不具有显著性相关关系,由此可以看出,硫化物含量和ORP是反映SRB数量的2个非常重要的指标。
表 3 SRB与其他因子的相关系数Table 3. Correlation between SRB and other chemical factors对子 parallelism 相关系数R correlation coefficient n number of sample P significance factor SRB-硫化物 SRB-sulfide 0.96 5 0.011 SRB-氧化还原电位 SRB-ORP -0.97 5 0.006 SRB-pH -2.71 5 0.611 -
表 1 大亚湾生态系统模型的基本输入和输出参数
Table 1 Basic input and output parameters of the trophic model of Daya Bay ecosystem
功能组
group name营养级
trophic level生物量/t·km-2
biomassP/B
(year-1)Q/B
(year-1)EE P/Q 海洋哺乳动物 marine Mammals (3.88) 0.008 0.045 14.768 (0.000) (0.003) 肉食性鱼类 carnivorous fish (3.05) 0.060 1.500 20.000 (0.144) (0.075) 底栖捕食鱼类 benthic-feeding fish (2.99) 0.100 1.000 2.700 (0.956) (0.370) 滤食性鱼类 zooplanktivorous fish (2.62) 0.230 1.900 5.000 (0.462) (0.380) 草食性鱼类 herbivorous fish (2.05) 0.147 2.300 6.500 (0.898) (0.354) 蟹类 crabs (2.77) 0.0208 1.200 11.600 (0.899) (0.103) 虾类 shrimps (2.36) 0.100 1.800 19.000 (0.806) (0.095) 头足类 cephalopods (3.08) 0.0184 3.100 11.970 (0.849) (0.259) 底栖动物 zoobenthos (2.15) 0.166 6.500 27.400 (0.800) (0.237) 水母 jellyfish (3.00) (0.0034) 5.011 25.050 0.95 (0.200) 浮游动物 zooplankton (2.00) 0.186 32.000 192.000 (0.351) (0.167) 珊瑚 corals (2.00) (0.0674) 1.090 9.000 0.95 (0.121) 沉水植物 submerged plant (1.00) 17.8 11.885 - 0.002 - 浮游植物 phytoplankton (1.00) 0.119 231.00 - (0.971) - 有机碎屑 detritus (1.00) 1.000 - - (0.071) - 注:括号内为输出参数
Note: Output parameters are bracketed.表 2 不同功能组的相对能量流动
Table 2 Relative flows of groups of trophic level
功能组
group name营养级 trophic level Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ Ⅷ 海洋哺乳动物 marine mammals 0.000 0.000 0.227 0.694 0.076 0.003 0.000 0.000 肉食性鱼类 carnivorous fish 0.000 0.110 0.748 0.126 0.016 0.000 0.000 0.000 底栖捕食鱼类 benthic-feeding fish 0.000 0.250 0.518 0.226 0.006 0.000 0.000 0.000 滤食性鱼类 zooplanktivorous fish 0.000 0.400 0.584 0.016 0.001 0.000 0.000 0.000 草食性鱼类 herbivorous fish 0.000 0.950 0.050 0.000 0.000 0.000 0.000 0.000 蟹类 crabs 0.000 0.250 0.728 0.022 0.000 0.000 0.000 0.000 虾类 shrimps 0.000 0.670 0.318 0.012 0.000 0.000 0.000 0.000 头足类 cephalopods 0.000 0.000 0.967 0.033 0.000 0.000 0.000 0.000 底栖动物 zoobenthos 0.000 0.938 0.062 0.000 0.000 0.000 0.000 0.000 水母 jellyfish 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 浮游动物 zooplankton 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 珊瑚 corals 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 沉水植物 submerged plant 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 浮游植物 phytoplankton 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 有机碎屑 detritus 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 表 3 不同营养级的能量转化效率
Table 3 Transfer efficiencies for each trophic level of the model in Daya Bay
来源 source 营养级 trophic level Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ 生产者 producer 6.8 10.9 8.8 3.3 碎屑 detritus 8.0 11.5 8.5 2.4 总能流 all flows 7.2 11.2 8.7 2.9 0.8 0.0 碎屑所占的能流比 proportion of total flow originating from detritus: 0.48 转换效率 transfer efficiencies 生产者转换效率 from primary producer: 8.7% 碎屑转换效率 from detritus: 9.3% 总转换效率 total: 8.9% -
[1] 广东省海岛资源综合调查大队等. 大亚湾海岛资源综合调查报告[R]. 广州: 广东科技出版社, 1992. https://www.las.ac.cn/front/book/detail?id=5997b7d1ba8db880c2b480fd2e65173d [2] 国家海洋局第三研究所. 大亚湾海洋生态文集(Ⅰ)[C]. 北京: 海洋出版社, 1989. https://max.book118.com/html/2017/0528/109990200.shtm [3] 国家海洋局第三研究所. 大亚湾海洋生态文集(Ⅱ)[C]. 北京: 海洋出版社, 1990. https://max.book118.com/html/2017/0528/109994712.shtm [4] 徐恭昭等. 大亚湾环境与资源[M]. 安徽: 安徽科学技术出版社, 1989. https://www.dushu.com/book/10300684/ [5] 贾晓平, 郭金富, 林钦, 等. 大型石化综合工程对大亚湾渔业生态环境与水生生物资源影响的预测评价[A]. 见: 贾晓平主编. 南海渔业生态环境与生物资源的污染效应研究[C]. 北京: 海洋出版社, 2004. [6] 张其永, 林秋眠, 林尤通, 等. 闽南-台湾浅滩渔场鱼类食物网研究[J]. 海洋学报, 1981, 3(2): 275-290. https://www.cnki.com.cn/Article/CJFDTotal-SEAC198102007.htm [7] 邓景耀, 孟田湘. 渤海鱼类食物关系[J]. 海洋水产研究, 1988, 9: 151-172. https://www.zhangqiaokeyan.com/academic-journal-cn_acta-ecologica-sinica_thesis/02012109842958.html [8] 韦晟, 姜为民. 黄海鱼类食物网研究[J]. 海洋与湖沼, 1992, 23(2): 182-191. https://www.cqvip.com/QK/90072X/19922/770477.html [9] GB 17378-1998, 海洋监测规范[S]. 北京: 标准出版社, 1998. [10] 国家技术监督局. GB 12763.6-91, 海洋调查规范[S]. 北京: 标准出版社, 1991. https://d.wanfangdata.com.cn/periodical/zgbzh-e201102008 [11] Polovina J J. Model of a coral reef ecosystem, PartⅠ: the ECOPATH model and its application to French Frigate Shoal[J]. Coral Reefs, 1984, 3(12): 1-11. doi: 10.1007/BF00306136
[12] Ulanowicz R E. Growth and development: Ecosystem phenomenology[M]. Springer-Verlag, New York, 1986. doi: 10.1007/978-1-4612-4916-0
[13] Christensen V, Pauly D. EcopathⅡ- a software for balancing steady-state ecosystem models and calculating network characteristics[J]. Ecol Modeling, 1992, 61(3-4): 169-185. doi: 10.1016/0304-3800(92)90016-8
[14] Christensen V, Pauly D. A guide to the ECOPATHⅡ program[M]. ICLARM software 6.72, 1992. http://www.mendeley.com/research/guide-ecopath-ii-software-system-version-21/
[15] Pauly D, Christensen V, Walters C. Ecopath, Ecosim, and ecospace as tools for evaluating ecosystem impact of fisheries[J]. ICES J Mar Sci, 2000, 57(3): 697-706. doi: 10.1006/jmsc.2000.0726
[16] Gulland J A. Fish stock assessment: a manual of basi methods[M]. John Wiley and Sons, New York, 1983. https://www.semanticscholar.org/paper/Fish-Stock-Assessment.-A-Manual-of-Basic-Methods-Elliott-Gulland/96fd586f5fbff273eeb3b471b15767443beccd07
[17] Pauly D. Fish population dynamics in tropical waters: a manual for use with programmable calculators[M]. ICLARM Stud Rew, 1984. 325. doi: 10.1016/0165-7836(86)90044-5
[18] Lafontaine J, Peters R H. Empirical relationships for marine primary production: the effect of environmental variables[J]. Oceanol Acta, 1986, 9(1): 65-72. https://www.researchgate.net/publication/230759530_Empirical_relationship_for_marine_primary_production_the_effect_of_environmental_variables
[19] Pauly D, Palomares M L. Shrimp comsumption by fish in Kuwait waters: a methodology, preliminary results and their implications for management and research[J]. Kuwait Bull Mar Sci, 1987, 9: 101-125. https://www.researchgate.net/publication/291754509_Shrimp_consumption_by_fish_in_Kuwait_waters_a_methodology_preliminary_results_and_their_implications_for_management_and_research
[20] Pitcher T J, Watson R, Courtney A, et al. Assessment of Hong Kong′s inshore fishery resources[R]. Fisheries Centre Research Reports, 1998, 6(1): 1-78. https://open.library.ubc.ca/collections/52383/items/1.0074776
[21] Pitcher T J, Buchary E A, Trujillo P. Spatial simulations of Hong Kong′s marine ecosystem: ecological and economic forecasting of marine protected areas with human-made reefs[R]. Fisheries Center Research Reports, 2002, 10(3): 1-72. https://www.researchgate.net/publication/276204667_Spatial_Simulations_of_Hong_Kong%27s_Marine_Ecosystem_Ecological_and_Economic_Forecasting_of_Marine_Protected_Areas_with_Human-Made_Reefs
[22] Lin H-J, Shao K-T, Kuo S-R, et al. A trophic model of a sandy barrier lagoon at Chiku in Southwestern Taiwan[J]. Est Coast Shelf Sci, 1999, 48(5): 575-588. doi: 10.1006/ecss.1998.0457
[23] Christensen V, Carl J, Pauly D. Ecopath with Ecosim: A use′s guide[M]. ICLARM, 2000. https://s3-us-west-2.amazonaws.com/legacy.seaaroundus/researcher/dpauly/PDF/2000/OtherItems/ECOPATH_WITH_ECOSIM_A_USERS_GUIDE.pdf
[24] Pauly D, Christensen V, Dalsgaard J, et al. Fishing down marine food webs[J]. Science, 1998, 279(5352): 860-863. doi: 10.1126/science.279.5352.860
[25] Tong L, Tang Q S. A preliminary approach on mass-balance ecopath model of the Bohai Sea[J]. Chinese J Appl Eco, 2000, 11(3): 435-440. https://pubmed.ncbi.nlm.nih.gov/11767650/
[26] 苏纪兰, 唐启升. 中国海洋生态系统动力学研究Ⅱ. 渤海生态系统动力学过程[M]. 北京: 科学出版社, 2002.212-312. https://www.zhangqiaokeyan.com/book-cn/081501282193.html