A preliminary study on community structure of mesopelagic fish in cold seep of Xisha Islands
-
摘要: 为了解冷泉生态系统鱼类种类组成，以及为相关区域的生物多样性保护和监测工作提供科学依据。采用中层拖网对西沙群岛冷泉海域内外区域中层鱼类资源进行调查，共采集鱼类12目31科62属106种，以灯笼鱼目和巨口鱼目为最主要类群。其中，冷泉海域11目30科54属84种；冷泉外海域6目13科32属60种。冷泉海域内外区域共有鱼类38种，共有优势种为斯氏蝰鱼 (Chauliodus sloani)、尾明角灯鱼 (Ceratoscopelus warmingii)、喀氏眶灯鱼 (Diaphus garmani)。Jaccard种类相似性指数 (I) 表明，冷泉海域内外区域中层鱼类种类相似度为中等不相似 (I=35.8%)。鱼类群落的平均分类差异指数 (Δ+) 冷泉海域显著大于冷泉外海域 (P<0.01)，但鱼类群落的分类差异变异指数 (Λ+) 冷泉海域小于冷泉外海域 (P<0.01)。从水层分布来看，冷泉海域内外中层鱼类群落物种组成存在明显差异 (P<0.01)；表层水域中层鱼类多样性冷泉区低于冷泉外部，而深层水域则相反。冷泉内外不同站位的群落组成之间在水层及时间上存在相关性，冷泉内部的白天深层和夜间深层相似性较高；冷泉内部夜间表层和冷泉外部夜间表层相似性较高。研究结果初步显示，冷泉内外中层鱼类群落组成差异显著，冷泉区中层鱼类种类更多，且深层水域的种类多样性较高；区分冷泉内外海域中层鱼类群落的关键种群为长鳍虹灯鱼 (Bolinichthys longipes) 和串灯鱼 (Vinciguerria nimbara) 等8种，累计贡献度达到56.06%。Abstract: To recognize the fish composition of cold seep ecosystems, and provide a scientific basis for biodiversity conservation and monitoring in the relevant regions, we investigated the mesopelagic fish resources inside and outside the cold seep in the Xisha Islands by a mesopelagic trawler. A total of 106 species of mesopelagic fish had been identified which belong to 62 Genera, 31 Families and 12 Orders, with the dominant groups of Myctophiformes and Stomiiformes. Among them, 84 species of fish belonging to 54 Genera, 30 Families, 11 Orders were collected inside the cold seep; while 60 species of fish belonging to 32 Genera, 13 Families and 6 Orders were collected outside the cold seep. There were 38 shared species of mesopelagic fish belonging to 22 Genera, 9 Families, 5 Orders inside and outside the cold seep, and the main shared species were Chauliodus sloani, Ceratoscopelus warmingii and Diaphus garmani. The Jaccard species similarity index shows that the species similarity of mesopelagic fish was moderately different inside and outside the cold seep (I=35.8%). The average taxonomic distinctness (Δ+) of mesopelagic fish was significantly higher inside the cold seep than outside (P<0.01), but the variation in the taxonomic distinctness (Λ+) of mesopelagic fish community was inverse (P<0.01). From the perspective of pelagic distribution, there were significant differences in the species composition of the mesopelagic fish communities inside and outside the cold seep (P<0.01). In the surface waters, the mesopelagic fish diversity inside the cold seep was lower than outside, while it's inverse in the deeper water layers. There was a correlation between the fish communities at different stations inside and outside the cold seep in terms of water stratigraphy and timing, with higher similarity between the deeper layers inside the cold seep; and higher similarity between the nocturnal surface layers inside the cold seep area and the nocturnal surface layers outside the cold seep. The results show that there are differences in the mesopelagic fish communities inside and outside the cold seep, with more species of mesopelagic fish inside the cold seep and a higher diversity of species in the deeper water layers. Eight key species that distinguish the mesopelagic fish communities in the waters inside and outside the cold seep include Bolinichthys longipes and Vinciguerria nimbara, having a cumulative contribution of 56.06%.
-
Keywords:
- Mesopelagic fish /
- Cold seep /
- Species composition /
- Dominant species /
- Xisha Islands
-
近几十年来，由于捕捞强度和环境污染等加剧，世界范围内海洋渔业资源的衰退现象日益严重[1]。据世界粮农组织统计，世界上1/4的鱼类种群处于过度捕捞、资源量剧减或恢复状态[2]。尽管物种灭绝的情况较少，但是一些鱼种出现了“经济价值灭绝(commercial extinction)”，即由经济鱼种衰竭为非经济鱼种[3]。实际上，对于大多数经济价值较高的海洋鱼类种群而言，捕捞可能是引起鱼类表型特征演化的潜在力量[4]；同时相对于海洋环境污染、水质恶化和人类活动引起的气候变化，捕捞对渔业资源造成的影响更大[5]。生物生长是基因与环境相互作用的结果，而鱼类的遗传力较低，其生长又受环境的强烈影响，研究中很难区分表型可塑性参数的变化和基因的影响，因而不同因子对鱼类种群的长远影响尚难以估量[6]。而生长、死亡和性成熟等生物学特征的演变将对种群的生物量、经济产量以及种群统计学产生重要作用，继而影响着渔业资源评估结果与管理目标的制定[7-8]。
北部湾是南海渔业资源生产力最高的海域之一，也是中国与邻国就渔业问题签订合作协定的区域之一。近几十年来，其海洋渔业资源发生了深刻的变化，重要经济种类的资源密度、捕捞渔船单产和渔获质量不断下降[9]。20世纪60年代的传统经济种类如笛鲷、石斑鱼等逐渐被发光鲷(Acropoma japonicum)、黄斑鲾(Leiognathus bindus)等[10]个体较小、营养级和经济价值较低的种类代替。2012年北部湾的捕捞产量约为85.7×104 t，其中中国渔船的产量为65.7×104 t[11]，依据袁蔚文[12]的资源量评估结果，当前北部湾渔业资源已承受过高的捕捞压力。北部湾经济鱼类如花斑蛇鲻(Saurida undosquamis)[13]、短尾大眼鲷(Priacanthus macracanthus)[14-15]、二长棘犁齿鲷(Evynnis cardinalis)[16-17]等出现了生活史特征的演化现象，如个体小型化、生长速度加快和性成熟提前等。这种演化现象会造成鱼类经济价值降低、渔业可持续产量减少等一系列不良结果[18]。因此，对北部湾蓝圆鲹群体生物学特征及开发状态进行评估，对评价渔业资源利用现状和制定科学的管理策略非常必要。
蓝圆鲹(Decapterus maruadsi)属近海暖水性中上层洄游鱼类[19]，广泛分布于东印度洋和中西太平洋海域，是中国东南沿海重要的经济鱼种，主要由拖网和围网捕获[20]。卢振彬等[21]对中国近海蓝圆鲹资源进行了分析，发现福建近海蓝圆鲹因捕捞压力加大，群体趋向简单化、低龄化和小型化，种群结构已遭到破坏；孙典荣和林昭进[22]发现北部湾海域蓝圆鲹资源量总体呈下降趋势，但近年来其资源密度在逐渐恢复。在高强度的捕捞压力下，北部湾蓝圆鲹表型特征是如何演变的、这种演变过程是否可逆，仍缺乏相关基础研究支撑。因此，本文根据1992—1993年(1992年)、1998—1999年(1998年)、2006年、2007年、2008—2009年(2009年)、2010年和2011—2012年(2012年)7个时期底拖网调查资料，对北部湾蓝圆鲹的群体组成、生长、死亡、性成熟和开发率等进行了分析，以期为蓝圆鲹群体的可持续发展提供科学依据和数据基础。
1. 材料与方法
1.1 实验材料
本研究所用样品主要来自中国水产科学研究院南海水产研究所在北部湾区域历年的底拖网调查资料。采样站点基本保持不变，均按照经纬度30′设置，所有年份的采样与生物学测定均依据《海洋调查规范GB 13763》进行，其中1992—1993年、1998—1999年、2008—2009年、2011—2012年的调查按季度进行，每相邻两年互补为春夏秋冬4个季度；2006年、2007年、2010年的调查为全年4个季度采样。生物学测定包括体长、体质量、性别以及性成熟度等。每站蓝圆鲹样本多于50尾时测定50尾，不足50尾时则全部测定。1992年测定635尾(缺乏性别数据)，1998年测定1 714尾，2006年测定1 049尾，2007年测定1 028尾，2009年测定1 226尾，2010年测定465尾，2012年测定1 318尾。
1.2 方法
1.2.1 群体结构分析
对蓝圆鲹7个时期的体长频次按月以10 mm为间距分别进行分组，其中占样品总数10%以上的体长组被定义为优势体长组。同一时期雌雄个体体长组成的显著性差异分析利用非参数Kolmogorov-Smirnov Z检验分析完成；通过非参数Kruskal-Wallis检验7个时期体长组成差异性是否显著；利用单因素方差分析检验7个时期体长均值的显著性差异；使用适合性卡方检验判断性成熟度数据较为充分的3个时期的雌雄比例是否为1∶1。
1.2.2 体长与体质量关系
体长体质量关系可以用以下公式表示：
$$ \quad\quad\quad\quad\quad\quad\quad\quad\quad\quad W = a{L^b} $$ (1) 其中W为体质量，L为体长。a为生长的条件因子(condition factor)，反映了资源群体所处的环境状况[23]。b为异速生长因子(allometric factor)，反映生长发育的不均匀性[24]；若b>3，则为正异速生长，体质量的增加快于体长的增长；若b=3，则为等速生长，体长和体质量增加等速；若b<3，则为负异速增长，体长的增长快于体质量的增加[25-26]。
肥满度(relative fatness)作为重要的渔业生态学指数，与体长和体质量关系密切，反映了个体的生理和营养状况[27]。在个体长度一定时，生理和营养状况越好的个体其体质量越大，肥满度值也就越高[28]。使用相对体质量Wr (%)的平均值表示肥满度[29]，计算公式为：
$$ \quad\quad\quad\quad\quad\quad\quad\quad{W_{\rm{r}}} = \frac{W}{{a{L^b}}} \times 100{\rm{\% }} $$ (2) 1.2.3 生长参数
蓝圆鲹的生长选用von Bertalanffy生长方程[30]拟合：
$$ \quad\quad\quad\quad\quad\quad{L_t} = L{}_\infty \{1 - \exp [ - K(t - {t_0})]\} $$ (3) 式中Lt为t龄时的体长，
$L_\infty$ 为渐近体长，K为生长系数，t0为理论生长起点年龄；其中$ L_\infty$ 、K可利用FAO开发的FiSATⅡ(V1.2.2)软件中的ELEFANⅠ(electronic length frequency analysisⅠ)技术根据体长频率的时间序列进行估算[31]，t0可以用以下经验公式求得:$$ \quad{\log _{10}}( - {t_0}) = - 0.392\;2 - 0.275\;1{\log _{10}}{L_\infty } - 1.038{\log _{10}}K $$ (4) 根据线性体长变换渔获曲线法可评估得到总死亡系数Z[30]。由Pauly[32]经验公式可求得自然死亡系数M：
$$ \begin{aligned} \quad\quad\quad\quad{\log _{10}}M =& - 0.006\;6 - 0.279{\log _{{\rm{1}}0}}{L_\infty } +\\ &0.654\;3{\log _{10}}K + 0.463\;4{\log _{10}}T \end{aligned} $$ (5) 其中T为该鱼种栖息水域的年平均水温(℃)。由于北部湾蓝圆鲹7个时期的年平均水温在实际测量中可能存在误差，因此采用这7个时期平均水温的均值计算，水温均值为26.20 ℃。
捕捞死亡系数(F)和开发率(E)分别通过以下公式计算：
$$ \quad\quad\quad\quad\quad\quad\quad\quad\quad F{\rm{ = }}Z - M $$ (6) $$ \quad\quad\quad\quad\quad\quad\quad\quad\quad E = F/Z $$ (7) 1.2.4 性成熟
蓝圆鲹的性成熟按照鱼类性腺成熟度Ⅵ期的标准划分，性腺Ⅲ期及以上个体认为达到性成熟。利用非线性回归拟合反正弦平方根形式的逻辑斯蒂(Logistic)模型可得到初次性成熟体长(L50)的95%置信区间[33]：
$$ \quad\quad\quad\quad\quad\quad{\rm{ASR}}({P_i}) = \frac{{{\rm{ASR\;(}}G{\rm{)}}}}{{1 + {{\rm e}^{ - \lambda \;({L_i} - {L_{50}})}}}} + {\varepsilon _i} $$ (8) 其中
${\rm{ASR}}\left( {{P_i}} \right) = \arcsin \sqrt {{P_i}} $ ；λ为瞬时成熟速率；G表示可得到的最高性成熟比例，在一定体长后所有蓝圆鲹个体均能达到性成熟，因此G取值1；Pi表示每个体长组对应的性成熟比例，本研究中，以5 mm的间距分组，其中每个体长组的组中值为Li；${\varepsilon _i}$ 为误差项。初次性成熟年龄(A50)以及95%置信区间由公式(3)转换得到：$$ \quad\quad\quad\quad\quad\quad{{{A}}_{50}} = {t_0} - \frac{1}{k}\ln \;(1 - \frac{{{L_{50}}}}{{{L_\infty }}}) $$ (9) 2. 结果
2.1 群体结构
由北部湾蓝圆鲹7个时期的体长频率分布(图1)和群体组成(表1)分析得到，1992年体长为90~270 mm，体长频率分布呈明显的单峰型，优势体长组为130~170 mm，占个体总数的64.2%；1998年体长为32~300 mm，体长频率分布呈较明显的单峰型，优势体长组为150~170 mm，占个体总数的18.9%；2006年体长为76~264 mm，体长频率分布呈较明显的单峰型，优势体长组为120~150 mm，占个体总数的47.8%，雌雄个体体长组成无显著性差异(z=0.86，P>0.01)；2007年体长为67~250 mm，体长频率分布呈较明显的单峰型，优势体长组为110~140 mm，占个体总数的46.8%，雌雄个体体长组成无显著性差异(z=1.883，P>0.01)；2009年体长为86~235 mm，体长频率分布呈较明显的单峰型，优势体长组为130~160 mm，占个体总数的61.5%，雌雄个体体长组成无显著性差异(z=1.185，P>0.01)；2010年体长为91~215 mm，体长频率分布呈较明显的单峰型，优势体长组为130~160 mm，占个体总数的65.4%，雌雄个体体长组成无显著性差异(z=0.435，P>0.01)；2012年体长为75~256 mm，体长频率分布呈较明显的双峰型，优势体长组为120~130 mm和140~170 mm，分别占个体总数的13.6%和44.1%，雌雄个体的体长组成无显著性差异(z=0.44，P>0.01)。
![]()
图 1 北部湾蓝圆鲹7个时期的体长频率分布Figure 1. Variation in length frequency of D. maruadsi in Beibu Gulf during seven periods表 1 北部湾蓝圆鲹7个时期的群体结构Table 1. Population structure of D. maruadsi in Beibu Gulf during seven periods年份
year尾数
ind.性比
sex ratio体长/mm body length 范围
distribution均值
mean优势组
dominant group1992 635 − 90~270 160.23 130~170 1998 1 714 0.48 32~300 155.29 160~170 2006 1 049 2.21 76~264 151 120~150 2007 1 028 2.6 67~250 139 110~140 2009 1 226 1.84 86~235 144 130~160 2010 465 0.656 96~215 144 130~160 2012 1 318 1.456 75~256 149 120~130；140~170 7个时期体长组的非参数Kruskal-Wallis检验表明其体长组成存在显著性差异(χ2=152.89，P<0.01)；7个时期体长的单因素方差分析表明其均值存在显著性差异(P<0.01)。
2.2 生长参数
7个时期蓝圆鲹的体长体质量关系分别为：
$$ 1992{\text{年}}, \;W = 5.27 \times {10^{ - 5}}{L^{2.786}}\;\;\;\big( {{R^2} = 0.917\;4, \;n = 635} \big) $$ (10) $$ 1998{\text{年}}, \;W = 1.43 \times {10^{ - 5}}{L^{2.993}}\;\;\big( {{R^2} = 0.957\;4,\;n = 1\;714} \big) $$ (11) $$ 2006{\text{年}},\;\;W = 9.77 \times {10^{ - 6}}{L^{3.048}}\;\big( {{R^2} = 0.948\;4,\;n = 1\;049} \big) $$ (12) $$ 2007{\text{年}},\;\;W = 3.79 \times {10^{ - 6}}{L^{3.24}}\;\big( {{R^2} = 0.962\;8,\;n = 1\;028} \big) $$ (13) $$ 2009{\text{年}},\;\;W = 6.3 \times {10^{ - 6}}{L^{3.139}}\;\big( {{R^2} = 0.964\;1,\;n = 1\;450} \big) $$ (14) $$ 2010{\text{年}},\;\;W = 9.57 \times {10^{ - 6}}{L^{3.04}}\;\big( {{R^2} = 0.864\;5,\;n = 459} \big) $$ (15) $$ 2012{\text{年}},\;\;W = 1.74 \times {10^{ - 6}}{L^{3.38}}\;\big( {{R^2} = 0.948\;5,\;n = 1\;318} \big) $$ (16) 蓝圆鲹肥满度的年际变化趋势见图2。北部湾蓝圆鲹群体2012年的肥满度(101.05)与1992年(101.86)相当，在1998年(103.20)略有升高，但在2006年(100.86)随即降低，此后呈连续波动状态。由公式(10)~(16)可得到，与1992年相比，2012年的条件因子b呈增大趋势，但b在7个时期中总体呈先增后减的趋势。
![]()
图 2 北部湾蓝圆鲹肥满度的年际变化Figure 2. Interannual variation in relative fatness of D. maruadsi in Beibu Gulf采用ELEFANⅠ方法估算出7个时期蓝圆鲹von Bertalanffy的生长方程分别为：
$$ \quad\quad\quad\quad 1992{\text{年}},\;\;{L_{\text{t}}} = 297 \times (1 - {{\rm e}^{ - 0.39 \times (t + 0.42)}}) $$ (17) $$ \quad\quad\quad\quad 1998{\text{年}},\;\;{L_{\rm{t}}} = 320 \times (1 - {{\rm e}^{ - 0.4 \times (t + 0.40)}}) $$ (18) $$\quad\quad\quad\quad 2006{\text{年}},\;\;{L_{\rm{t}}} = 304.5 \times (1 - {{\rm e}^{ - 0.33 \times (t + 0.22)}}) $$ (19) $$ \quad\quad\quad\quad 2007{\text{年}},\;\;{L_{\rm{t}}} = 329 \times (1 - {{\rm e}^{ - 0.76 \times (t + 0.21)}}) $$ (20) $$ \quad\quad\quad\quad 2009{\text{年}},\;\;{L_{\rm{t}}} = 267.5 \times (1 - {{\rm e}^{ - 0.52 \times (t + 0.32)}}) $$ (21) $$ \quad\quad\quad\quad 2010{\text{年}},\;\;{L_{\rm{t}}} = 255 \times (1 - {{\rm {e}}^{ - 0.53 \times (t + 0.23)}}) $$ (22) $$ \quad\quad\quad\quad 2012{\text{年}},\;\;{L_{\rm{t}}} = 239 \times (1 - {{\rm {e}}^{ - 0.5 \times (t + 0.35)}}) $$ (23) 蓝圆鲹极限体长(
$ L_\infty$ )和生长系数(K)的年际变化见图3。$ L_\infty$ 和K分别从1992年的297 mm和0.39上升到2007的329 mm和0.76，之后逐年下降，2012年为239 mm和0.5，其中$ L_\infty$ 的总体趋势为降低，K的总体趋势为增加。![]()
图 3 北部湾蓝圆鲹极限体长 ($L_\infty $ ) 和生长系数 (K) 的年际变化Figure 3. Interannual variation in asymptotic length ($L_\infty $ ) and growth coefficient (K) of D. maruadsi in Beibu Gulf2.3 死亡参数
运用渔获量曲线法估算7个时期蓝圆鲹总死亡系数的回归线(图4)。其中以“未全面补充年龄段个体体长接近渐进体长的年龄段不能用来回归”为回归点的选取依据[30]。7条直线中前5条直线的斜率不断增大，后2条直线斜率逐渐减小，即1992年到2009年5个世代的总死亡系数Z不断增大，从1992年的1.66持续增加到2009年的2.77，此后开始逐年减小，2010年、2012年的分别为2.17、1.68。这7个世代中，捕捞死亡系数F也呈明显的先增后减趋势，从1992年的1.17逐年增加到2009年的2.16，随后逐年降低，2010年、2012年分别为1.54、1.06 (表2)。开发率(E)、自然死亡系数(M)则呈波动状态，但整体趋势为先增后减。
表 2 北部湾蓝圆鲹7个时期死亡参数的变化Table 2. Total mortality coefficient, natural mortality coefficient, fishing mortality coefficient and exploitation of D. maruadsi in Beibu Gulf during seven periods项目
item年份 year 1992 1998 2006 2007 2009 2010 2012 总死亡系数 (Z) total mortality coefficient 1.66 1.74 1.81 2.71 2.77 2.17 1.68 自然死亡系数 (M) natural mortality coefficient 0.49 0.49 0.44 0.74 0.61 0.63 0.62 捕捞死亡系数 (F) fishing mortality coefficient 1.17 1.25 1.37 1.97 2.16 1.54 1.06 开发率 (E) exploitation rate 0.7 0.72 0.76 0.73 0.78 0.71 0.63 ![]()
图 4 根据渔获量曲线方法估算的北部湾蓝圆鲹总死亡系数Figure 4. Length converted catch curve for D. maruadsi in Beibu Gulf2.4 性成熟
由于部分时期缺乏性成熟度数据，现仅对数据较为充分的3个时期(1998年、2007年、2012年)进行性成熟分析，卡方检验表明北部湾蓝圆鲹这3个时期的雌雄比与1∶1均有显著性差异，1998年(♀∶♂=0.66，χ2=688.74，P<0.01)为雌性少于雄性，2007年(♀∶♂=2.67，χ2=11.56，P<0.01)和2012年(♀∶♂=1.46，χ2=50.02，P<0.01)为雌性多于雄性。
通过非线性回归方法拟合的ASR逻辑斯蒂曲线可以得出这3个时期的L50分别为180.16 mm (177.12~183.21为95%置信区间，下同)、174.91 mm (170.84~178.97)和180.87 mm (177.96~183.78) (图5)。通过von Bertalanffy生长方程逆算的A50分别为1.860 (1.806~1.915)、0.888 (0.854~0.923)和2.538 (2.439~2.640) (图6)。
![]()
图 5 北部湾蓝圆鲹雌性个体性成熟比例观测值与ASR模型拟合的逻辑斯蒂曲线Figure 5. Observed proportion of mature individuals and logistic curve fitted by ASR model of D. maruadsi in Beibu Gulf during three periods![]()
图 6 蓝圆鲹雌性个体的初次性成熟体长 (L50) 与初次性成熟年龄 (A50) 的演变趋势误差线为95%置信区间Figure 6. Fluctuation in L50 and A50 of female D. maruadsiError bars show the 95% confidence interval.3. 讨论
依据体长变换渔获量曲线法估算死亡参数，本研究发现北部湾蓝圆鲹的总死亡系数和捕捞死亡系数自20世纪90年代起开始逐年升高，于2009年达到最高值，分别为2.77和2.16，随后逐年降低，2012年下降至1.68和1.06。其开发率也呈现出相同的变化趋势，从1992年的0.7上升至2009年的0.78，之后降低至2012年的0.63，其中2010年北部湾的蓝圆鲹开发率与20世纪90年代早期相当。而依据Gulland[34]提出的一般鱼类最适利用率为0.5来分析，结合表2可发现，自20世纪90年代以来，北部湾蓝圆鲹群体一直处于过度捕捞状态，进入21世纪后虽仍处于过度开发状态，但这一状态正在逐年改善。这一结果与张魁等[35]等利用Catch-MSY模型对蓝圆鲹的资源评估结果相似。
在经过较大强度的捕捞胁迫后，大部分鱼类的表型性状会产生适应性响应[36-37]，其中以体长和性成熟年龄的变化较为明显。本研究表明，北部湾蓝圆鲹的资源群体在长期的高捕捞压力下，其结构组成发生了显著变化，其优势体长组由1992年的130~170 mm，减小为2007年的110~140 mm，随后逐年回升，2012年为120~130 mm和140~170 mm。从生长参数来看，蓝圆鲹生长系数(K)则从1992年的0.39上升到2007的0.76，之后逐年下降，2012年为0.5，但总体趋势为波动上升。而平均体长相比1992年的160.23 mm，2012年只有149 mm，呈下降趋势；对7个世代的对比分析发现，从1992年到2012年，北部湾蓝圆鲹平均体长呈先减后增的变化趋势，拐点时间为2007年。作为渔业资源管理重要参考点的L50和A50[38]也呈现出与平均体长相同的变化趋势，L50和A50分别由1998年的180.16 mm和1.860减小至2007年的174.91 mm和0.888，之后回升至2012年的180.87 mm和2.538。表明群体虽然有逐年恢复的趋势，但其小型化的现象依然存在，且个体生长速度加快。而张静等[39]2012—2013年在北部湾口和南沙群岛西南部海域进行的底拖网调查资料表明，与1997—1999年的南海北部调查数据对比，蓝圆鲹体长下降明显，优势体长组降低近40~50 mm，其资源结构也向着低龄化、小型化转变，同为中上层鱼类的竹䇲鱼也出现与此相似的转变。
造成鱼类个体的表型特征对捕捞压力响应的原因可能是由于最小网目尺寸(minimum mesh size)和最小捕捞尺寸(minimum landing size)的影响，使捕捞具有选择性，这种人为的选择性会导致开发群体的进化过程受到影响，最终造成个体的生物学特性发生改变[40]。捕捞过程中对鱼类的表型特征(体长)会有一定的选择性[41]；捕捞死亡系数经常数倍高于自然死亡系数，这也增强了资源群体的选择性[42]；高强度的捕捞胁迫使得资源群体密度变得稀疏，降低了种间竞争及捕食压力，因此小个体可以快速补充进入渔业[43]。如Neja[44]发现大西洋鲭(Scomber scombrus)幼体的生长与其群体密度强相关，Jensan和Burns[45]认为这种强相关关系与大西洋鲭的摄食压力变小有关。尽管捕捞驱动进化(fisheries-induced evolution，FIE)的现象已在众多鱼类群体中得到证实，如孔雀鱼(Poecilia reticulata)[46]性成熟年龄降低，大西洋银汉鱼(Menidia menidia)[36]个体的平均体长减小，但当终止捕捞活动时，鱼类个体的表型特征能否发生可逆性恢复尚未有统一的结论。如de Roos等[47]通过食物网模型，发现鱼类资源的群体结构以及表型特征容易对捕捞压力产生响应，且这种响应是不可逆的。而Fukuwaka和Morita[48]分析了日本北海道区域的大麻哈鱼(Oncorhynchus keta)，认为在降低捕捞压力后，已减小的L50能够有所恢复。
尽管当前北部湾渔业资源的捕捞压力过高，但自1999年起，中国政府开始对包括北部湾在内的南海区域实施捕捞渔船、渔网指标“双控”、海洋零捕捞增长、休渔期、休渔区以及转产转业等减轻捕捞压力的举措，而同时由于上涨的燃油价格，使得捕捞成本增加，造成北部湾的开发强度呈下降趋势[49]。本研究发现，北部湾蓝圆鲹2012年的开发率已低于1992年，为7个时期中最低，这一现象与陈作志等[50]对北部湾深水金线鱼(Nemipterus bathybius) 1992年、1997年、2012年3个世代的渔业调查分析结果相似，表明在捕捞强度降低后，蓝圆鲹等鱼类个体的小型化现象趋于停止甚至好转。北部湾海域是中国与越南实施海洋渔业共同开发和管理的海域，但越南尚未有完善的休渔措施[11]。1986—2006年，越南的渔船功率增加了11.6倍，捕捞产量增加了2.9倍，激增的渔船量使得从1995年起北部湾近岸就出现过度开发现象，且这一形势已蔓延至整个北部湾海域[51]。北部湾渔业资源的养护，越方也负有重大责任和义务。因此，中国在严控北部湾捕捞产量增长的同时，也需加强与越南的沟通合作，以期尽可能加大渔业资源的修复力度。
-
![]()
图 1 南海西部西沙海域中层鱼类采样站位图
Figure 1. Sampling station for mesopelagic fish in Xisha Islands of western South China Sea
![]()
图 2 南海西沙冷泉内外区域站位中层鱼类各目种类数占比
Figure 2. Percentage of number of mesopelagic fish species by Order in sea area inside and outside cold seep in Xisha Islands of South China Sea
![]()
图 3 南海西沙冷泉内外区域中层鱼类站位聚类分析
Figure 3. CLUSTER results of mesopelagic fish species in sea area inside and outside cold seep in Xisha Islands of South China Sea
![]()
图 4 南海西沙冷泉内外区域中层鱼类PCoA主坐标排序分析
Figure 4. PCoA (Principal co-ordinates analysis) ordination of merger of mesopelagic fish in sea area inside and outside cold seep in Xisha Islands of South China Sea
表 1 南海西沙冷泉内外区域中层鱼类拖网采样站位信息
Table 1 Information of sampling stations for mesopelagic trawl in sea area inside and outside cold seep in Xisha Islands of South China Sea
站位编号
Station No.采样水层
Sampling water layer站位水深
Station depth/m平均拖速
Mean towing speed/(m·s–1)放网时段
Trawl time拖网时长
Trawl duration/hA1 深层 400~600 2.06 白天 1 A2 深层 400~600 1.90 夜间 1 A3 表层 75 2.26 夜间 1 B1 深层 400~600 1.85 白天 1 B2 深层 400~600 1.90 夜间 1 B3 表层 75 1.90 夜间 1 
下载: 导出CSV
表 2 南海西沙冷泉内外区域站位种类组成差异
Table 2 Difference of species composition in each station in sea area inside and outside cold seep in Xisha Islands of South China Sea
站位 Station A1 A2 A3 B1 B2 B3 A1 31 17 10 A2 39.24% 16 16 A3 29.31% 22.54% 15 B1 18.18% 15 6 B2 30.95% 30.0% 16 B3 40.54% 16.67% 27.12% 种类数
Species number49 61 26 16 49 26 注：对角线右上方为站位间共有种类数，对角线左下方为Jaccard相似性系数。 Note: The number of shared species between stations is above the diagonal line, and the Jaccard similarity coefficient is below the diagonal line.
下载: 导出CSV
表 3 南海西沙冷泉内外区域不同站位鱼类平均分类多样性指数及分类差异变异指数
Table 3 Average taxonomic distinctness (Δ+) and variation in taxonomic distinctness (Λ+) in each station in sea area inside and outside cold seep in Xisha Islands of South China Sea
站位
Station平均分类差异指数
Δ+分类差异变异指数
Λ+A1 70.24 275.45 A2 70.74 334.20 A3 70.16 324.31 B1 63.17 453.31 B2 64.20 454.80 B3 55.82 531.40
下载: 导出CSV
表 4 南海西沙冷泉内外区域各站位中层鱼类优势种
Table 4 Dominant species of mesopelagic fish in each station in sea area inside and outside cold seep in Xisha Islands of South China Sea
站位
Station优势种
Dominant species相对重要性
指数
IRI优势度
百分比
IRI%A1 斯氏蝰鱼 Chauliodus sloani 2 957.81 14.79 蛇口异星杉鱼 Heterophotus ophistoma 2 870.73 14.35 长钻光鱼 Gonostoma elongatus 1 335.23 6.68 天纽珍灯鱼 Lampanyctus tenuiformis 1 320.00 6.60 褶胸鱼 Sternoptyx diaphana 1 102.62 5.51 A2 斯氏蝰鱼 C. sloani 2 928.60 14.64 天纽珍灯鱼 Lampanyctus tenuiformis 1 713.85 8.57 A3 尾明角灯鱼 Ceratoscopelus warmingii 5 206.78 26.03 喀氏眶灯鱼 Diaphus garmani 4 072.03 20.36 前鳍盗目鱼 Lestidiops jayakari 2 427.66 12.14 马来眶灯鱼 D. malayanus 1 720.47 8.60 瓦氏眶灯鱼 D. watasei 1 656.03 8.28 B1 长银斧鱼 Argyropelecus affinis 5 311.29 26.56 西钻光鱼 G. atlanticum 3 272.97 16.36 灿烂眶灯鱼 D. fulgens 2 186.23 10.93 银斧鱼 A. hemigymnus 1 893.63 9.47 串灯鱼 Vinciguerria nimbara 1 611.99 8.06 黑柔骨鱼 Malacosteus niger 1 578.59 7.89 B2 斯氏蝰鱼 C. sloani 3 123.79 15.62 近壮灯鱼 Hygophum proximum 1 439.34 7.20 西钻光鱼 G. atlanticum 1 399.58 7.00 B3 尾明角灯鱼 C. warmingii 4 494.93 22.47 长鳍虹灯鱼 Bolinichthys longipes 3 239.39 16.20 喀氏眶灯鱼 D. garmani 2 692.64 13.46 串灯鱼 Vinciguerria nimbara 2 131.97 10.66 条带眶灯鱼 D. brachycephalus 1 785.02 8.93
下载: 导出CSV
表 5 南海西沙冷泉内外区域中层鱼类多样性指数
Table 5 Diversity index of mesopelagic fish in sea area inside and outside cold seep in Xisha Islands of South China Sea
站位
StationShannon-Wiener
多样性指数
H'Pielou 均匀度
指数
J'Margalef 丰富度
指数
DA1 3.47 0.89 6.45 A2 3.54 0.86 7.27 A3 2.04 0.63 2.30 B1 2.34 0.84 2.67 B2 3.32 0.85 4.37 B3 2.23 0.68 2.58
下载: 导出CSV
表 6 南海西沙冷泉内外区域中层鱼类百分比相似性分析
Table 6 Results of similarity percentage analysis (SIMPER) of mesopelagic fish in sea area inside and outside cold seep in Xisha Islands of South China Sea
物种
Species平均不相似度
Average dissimilarity/%累计贡献度
Cumulative contribution/%数量百分比 Quantity percentage/% 冷泉内
Inside cold seep冷泉外
Outside cold seep长鳍虹灯鱼 Bolinichthys longipes 5.31 10.40 0.64 8.87 串灯鱼 Vinciguerria nimbara 5.12 20.42 0.17 8.28 喀氏眶灯鱼 Diaphus garmani 3.67 27.61 27.46 10.52 近壮灯鱼 Hygophum proximum 3.48 34.42 0.60 5.91 条带眶灯鱼 Diaphus brachycephalus 3.18 40.63 0.04 5.10 西钻光鱼 Gonostoma atlanticum 3.10 46.69 0.56 5.28 马来眶灯鱼 Diaphus malayanus 2.55 51.69 12.81 3.57 多鳞孔头鲷 Melamphaes polylepis 2.23 56.06 0.68 3.98
下载: 导出CSV
-
[1] SUESS E. Marine cold seeps and their manifestations: geological control, biogeochemical criteria and environmental conditions[J]. Int J Earth Sci, 2014, 103(7): 1889-1916. doi: 10.1007/s00531-014-1010-0
[2] WANG J L, WU S G, KONG X, et al. Subsurface fluid flow at an active cold seep area in the Qiongdongnan Basin, northern South China Sea[J]. J Asian Earth Sci, 2018, 168: 17-26. doi: 10.1016/j.jseaes.2018.06.001
[3] MIRONOV A N. New taxa of stalked crinoids from the suborder Bourgueticrinina (Echinodermata, Crinoidea)[J]. Zool Stud, 2000, 79(6): 712-728.
[4] KOVA T Y, MEKHANIKOVA I V, SIDELEVA V G, et al. Trophic relationships between macroinvertebrates and fish in St. Petersburg methane seep community in abyssal zone of Lake Baikal[J]. Contemp Probl Ecol, 2017, 10(2): 147-156. doi: 10.1134/S1995425517020123
[5] CRUAUD P, DECKER C, OLU K, et al. Ecophysiological differences between vesicomyid species and metabolic capabilities of their symbionts influence distribution patterns of the deep-sea clams[J]. Mar Ecol, 2019, 40(3): e12541.
[6] MACDONALD I R, LEIFER I, SASSEN R, et al. Transfer of hydrocarbons from natural seeps to the water column and atmosphere[J]. Geofluids, 2002, 2(2): 95-107. doi: 10.1046/j.1468-8123.2002.00023.x
[7] 蒲燕萍, 孙春岩, 陈世成, 等. 南海琼东南盆地-西沙海槽天然气水合物地球化学勘探与资源远景评价[J]. 地质通报, 2009, 28(11): 1656-1661. doi: 10.3969/j.issn.1671-2552.2009.11.015 [8] 万志峰, 张伟, 陈崇敏, 等. 琼东南盆地冷泉差异发育特征及其深部控制机理[J]. 海洋地质前沿, 2021, 37(7): 1-10. [9] OLIVAR M P, BERNAL A, MOLI B, et al. Vertical distribution, diversity and assemblages of mesopelagic fishes in the western Mediterranean[J]. Deep-Sea Res I, 2013, 62: 53-69.
[10] CLAVEL-HENRY M, PIRODDI C, QUATTROCCHI F, et al. Spatial distribution and abundance of mesopelagic fish biomass in the Mediterranean Sea[J]. Front Mar Sci, 2020, 7: 573986. doi: 10.3389/fmars.2020.573986
[11] SOBRADILLO B, BOYRA G, MARTINEZ U, et al. Target strength and swimbladder morphology of Mueller's pearlside (Maurolicus muelleri)[J]. Sci Rep, 2019, 9(1): 17311. doi: 10.1038/s41598-019-53819-6
[12] BADCOCK J. The vertical distribution of mesopelagic fishes collected on the sond cruise[J]. J Mar Biolog Assoc UK, 1970, 50(4): 1001-1044. doi: 10.1017/S0025315400005920
[13] EDUARDO L N, LUCENA-FREDOU F, MINCARONE M M, et al. Trophic ecology, habitat, and migratory behaviour of the viperfish Chauliodus sloani reveal a key mesopelagic player[J]. Sci Rep, 2020, 10(1): 20996. doi: 10.1038/s41598-020-77222-8
[14] BOYD P W, CLAUSTRE H, LEVY M, et al. Multi-faceted particle pumps drive carbon sequestration in the ocean[J]. Nature, 2019, 568(7752): 327-335. doi: 10.1038/s41586-019-1098-2
[15] CAVAN E L, LAURENCEAU-CORNRC E C, BRESSAC M, et al. Exploring the ecology of the mesopelagic biological pump[J]. Prog Oceanogr, 2019, 176: 102125. doi: 10.1016/j.pocean.2019.102125
[16] BACKUS R H, CRADDOCK J E, HAEDRICH R L, et al. The distribution of mesopelagic fishes in the Equatorial and Western North Atlantic Ocean[J]. J Mar Res, 1970, 28: 179-201.
[17] DOYA C, CHATZIEVANGELOU D, BAHAMON N, et al. Seasonal monitoring of deep-sea megabenthos in Barkley Canyon cold seep by internet operated vehicle (IOV)[J]. PLOS ONE, 2017, 12(5): e0176917. doi: 10.1371/journal.pone.0176917
[18] ROSS S W, QUATTRINI A M, ROA-VARON A Y, et al. Species composition and distributions of mesopelagic fishes over the slope of the north-central Gulf of Mexico[J]. Deep-Sea Res II, 2010, 57(21/22/23): 1926-1956.
[19] 徐姗楠, 郭建忠, 范江涛, 等. 夏季大亚湾鱼类群落结构与多样性[J]. 生态学杂志, 2020, 39(4): 1254-1264. [20] SIMPSON E H. Measurement of diversity[J]. Nature, 1949, 163(4148): 688. doi: 10.1038/163688a0
[21] CLARKE K R, WARWICK R. A further biodiversity index applicable to species lists: variation in taxonomic distinctness[J]. Mar Ecol Prog Ser, 2001, 216: 265-278. doi: 10.3354/meps216265
[22] PIANKA E R. Ecology of the agamid lizard Amphibolurus isolepis in Western Australia[J]. Copeia, 1971, 3(3): 527-536.
[23] KREBS C J. Ecological methodology[M]. New York: Harper Collins Publishers, 1989: 1-624.
[24] PIELOU E C. The use of information theory in the study of ecological succession[J]. J Theor Biol, 1966, 10: 370-383. doi: 10.1016/0022-5193(66)90133-0
[25] MARGALEF R. Information theory in ecology[J]. Gen Syst, 1958, 3: 36-71.
[26] 朱晓芬, 陈彬, 俞炜炜, 等. 厦门湾大型底栖动物分类学多样性指数及分类充分性[J]. 生态学报, 2018, 38(15): 5554-5565. [27] 舒璐, 林佳艳, 徐源, 等. 基于环境DNA宏条形码的洱海鱼类多样性研究[J]. 水生生物学报, 2020, 44(5): 1080-1086. doi: 10.7541/2020.125 [28] 舒卫先, 胡菊香, 陈胜, 等. 流域治理后生态恢复期的沙颍河干流底栖动物群落结构分析[J]. 水生态学杂志, 2014, 35(4): 35-42. doi: 10.3969/j.issn.1674-3075.2014.04.006 [29] 王雪辉, 杜飞雁, 邱永松, 等. 1980—2007年大亚湾鱼类物种多样性、区系特征和数量变化[J]. 应用生态学报, 2010, 21(9): 2403-2410. [30] 许友伟, 江艳娥, 范江涛, 等. 南海中部海域秋冬季中层渔业生物群落结构的初步研究[J]. 南方水产科学, 2016, 12(4): 49-56. doi: 10.3969/j.issn.2095-0780.2016.04.006 [31] 袁梦, 陈作志, 张俊, 等. 南海北部陆坡海域中层渔业生物群落结构特征[J]. 南方水产科学, 2018, 14(1): 85-91. doi: 10.3969/j.issn.20950780.2018.01.011 [32] ZHANG J, WANG X L, JIANG Y, et al. Species composition and biomass density of mesopelagic nekton of the South China Sea continental slope[J]. Deep-Sea Res II, 2019, 167(C): 105-120.
[33] BARJAU-GONZALEZ E, RODRIGUEZ-ROMREO J, GALVAN-MAGANA F, et al. Changes in the taxonomic diversity of the reef fish community of San José Island, Gulf of California, Mexico[J]. Biodivers Conserv, 2012, 21(14): 3543-3554. doi: 10.1007/s10531-012-0378-z
[34] 李媛洁, 张俊, 陈作志, 等. 南沙群岛渚碧礁鱼类分类多样性研究[J]. 南方水产科学, 2020, 16(1): 36-41. doi: 10.12131/20190159 [35] 李圣法. 东海大陆架鱼类群落生态学研究-空间格局及其多样性[D]. 上海: 华东师范大学, 2005: 1-154. [36] 曾晓光, 李娜娜, 杨权, 等. 南沙群岛西南部陆架海域鱼类分类的多样性[J]. 水产学报, 2012, 36(4): 592-600. [37] 张衡. 鱼类分类多样性估算方法在长江河口区的应用[J]. 华东师范大学学报 (自然科学版), 2007(2): 11-22. [38] 史赟荣, 李永振, 卢伟华, 等. 东沙群岛珊瑚礁海域鱼类物种分类多样性研究[J]. 南方水产, 2009, 5(2): 10-16. [39] WANG X L, ZHANG J, ZHAO X Y, et al. Vertical distribution and diel migration of mesopelagic fishes on the northern slope of the South China Sea[J]. Deep-Sea Res II, 2019, 167(C): 128-141.
-
期刊类型引用(4)
1. 梁前才，车南青，李玲，黎涵，曹际振，李安兴. 鱼类链球菌病的研究进展及防控策略. 广东畜牧兽医科技. 2023(04): 54-59 . 
百度学术
2. 吴凡，谢彩霞，汪志文，夏洪丽，鲁义善，简纪常. 尼罗罗非鱼Crtam基因的克隆鉴定与mRNA表达分析. 基因组学与应用生物学. 2021(Z1): 1943-1951 . 百度学术
3. 吴甘林，邓玉婷，张瑞泉，谭爱萍，赵飞，姜兰. 罗非鱼不同养殖阶段喹诺酮类耐药菌群数量与组成分析. 南方农业学报. 2020(12): 3073-3082 . 百度学术
4. 袁玉梅，石存斌，陶家发，张德锋，孙承文，巩华，黄志斌，赖迎迢. 罗非鱼无乳链球菌Sip-Pgk融合蛋白乳酸菌口服疫苗制备及其免疫效果的研究. 南方水产科学. 2019(06): 9-18 . 本站查看
其他类型引用(3)
计量
- 文章访问数: 926
- HTML全文浏览量: 190
- PDF下载量: 94
- 被引次数: 7
粤公网安备 44010502001741号
