Evaluation on carrying capacity of fishery resources in coastal waters of Dapeng Peninsula, Shenzhen
-
摘要:
为探究深圳大鹏半岛沿岸海域渔业资源承载力情况,文章参考《资源环境承载能力监测预警技术方法(试行)》和国内外渔业资源承载力评价方法,结合大鹏半岛沿岸海域实际情况,构建了大鹏半岛沿岸海域渔业资源承载力评价指标体系。基于2012—2016年每年秋季大鹏半岛沿岸海域渔业资源监测调查数据,采用2016年数据与近3年(2013—2015年)数据比较、2016年数据与近4年(2012—2015年)数据比较两种方法,对大鹏半岛沿岸海域渔业资源承载力进行评估。结果表明,两种方法所得结果一致,均为大鹏半岛沿岸海域游泳动物指数呈下降趋势、鱼卵仔稚鱼指数显著下降,渔业资源承载力已临界超载。建议进一步加强海域渔业资源监测调查和渔业资源承载力评价预警,为推进渔业资源修复保护和增殖养护提供技术支持。
Abstract:To explore the carrying capacity of fishery resources in the coastal waters of Dapeng Peninsula, Shenzhen, considering the actual situation of that sea area, we established an evaluation index system by referring to the Resource and Environmental Carrying Capacity Monitoring and Early Warning Technology Method (Trial) and the evaluation methods of fishery resources carrying capacity at home and abroad. Based on the survey data of fishery resources in that waters every autumn in 2012 to 2016, we applied two methods [comparing the data between 2016 and the previous three years (2013 to 2015) and comparing the data between 2016 and the previous four years (2012 to 2015)] to evaluate the carrying capacity. The results of both methods show that the swimming animal index had a downward trend. The fish spawn and larvae index decreased significantly. The fishery resources carrying capacity was almost overloaded. Therefore, it is recommended to strengthen the monitoring and investigation of fishery resources as well as the evaluation and early warning of fishery resource carrying capacity, so as to provide technical support for promoting the restoration and protection of fishery resources as well as proliferative maintenance.
-
鱼类运动轨迹监测技术主要有视觉监测技术和声波监测技术。视觉监测技术以视频、图片观察和手工记录为主[1-7],可以在鱼类不受干扰的情况下较好地对其行为进行观察、记录和分析,但后期数据处理难度较大、误差大且监测数据连续性差。早期的研究主要是对鱼类的行为特征进行观察和描述,基本处于“定性”分析阶段,很难得到鱼类运动行为响应的量化指标。声波监测技术分为被动声学法和主动声学法两种[8]:被动声学法是通过鱼类发出的声波特征进行探测和识别[9];主动声学法是通过声波设备接收鱼类的回波信号并进行分析和识别,具有受水域环境影响小、应用范围广泛等优势,是当今国内外研究和应用最为广泛的一种方法[10]。目前,声波监测技术应用最广泛的有鱼探仪,采用鱼探仪对水域鱼类资源的丰度[11-16]、时空分布情况 [17-21]和活动规律[22-23]等方面的研究已取得了很多成果,但该技术仅对鱼群的种类、大小和距离进行探测,无法准确定位到鱼类的位置及其运动轨迹。
声学标签系统 (Acoustic Tag System,ATS) 是声波监测技术的一种被动声学法,根据鱼类个体大小和研究周期选择合适的声学标签 (Acoustic tag,也称声学信号发射器) 类型和参数对鱼类运动轨迹进行监测[24],已逐步成为鱼类行为学研究最主要的手段。目前,国外对声学标签系统的应用研究成果主要有鱼类资源的丰度评估[25]、鱼类的游泳模式[26]、栖息地特征评价[27-29]、鱼类的产卵场地[30]、鱼类的生存情况[31]、鱼类的行为差异[32]、鱼类行为模型[33]等。我国仅有环境变化对鱼类行为的影响[34-35]和水利工程建设对鱼类洄游能力[36-37]等相关方向的研究成果报道。
声学标签系统主要用于江河湖库、海洋、大坝、河口海岸、船闸码头等的鱼类运动轨迹监测,并通过鱼类运动行为响应情况评估鱼类的丰度、种群结构、生理行为、迁徙及栖息地变化等,但应用在水生态环境对鱼类行为产生影响等方面的研究成果相对较少。因此,可将声学标签监测技术在鱼类运动行为学的应用与养殖、生态学、声学、数学模型以及计算机仿真技术等学科相结合,进行多学科的交叉利用与研究。
1. 声学标签监测技术
1.1 声学标签系统的组成
声学标签系统由硬件和软件2个部分组成。硬件部分由声学标签信号接收器、水听器、标签激活器、标签、标签检测器和电脑组成,软件部分包括Tag Programmer (用于激活及休眠标签)、Acoustic Tag (用于原始数据收集及存储) 和Mark Tags (用于原始数据文件中的环境噪音和标签回声信息处理) 3个模块。
1.2 声学标签系统的工作原理
声学标签系统采用4个水听器 (Hydrophone,简写为H1、H2、H3和H4) 接收移植或捆绑于鱼类身上的声学标签发射出来的声波,通过数据线传输到信号终端处理器进行信号处理,最后经计算机终端去噪处理后即可得到鱼类的二维、三维行为轨迹坐标。
1.3 声学标签系统的操作方法
声学标签系统的具体操作方法为:
1) 水听器布设。将水听器下端固定于2个不同平面 (图1),并设置水听器底部的三维坐标及相关参数。水听器的坐标可根据监测目的选择大地坐标或相对坐标进行设置。4个水听器监测的水域范围一般为0.5 km。如监测范围为河流等线状水域可进行多个水听器线状布设;如监测范围为大水域可根据水域情况布设多个水听器,最多可布设16个,或将水听器固定于监测船进行监测,如果监测船更换位置需要重新设置水听器坐标。
2) 采用标签激活器对声学标签进行唯一编码 (即发射频率) 激活,并用标签检测器检测是否激活成功。编码根据声波信号的发射频率进行设置,设置范围一般为1 000~6 000 Hz。可根据实际情况需要选择编码大小,如设置的发射频率为3 000 Hz即每3 s可获得1个声波信号,1个信号代表1个轨迹点。每1 h为一组数据计,一组数据有1 200个轨迹点,一年可监测到的轨迹点约1 051.2 × 104个。
3) 将已激活的声学标签捆绑或移植于试验鱼身上,水听器将接收到声学标签发射出来的声波信号并通过信号线传输至信号处理器,在PC终端便可实时监控鱼类的二维、三维运动轨迹 (图2)。
4) 采用数据处理软件Mark Tags模块对采集的原始信号进行去噪处理即可获得鱼类二维、三维运动轨迹坐标。
2. 鱼类运动轨迹数据处理方法
在数据监测过程中,除了接收到这些具有特定发射频率的标签发射的声波信号之外,还会接收到其他噪声的声波信号 (图3),或者接收的信号不连续,这些信号数据通常被称为“异常数据”。由于异常数据的存在导致监测到的原始轨迹数据杂乱无序 (图2),很难从这些庞大的试验数据中获取隐藏于其中的信息和数据变化规律。如果对这些原始数据直接进行分析,得到的结果将是不准确甚至是错误的。因此,在数据分析,首先要对这些杂乱无序的大数据进行处理,对原始的运动轨迹数据进行去噪、清洗,转换成简洁、高效的鱼类二维、三维行为轨迹数据 (图4),最后加载到数据库中。
监测到的异常数据类型主要有4类:噪声值、缺失值、异常值 (离群值) 和重复值。数据清洗过程通常是将重复和多余的数据筛选后清除、将缺失的数据补充完整、将错误的数据纠正或删除,最后整理可用的数据库,为数据挖掘和分析的准确性奠定基础。以下是根据不同异常数据类型采用不同的数据清洗方法。
2.1 噪声值
在采集声学标签声波信号的同时也会接收到其他噪声声波的信号 (图5),通过数据处理软件Mark Tags模块设置相关的去噪参数,或者添加噪声过滤器即可对采集到的原始信号进行去噪和过滤处理,得到声学标签的声波信号 (图6)。
2.2 缺失值
缺失值的清洗方法有删除法、均值插补法、热卡填补法、最近距离决定填补法、回归填补法、多重填补方法、K-最近邻法、有序最近邻法和基于贝叶斯的方法等。根据数据的缺失程度选择不同的清洗方法,其中删除法方法适用于数据比较多且缺失值的数量占整个数据的比例相对较小的情况,可以直接将缺失值删除掉,是一种最简单有效的方法。
在鱼类运动轨迹数据采集过程中,会出现由于环境客观条件造成信号无法获取而导致轨迹缺失的情况,可通过均值插补法进行插补 (图7),但某个时间段的缺失值较多则不再进行插补。如设置的声波信号发射频率为3 s,若信号中断时间超过10 s时将不再进行插补 (图8)。采用数据处理软件的Mark Tags模块过程中即可实现缺失值的自动插补。
2.3 异常值 (离群点)
异常值清洗方法主要有统计分析法、3∂原则、箱型图分析等,其中统计分析法是一种常用的简单方法。对数据库进行简单的筛选和统计,分别对鱼类运动轨迹三维坐标的x、y、z值进行排序,取数据的最大值和最小值来判断是否超出取值范围,若超出该范围则作删除处理。鱼类运动轨迹坐标应出现在水域内如a点,如果轨迹点如b点出现在水域外则显然不合常理,此为异常值,b点坐标应被删除 (图9)。
2.4 重复值
对于数据重复值的判断主要采用排序与合并的方法。重复值的检测方法主要有2个步骤:首先将数据库的数据按一定规律排序,然后通过比较相邻数据的相似情况来检测数据是否重复。对鱼类的三维运动轨迹重复值进行统计,通过判断鱼类在某处的重现性来获得鱼类的运动行为规律以及对栖息场所的偏好等信息。
要从这些庞大的数据中获得鱼类重现频率较高的轨迹坐标,需要采用数据查询语言对数据进行查询以获取有用的数据信息。数据查询语言Frequency函数是统计各区间段数值频率的一个函数,其结构为Frequency (Data_array和bins_array),其中Data_array是用于判断的数组或者数据区域,而bins_array是用于输出结果数据的分割点。其查询步骤为:输入函数Frequency→设定参数data_array和bins_array→选择输出结果的单元格区域→按F2键输入函数公式→按“Ctrl+Shift+Enter”可返回一个数组。
为获得鱼类的运动轨迹分布情况,可将轨迹坐标分别投影在xy平面和xz平面上即可得到轨迹散点图。散点图在大数据分析中的作用尤为明显,其可以展示数据的分布和聚合情况。通过散点图中散点的疏密程度和变化趋势来获得鱼类运动行为规律,轨迹越密集说明鱼类出现的频率越高。
3. 声学标签系统的应用实例
实测水域为某高校校园的景观湖,水域面积为5.3 × 103 m2,水深为1.6~2.5 m,湖内修建有景观亭 (图10)。湖内水源的主要补给来源为天然降水和校内污水处理站的中水补给,湖中有鲤 (Cyprinus carpio)、鲫 (Carassius auratus) 等湖泊常见鱼种。景观湖为封闭水域,自净能力较差,水中的氮 (N)、磷 (P) 营养元素长期积累容易导致水体的富营养化,尤其是在夏季,气温高容易导致水质恶化,严重影响景观湖的生态功能和景观效果。
由于春、夏季水域生态环境变化较大,因此对鱼类在春、夏季的运动行为轨迹进行了监测。发现春、夏两季鱼类的运动轨迹变化明显 (图11,图12):春季到夏季,在水平方向上鱼类离污染源的距离越大,垂直方向上鱼类由底层往表层迁移。引起鱼类运动行为发生变化的主要原因除了水温外,还与水中溶解氧 (Dissolved oxygen, DO) 含量和藻类死亡释放的毒素如微囊藻毒素的分布有关。夏季水温高,日照充足,加上N、P营养物质的不断输入,湖内藻类及其他浮游生物迅速繁殖并消耗水中大量的氧气 (O2) 导致水中DO含量不断下降,且得不到及时补充,致使O2收支不平衡,藻类死亡并释放衍生污染物。说明鱼类运动行为发生改变主要与水环境因子有关,通过评价鱼类运动行为与水环境因子的相关性,可为养殖水质及水生态健康评价等提供参考依据。
4. 应用前景
目前,鱼类声学标签监测技术多应用于江河湖库和海洋等鱼类丰度、种群结构等方面的研究。采用该技术对自然水域环境中鱼类的运动行为轨迹进行实时监测,可为研究水利工程建设对鱼类洄游能力、产卵以及栖息地的影响,鱼类在生境变化过程的行为响应,鱼类毒理性行为响应,鱼类行为水质监测系统建设,水生态环境健康评价以及水生态修复效果评价等提供科学的依据,是一种实时有效、快速的先进技术。
此外,根据鱼类对不同养殖水质环境的行为响应可为水产养殖提供重要的指导信息[38-39]。目前,该声学标签监测技术在水产养殖业的应用较少,具有广泛的应用前景。如通过研究鱼类的运动轨迹变化规律,掌握鱼类的栖息场所,指导人工繁殖和幼鱼培育,可以提高幼鱼的成活率和质量;此外,鱼食投放和污染物的排入会引起水质的改变,鱼类如果表现出逃避行为,且大多数都集中在洁净水的一端,表明水质遭受到污染。因此,通过连续测定鱼类的运动轨迹,对比鱼类当前运动轨迹与历史运动轨迹的变化,可实现渔业养殖水质环境的实时在线监测和预警,为提高水产养殖的产量和质量提供可靠的科学依据。
利用声学标签监测技术对鱼类运动行为轨迹进行24 h监测,获得鱼类的实时三维运动轨迹坐标并进行相关的数据分析及应用,是一种先进的技术手段。与其他监测技术相比,声学标签监测技术具有原位观察、数据处理方法简单、可24 h实时监控鱼类的二维、三维运动轨迹等优势。但该技术是通过接收鱼类身上声学标签发送的声波信号来确定鱼类的位置,在数据监测过程中可能会出现鱼类死亡、声学标签丢失和标签电量不足等情况,因此需要技术人员及时对监测数据进行实时监控和数据处理分析,以保证监测数据的连续性和准确性。
-
表 1 深圳大鹏半岛沿岸海域渔业资源承载力评价指标体系
Table 1 Evaluation system of fishery resources carrying capacity in coastal waters Dapeng Peninsula, Shenzhen
评价目标
evaluation target一级指标
primary indicator二级指数及其赋值 secondary indicator and assignment 二级指标
secondary indicator阈值范围
threshold range指示结果
indicating result赋值
assignmentF F1 ES >10% 显著下降 1 5%~10% 下降 2 ≤5% 基本稳定 3 TCPUE >30% 显著下降 1 10%~30% 下降 2 ≤10% 基本稳定 3 CTL >5% 显著下降 1 3%~5% 下降 2 ≤3% 基本稳定 3 SDI >50% 显著下降 1 25%~50% 下降 2 ≤25% 基本稳定 3 F2 FE >30% 显著下降 1 10%~30% 下降 2 ≤10% 基本稳定 3 FL >30% 显著下降 1 10%~30% 下降 2 ≤10% 基本稳定 3 表 2 深圳大鹏半岛沿岸海域游泳动物指数和鱼卵仔稚鱼指数分级评价标准
Table 2 Graded evaluation criteria for swimming animal index, and fish spawn and larvae index in coastal waters of Dapeng Peninsula, Shenzhen
一级指标
primary indicator评价阈值范围
evaluation threshold range评价结果
evaluation resultF1 F1>2.5 基本稳定 1.5<F1≤2.5 呈下降趋势 F1≤1.5 显著下降 F2 F2>2.5 基本稳定 1.5<F2≤2.5 呈下降趋势 F2≤1.5 显著下降 表 3 深圳大鹏半岛沿岸海域渔业资源综合承载指数分级评价标准
Table 3 Graded evaluation criteria for comprehensive carrying index of fishery resources in coastal waters of Dapeng Peninsula, Shenzhen
评价阈值范围
evaluation threshold range评价结果
evaluation resultF>2.5 海洋渔业资源可载 1.5<F≤2.5 海洋渔业资源临界超载 F≤1.5 海洋渔业资源超载 -
[1] 张露, 袁华荣, 佟飞, 等. 广东深圳大鹏半岛海域秋季渔业资源群落结构分析[J]. 安徽农业科学, 2018, 46(13): 89-93. [2] 陈涛, 林金錶, 郭金富, 等. 大亚湾真鲷资源状况研究[J]. 热带海洋学报, 2003, 22(3): 30-35. doi: 10.3969/j.issn.1009-5470.2003.03.005 [3] 霍军. 近海渔业资源承载力的特征及评估指标体系初探[C]//2009中国海洋论坛论文集. 青岛: 中国海洋大学出版社, 2009: 251-258. [4] JENSEN A L. Maximum harvest of a fish population that has the smallest impact on population biomass[J]. Fish Res, 2002, 57(1): 89-91. doi: 10.1016/S0165-7836(01)00337-X
[5] SHEPHARD S, BROPHY D, REID D G. Can bottom trawling indirectly diminish carrying capacity in a marine ecosystem?[J]. Mar Biol, 2010, 157(11): 2375-2381. doi: 10.1007/s00227-010-1502-9
[6] 韩增林, 狄乾斌, 刘锴. 辽宁省海洋水产资源承载力与可持续发展探讨[J]. 海洋开发与管理, 2003, 20(2): 52-57. doi: 10.3969/j.issn.1005-9857.2003.02.017 [7] 杨洋, 刘志国, 何彦龙, 等. 基于非平衡产量模型的海洋渔业资源承载力评估——以浙江省为例[J]. 海洋环境科学, 2016, 35(4): 534-539. [8] PERRY R I, SCHWEIGERT J F. Primary productivity and the carrying capacity for herring in NE Pacific marine ecosystems[J]. Prog Oceanogr, 2008, 77(2/3): 241-251.
[9] 叶孙忠, 罗冬莲, 杨芳, 等. 东山湾渔业资源承载力评价指标体系构建及评估[J]. 海洋环境科学, 2018, 37(4): 493-498. [10] 刘尊雷, 袁兴伟, 杨林林, 等. 基于交通灯方法的小黄鱼资源承载力诊断与评价[J]. 海洋环境科学, 2018, 37(4): 521-527. [11] 陈丕茂, 袁华荣, 贾晓平, 等. 大亚湾杨梅坑人工鱼礁区渔业资源变动初步研究[J]. 南方水产科学, 2013, 9(5): 100-108. doi: 10.3969/j.issn.2095-0780.2013.05.016 [12] 贾春斌, 庄世鹏. 深圳市东部沿海渔业资源现状与评估[J]. 中国水产, 2009(8): 21-23. doi: 10.3969/j.issn.1002-6681.2009.08.013 [13] 余景, 毛江美, 袁华荣, 等. 深圳鹅公湾水域渔业资源季节变动分析[J]. 南方农业学报, 2017, 48(6): 1106-1112. doi: 10.3969/j.issn.2095-1191.2017.06.27 [14] 赵漫, 余景, 陈丕茂, 等. 深圳鹅公湾渔业水域生态系统健康状况评价[J]. 南方农业学报, 2016, 47(6): 1025-1031. doi: 10.3969/j:issn.2095-1191.2016.06.1025 [15] 国家发展和改革委员会.资源环境承载能力监测预警技术方法(试行)[Z]. 发展规划[2016] 2043号, 2016: 23-25. [16] PAULY D, CHRISTENSEN V V, DALSGAARD J, et al. Fishing down marine food webs[J]. Science, 1998, 279(5352): 860-863. doi: 10.1126/science.279.5352.860
[17] 张其永, 林秋眠, 林尤通, 等. 闽南-台湾浅滩渔场鱼类食物网研究[J]. 海洋学报, 1981, 3(2): 275-290. [18] REYNOLDS J F, LUDWIG J A. Statistical ecology: a primer on methods and computing[M]. New York: John Wiley&Sons, 1988: 1-337.
[19] 夏一璐, 陈琼, 赵荣磊, 等. 伏休前后舟山渔场单拖船低值杂鱼渔获物组成比较分析[J]. 浙江海洋学院学报(自然科学版), 2015, 34(6): 520-525, 558. doi: 10.3969/j.issn.1008-830X.2015.06.004 [20] 朱国平, 张衡, 王家樵, 等. 大西洋金枪鱼渔业平均营养级的长期变动[J]. 生态科学, 2009, 28(2): 97-101. doi: 10.3969/j.issn.1008-8873.2009.02.001 [21] 万瑞景, 孙珊. 黄、东海生态系统中鱼卵、仔稚幼鱼种类组成与数量分布[J]. 动物学报, 2006, 52(1): 28-44. doi: 10.3969/j.issn.1674-5507.2006.01.004 [22] 赵静, 张秀梅, 卞晓东, 等. 年葫芦岛附近海域鱼卵仔稚鱼种类组成与数量分布[J]. 中国海洋大学学报(自然科学版), 2011, 41(11): 34-42. [23] MOUILLOT D, LAUNE J, TOMASINI J A, et al. Assessment of coastal lagoon quality with taxonomic diversity indices of fish, zoobenthos and macrophyte communities[J]. Hydrobiologia, 2005, 550(1): 121-130. doi: 10.1007/s10750-005-4368-y
[24] 乔延龙, 林昭进, 邱永松. 北部湾秋、冬季渔业生物群落结构特征的变化[J]. 广西师范大学学报(自然科学版), 2008, 26(1): 100-104. doi: 10.3969/j.issn.1001-6600.2008.01.025 [25] 杨正先, 张志锋, 韩建波, 等. 海洋资源环境承载能力超载阈值确定方法探讨[J]. 地理科学进展, 2017, 36(3): 313-319. [26] 郭建忠, 陈作志, 许友伟, 等. 大亚湾鱼类资源数量的时空分布特征[J]. 中国海洋大学学报: 自然科学版, 2018, 48(12): 47-55. [27] 农业部渔业局. 中国渔业统计年鉴[M]. 北京: 中国农业出版社, 2016: 44-47. [28] 孙才志, 曹威威, 肖春柳. 长山群岛海洋渔业资源的平均营养级变化特征[J]. 海洋通报, 2019, 38(1): 87-95. [29] 高彦洁, 吕振波, 杨艳艳, 等. 莱州湾春季鱼卵仔稚鱼群落结构和物种多样性[J]. 生态学报, 2016, 36(20): 6565-6573. [30] 罗冬莲. 悬浮物对鱼卵仔稚鱼的影响分析及其损失评估——以厦漳跨海大桥工程为例[J]. 海洋通报, 2010, 29(4): 439-443. doi: 10.3969/j.issn.1001-6392.2010.04.015