Analysis of grey relationship between morphological traits and body weight of pearl oyster (Pinctada fucata)
-
摘要:
文章采用灰色关联分析方法对合浦珠母贝(Pinctada fucata)12个家系的壳长、壳高、壳宽和绞合线长4个形态性状与体质量的关联度进行了分析和评价。结果显示，12个家系的平均体质量、壳长、壳高、壳宽和绞合线长分别为35.96 g、57.36 mm、59.24 mm、23.55 mm和56.93 mm。各个性状的变异系数较大(8.30%~21.30%)，其中壳高和壳宽的变异系数相对较小，而体质量的变异系数相对较大。体质量与各形态性状灰色关联系数之和以壳高为最大(8.950 8)，且家系间较均一，壳长为最小(8.189 8)，但家系间变化幅度较大。体质量与4个形态性状的关联度大小趋势为壳高(0.745 9)＞绞合线长(0.714 8)＞壳宽(0.705 9)＞壳长(0.682 5)。根据上述结果，可将壳高作为间接选育合浦珠母贝体质量性状的重要评价指标，以提高体质量的选育效率。
Abstract:We conducted a grey relational analysis to examine the relationship between shell traits and body weight of 12 selective breeding families of Pinctada fucata. The results show that the average body weight, shell length, shell height, shell width and hinge length were 35.96 g, 57.36 mm, 59.24 mm, 23.55 mm and 56.93 mm, respectively. The coefficients of variation in body weight and shell traits were ranked as follows: body weight (21.30%) > shell length (9.91%) > hinge length (9.12%) > shell height (8.85%) > shell width (8.30%). The grey relational degrees between each shell trait and body weight were ranked as follows: shell height (0.745 9) > hinge length (0.714 8) > shell width (0.705 9) > shell length (0.682 5). It is suggested that shell height which had the closest relationship with body weight could be regarded as an important factor for evaluation of P.fucata breeding practice for high growth lines selection and improvement of selection efficiency.
-
南澳岛位于116°53′E~117°19′E，23°11′N~23°32′N，地理位置优越，海洋生物资源丰富。沿岸海域往往是生态系统中的敏感区域，其资源环境受诸多因素的影响。近年来，由于南澳岛近岸资源的不合理开发利用，生态环境已遭受一定程度的破坏，渔业资源明显衰退。
准确高效的渔业资源评估方法是制定科学长效的生态渔业管理机制的先决条件。传统的渔业资源评估主要依托于网捕采样调查，而任何网具对于捕捞对象都具有一定选择性，并且网捕采样通常会对调查水域渔业资源及鱼类栖息环境造成一定的破坏，其应用范围通常受水域生态环境和海底地形等条件的约束[1-2]。近年来，水下记录、摄影等技术成为获取鱼类群落信息，进行鱼类资源生态监管常用的一类方法。由于此类方法易受光照、悬浮颗粒等水域环境条件的限制，导致其评估结果可信度降低，应用范围受限[3-4]。声学评估作为一种后起的生态渔业评估方法正不断推广应用并取得了较好的成效；与传统渔业调查法相比，该方法便捷、高效、采样面积大、对渔业资源和水域生境破坏力度小，更能适应当前生态渔业发展的需求，应用前景十分广阔[5-8]。声学方法自20世纪80年代引入中国，并逐渐成为中国海洋生物资源调查研究的重要手段之一[9-10]。2007—2009年，李娜娜等[11]采用水声学方法对大亚湾杨梅坑人工鱼礁水域资源声量进行了评估；张俊等[12]针对南沙南部陆架海域渔业资源于2013年进行了4个航次的声学调查，分析了鱼类资源量空间分布特征和季节变化趋势；2013年秋季，曾雷等[13]采用声学的方法评估了百色水库渔业资源状况，并分析了鱼类昼夜行为特征。声学评估法与直接捕捞法、水下摄影等方法的有机结合被认为是未来渔业资源评估与生态监测的最佳选择之一[2]。
本研究于2014年秋季和2015年春季，综合采用声学探测和拖网调查对南澳岛部分海域渔业资源进行探查，了解不同季节该海域渔业资源动态、结构组成、空间分布等特征，并分析探讨影响渔业资源评估结果的主要因素，以期为该海域渔业资源的合理开发利用及管理策略的制定提供科学依据。
1. 材料与方法
1.1 调查海域概况与调查设置
研究海域地理空间位置见图1，大致范围为117.021 9°E~117.427 8°E，23.274 3°N~23.635 8°N。累积进行了3个航次的声学调查，调查时间分别为2015年4月25日(春季)、2014年9月10日(秋季)和2014年12月14日(冬季)。由于设备故障，冬季航次声学调查结果不能用作渔业资源评估。本文基于2014年秋季与2015年春季航次声学探测数据对南澳岛海域渔业资源进行评估分析。南澳岛海域秋季声学探测区域水深范围为1.84~38.4 m，平均水深为15.61 m；春季声学走航区域水深范围为2.33~40.2 m，平均水深为24.59 m。
![]()
图 1 南澳岛海域声学调查区域与拖网站位分布图A3~A14. 秋季拖网位点；S3~S14. 春季拖网位点Figure 1. Survey area of acoustic assessment and distribution of trawl sampling sites in Nan'ao Island areaA3–A14. trawl sampling sites in autumn; S3–S14. trawl sampling sites in spring1.2 声学数据采集
声学数据采集借助挪威Simrad EY60分裂波束科学鱼探仪，工作频率分别为120 kHz (2014年秋季)和70 kHz (2015年春季)。实时经纬度信息通过GPS (Gamin GPSCSx，美国)与声学数据同步采集，各次调查科学鱼探仪的主要技术参数见表1。由于Simrad EY60声学评估硬件系统缺乏长期的稳定性，且换能器性能受不同环境介质的影响[14]，故各季度调查前按照国际通用方法对科学鱼探仪系统的收发增益系数进行现场校正[15]。为防震减噪，换能器通过导流罩内固定于船体右舷外侧且不与船体接触，吃水0.8 m，走航速度约5~7 kn。通常声学调查航线依据调查区域的特征设计成平行式或“之”字式，受调查海域自然地理状况及其他人为因素的限制，声学调查航迹并不规则(图1)。参照Aglen[16]覆盖率计算公式
$D{\rm{ = }}L/\sqrt A $ ，式中L为声学航程，A为调查海域面积。2次调查声学评估海域覆盖率D分别为7.34 (2014年秋季)和8.07 (2015年春季)。表 1 EY60科学鱼探仪主要技术参数设定Table 1. Setting of main technical parameters of scientific echosounder EY60技术参数
technical parameter2014年秋季
autumn 20142015年春季
spring 2015换能器频率/kHz transducer frequency 120 70 发射功率/W transmitting power 200 300 脉冲宽度/μs pulse duration 128 512 等效波束角/dB equivalent beam angle –21 –21 换能器增益/dB transducer gain 26.8 25.69 横向波束宽度/(°) athwardship 3 dB beam width 7 6.53 纵向波束宽度/(°) alongship 3 dB beam width 7 6.43 吸收系数/dB·km–1 absorption coefficient 45.3 18.99 声速/m·s–1 sound speed 1 531.82 1 536.75 1.3 生物学数据获取
调查海域渔业资源生物学信息的获取主要通过底拖网进行，用以辅助渔业资源声学回波映像的识别与积分分配。共设置8个站位(图1)，底拖网信息见表2。各站位渔获样品均进行现场分类并计数取样，数量少于50尾的物种全部取样；若单一物种渔获数量大于50尾，则根据其大小组成特点按比例随机取样，并测量记录每种渔获的长度(体长/叉长/胴长)和体质量，体长精确至1 cm，体质量精确至1 g。
表 2 主要生物学取样设置Table 2. Main setting of biological sampling主要参数
main parameter2014年秋季
autumn 20142015年春季
spring 2015采样类型 sample type 底拖网 底拖网 发动机功率/kW engine power 280 280 平均拖网速度/km·h–1 average trawling speed 5.56 4.93 平均采样时长/h average trawling time 0.55 0.46 网口宽度/m net width 6 4 网口高度/m net height 2 2 网目/cm net size 3×3 3×3 1.4 数据处理与分析
声学数据的处理主要借助专业分析软件Echoview 6.1 (Myriax)，根据回波映像特征，剔除明显的环境噪声数据。利用海底检测功能板块生成探测区域海底轮廓并进行手动编辑。重新设置积分起始水层为1.5 m以屏蔽海表航行噪声干扰，积分终止水层设置为海底之上0.5 m以排除海底回波信号干扰(图2)。参照多种类海洋渔业资源声学评估方法[5,17]，基本积分航程单元(elementary distance sampling unit，EDSU，n mile)设置为0.5 n mile，以底拖网采样的渔获物组成信息作为回波积分分配的主要依据。通过阈值响应(threshold response)模块分析Sv阈值的变化对回波面积积分(NASC)的影响，从而大致判断目标群体声学回波的Sv可信阈值区间。目标群体与非目标群体之间Sv并无明显界线，一般为NASC随Sv加速变化的拐点[17]。Sv积分阈值设置为–70 dB，以减小非目标信号对回波积分的影响。同时根据单体回波映像，利用Echoview软件中单体检测和轨迹追踪模块[18]，进而获得回波单体大小、频率组成与垂直空间分布特征。
1.5 鱼类资源密度评估方法
1.5.1 回波积分法
1) 在不考虑渔获种类组成及体长结构分布的情况下，根据目标强度现场测定法，各基本积分航程单元内渔业资源数量密度(number density，尾·km–2)计算式为
$$ \quad\quad\quad\quad\quad\quad\quad\overline \rho = \frac{{\rm NASC}}{{4\pi \cdot {\sigma _{bs}} \cdot {{1.852}^2}}} $$ (1) 式中σbs为后向散射截面(backscattering cross-section，m2)，它与目标强度(target strength，TS)的关系为
$$ \quad\quad\quad\quad\quad\quad\quad TS = 10\log {\sigma _{bs}} $$ (2) 2) 参照多种类海洋渔业资源声学评估方法，以拖网采样的渔获物信息作为积分值分配的主要依据，分析区域内第i物种的数量密度
${\rho _{i,a}}$ (number density，尾·km–2)和资源量密度${\rho _{i,b}}$ (biomass density，t·km–2)分别为$$ \quad\quad\quad\quad\quad\quad\quad {\rho _{i,a}} = {c_i}\frac{{\rm NASC}}{{4\pi \overline \sigma \cdot {{1.852}^2}}} $$ (3) $$ \quad\quad\quad\quad\quad\quad\quad {\rho _{i,b}} = {\rho _{i,a}}\overline {{w_i}} {10^{ - 6}} $$ (4) 式中ci (%)为分析海域内第i物种渔获物数量百分比，
$\overline \sigma $ 为分析海域内所有声学评估种类的平均后向散射截面，$\overline {{w_i}} $ 为第i物种的平均体质量(body mass，g)，且有$$ \quad\quad\quad\quad\quad\quad\quad \overline \sigma = \sum\limits_{i = 1}^n {{c_i}{{10}^{\frac{{T{S_i}}}{{10}}}}} $$ (5) $$ \quad\quad\quad\quad\quad\quad\quad T{S_i} = 20\log {L_i} + {b_{20,i}} $$ (6) 式中TSi为第i物种的目标强度，n为参与声学评估的物种数，Li为第i物种平均体长(body length，cm)，
${b_{20,i}}$ 为第i物种的参考目标强度。1.5.2 扫海面积法
游泳生物资源量密度B (kg·km–2) 根据扫海面积法估算，公式为
$$ \quad\quad\quad\quad\quad\quad\quad B = Y/A(1-E) $$ (7) 式中Y为平均渔获率 (kg·h–1)；B为资源量密度 (kg·km–2)；A为每小时扫海面积(km2·h–1)；E为逃逸率，取0.5。
2. 结果
2.1 渔获组成特征
2014年秋季，捕获游泳生物和底栖无脊椎动物共73种，其中鱼类40种，头足类5种，虾类6种，蟹类16种，虾蛄类6种，总渔获数量为2 598尾，总渔获量为51.19 kg。2015年春季共捕获游泳生物和底栖无脊椎动物共54种，其中鱼类26种，头足类3种，虾类6种，虾蛄类4种，蟹类15种，总渔获数量为6 251尾，总渔获量为28.33 kg。
为排除海底回波信号干扰，历次调查海底之上0.5 m内均被视为声学探测的盲区，故底栖的鲆鲽类、鰕虎鱼类、蛸类、虾蟹类等非常贴底的生物均不参与声学评估。
根据生物学拖网采样结果，声学评估种类中参与积分分配的种类主要为相对重要性指数(IRI)大于100的常见种和优势种[19]。不同季节，声学评估鱼类组成存在较大差异，2014年秋季参与声学评估的种类主要有39种；2015年春季参与评估的种类主要有18种(表3)。
表 3 声学评估优势种类生物学信息Table 3. Biological information of dominant species in acoustic assessment季节(年份)
season (year)物种
species数量百分比/%
quantity percentage生物量百分比/%
biomass percentage相对重要性指数
IRI体长范围/mm
body length range体质量范围/g
body mass range秋季 (2014)
autumn短尾大眼鲷
Priacanthus macracanthus45.3 34.8 6 678 72~88 8~21 丝背细鳞鲀
Stephanolepis cirrhifer12.0 20.8 2 187 68~107 6~55 丽叶鲹
Caranx kalla9.3 6.5 1 322 60~135 4~85 杜氏枪乌贼
Loligo duvaucelii4.9 6.5 1 134 50~123* 6~58 鹿斑鲾
Secutor ruconius7.4 2.3 647 27~70 2~12 二长棘鲷
Parargyrops edita3.5 9.0 624 95~108 32~56 斑鳍白姑鱼
Argyrosomus pawak6.2 2.2 280 54~105 5.3~30 黄斑篮子鱼
Siganus canaliculatus0.6 0.9 129 74~126 5~74 春季 (2015)
spring鹿斑鲾
Secutor ruconius86.9 68.0 11 613 36~63 1.6~10 二长棘鲷
Parargyrops edita5.9 5.4 845 35~62 1.7~58 杜氏枪乌贼
Loligo duvaucelii1.0 6.9 689 42~182* 4.9~140 竹䇲鱼
Trachurus japonicus3.2 4.1 458 60~90 2.4~11 刺鲳
Psenopsis anomala0.5 2.6 229 65~110 8.5~45 短吻鲾
Leiognathus brevirostris1.1 3.6 179 58~94 4.2~24 皮氏叫姑鱼
Johnius belangerii0.3 2.6 145 92~160 14~75 沙带鱼
Lepturacanthus savala0.5 1.8 114 66~150 3.6~45 注：*. 胴长
Note: *. mantle length参与声学评估种类渔业资源密度扫海面积法计算结果如图3所示：2014年秋季平均数量密度为63 246尾·km–2，生物量密度为1 075.98 kg·km–2；2015年春季平均数量密度为221 984尾·km–2，生物量密度为949.74 kg·km–2。
![]()
图 3 调查海域各拖网位点参与声学评估种类数量与生物量密度Figure 3. Quantity and biomass density of acoustic assessed species in each trawling site of survey area2.2 声学评估种类资源密度
1) 假定调查海域渔获种类组成单一，根据目标强度现场测定法，2014年秋季和2015年春季渔业资源数量密度分别为(64 542±132 583)尾·n mile–2和(25 472±45 814)尾·n mile–2；生物量密度分别为(1 053.32±2 163.75) kg·n mile–2和(102.4±184.17) kg·n mile–2。
2) 依据多种类渔业资源积分分配原则，调查海域2014年秋季和2015年春季参与声学积分分配的种类，其平均数量密度分别为(26 806±39 968)尾·n mile–2和(34 416±384 97)尾·n mile–2；平均生物量密度分别为(403.68±601.9) kg·n mile–2和(124.5±139.26) kg·n mile–2。
在假定种类组成单一和混合2种不同条件下，单因素方差分析结果见图4。2015年春季渔业资源数量密度(P=0.113>0.05)与生物量密度(P=0.310>0.05)均未表现出显著性差异。而2014年秋季资源数量密度(P=0.007<0.05)和生物量密度(P=0.004<0.05)在2种不同条件下均表现出显著性差异。
![]()
图 4 种类单一和混合2种不同条件下调查海域渔业资源声学评估误差线. 地理空间差异；**. α=0.05水平下差异显著Figure 4. Acoustic assessment of fishery resources in survey area in single or multi-species assumptionsError bar. geographical spatial difference; **. significant difference at α=0.05 level根据多种类渔业资源声学评估积分分配原则，调查海域各季度不同种类数量密度与生物量密度组成情况见表4。2014年秋季调查海域优势种类主要为短尾大眼鲷(Priacanthus macracanthus)，其数量密度和生物量密度分别为12 142.92 尾·n mile–2和150.22 kg·n mile–2。2015年春季调查海域优势种类主要为鹿斑鲾(Secutor ruconius)，其数量密度和生物量密度分别为27 791.47尾·n mile–2和87.42 kg·n mile–2。整体上调查海域渔业资源较为贫乏，且个体偏小，应多予保护。
表 4 调查海域声学评估种类丰度与生物量密度Table 4. Abundance and biomass density of acoustic assessed species in survey area季节(年份)
season (year)种类
speciesb20/dB 丰度密度/尾·n mile–2
abundance density生物量密度/kg·n mile–2
biomass density秋季 (2014)
autumn短尾大眼鲷 Priacanthus macracanthus –68.0 12 142.92 150.22 丝背细鳞鲀 Stephanolepis cirrhifer –76.0 3 216.67 87.92 丽叶鲹 Caranx kalla –72.5 2 492.92 27.80 鹿斑鲾 Secutor ruconius –72.5 1 983.61 9.67 斑鳍白姑鱼 Argyrosomus pawak –68.0 1 661.95 9.46 杜氏枪乌贼 Loligo duvaucelii –80.0 1 313.47 28.32 二长棘鲷 Parargyrops edita –68.0 938.19 39.55 黄斑篮子鱼 Siganus canaliculatus –72.9 160.83 3.50 春季 (2015)
spring鹿斑鲾 Secutor ruconius –72.5 27 791.47 87.42 二长棘鲷 Parargyrops edita –68.0 1 886.88 6.94 竹䇲鱼 Trachurus japonicus –72.5 1 023.39 5.37 杜氏枪乌贼 Loligo duvaucelii –80.0 319.81 9.62 短吻鲾 Leiognathus brevirostris –72.5 351.79 4.58 刺鲳 Psenopsis anomala –80.0 159.90 3.53 沙带鱼 Lepturacanthus savala –66.1 159.90 2.55 皮氏叫姑鱼 Johnius belangerii –68.0 95.94 3.66 注：b20值来自文献记载
Note: The b20 values come from documentary records.调查海域渔业资源数量在地理空间上的分布见图5，其结果表明：1)调查海域渔业资源基本处于较贫乏水平。2)不同季节，声学评估种类存在较大差异，仅鲾类在2个季度调查中均占有较高的比重，枪乌贼类为2个季度常见种。3)不同季节，渔业资源在空间尺度上分布极不均匀，且无明显的空间梯度变化趋势，仅在临近南澳岛东北部(红色标记区域)存在小范围的渔业资源密集区。
2.3 回波单体TS频率组成与垂直分布
据单回波检测结果(图6)，以2 dB为分组单元，调查海域2014年秋季回波单体TS变化区间为–64~–33 dB，其中–54 dB以下个体累积数量百分比高达90.37%。此外，根据回波单体在垂直空间上的分布特征，以15 m水深为界可将调查水域大致分为2~15 m (中上层)和15~35 m (中下层)两个水层。中下层水域以小型个体为主，回波单体TS基本在–54 dB以下；中上层个体大小组成相对均匀。
![]()
图 6 调查海域回波单体频率组成与垂直分布Figure 6. Frequency and vertical distribution of single echo detection in survey area2015年春季调查海域回波单体TS变化范围为–66~–33 dB，各分组单元TS频率分布大致呈逐级下降的变化趋势。在垂直空间上，以15 m水深为界可将调查海域大致分为2~15 m (中上层)和15~35 m (中下层) 2个水层。不同水层，回波单体TS大小组成表现出相似的特征，随回波单体TS的增大其所占数量百分比逐渐下降。
3. 讨论
3.1 渔业资源时空分布特征
各站位拖网渔获分析结果显示，南澳岛海域2014年秋季和2015年春季参与声学评估种类个体均质量分别为16.32 g·尾–1和4.02 g·尾–1。不同季节，渔业资源声学评估呈现出相似的特征，秋季鱼类个体相对较大。此外，南澳岛海域2次调查回波单体TS频率组成亦表现出明显的季节差异。较2014年秋季，2015年春季调查海域回波单体TS变化幅度更大，且小型个体所占比重更高。张俊等[5]基于南海拓林湾渔业资源研究指出，产卵引起的幼鱼增加及休渔结束后对大个体鱼类的大规模捕捞是导致鱼类个体均质量季节差异的主要原因。连玉喜等[20]研究显示，春季是大多数鱼类的繁殖期，通过增殖，水体中仔稚鱼密度大大增加，经过夏季和秋季的生长，鱼类平均个体大小明显大于春季。上述研究结果表明，调查海域鱼类资源的季节动态规律受鱼类繁殖、生长、死亡等因素的共同调节，结果符合生态系统中种群动态理论与能量传递规律[21-24]。
基于扫海面积法，不同站位渔业资源密度在水平空间上存在极大差异，2次调查均有极大异常值出现。结果可能与特定鱼类生理生活习性相关，鹿斑鲾在S14站位及短尾大眼鲷在A3站位的大规模出现，对调查海域不同季节渔业资源水平空间分布差异产生了较大影响。根据回波单体在垂直方向上的大小组成分析可知，以15 m水深为分水岭，秋季中小型个体有明显向中下层水域迁移的趋势。相关研究表明，鱼类在时空尺度上的分布受气候、水文、水质、饵料及人类活动等诸多因素的影响[25-28]。由于缺乏相应的理化环境数据，对于影响调查水域鱼类资源时空变动的主要生物/非生物因子还有待进一步研究。
3.2 声学评估误差分析
声学评估作为海洋渔业资源评估与鱼类行为分析的一种有效方法，其评估结果受诸多因素的影响。传统的声学调查结合调查区域具体的地形地貌特征将调查航线设计为“之”字型或平行断面型2种，由于调查海域地理环境复杂且海上阻隔较多，导致声学调查航迹不规则，增加了声学取样产生随机误差的可能性[11]。受声学近场效应及海底声学探测盲区的影响，调查海域声学积分起始水层设置为海表1.5 m以下至海底0.5 m之上，故临近海表浮游性鱼类和底栖种类不在回波积分范围之内，从而给声学评估结果带来一定的偏差[29-30]。国内外相关研究结果表明，部分鱼类对船舶噪声存在明显的逃避行为，从而使声学评估结果偏小[31-32]。调查海域春季优势群体鹿斑鲾受船舶噪声的影响有明显的跃水现象，可能对声学评估结果带来一定影响。此外，声学回波映像处理结果显示，调查海域各季度均有小型鱼群聚集的现象，2014年秋季在中下层这类集群现象更为显著。据回波积分原理，当研究水域渔业资源存在集群现象时，遮蔽效应往往使资源数量密度评估结果偏低[33]。此外，受调查水域鱼类生理生活习性及理化环境因子的调节，鱼类资源在垂直方向上的组成与分布可能存在一定差异[34-35]。因此，渔业资源底层拖网数据用于整个水域渔业资源回波积分分配是影响渔业资源声学评估结果的另一重要因素。
3.3 不同评估方法的比较
基于回波积分原理分析比较了在鱼种单一和混合2种不同条件下，渔业资源声学评估结果的差异。研究显示，2015年春季渔获种类组成单一和混合2种不同条件下，调查水域渔业资源密度差异不显著；而2014年秋季不同条件下渔业资源声学评估结果存在显著性差异，这可能与调查海域不同季节参与声学评估种类优势集中程度相关。2015年春季调查水域用于声学评估种类优势集中程度极高，鹿斑鲾在渔获数量中所占比重高达86.87%，该时期很可能为鹿斑鲾繁殖盛期。在渔获种类组成单一和混合2种不同的假设前提下，调查海域该季度参与声学评估回波积分分配的种类主要为鹿斑鲾，因而评估结果基本一致。由此可知，在渔获优势集中度较高的特定水域，使用基本积分单元内平均目标强度现场测定法可作为简化渔业资源声学评估方法的一种有效技术手段。而2014年秋季调查水域参与声学评估种类优势集中程度相对较低，第一优势群体所占比重仅45.31%。单因素方差分析结果显示，在2种不同条件下渔业资源声学评估结果存在显著性差异。参照回波积分分配原则，调查海域渔获优势集中程度通过基本声学航程单元内平均声学反射系数(
$\overline \sigma $ )的估算 [公式(5)]，对渔业资源声学评估结果产生了影响。此外，本研究分别利用声学评估法和扫海面积法对目标海域不同季节渔业资源密度进行了评估与分析，结果表明，各季度渔业资源声学评估结果明显低于扫海面积法计算结果。扫海面积法作为传统的渔业资源评估方法存在诸多局限。理论上，参照扫海面积法计算公式中，逃逸率E是影响评估结果的重要因子之一。国内外相关研究表明，渔业资源捕捞效率受网具类型、规格、地理形态特征、鱼类行为等众多因素的影响[36-37]。任何调查网具对于渔业资源生物学信息的获取均具有一定的地域局限性和选择性[38-39]。拖网捕捞作为一种主动的渔捞作业方式，其捕捞效率往往高于其他定置渔捞方法。复杂的地理环境往往能为游泳生物提供天然的庇护场所，因此，复杂的水域环境逃逸率往往较高[5,40]。受鱼类昼夜和季节迁移习性的调节，鱼类资源时空分布特征可能存在明显的差异，从而对渔业捕捞效率产生影响。相关研究显示，鱼类在夜间活动较为频繁，且有向表层迁移的现象；而在白天鱼类有向底层迁移的规律，且活动较为稀少[25,37]。此外，鱼类随季节水温的变化有从近岸向深水水域迁移的趋势[34]。因此，针对不同季节，底拖网地理位置的规划设计亦会对捕捞效率造成一定的影响。
扫海面积法采用站点抽样调查的方式，其覆盖率明显低于声学评估方法。渔业资源空间分布存在较大差异时，应尽量增加抽样站位数以削弱随机抽样产生的误差。本研究扫海面积法及声学评估结果均表明，调查海域渔业资源空间分布差异极大。水平方向2014年秋季A3站位短尾大眼鲷资源密度明显高于其他站位，2015年春季S14位点鹿斑鲾资源密度显著高于其他站位。而根据回波单体垂直分布特征，2014年秋季小型个体有向中下层聚集的趋势，故仅基于底拖网渔获数据分析整个水域渔业资源密度可能高估了该水域渔业资源状况。因此，在渔业资源垂直空间分布差异较大时，建议采用分层拖网的取样方法。综上分析可知，渔业资源空间分布特征是影响扫海面积法渔业资源评估结果的另一重要因子。而有关研究结果显示，生物/非生物环境参数是影响渔业资源变动与空间分布的关键要素[26,28,35,41]。因此，在今后的渔业资源研究工作中应同步结合理化环境因子时空梯度变化等方面的研究，以期更加科学地开展渔业资源评估工作。
声学评估法作为一种新兴的生态渔业评估手段，与传统渔业调查法相比，具有科学、便捷、高效且对渔业资源与水域生境的破坏力度小等特点[5-6]，因而更能适应现时代渔业调查的需要，其应用前景十分广阔。然而，该方法在种类鉴定、映像识别、噪声消除等方面还有待进一步研究。因此，以声学评估为主体，综合采用传统生物学取样技术是当前渔业资源评估的有效方法，能为渔业资源科学管理及合理开发利用提供重要理论支撑。
-
![]()
图 1 合浦珠母贝不同家系5个性状的箱线图
盒子由中位数、第一四分位数和第三四分位数组成；温和异常值和极端异常值分别用“○”和“*”表示
Figure 1. Box plots of five traits of different families of P.fucata
The medium, 25th and 75th percentiles are plotted as boxes; mild outliers and extreme outliers are shown as ○ and *, respectively.
表 1 合浦珠母贝不同家系生长性状无量纲化结果
Table 1 Nondimensionalization of growth traits of different P.fucata families
性状trait 变量编号No. of variables 家系编号No.of family 1 2 3 4 5 6 7 8 9 10 11 12 体质量/g body weight X′0 -0.261 4 0.443 2 1.005 9 -0.901 4 -1.628 4 -0.089 6 -1.048 3 -0.704 7 -0.107 1 1.777 8 1.167 8 0.346 1 壳长/mm shell length X′1 0.310 3 0.742 5 0.7721 -1.009 9 -0.702 0 -0.382 3 -1.385 8 0.298 5 1.381 9 1.660 1 -0.873 7 -0.811 6 壳高/mm shell height X′2 -0.002 8 0.239 1 0.580 2 -1.049 5 -1.241 7 -1.109 2 -1.341 0 -0.042 5 1.388 5 1.584 0 0.477 6 0.517 3 壳宽/mm shell width X′3 0.670 5 0.345 0 1.189 4 -0.143 3 -1.221 6 0.741 8 -0.489 1 -0.946 9 -1.964 2 1.362 3 0.223 0 0.233 1 绞合线长/mm hinge length X′4 -0.313 2 0.728 9 0.403 5 -1.214 5 -0.763 9 -1.189 5 -1.311 6 0.553 7 0.547 4 2.059 0 0.397 2 0.103 0 
下载: 导出CSV
表 2 合浦珠母贝不同家系关联系数母序列和子序列的绝对差值
Table 2 Absolute difference between reference array and comparative array of different P.fucata families
绝对差值absolute difference (Δi) 家系编号No.of family 1 2 3 4 5 6 7 8 9 10 11 12 Δ1 0.571 7 0.299 3 0.233 8 0.108 5 0.926 4 0.292 7 0.337 5 1.003 2 1.489 0 0.117 7 2.041 5 1.157 7 Δ2 0.258 6 0.204 1 0.425 7 0.148 1 0.386 7 1.019 6 0.292 7 0.662 2 1.495 6 0.193 8 0.690 2 0.171 2 Δ3 0.931 9 0.098 2 0.183 5 0.758 1 0.406 8 0.831 4 0.559 2 0.242 2 1.857 1 0.415 5 0.944 8 0.113 0 Δ4 0.051 8 0.285 7 0.602 4 0.313 1 0.864 5 1.099 9 0.263 3 1.258 4 0.654 5 0.281 2 0.770 6 0.243 1
下载: 导出CSV
表 3 合浦珠母贝5个性状的描述性统计量
Table 3 Descriptive statistics of five traits of P.fucata
性状trait 平均值mean 最大值maximum 最小值minimum 极差ranges 标准差SD 变异系数/% CV 体质量/g body weight 35.96 68.81 17.27 51.54 7.66 21.30 壳长/mm shell length 57.36 80.84 43.82 37.02 5.69 9.91 壳高/mm shell height 59.24 81.58 45.07 36.51 5.24 8.85 壳宽/mm shell width 23.55 31.50 17.25 14.25 1.96 8.30 绞合线长/mm hinge length 56.93 78.68 42.36 36.32 5.19 9.12
下载: 导出CSV
表 4 合浦珠母贝不同家系体质量与各形态性状的灰色关联系数
Table 4 Grey relational coefficient between body weight and each morphological trait of different P.fucata families
家系编号No. of family 关联系数(ξi) relational coefficient ξ1 ξ2 ξ3 ξ4 1 0.673 5 0.838 4 0.549 3 1.000 0 2 0.812 5 0.875 7 0.958 5 0.821 0 3 0.854 9 0.741 5 0.890 6 0.660 8 4 0.949 8 0.917 6 0.602 9 0.804 1 5 0.550 8 0.762 1 0.751 3 0.568 9 6 0.816 6 0.525 7 0.579 1 0.505 8 7 0.789 7 0.816 6 0.678 9 0.835 3 8 0.529 9 0.637 3 0.849 2 0.470 6 9 0.427 4 0.426 2 0.372 7 0.640 2 10 0.942 1 0.883 1 0.746 8 0.823 8 11 0.350 2 0.626 9 0.545 7 0.598 7 12 0.492 3 0.899 8 0.946 0 0.848 6 合计total 8.189 7 8.950 9 8.471 0 8.577 8
下载: 导出CSV
表 5 合浦珠母贝不同家系各形态性状对体质量的灰色关联度
Table 5 Grey relational degree of each shell trait of different P.fucata families
性状trait 关联度relational degree 关联序relational order 壳长/mm shell length 0.682 5 4 壳高/mm shell height 0.745 9 1 壳宽/mm shell width 0.705 9 3 绞合线长/mm hinge length 0.714 8 2
下载: 导出CSV
-
[1] 郭华阳, 陈明强, 王雨, 等. 黄边糙鸟蛤野生群体主要经济性状间的相关性及通径分析[J]. 南方水产科学, 2013, 9(2): 1-8. doi: 10.3969/j.issn.2095-0780.2013.02.001 [2] 刘伟, 苏胜彦, 董在杰, 等. 3个鲤群体的微卫星标记与生长性状相关性分析[J]. 南方水产科学, 2012, 8(3): 17-24. doi: 10.3969/j.issn.2095-0780.2012.03.003 [3] 肖述, 喻子牛, 符政君. 香港巨牡蛎雌雄群体的数量性状通径分析[J]. 南方水产科学, 2011, 7(4): 1-9. doi: 10.3969/j.issn.2095-0780.2011.04.001 [4] 刘永新, 刘英杰, 周勤, 等. 牙鲆主要生长性状与体质量的灰色关联度分析[J]. 中国水产科学, 2014, 21(2): 205-213. doi: 10.3724/SP.J.1118.2014.00205 [5] 王玉涛, 郭丽君, 邓勇, 等. 利用皮特兰作父系商品猪肉品质的灰色关联度分析[J]. 上海交通大学学报: 农业科学版, 2008, 26(6): 524-527. doi: 10.3969/j.issn.1671-9964.2008.06.005 [6] 李世凯, 江敏, 潘若雷, 等. 罗氏沼虾养殖塘叶绿素a与水质因子的灰色关联分析[J]. 生态与农村环境学报, 2014, 30(2): 257-261. doi: 10.3969/j.issn.1673-4831.2014.02.019 [7] 权月伟, 米换房, 翟雷霞, 等. 棉花品种主要数量性状与皮棉产量的灰色关联分析[J]. 棉花科学, 2013, 35(4): 33-36. doi: 10.3969/j.issn.2095-3143.2013.04.006 [8] 马文清, 郭强, 邱文武, 等. 利用灰色关联分析法评价甘蔗自育新品系[J]. 广东农业科学, 2014, 41(5): 71-78. doi: 10.3969/j.issn.1004-874X.2014.05.017 [9] 吕耀平. 影响越冬罗非鱼死亡率环境因子的灰色关联分析[J]. 水产科学, 2005, 24(8): 14-16. doi: 10.3969/j.issn.1003-1111.2005.08.005 [10] 苏胜彦, 董在杰, 曲疆奇, 等. 3个鲤群体杂交后代生长性状的灰色关联及复合杂交后代的体重预测分析[J]. 水产学报, 2011, 35(1): 20-26. doi: 10.3724/SP.J.1231.2011.17109 [11] 龙敏明, 黄桂菊, 邹记兴, 等. 育珠对合浦珠母贝N19和Prismalin-14基因表达水平的影响[J]. 南方水产科学, 2013, 9(5): 58-63. doi: 10.3969/j.issn.2095-0780.2013.05.010 [12] 罗会, 刘宝锁, 黎火金, 等. 合浦珠母贝不同地理群体的形态差异和判别分析[J]. 广东农业科学, 2013, 40 (12): 171-174. doi: 10.3969/j.issn.1004-874X.2013.12.052 [13] TAKASHI A, TAKASHI I, NARIAKI I, et al. Post-operative care of implanted pearl oysters Pinctada fucata in low salinity seawater improves the quality of pearls[J]. Aquaculture, 2014, 422/423: 232-238. doi: 10.1016/j.aquaculture.2013.12.022
[14] OKAMOTO C, KOMARU A, HAYASHI M, et al. Variation of shell-closing strength among several families in pearl oyster, Pinctada fucata martensii[J]. Aquac Sci, 2006, 54(4): 525-529. doi: 10.11233/aquaculturesci1953.54.525
[15] ISHIKAWA T, OKAMOTO C, HAYASHI M, et al. The heritability of the shell-closing strength in the Japanese pearl oyster, Pinctada fucata martensii[J]. Aquac Sci, 2009, 57(1): 77-82. doi: 10.11233/aquaculturesci.57.77
[16] 黎火金, 刘宝锁, 罗会, 等. 合浦珠母贝幼贝生长性状的遗传参数估计[J]. 中国水产科学, 2013, 20(6): 1-6. doi: 10.3724/SP.J.1118.2013.01182 [17] HE M X, GUAN Y Y, YUAN T, et al. Realized heritability and response to selection for shell height in the pearl oyster Pinctada fucata (Gould)[J]. Aquac Res, 2008, 39(8): 801-805. doi: 10.1111/j.1365-2109.2008.01889.x
[18] DENG YW, FU S, DU X D, et al. Realized heritability and genetic gain estimates of larval shell length in the Chinese pearl oyster Pinctada martensii at three different salinities[J]. N Am J Aquac, 2009, 71(4): 302-306. doi: 10.1577/A08-024.1
[19] YU S, HE M X. Differential gene expression identified by RNA-Seq and qPCR in two sizes of pearl oyster (Pinctada fucata)[J]. Gene, 2014, 538(2): 313-322. doi: 10.1016/j.gene.2014.01.031
[20] LIU W G, HUANG X D, LIN J S, et al. Effect of temperature on gene expression in the pearl oyster Pinctada fucata[J]. J Ocean Univ China, 2014, 13(3): 509-515. doi: 10.1007/s11802-014-2147-y
[21] 喻达辉, 王小玉, 黄桂菊, 等. 合浦珠母贝遗传连锁图谱的构建[J]. 中国水产科学, 2007, 14(3): 361-368. https://www.nstl.gov.cn/paper_detail.html?id=f3613d6f72cb32b3e95aeffe0a13e956 [22] 何毛贤, 史兼华, 林岳光, 等. 马氏珠母贝生长性状的相关分析[J]. 海洋科学, 2006, 30(11): 1-4. https://cstj.cqvip.com/Qikan/Article/Detail?id=23295382&from=Qikan_Article_Detail [23] DENG Y W, DU X D, WANG Q H, et al. Correlation and path analysis for growth traits in F1 population of pearl oyster Pinctada martensii[J]. Mar Sci Bull, 2008, 10(2): 68-73. https://www.researchgate.net/publication/267423852_Correlation_and_Path_Analysis_for_Growth_Traits_in_F_1_Population_of_Pearl_Oyster_Pinctada_martensii
[24] 陈飞飞, 黄桂菊, 陈明强, 等. 合浦珠母贝三亚养殖群体生长性状的相关与通径分析[J]. 广东农业科学, 2012, 39(9): 122-125. doi: 10.3969/j.issn.1004-874X.2012.09.040 [25] 谷龙春, 李金碧, 喻达辉, 等. 合浦珠母贝双列杂交家系的建立与遗传分析[J]. 水产学报, 2010, 34(1): 26-31. doi: 10.3724/SP.J.1231.2010.06498 [26] 李建立, 王春德, 李朝霞, 等. 紫扇贝和海湾扇贝杂交家系的生长和通径分析[J]. 海洋科学, 2012(8): 15-20. https://www.nstl.gov.cn/paper_detail.html?id=13d0b46ad914629fc344d8e145bcf16c [27] 孙泽伟, 郑怀平, 杨彦鸿, 等. 近江牡蛎养殖群体数量性状间的相关及通径分析[J]. 中国农业通报, 2010, 26(6): 332-336. doi: 10.11924/j.issn.1000-6850.2009-2269 [28] 刘广文, 林坚士, 何毛贤. 不同贝龄华贵栉孔扇贝数量性状的通径分析[J]. 南方水产科学, 2012, 8(1): 43-48. doi: 10.3969/j.issn.2095-0780.2012.01.007 [29] 严福升, 王志刚, 刘旭东, 等. 3月龄牙鲆形态性状对体质量的通径分析[J]. 渔业科学进展, 2010, 31(2): 45-50. http://journal.yykxjz.cn/yykxjz/ch/reader/view_abstract.aspx?file_no=20100207&flag=1 [30] 王凯, 刘海金, 刘永新, 等. 牙鲆形态性状对体重的影响效果分析[J]. 上海水产大学学报, 2008, 17(6): 655-660. http://shhydxxb.ijournals.cn/shhy/article/abstract/200806161 -
期刊类型引用(5)
1. 石娟，刘永，李纯厚，王腾，赵金发，宋晓宇，谢宏宇. 珠江口棘头梅童鱼营养生态位研究. 南方水产科学. 2024(03): 56-65 . 
本站查看
2. 郑秋实，刘永，李纯厚，吴鹏，肖雅元，林琳，刘玉，邹剑. 季节变化对珠江口棘头梅童鱼肠道和水体中细菌群落结构的影响. 南方水产科学. 2024(04): 77-87 . 本站查看
3. 马菁菁，陈海刚，张喆，田斐，唐振朝，熊倩，张林宝. 2018—2023年珠江口鱼类群落结构变化及其与环境因子的关系. 南方水产科学. 2024(06): 62-73 . 本站查看
4. 熊朋莉，徐姗楠，陈作志，张帅，蒋佩文，范江涛. 珠江河口棘头梅童鱼时空分布以及影响因素的研究进展. 海洋科学. 2022(08): 79-87 . 百度学术
5. 汪润芝，原峰，李崇淑，龙焰，许清燕，奚冰柔，冼嘉俊，叶锦韶. 海上风电场建设与运行对渔业资源群落结构及空间分布的影响研究. 海洋开发与管理. 2022(12): 98-106 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 3625
- HTML全文浏览量: 191
- PDF下载量: 1639
- 被引次数: 8
粤公网安备 44010502001741号
