Pattern of fish community and its relationship with environmental factors in Fangchenggang (Qinzhou coastal area of Beibu Gulf)
-
摘要: 依据2016年北部湾防城港—钦州近岸海域的底拖网调查数据，分析了鱼类群落组成和资源的时空分布；同时结合温度、盐度、溶解氧、pH、深度5种水环境要素，探讨了该海域鱼类群落与环境因子间的关系。结果表明，该海域共捕获鱼类152种，隶属于15目55科96属，二长棘鲷 (Parargyrops edita)、鹿斑鲾 (Secutor ruconius)、多齿蛇鲻 (Saurida tumbil)、褐菖鲉 (Sebastiscus marmoratus) 等为优势种。资源量以春季最丰富，多样性指标春、夏季较高，冬季最低。聚类和NMDS分析得出，该海域可划分为4个鱼类组群，组内平均相似度介于18.27% (组群I) 至32.49% (组群II)。ANOSIM分析表明各组群间鱼类组成差异极显著 (R=0.703~0.982, P<0.01)，组群I与IV间存在最大相异性 (98.02%)。冗余分析 (RDA) 表明，底层温度、溶解氧、深度和底层盐度是影响北部湾防城港—钦州近岸海域鱼类分布的主要环境因子。Abstract: We analyzed the fish assemblage composition and spatio-temporal distribution of fish stocks based on the bottom trawl survey data in Fangchenggang (Qinzhou coastal area of Beibu Gulf) in 2016. Combining with the information of temperature, salinity, dissolved oxygen, pH and depth, we investigated the relationship between fish community structure and environmental factors. A total of 152 fish species belonging to 15 orders, 55 families and 96 genera had been collected, and Parargyrops edita, Secutor Ruconius, Saurida tumbil and Sebastiscus marmoratus were the dominant species. The abundance of fish stocks was the highest in spring. The diversity index was higher in spring and summer, but lower in winter. Four distinct groups were defined by Cluster and non-metric multidimensional scaling (NMDS) analysis, and the average similarity of the groups ranged from 18.27% (Group I) to 32.49% (Group II). The analysis of similarity (ANOSIM) shows that the differences of fish composition among different groups were extremely significant (R=0.703−0.982, P<0.01), with the most sinificant difference of 98.02% between Group I and Group IV. The results of redundancy analysis (RDA) show that the bottom temperature, dissolved oxygen, depth and bottom salinity of bottom were the main environmental factors that affected fish distribution in Fangchenggang.
-
Keywords:
- Fish composition /
- Community structure /
- Species diversity /
- Environmental factors /
- Beibu Gulf
-
南太平洋长鳍金枪鱼 (Thunnus alalunga) 主要生活在10°S—30°S之间的热带和亚热带水域，渔业资源丰富。长鳍金枪鱼是延绳钓渔业主要捕捞目标物种，占金枪鱼每年渔获量的1/3，近年来产量仍在增加[1-2]。目前国内南太平洋长鳍金枪鱼渔业资源研究多基于我国渔捞日志数据处理得到的多尺度网格聚合数据，或中西太平洋渔业委员会提供的5°方格数据[3-4]，且多使用单一尺度数据进行研究[5-6]。鱼类以个体、种群、群落的形式分布于特定空间，具有高度的空间自相关性，不同研究尺度下，其空间格局差异较大，而这种不同空间尺度下空间自相关性得出的结果往往会自相矛盾，甚至影响到渔业资源评估[3]。因此，研究不同空间尺度在空间格局上的差异，以及空间数据的收集差异导致我国大陆数据和中西太平洋渔业委员会数据的空间格局差异，对长鳍金枪鱼渔业资源的合理利用和保护有重要意义。
本研究选取2015—2019年中国大陆地区 (Chinese Flag, CN) 在南太平洋的渔捞日志数据，并按照不同空间尺度聚合。此外，对比了中西太平洋渔业委员会 (Western Center Pacific Fisheries Commission, WCPFC) 下载的长鳍金枪鱼延绳钓渔业数据，运用标准差椭圆、莫兰指数和热点分析等方法分析了长鳍金枪鱼的空间分布格局，并对结果进行了对比，以探明南太平洋长鳍金枪鱼渔业资源的空间格局特征，确定是否存在空间尺度差异，以期为南太平洋长鳍金枪鱼渔业资源管理提供科学的理论依据。
1. 材料与方法
1.1 数据来源
本研究数据为2015—2019年每年1—12月中国大陆地区所有延绳钓渔捞日志数据，及同时期和同区域WCPFC发布的5°网格延绳钓渔业数据 (表1，由中国远洋渔业数据中心提供)。包括作业位置 (经纬度)、投钩数量、渔获量和作业时间等。本文选择的研究区域为南太平洋150°E—130°W，0—50°S，见图1中阴影区范围，是南太平洋长鳍金枪鱼主要出现海域[7]。
表 1 南太平洋长鳍金枪鱼延绳钓渔业船次Table 1. Number of vessels of longline fishery of T. alalunga in South Pacific年份Year 船次 (CN)Number of vessels 船次 (WCPFC)Number of vessels 2015 189 1 752 2016 138 1 695 2017 159 1 764 2018 150 1 764 2019 156 1 817 将中国大陆的渔业数据按季度处理成空间分辨率为1°×1° (后文简称为CN-1)、2°×2° (后文简称为CN-2) 和5°×5° (后文简称为CN-5)，中西太平洋渔业委员会数据按季度处理成空间分辨率为5°×5° (后文简称为WCPFC-5)。用单位捕捞努力量渔获量 (CPUE) 表征长鳍金枪鱼渔业资源分布状况：
$$Y_{\rm{CPUE}} =\dfrac{{U}_{{\rm{catch}}}}{{f}_{{\rm{hooks}}}} $$ (1) 式中：
$Y_{\rm{CPUE}} $ 表示单位捕捞努力量渔获量；$ {U}_{\mathrm{c}\mathrm{a}\mathrm{t}\mathrm{c}\mathrm{h}} $ 表示该作业渔区的累计渔获量 (kg)；$ {f}_{\mathrm{h}\mathrm{o}\mathrm{o}\mathrm{k}\mathrm{s}} $ 为该作业渔区的累计投钩数量 (千钩)；时间尺度为季度。1.2 数据分析方法
1.2.1 标准差椭圆
空间格局研究中常采用标准差椭圆 (Standard deviational ellipse, SDE) 方法，从而获得空间格局特征的中心性、展布范围、方向和形状等。本文采用加权标准差椭圆方法，基于南太平洋长鳍金枪鱼的空间位置，用其CPUE表示相应的权重，计算长鳍金枪鱼空间分布标准差椭圆各项指标。第一级标准差椭圆 (后文简称标准差椭圆)，表示68%的资源量分布在椭圆内。获得的空间格局特征主要指标物理意义如下：1) 标准差椭圆的中心，指示空间格局的分布中心 (渔场中心)；2) 标准差椭圆长轴，反映空间格局的分布范围；3) 椭圆扁率，短轴与长轴标准差比值，反映空间格局的分布形状；4) 方位角，表示长轴与正北向的夹角，反映空间格局的主要趋势方向。
1.2.2 空间格局的差异指数
空间格局的差异指数能够比较同一区域，采用不同计算方法获得的空间格局之间差异程度的指标。空间格局差异指数是在标准差椭圆基础上，计算不同标准差椭圆在空间上的重叠水平，反映不同空间格局的整体差异。具体计算公式如下：
$$ {\mathrm{S}\mathrm{D}\mathrm{I}}_{\mathrm{1,2}}=1-\dfrac{\mathrm{A}\mathrm{r}\mathrm{e}\mathrm{a}({{\rm{SDE}}}_{1}\cap {{\rm{SDE}}}_{2})}{\mathrm{A}\mathrm{r}\mathrm{e}\mathrm{a}({{\rm{SDE}}}_{1}\cup {{\rm{SDE}}}_{2})} $$ (2) 式中：
$ {\mathrm{S}\mathrm{D}\mathrm{E}}_{1} $ 、$ {\mathrm{S}\mathrm{D}\mathrm{E}}_{2} $ 表示两种空间分布标准差椭圆，$ \mathrm{A}\mathrm{r}\mathrm{e}\mathrm{a} $ 为面积；$ {\mathrm{S}\mathrm{D}\mathrm{I}}_{\mathrm{1,2}} $ 为空间差异指数，介于0~1，值越大表明空间格局整体差异越大[8]。1.2.3 空间格局的全局分析
全局莫兰指数 (Moran's I) 能够从总体上判定资源空间格局是否存在分布聚集性。值域为−1~1，若Moran's I为正，则指示有聚集分布趋势；反之，则指示离散分布趋势[9]。在ArcGIS 10.2软件环境运行下，Moran's I工具还会得到Z得分——标准差的倍数。采用随机条件下的近似正态分布的统计量Z进行检验，当I>0，
$ \left|Z\right| $ >1.65时，则表明空间分布为正相关性，即长鳍金枪鱼呈聚集分布，Z得分越高则相关性越明显。具体的I值和Z得分计算可参考ArcGIS在线文档[10]。1.2.4 空间格局的局部分析
由于Moran's I无法准确识别资源空间格局的局部分布特征[11]。在ArcGIS 10.2中采用热点分析工具实现，能够获得具有统计显著性的资源分布热点和冷点。热点分析对数据集中的每一个要素计算Getis-Ord Gi* 统计。通过得到的Z得分和P值，可以获得高值或低值要素在空间上发生聚类的位置，即热点和冷点，具体计算公式和原理可参考ArcGIS在线文档[12]。
2. 结果
2.1 南太平洋地区长鳍金枪鱼不同空间尺度空间格局特征
2.1.1 空间格局特征
2015—2019年中国大陆和WCPFC的长鳍金枪鱼空间格局分布形状呈现东-西分布，其资源标准差椭圆分布如图2所示。其中Q1—Q4分别代表第一到第四季度。
![]()
图 2 2015—2019年不同空间尺度的长鳍金枪鱼资源第一级标准差椭圆Figure 2. First order standard deviational ellipse spatial distribution of T. alalunga at different scales from 2015 to 20192015—2019年CN-1、CN-2和CN-5渔业资源的空间分布范围显示，三者的分布椭圆在第一季度变化趋势基本一致，资源中心的移动趋势均大致相同，总体上CN-1、CN-2和CN-5渔业资源的空间分布范围覆盖面积有很大的重合 (图2)。
2015—2019年CN-5和WCPFC-5渔业资源的空间分布范围显示 (图2)：1) CN-5较WCPFC-5的资源中心偏东；2) CN-5的移动范围较WCPFC-5的大且分布范围较广；3) CN-5的分布椭圆较为偏东；4) CN-5的渔业资源的空间分布范围呈东西方向移动，WCPFC-5的分布椭圆呈南北方向移动。
由图3可见，CN-5和WCPFC-5空间格局差异指数在0.22~0.51之间波动，且多数大于0.3。空间格局差异指数小于0.3则认为两者差异性较小，所以CN-1和CN-2、CN-2和CN-5的空间差异性较小，CN-5和WCPFC-5有较大的空间差异性。
![]()
图 3 2015—2019年不同尺度的长鳍金枪鱼资源空间格局差异指数Figure 3. Spatial pattern difference index of T. alalunga resources at different scales from 2015 to 20192.1.2 空间格局的资源中心变化特征
2015-Q1表示2015年第一季度。2015—2019年南太平洋地区长鳍金枪鱼资源中心移动显示，中国大陆数据3种空间尺度的资源中心移动趋势大致相同，CN-2的资源中心整体上跨度较其他尺度数据的大 (图2)，可见三者的差异较小；同一空间尺度的中国大陆数据和WCPFC-5数据的资源中心移动方向存在明显差异，整体上中国大陆数据的资源中心较为偏东 (图2)，说明两者差异较大。
2.1.3 资源空间格局的分布范围
标准差椭圆的长轴反映长鳍金枪鱼空间分布的范围，椭圆内部的区域是南太平洋长鳍金枪鱼分布的主体区域。2015—2019年，CN-1、CN-2和CN-5的长轴标准差有大致相同的变化趋势，分布范围差异性较小 (图4)，资源分布呈东西方向变动且波动范围较大；CN-5和WCPFC-5的长轴标准差变化趋势大致相同 (图4)，两者分布范围差异性较小，且后者波动较小，说明其资源分布范围相对稳定。
![]()
图 4 2015—2019年不同尺度的长鳍金枪鱼资源空间格局的分布范围Figure 4. Distribution range of spatial pattern of T. alalunga resources at different scales from 2015 to 20192.1.4 椭圆扁率
长短半轴的值差距越大 (扁率越大)，表示数据的方向性越明显，反之则方向性越不明显。2015—2019年，CN-1、CN-2和CN-5的空间分布形状 (短轴/长轴) 基本呈现相同的变化趋势 (图5)，分布形状差异不明显；CN-5和WCPFC-5分布形状变化趋势大致相同 (图5)，整体上两者差异较小。CN和WCPFC的扁率有明显的季节性规律，整体差异较小。扁平化表现为捕捞范围在其长轴 (东西) 方向上增长明显或短轴 (南北) 方向上收缩明显；反之则表现为捕捞范围在其长轴 (东西) 方向上收缩明显或短轴 (南北) 方向上增长明显。
![]()
图 5 2015—2019年不同尺度的长鳍金枪鱼资源标准差椭圆扁率分布Figure 5. Oblateness of standard deviation ellipticity of T. alalunga resources at different scales from 2015 to 20192.1.5 方位角变化
标准差椭圆方位角为90° 时，表示空间格局分布趋势为水平东西方向，大于90° 可以看作东侧资源向南移动幅度大于西侧，小于90° 则相反。2015—2019年，CN-1、CN-2和CN-5的方位角变化趋势基本一致 (图6)，总体来看不同空间尺度资源分布的方位角差异性较小；CN-5和WCPFC-5的方位角变化趋势大致相同，总体来看两者的方位角变化趋势差异性较小。变化趋势具体变现为，当方位角增加，空间分布椭圆呈顺时针旋转，表明长鳍金枪鱼的捕捞向椭圆西北或东南部的增长速度相对较快；当方位角减小，空间分布椭圆呈逆时针旋转，表明长鳍金枪鱼的捕捞向椭圆西南部或东北部的增长速度相对较快。
![]()
图 6 2015—2019年不同尺度的长鳍金枪鱼资源标准差椭圆方位角分布Figure 6. Azimuth of standard deviation ellipticity of T. alalunga resources at different scales from 2015 to 20192.2 资源空间格局全局特征
本文把CPUE作为观测变量，根据全局自相关方法获得CN-1、CN-2、CN-5以及WCPFC-5下的南太平洋长鳍金枪鱼渔业资源空间分布特征 (表2)。整体上看，中国大陆三种空间尺度数据和WCPFC数据下的全局空间自相关莫兰指数I值均为正值，并且
$ \left|Z\right| $ >1.96，表明南太平洋长鳍金枪鱼渔业资源的空间分布不是随机的，而是呈现出较强的空间聚集性，第三、第四季度Z得分高于第一、第二季度，说明第三、第四季度渔业资源的聚集程度相对较高。CN-1、CN-2和CN-5得出的Z得分随空间尺度变大而提高，CN-2和CN-5的Z得分差距不明显，说明两者存在微小差异；CN-5和WCPFC-5得出的莫兰指数和Z得分有很大差距，尤其是Z得分，说明两者差异较大。表 2 南太平洋长鳍金枪鱼资源空间格局的全局分布聚集性显著性检验Table 2. Significance test of global distribution and aggregation of spatial pattern of T. alalunga in South Pacific统计时间
Statistical time中国大陆1°
CN-1中国大陆2°
CN-2中国大陆5°
CN-5中西太平洋渔业委员会5°
WCPFC-5I Z_score I Z_score I Z_score I Z_score 第一季度
Quarter 10.032 3 3.975 1 0.236 6 15.015 6 0.337 4 18.301 4 0.628 7 43.204 7 第二季度
Quarter 20.032 1 2.184 4 0.047 3 10.867 5 0.371 1 17.594 5 0.563 7 47.969 7 第三季度
Quarter 30.042 2 6.137 2 0.064 5 20.338 6 0.487 0 30.940 6 0.614 7 50.730 2 第四季度
Quarter 40.038 1 6.675 1 0.046 4 20.380 1 0.532 6 29.146 5 0.637 4 49.767 8 注：Z得分的临界值为1.65、1.96和2.58分别对应的置信度为90%、95%和99%。Note: The critical values of Z score are 1.65, 1.96 and 2.58, respectively, with confidence levels of 90%, 95% and 99%. 2.3 资源空间格局的局部特征
对局部空间格局的聚集特征使用热点分析方法，对CN-1、CN-2、CN-5及WCPFC-5的南太平洋长鳍金枪鱼渔业资源进行计算。结果表明，CN-1第三季度有2个渔业资源热点区域，分别为130°W—140°W、32°S—34°S，158°W—162W°、24°S—27°S (图7)，CN-2第三季度有3个渔业资源热点区域，大概位置分别为131°W—143°W、29°S—35°S，151°W—165°W、23°S—27°S和159°E—165°E、25°S—29°S (图7)，CN-5数据下第三季度的渔业资源热点区域大概为132.5°W—157.5°E、27.5°S—32.5°S的狭长地带 (图7)，第三季度WCPFC-5的渔业资源热点区域大致为132.5°W—172.5°W、27.5°S—42.5°S以及157.5°E—167.5°E、22.5°S—32.5°S (图7)。整体来看，CN-1、CN-2和CN-5有相似的渔业资源热点区域，尺度越大，热点区域分布越广，CN-5和WCPFC-5下的渔业资源热点区域亦有大面积重合，其中WCPFC-5的热点区域范围更广。
![]()
图 7 南太平洋长鳍金枪鱼渔业资源第三季度不同尺度的空间格局的局部聚集特征Figure 7. Aggregation characteristics of spatial pattern of T. alalunga fishery resources at different scales in South Pacific in third quarter3. 讨论
3.1 长鳍金枪鱼资源的空间格局特征分析
本研究发现，长鳍金枪鱼资源空间格局呈现明显的东西分布特征，且资源中心随着季节和年际变化十分显著，季节变化空间格局趋势整体上均呈东西方向移动范围，并大于南北方向。通过该海域的环境分布研究发现海表温度、叶绿素浓度、洋流和海表风场均表现出显著的东西向分布特征[13-16]。这些研究均指出了空间格局的纬向分布特征。从格局分布形态看 (扁率) 其季节性变化特征明显，在第二、第三季度扁率低，说明短轴与长轴差距小，方向性不明显；因为第二、第三季度南太平洋长鳍金枪鱼属于盛渔期 (5—8月)，渔业资源丰富，纬向分布范围广。张嘉容等[17]研究发现的盛渔期与本文研究结果一致。渔业资源的方位角变化有明显的季节规律性，整体上在第三季度达到最大值，原因是渔业资源在东南和西北海域聚集性明显。
3.2 不同空间尺度空间格局特征分析
本研究发现，近5年CN-1和CN-2、CN-2和CN-5空间尺度数据的空间格局差异指数较小，三者的资源分布移动方向基本一致，分布中心、方位角、长轴和短轴比值变化总体上也呈大致相同的趋势；大尺度数据空间分布椭圆范围略大，两者总体差异性较小。造成这种差异性的主要原因是数据综合导致渔区内的CPUE不同从而影响权重，而分辨率的不同并未影响其移动规律，两者均能反映长鳍金枪鱼的空间分布规律。魏广恩和陈新军[18]的研究结果显示，随着空间尺度的变大对应小渔区的CPUE也会改变，作业渔场覆盖海域明显变大，与本研究结果一致。分辨率会影响变化规律的精准性，小尺度数据更能准确真实地反映长鳍金枪鱼渔业资源的变化规律，更具参考价值，而CN-2、CN-5则不能准确反映长鳍金枪鱼渔业资源的变化规律，误差较大。
渔场变动会影响标准差椭圆的长轴、方位角和分布形状等，影响渔场变动的因素较多，其中最主要的是温度、盐度和海表面风场等海洋环境因子[19-21]。整体来看，CN-1、CN-2和CN-5空间尺度的空间格局差异性相对较小，由于是同源的中国大陆数据，所以造成这种差异性的主要原因在于空间分辨率的不同而导致渔区内CPUE的不同。
3.3 相同空间尺度的空间格局特征分析
CN-5和WCPFC-5的空间差异指数相对较大。从资源分布范围来看，其变化趋势有一定差异，WCPFC的数据综合了各个成员国几乎所有的渔业数据，数据量大，而中国大陆的渔业数据量少，且WCPFC-5的CPUE大于CN-5的，CPUE作为标准差椭圆的权重会造成这种差异性。已有研究表明，CPUE与资源量之间的正比关系主要受捕捞能力 (渔船吨位、马力、渔具和助渔设备等)、时空要素 (年、月、经度、纬度) 和环境条件[22-24] 。在时间分辨率、空间分辨率和海洋环境相同条件下，影响CPUE的主要原因是捕捞能力。目前，中国远洋渔业仍处于大而不强的发展阶段，远洋渔业存在设施陈旧，渔业捕捞勘探等技术缺乏创新，经营管理水平较低等问题，因此中国大陆的捕捞能力低于WCPFC其他成员国，导致CPUE较低[25-27]。从空间分布中心来看，两者渔场中心差异性较大，移动趋势整体上不同，WCPFC-5的渔场中心偏西，结合其CPUE，可能是西部海域更适合长鳍金枪鱼生存，渔业资源丰富、渔获量大，另外为了保护渔业资源、限制资源的过度开发，区域渔业管理组织对各国渔船总数进行了限制，根据WCPFC 2019年统计年鉴，2015—2019年我国在中西太平洋区域延绳钓渔船数量从429艘下降到370艘。两者长轴变化、方位角移动和空间分布形状变化趋势大致相同。长鳍金枪鱼是一种高度洄游的中上层鱼类，研究表明环境变化对其分布和渔场变动有显著影响[28-29]，渔场的分布范围、渔场中心的移动及其分布形状等的变化与海洋环境有着密切联系，在以后的研究中应考虑环境因子的影响。综合来看，中国大陆大尺度数据和WCPFC-5的差异性主要由数据源和数据量的不同而导致渔区内CPUE的不同。
3.4 中国大陆数据和WCPFC空间格局的全局和局部聚集特征分析
CN-5和WCPFC-5的莫兰指数I值和Z得分有很大差距，产生这种现象的原因可能是中国大陆渔船数量较少 (表1)，分布较为分散，且CPUE差距较大，但四者均呈现聚集性分布。原作辉等[30]对中西太平洋长鳍金枪鱼延绳钓渔场分布的研究显示，各月捕捞努力量存在空间自相关并呈现显著聚集分布模式，与本研究结果一致。刘禹希等[31]用全局自相关研究了南海北部海域大眼鲷 (Priacanthus tayenus)，结果表明年际间资源密度普遍为正相关。
第三季度CN-1有2个大的热点区域，CN-2有3个大的热点区域，CN-5有1个经度方向的热点地带，三者的覆盖范围有很大重叠，但是仍有一定的差异性，造成这种差异性的原因主要是分辨率的不同，整体上看这种差异性较小，受分辨率的影响较少，在可接受范围内。WCPFC-5有两个成片的热点区域，相对于CN-5的热点地带覆盖范围更广，热点区域更多，但WCPFC-5的热点区域可覆盖CN-5的热点区域，原因是两种数据的覆盖范围和CPUE不同，WCPFC-5的覆盖范围更广、CPUE值更高，源数据高值较为连续，足以形成连续的热点地带，两者的差异性总体上在可接受范围内。形成热点的原因是该区域的渔业资源密度被高值包围，长鳍金枪鱼渔业资源丰富度高，空间自相关呈正相关并且自相关性较强。整体上看，四者热点区域差异性不大，受分辨率和数据的影响较小，未在文中展现的其他季度的热点区域也有相同结论。江承旭[32]分析斐济专属经济区2013—2015年长鳍金枪鱼的渔场，侯娟等[5]研究南太平洋长鳍金枪鱼渔场，两者得出的高CPUE渔区与本文结论有一定的差异，原因是其数据来源于中国某一渔业公司，作业区域覆盖范围较小。
3.5 展望
本研究选取2015—2019年的数据，数据量有限，得出的结果可能会存在局限性，今后研究应选取更长时间尺度的数据。本研究中，中国大陆数据按照空间分辨率为1°×1°、2°×2°、5°×5° 聚合，中西太平洋渔业委员会的数据空间分辨率为5°×5°，这种大分辨率在实际中对应着很大的空间范围，会与实际热点区域有差距，如果加以应用，可能会有误差。今后研究应选用如分辨率为0.5°×0.5°、0.25°×0.25° 等较小尺度的数据，提高数据精度。长鳍金枪鱼渔业资源的分布和洄游受环境影响，许多已有研究认为温度和厄尔尼诺等环境因子会影响长鳍金枪鱼的分布[28,33-34]，杨晓明等[35]在中西太平洋围网鲣鱼的热点分析中加入了环境因子，对鲣 (Katsuwonus pelamis) 源的热点与海洋环境的关系进行了分析，在今后的长鳍金枪鱼渔业资源的热点分析中也应加入适当的环境因子，结合其他模型开展进一步研究。
-
![]()
图 1 北部湾防城港—钦州近岸海域鱼类资源调查站位图
Figure 1. Sampling station of fishery resources in Fangchenggang (Qinzhou coastal area of Beibu Gulf)
![]()
图 5 北部湾防城港—钦州近岸各季节鱼类群落多样性指数
Figure 5. Biodiversity index of fish community in Fangchenggang (Qinzhou coastal area of Beibu Gulf) in four seasons
![]()
图 6 2016年北部湾防城港—钦州近岸鱼类聚类图
Figure 6. Cluster analysis of fish species collected from Fangchenggang (Qinzhou coastal area of Beibu Gulf) in 2016
![]()
图 7 鱼类群落结构NMDS分析
Figure 7. Non-metric multidimensional scalinganalysis of fish community
![]()
图 8 北部湾防城港—钦州近岸四季环境因子
注：不同的大写字母表示差异极显著 (P<0.01)；不同的小写字母表示差异显著 (P<0.05)。
Figure 8. Environmental factors in Fangchenggang (Qinzhou coastal area of Beibu Gulf) in four seasons
Note: Different uppercase and lowercase letters on the columns indicate extremely significant difference (P<0.01) and significant difference (P<0.05).
![]()
图 9 各季节环境因子相关性热图
Figure 9. Correlation heat maps of environmental factors in four seasons
![]()
图 10 基于冗余分析的鱼类群落和环境因子排序图
Figure 10. Ordination bi-plot of fish communities and environmental factors based on redundancy analysis
表 1 北部湾防城港—钦州近岸渔获鱼类种类组成
Table 1 Species composition of fish in Fangchenggang (Qinzhou coastal area of Beibu Gulf)
纲 Class 目 Order 科 Family 属 Genus 种 Species 软骨鱼纲 Chondrichthyes 真鲨目 Carcharhiniformes 2 (3.64%) 2 (2.08%) 2 (1.32%) 鲼形目 Myliobatiformes 1 (1.82%) 1 (1.04%) 1 (0.66%) 硬骨鱼纲 Osteichthyes 银汉鱼目 Atheriniformes 1 (1.82%) 1 (1.04%) 1 (0.66%) 鳗鲡目 Anguilliformes 2 (3.64%) 4 (4.17%) 6 (3.95%) 鳕形目 Gadiformes 1 (1.82%) 1 (1.04%) 1 (0.66%) 鲱形目 Clupeiformes 3 (5.45%) 6 (6.25%) 11 (7.24%) 鲇形目 Siluriformes 2 (3.64%) 3 (3.13%) 3 (1.97%) 鲻形目 Mugiliformes 2 (3.64%) 3 (3.13%) 3 (1.97%) 灯笼鱼目 Myctophiformes 2 (3.64%) 3 (3.13%) 4 (2.63%) 颌针鱼目 Beloniformes 1 (1.82%) 1 (1.04%) 1 (0.66%) 刺鱼目 Gasterosteiformes 1 (1.82%) 1 (1.04%) 1 (0.66%) 鲈形目 Perciformes 26 (47.27%) 52 (54.17%) 91 (59.87%) 鲉形目 Scorpaeniformes 4 (7.27%) 5 (5.21%) 5 (3.29%) 鲽形目 Pleuronectiforme 5 (9.09%) 8 (8.33%) 14 (9.21%) 鲀形目 Tetraodontiforme 2 (3.64%) 5 (5.21%) 8 (5.26%) 
下载: 导出CSV
表 2 各季节相对重要性指数前5的鱼类优势种
Table 2 Importance relative index of five dominated fish species in four seasons
季节
Season种类
Species丰度
RW/%生物量占比
RN/%出现频率
RFO/%相对重要性指数
IRI春季
Spring矛尾虾虎鱼 C. stigmatias 17.76 47.27 96.55 6 279.21 二长棘鲷 P. edita 19.26 10.12 93.10 2 735.73 褐菖鲉 S. marmoratus 9.24 11.84 79.31 1 672.29 竹䇲鱼 T. japonicus 9.36 6.10 51.72 799.41 短吻鲾 L. brevirostris 6.62 2.75 65.52 613.76 夏季
Summer鹿斑鲾 S. ruconius 10.06 31.22 96.55 3 985.99 长鲾 L. elongatus 13.28 24.62 82.76 3 137.32 蓝圆鲹 D. maruadsi 11.69 8.42 100.00 2 011.67 二长棘鲷 P. edita 14.94 3.50 93.10 1 716.24 康氏小公鱼 S. commersonii 3.41 11.42 75.86 1 124.93 秋季
Autumn多齿蛇鲻 S. tumbil 9.71 41.07 86.2 4 377.55 鹿斑鲾 S. ruconius 7.02 19.99 93.10 2 515.05 棕腹刺鲀 G. spadiceus 16.25 4.72 96.55 2 024.46 二长棘鲷 P. edita 10.49 2.52 79.31 1 031.89 长丝虾虎鱼 M. cryptocentrus 3.87 5.76 86.21 829.84 冬季
Winter褐菖鲉 S. marmoratus 12.20 31.01 58.62 2 533.23 二长棘鲷 P. edita 9.87 38.85 34.48 1 680.16 鹿斑鲾 S. ruconius 4.31 4.11 82.76 696.59 多齿蛇鲻 S. tumbil 11.46 5.82 34.48 595.70 花斑蛇鲻 S. undosquamis 8.79 4.60 31.03 415.12
下载: 导出CSV
-
[1] ROZAS L P, ODUM W E. The role of submerged aquatic vegetation in influencing the abundance of nekton on contiguous tidal fresh-water marshes[J]. J Exp Mar Biol Ecol, 1988, 114(2): 289-300.
[2] ZIJLSTRA J J. On the importance of the Wadden Sea as a nusery area in relation to the conservation of the southern North Sea fishery resources[J]. Symp Zool Soc Lond, 1972, 29(2): 233-258.
[3] PIHL L, WENNHAGE H. Structure and diversity of fish assemblages on rocky and bottom shores on the Swedish west coast[J]. J Fish Biol, 2002, 61(2): 148-166.
[4] BLABER S J M, BREWER D T, SALINI J P. Fish communities and the nursery role of the shallow inshore waters of a tropical bay in the Gulf of Carpentaria, Australia[J]. Estuar Coast Shelf Sci, 1995, 40(2): 177-193. doi: 10.1016/S0272-7714(05)80004-6
[5] SELLESLAGH J, AMARA R. Environmental factors structuring fish composition and assemblages in a small macrotidal estuary (eastern English Channel)[J]. Estuar Coast Shelf Sci, 2008, 79(3): 507-517. doi: 10.1016/j.ecss.2008.05.006
[6] BARLETTA M, BARLETTA-BERGAN A, SAINT-PAUL U, et al. The role of salinity in structuring the fish assemblages in a tropical estuary[J]. J Fish Biol, 2005, 66(1): 45-72. doi: 10.1111/j.0022-1112.2005.00582.x
[7] FRANCO A C S, SANTOS L N. Habitat-dependent responses of tropical fish assemblages to environmental variables in a marine-estuarine transitional system[J]. Estuar Coast Shelf Sci, 2018, 211(2): 110-117.
[8] DIRECTIVE W F. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 establishing a framework for community action in the field of water policy[J]. Off J Eur Comm, 2000, 22(12): 20-35.
[9] CABRAL H N. Distribution and abundance patterns of flatfishes in the Sado estuary, Portugal[J]. Estuaries, 2000, 23(1): 351-358.
[10] THIEL R, SEPÚLVEDA A, KAFEMANN R, et al. Environmental factors as forces structuring the fish community of the Elbe estury[J]. J Fish Biol, 1995, 46(1): 47-69. doi: 10.1111/j.1095-8649.1995.tb05946.x
[11] MARSHALL S, ELLIOTT M. Environmental influences on the fish assemblage of the Humber estuary, U.K.[J]. Estuar Coast Shelf Sci, 1998, 46(2): 175-184. doi: 10.1006/ecss.1997.0268
[12] 王雪辉, 邱永松, 杜飞雁, 等. 北部湾鱼类群落格局及其与环境因子的关系[J]. 水产学报, 2010, 34(10): 1579-1586. [13] 王理想. 北部湾海域春、秋季鱼类群落结构初步研究[D]. 青岛: 中国海洋大学, 2009: 35-37. [14] 罗春业, 李英, 朱瑜, 等. 广西北部湾鱼类区系的再研究[J]. 广西师范大学学报(自然科学版), 1999, 22(2): 85-89. [15] 杨璐, 曹文清, 林元烧, 等. 夏季北部湾九种经济鱼类的食性类型及营养生态位初步研究[J]. 热带海洋学报, 2016, 35(2): 66-75. doi: 10.11978/2014133 [16] 李渊, 王燕平, 张静, 等. 北部湾口海域鱼类分类多样性的初步探讨[J]. 应用海洋学学报, 2016, 35(2): 229-235. doi: 10.3969/J.ISSN.2095-4972.2016.02.012 [17] 王淼娣, 王雪辉, 孙典荣, 等. 基于长度贝叶斯生物量估算法评估北部湾大头白姑鱼资源状况[J]. 南方水产科学, 2021, 17(2): 20-27. doi: 10.12131/20200202 [18] 孙典荣, 林昭进. 北部湾主要经济鱼类资源变动分析及保护对策探讨[J]. 热带海洋学报, 2004, 27(2): 62-68. doi: 10.3969/j.issn.1009-5470.2004.02.008 [19] 傅昕龙. 北部湾西北部近海鱼类资源分布和群落结构特征及其影响因素. [D]. 上海: 上海海洋大学, 2018: 29-33. [20] 何思璇, 何斌源. 防城河口湾鱼类群落结构及其与环境因子关系研究[J]. 热带海洋学报, 2019, 38(5): 86-97. [21] 袁华荣, 陈丕茂, 贾晓平, 等. 北部湾东北部游泳生物资源现状[J]. 南方水产科学, 2011, 7(3): 31-38. doi: 10.3969/j.issn.2095-0780.2011.03.006 [22] 乔延龙, 林昭进. 北部湾地形、底质特征与渔场分布的关系[J]. 海洋湖沼通报, 2007(1): 232-238. [23] 陈波, 侍茂崇. 北部湾海洋环流研究进展[J]. 广西科学, 2019, 26(6): 595-603. [24] 杨洁清, 史赟荣, 沈新强, 等. 湄洲湾海域甲壳动物群落多样性分析[J]. 南方水产科学, 2017, 13(5): 16-24. doi: 10.3969/j.issn.2095-0780.2017.05.003 [25] 劳齐斌, 刘国强, 申友利, 等. 北部湾入海河流营养盐的分布特征及入海通量研究[J]. 海洋学报, 2020, 42(12): 93-100. [26] PINKAS L, OLIPHANT M S, IVERSON I L K. Food habits of albacore, bluefin tuna, and bonito in California waters[J]. Fish Bull, 1971, 152(7): l-105.
[27] CLARKE K R, GORLEY R N, SOMERFIELD P J, et al. Change in marine communities: an approach to statistical analysis and interpretation[M]. Plymouth: Primer-E Ltd. , 1994: 61-77.
[28] BRAAK C J F. CANOCO: a FORTRAN program for canonical community ordination by partial etrended canonical correspondence analysis, principal components analysis and redundancy analysis (version 2.1) [M]. Wageningen: Agricultural Mathematics Group, 1988: 45-50.
[29] KNEIB R T. Early life stages of resident nekton in intertidal marshes[J]. Estuaries, 1997, 20(4): 214-230.
[30] ELLIOTT M, DEWAILLY F. The structure and components of European estuarine fish assemblages[J]. Aquat Ecol, 1995, 29(6): 397-417.
[31] 孙典荣. 北部湾渔业资源与渔业可持续发展研究[D]. 青岛: 中国海洋大学, 2008: 40-44. [32] 陈作志, 邱永松. 北部湾二长棘鲷的生态分布[J]. 海洋水产研究, 2005, 26(3): 16-21. [33] PESSANHA A L M, ARAÚJO F G. Spatial, temporal and diel variations of fish assemblages at two sandy beaches in the Sepetiba Bay[J]. Estuar Coast Shelf Sci, 2003, 57(5/6): 817-828.
[34] SAYER M D J, CAMERON K S, WILKINGSON G. Fish species found in the rocky sublittoral during winter months as revealed by the underwater application of the anaesthetic quinaldine[J]. J Fish Biol, 1994, 44(4): 351-353.
[35] SAYER M D J, GIBSON R N, ATKINSON R J A. Growth, diet and condition of goldsinny on the west coast of Scotland[J]. J Fish Biol, 1995, 43(2): 157-167.
[36] VENDEL A L, LOPES S G, SANTOS C, et al. Fish assemblages in a tidal flat[J]. Braz Arch Biol Technol, 2003, 46(5): 233-242.
[37] 王言丰, 余景, 陈丕茂, 等. 北部湾灯光罩网渔场时空分布与海洋环境关系分析[J]. 热带海洋学报, 2019, 38(5): 68-76. [38] 王雪辉, 邱永松, 杜飞雁, 等. 北部湾鱼类多样性及优势种的时空变化[J]. 中国水产科学, 2011, 18(2): 427-436. [39] 王雪辉, 邱永松, 杜飞雁, 等. 北部湾秋季底层鱼类多样性和优势种数量的变动趋势[J]. 生态学报, 2012, 32(2): 333-342. [40] 袁蔚文. 北部湾底层渔业资源的数量变动和种类更替[J]. 中国水产科学, 1995, 2(2): 57-65. doi: 10.3321/j.issn:1005-8737.1995.02.006 [41] PRISTA N, VASCONCELOS R P, COSTA M J, et al. The demersal fish assemblage of the coastal area adjacent to the Tagus estuary (Portugal): relationships with environmental conditions[J]. Oceanologica Acta, 2003, 26(5/6): 525-536.
[42] MOLINA A, DUQUE G, COGUA P. Influences of environmental conditions in the fish assemblage structure of a tropical estuary[J]. Mar Biodivers, 2020, 50(1): 1-13. doi: 10.1007/s12526-019-01025-y
[43] 王小林. 海州湾及邻近海域鱼类群落结构的时空变化[D]. 青岛: 中国海洋大学, 2013: 42-47. [44] WHITFIELD A K. Ichthyofaunal assemblages in estuaries: a South African case study[J]. Rev Fish Biol Fish, 1999, 9(2): 151-186. doi: 10.1023/A:1008994405375
[45] FONDS M. The seasonal distribution of some fish species in the Western Dutch. Wadden Sea[J]. J Fish Biol, 1978, 32(4): 223-235.
[46] MEYNECKE J O, LEE S Y, DUKE N C, et al. Effect of rainfall as a component of climate change on estuarine fish production in Queensland, Australia[J]. Estuar Coast Shelf Sci, 2006, 69(3/4): 491-504.
[47] LONERAGAN N R. River flows and estuarine ecosystems: implications for coastal fisheries from a review and a case study of the Logan River, southeast Queensland[J]. Aust J Ecol, 1999, 24(4): 431-440. doi: 10.1046/j.1442-9993.1999.00975.x
[48] GILLSON J, SCANDOL J, SUTHERS I. Estuarine gillnet fishery catch rates decline during drought in eastern Australia[J]. Fish Res, 2009, 99(1): 26-37. doi: 10.1016/j.fishres.2009.04.007
[49] POXTON M G, ALLOUSE S B. Water quality criteria for marine fisheries[J]. Aquacult Eng, 1982, 1(3): 153-191. doi: 10.1016/0144-8609(82)90026-7
[50] POMFRET J R, ELLIOTT M, REILLY M G. Spatial and temporal patterns in the fish communities in two UK North Sea estuaries[J]. Estuar Coast, 1991, 11(3): 227-284.
[51] MAES J, STEVENS M, BREINE J. Modelling the migration opportunities of diadromous fish species along a gradient of dissolved oxygen concentration in a European[J]. Estuar Coast Shelf Sci, 2007, 75(5): 151-162.
[52] 胡翠林, 张洪亮, 张亚洲, 等. 带鱼保护区春秋季鱼类群落特征及与环境因子的关系[J]. 水产学报, 2018, 42(5): 694-703. [53] BARROS F, UNDERWOOD A J, LINDEGARTH M. The influence of rocky reefs on structure of benthic macrofauna in the nearby soft-sediments[J]. Estuar Coast Shelf Sci, 2001, 52(2): 191-199. doi: 10.1006/ecss.2000.0734
[54] 邱永松. 南海北部大陆架鱼类群落的区域性变化[J]. 水产学报, 1988, 10(4): 303-313. [55] GELWICK F P, AKIN S, ARRINGTON A D. Fish assemblage structure in relation to environmental variation in a Texas gulf coastal Wetland[J]. Estuaries, 2001, 24(3): 285-296.
[56] RENFRO W C. Survival and migration of fresh-water fishes in salt water[J]. Tex J Sci, 1959, 11(5): 172-180.
-
期刊类型引用(0)
其他类型引用(1)
-
其他相关附件
-
PDF格式
20210255-附录A 点击下载(174KB)
-
计量
- 文章访问数: 974
- HTML全文浏览量: 296
- PDF下载量: 84
- 被引次数: 1
粤公网安备 44010502001741号
