Environmental impact mechanism of skipjack tuna fishery in Western and Central Pacific Ocean based on Multi-scale Geographical Weighted Regression Model (MGWR)
-
摘要:
鲣 (Katsuwonus pelamis) 是中西太平洋金枪鱼围网捕捞的重要资源，其资源分布受环境影响明显。为探索环境对鲣渔获率影响的空间异质性特征，利用中西太平洋渔业委员会 (Western and Central Pacific Fisheries Commission, WCPFC) 所公布的2005—2019年中西太平洋金枪鱼围网综合的1°×1°渔业及海洋环境数据，对标准化后的环境因子及渔获率选用多尺度地理加权回归 (Multi-scale Geographically Weighted Regression, MGWR) 方法进行研究。结果表明：1) 与传统广义加性模型 (Generalized Additive Model, GAM) 相比，考虑环境影响空间异质性问题的地理加权回归模型 (Geographically Weighted Regression, GWR) 和MGWR拟合优度 (R2) 有明显提升，校正后拟合优度 (Adjusted R2) 分别为0.273、0.846和0.871，且拟合结果的空间分布形态更符合真实情况。2) 各环境因子对鲣资源分布存在显著的空间非平稳性影响。各海洋环境因子对鲣渔获率分布影响的空间异质性程度 (各环境变量变异系数大小) 依次为水下55 m东西向海流速度 (Sea water X velocity at 55 m depth, U55)>海表面温度 (Sea surface temperature, SST)>净初级生产力 (Net primary productivity, NPP)>100 m盐度 (Sea water salinity at 100 m depth, S100)>55 m南北向海流速度 (Sea water Y velocity at 55 m depth, V55)。3) 各环境因子的影响存在明显尺度效应差异，NPP的作用尺度为44，其次为S100和U55 (均为48)，SST的为54，V55为全局尺度。4) 总体上，S100、NPP、SST、V55和U55对鲣渔获率正向影响比例依次为73.5%、64.8%、66.8%、80.8%和32.3%；其中S100、NPP和SST对鲣渔获率空间分布的影响相似，具体表现为东西向差异，170°E以西主要为正向影响，170°E以东为负向影响；U55为负向影响为主的因子。
-
关键词:
- 鲣 /
- 多尺度地理加权回归模型 /
- 空间异质性 /
- 中西太平洋
Abstract:Katsuwonus pelamis is an important resource for tuna purse seine fishing in the Western and Central Pacific Ocean, and its resource distribution is significantly affected by environment. In order to explore the characteristics of spatial heterogeneity of environmental impact on tuna catch rate, we used the 1°×1° fishery and marine environmental data of the Western and Central Pacific Ocean tuna purse-seine published by the Western and Central Pacific Fisheries Commission (WCPFC) from 2005 to 2019, and investigated the standardized environmental factors and catch rates by using Multi-scale Geographically Weighted Regression (MGWR) method. The results show that: 1) Compared with the traditional Generalized Additive Model (GAM), the Geographically Weighted Regression (GWR) and MGWR with spatial heterogeneity of environmental impacts improved the fit performance significantly. 2) Significant spatial non-stationarity was found for each environmental factor on the distribution of tuna resources. The degree of spatial heterogeneity (The magnitude of the coefficient of variation) of each environmental factor on the distribution of tuna catch rate followed a descending order of Sea water X velocity at 55 m depth (U55) > Sea surface temperature (SST) > Net primary productivity (NPP) >Sea water salinity at 100 m depth (S100) > Sea water Y velocity at 55 m depth (V55). 3) The effects of the environmental factors were found to have significant scale effects. 4) Overall, the positive effects of S100, NPP, SST, V55 and U55 on the catch rate of tuna were 73.5%, 64.8%, 66.8%, 80.8% and 32.3%, respectively. The effects of S100, NPP and SST on the spatial distribution of bonito catch rate were similar, specifically in terms of east-west differences, with positive effects mainly west of 170°E and negative effects east of 170°E. U55 was the main factor with negative effects.
-
长鳍光唇鱼 (Acrossocheilus longipinnis) 隶属鲤形目、鲤科、鲃亚科、光唇鱼属，俗称花榄鱼、五杠鲮、火烧鲮、阿肚鲮，为红水河珍稀特有鱼类，也是当地的重要经济物种。长鳍光唇鱼体长而侧扁，体色为银灰色，鳞片密集，侧线明显，体侧有5条淡黄色条带环绕，雄鱼背鳍末根分枝鳍条和第1根不分枝鳍条延长成丝状，形态独特，色彩美丽，在观赏鱼市场广受欢迎，被誉为红水河中的瑰宝[1]。然而，受梯级水电开发、采砂活动、过度捕捞、环境污染及外来鱼类入侵等因素的影响，近年来渔业资源调查显示长鳍光唇鱼野外种群资源下降明显[2-3]，已被列为易危种 (VU)[4]，开展人工驯养及繁育技术研究并实施放流以恢复其野外种群资源势在必行[5]。
鱼类早期发育研究对人工繁殖、苗种培育、资源保护具有重要指导意义，是开展人工繁育工作的基础[6]。一些学者对光唇鱼属的宽口光唇鱼 (A. monticola)[7]、光唇鱼 (A. fasciatus)[8-11]、云南光唇鱼 (A. yunanersis)[12-15] 等鱼类开展了人工繁殖及胚胎发育观察研究。国内外关于长鳍光唇鱼的研究极少，未见其人工繁育技术方面的相关报道。2023年4月，本研究团队成功完成了4次长鳍光唇鱼的人工繁育实验，并对其胚胎和仔鱼发育过程进行了系统观察，详细记录了各发育阶段的形态特征，为开展长鳍光唇鱼苗种规模化生产及增殖放流提供技术指导，对保护红水河鱼类种群资源及多样性具有重要意义。
1. 材料与方法
1.1 亲鱼培育
实验用长鳍光唇鱼亲鱼于2020年采集自红水河来宾市珍稀鱼类自然保护区，于来宾市红水河珍稀鱼类增殖保护站驯养培育。2023年3月1日，挑选3龄以上亲鱼250尾 (雌雄比约1∶1)，置于室内培育缸 (直径3 m，水深1.2 m) 开展强化培育，每缸50尾，溶解氧9.23 mg·L−1，pH 7.87，温度18~20 ℃。每天8:00和19:00各投喂1次拌合饲料 [硅藻和鲤 (Cyprinus carpio) 配合饲料]，每次拌合饲料投喂量为鱼总质量的2%。每天下午冲水1次，每次2 h。
1.2 催产和受精
挑选体质量大于0.3 kg的亲鱼用于人工繁殖。其中，雄鱼腹部松软，轻压有乳白色精液流出；雌鱼腹部明显隆起，生殖孔红肿。
催产剂由促黄体素释放激素A3 (LHRH-A3)、地欧酮 (DOM) 和绒毛膜促性腺激素 (HCG) 组成，采用2次注射法：第一针雌鱼注射LHRH-A3 2 μg·kg−1，雄鱼不注射；12 h后注射第二针，雌鱼LHRH-A3 15 μg·kg−1+DOM 5 mg·kg−1+HCG 1 000 IU·kg−1，雄鱼注射雌鱼50%剂量。
第二针催产剂注射24 h 后进行干法授精，采用人工挤压鱼体腹部的方式收集卵子和精液：亲鱼捞出后擦干鱼体表面水分，用毛巾固定好鱼体后，手从前向后挤压，将卵子挤出置于洗净擦干的塑料盆中，待所有雌鱼排卵完毕后，将雄鱼精液挤至卵子上，要求精液呈乳白色，无水、尿、血和粪便，入水后立即散开，显微镜下可见其快速运动。雌雄比约1∶1，加入适量9 g·L−1生理盐水稀释后用羽毛搅拌均匀，加水激活精子，静置2 min后冲洗数次，洗去黏液和杂质后将受精卵置于尤先科孵化器 (XYK-3260，水工程生态研究所监制) 中流水孵化，每平方米水面布置约1×104粒受精卵。孵化水温22~24 ℃，溶解氧质量浓度≥7 mg·L−1，每小时搅动受精卵12~15次，每天早晚用塑料吸管挑出死卵。
1.3 胚胎发育观察
对发育的受精卵及仔鱼进行持续观察和拍照，每次从孵化器中随机取30~50粒受精卵，在ZEISS Stemi 305体视显微镜 (北京普瑞赛司仪器有限公司) 下观察，当超过50%的受精卵发生形态变化时则判定为进入下一发育时期，并拍照记录和进行后续测量。原肠期前每次取受精卵间隔不超过10 min，之后每次取受精卵间隔不超过2 h。参考已有研究中的方法[6-7,11-12,16]，并根据显微镜下对长鳍光唇鱼胚胎发育的观察，对其胚胎发育时期进行划分。
1.4 数据采集和分析
人工繁殖效率、胚胎发育和仔稚鱼发育过程相关指标计算公式为
催产率=产卵鱼数/催产鱼总数×100%
受精率=原肠中期活卵数/产卵总数×100%
孵化率=鱼苗水花数/受精卵总数×100%
K= N×T
AGR=(Lt2−Lt1)/(t2−t1)
SGR=[(lnLt2−lnLt1)/(t2−t1)]×100%
式中：K为积温 (℃·h)；N为发育所需时间 (h)；T为平均水温 (℃)；AGR为绝对生长率 (mm·d−1)；SGR为特定生长率 (%·d−1)；Lt2和Lt1分别为发育到t2和t1时的全长[17]。
采用Excel 2021、Graphpad Prism 9.1和Adobe Photoshop CC 2023软件整理实验数据并统计作图，数据以“平均值±标准差 (
$\overline { x}\pm s $ )”表示。2. 结果
2.1 人工繁殖情况
2023年4月14—28日，共进行了4次人工繁殖实验 (表1) ，共催产雌鱼113尾，产卵78尾，催产率为 (68.95±13.26)%、受精率为 (68.50±19.82)%、孵化率为 (80.00±10.23)%，共获得初孵仔鱼3.80×104尾。
表 1 长鳍光唇鱼人工繁殖情况Table 1. Artificial propagation of A. longipinnis项目Item 2023-04-14 2023-04-17 2023-04-25 2023-04-28 雌雄亲鱼数量比
Male-to-female parental ratio9∶9 51∶44 30∶14 23∶15 催产率 Oxytocic rate/% 77.78 78.43 50.00 69.57 产卵量 Spawning capacity/104粒 1.30 3.20 0.90 1.10 受精率 Fertilization rate/% 70.00 72.00 42.00 90.00 孵化率 Hatchability/% 66.00 89.00 79.00 86.00 出苗量Seedling output/104尾 0.60 2.05 0.30 0.85 2.2 胚胎发育
长鳍光唇鱼受精卵为沉性卵、具弱黏性，呈金黄色 (图1-a)，卵径(2.42±0.06) mm，吸水膨胀后(3.03±0.08) mm，增大了约26%。长鳍光唇鱼胚胎发育过程历经胚盘形成期、卵裂期、囊胚期、原肠期、神经胚期、器官形成期和孵化期，共计7个阶段、28个时期，经66~74 h孵化出膜，详见图1和表2。
![]()
图 1 长鳍光唇鱼胚胎发育时序注：a. 受精卵；b. 胚盘期；c. 2细胞期；d. 4细胞期；e. 8细胞期；f. 16细胞期；g. 32细胞期；h. 64细胞期；i. 多细胞期；j. 囊胚早期；k. 囊胚中期；l. 囊胚晚期；m. 原肠早期；n. 原肠中期；o. 原肠晚期；p. 神经胚期；q. 胚孔封闭期；r. 肌节出现期；s. 眼基期；t. 眼囊期；u. 尾芽期；v. 尾泡期；w. 眼晶体形成期；x. 尾鳍出现；y. 肌肉效应期；z. 耳石出现期；aa. 心脏搏动期；ab. 孵化期；ac. 破膜期；ad. 初孵仔鱼。标尺=1 mm。Figure 1. Embryonic development sequence of A. longipinnisNote: a. Fertilized egg; b. Blastodem stage; c. 2-cell stage; d. 4-cell stage; e. 8-cell stage; f. 16-cell stage; g. 32-cell stage; h. 64-cell stage; i. Multicellular stage; j. Early-blastula stage; k. Mid-blastula stage; l. Late-blastula stage; m. Early-gastrula stage; n. Mid-gastrula stage; o. Late-gastrula stage; p. Neural stage; q. Closure of blastopore stage; r. Appearance of myomere; s. Optic rudiment stage; t. Eye vesicle formation period; u. Tail bud stage; v. Tail vesicle period; w. Crystal stage of eyes; x. Rudiment of tail fin; y. Muscular contraction period; z. Appearance of statolith; aa. Heart beating period; ab. Hatching stage; ac. Embryo hatching; ad. Newly hatched larvae. Bar=1 mm.表 2 长鳍光唇鱼胚胎发育时序Table 2. Timing of embryonic development sequence of A. longipinnis发育时期
Developmental period受精后时间
Time after fertilization持续时间
Duration积温
Accumulative temperature/(℃·h)图序
Plate受精卵 Fertilized egg 0 57 min 20.90 1-a 胚盘形成期 Blastodem stage 57 min 29 min 10.63 1-b 2细胞期 2-cell stage 1 h 26 min 36 min 13.20 1-c 4细胞期 4-cell stage 2 h 2 min 14 min 5.13 1-d 8细胞期 8-cell stage 2 h 16 min 12 min 4.40 1-e 16细胞期 16-cell stage 2 h 28 min 15 min 5.50 1-f 32细胞期 32-cell stage 2 h 43 min 28 min 10.27 1-g 64细胞期 64-cell stage 3 h 11 min 30 min 11.50 1-h 多细胞期 Multicellular stage 3 h 41 min 1 h 28 min 33.73 1-i 囊胚早期 Early-blastula stage 5 h 9 min 1 h 5 min 26.00 1-j 囊胚中期 Mid-blastula stage 6 h 14 min 53 min 21.20 1-k 囊胚晚期 Late-blastula stage 7 h 7 min 1 h 27 min 34.80 1-l 原肠早期 Early-gastrula stage 8 h 34 min 2 h 37 min 62.80 1-m 原肠中期 Mid-gastrula stage 11 h 11 min 1 h 59 min 45.62 1-n 原肠晚期 Late-gastrula stage 13 h 10 min 1 h 51 min 42.55 1-o 神经胚期 Neural stage 15 h 1 min 1 h 36 min 36.80 1-p 胚孔封闭期 Closure of blastopore stage 16 h 37 min 5 h 6 min 112.20 1-q 肌节出现期 Appearance of myomere 21 h 43 min 6 h 54 min 151.80 1-r 眼基期 Optic rudiment stage 28 h 35 min 1 h 19 min 28.97 1-s 眼囊期 Eye vesicle formation period 29 h 54 min 5 h 40 min 124.67 1-t 尾芽期 Tail bud stage 35 h 34 min 31 min 11.37 1-u 尾泡期 Tail vesicle period 36 h 5 min 1 h 46 min 38.87 1-v 眼晶体形成期 Crystal stage of eyes 37 h 51 min 51 min 18.70 1-w 尾鳍出现 Rudiment of tail fin 38 h 42 min 1 h 16 min 27.87 1-x 肌肉效应期 Muscular contraction period 39 h 58 min 6 h 4 min 139.53 1-y 耳石形成期 Appearance of statolith 46 h 2 min 7 h 168.00 1-z 心跳期 Heart beating period 53 h 2 min 7 h 23 min 169.82 1-aa 孵化期 Hatching stage 60 h 25 min 5 h 39 min 124.30 1-ab 破膜期 Embryo hatching 66 h 4 min 3~5 min 1-ac 初孵仔鱼 Newly hatched larvae 66 h 8 min 1-ad 2.2.1 胚盘形成期
受精后，长鳍光唇鱼受精卵内原生质向动物极汇聚，逐渐隆起形成胚盘，在受精57 min后胚盘隆起最高，约达卵黄高度的1/4，在黄色卵黄顶端形成1个淡乳白色球状结构 (图1-b)。
2.2.2 卵裂期
长鳍光唇鱼受精卵属于盘状卵裂，胚盘部分进行分裂而卵黄部分不参与分裂。在受精1 h 26 min后，受精卵进行了第1次分裂，胚盘中央出现1条分裂沟，由胚盘处的1个细胞分裂成2个同样大小的细胞，到达2细胞期 (图1-c)；受精卵卵裂速度较快，在受精2 h 2 min后，经历第2次卵裂，出现1条与第1次分裂相垂直的分裂沟，将2细胞分裂成为4细胞，到达4细胞期 (图1-d)；在进入4细胞期14 min后，出现2条与初次分裂平行的分裂沟，将4细胞分裂成8细胞，到达8细胞期 (图1-e)；又经过12 min，受精卵经历第4次卵裂，到达16细胞期 (图1-f)；再经过15 min，受精卵经历第5次卵裂，到达32细胞期，细胞大小不再均一，排列也变得不规则 (图1-g)；再经过28 min，受精卵经历第6次分裂，形成细胞密集排列的细胞团，到达64细胞期 (图1-h)；在受精3 h 41 min后，到达多细胞期，随着细胞团不断分裂，细胞体积逐渐变小，开始重叠排列，但细胞边缘仍清晰可见 (图1-i)。
2.2.3 囊胚期
在受精5 h 9 min后，随着细胞团分裂，细胞间界限不清晰，边缘逐渐平滑，形成中部向上隆起的帽状半透明细胞团，到达囊胚早期 (图1-j)；囊胚早期持续1 h 5 min后，囊胚层高度逐渐降低，颜色逐渐加深，细胞不断下包，向植物极延伸，但胚层与卵黄之间的过渡部分凹痕仍然明显，到达囊胚中期 (图1-k)；又经过53 min，囊胚层细胞不断下包扩散，形成紧贴卵黄边缘平滑的囊胚层，进入囊胚晚期 (图1-l)。
2.2.4 原肠期
在受精8 h 34 min后，胚层细胞下包至卵黄1/3处，在卵黄顶部形成帽状“保护套”，出现了清晰可见的胚环，进入原肠早期 (图1-m)；又经过2 h 37 min，胚层已包裹至卵黄的一半，此时胚环背侧开始出现增厚并形成胚盾，标志着进入原肠中期 (图1-n)；再经过1 h 59 min，胚层继续下包至卵黄的4/5，进入原肠晚期 (图1-o)。
2.2.5 神经胚期
在受精15 h 1 min后，胚层一侧从顶端开始增厚隆起形成脑泡，出现了胚胎神经板雏形，胚盾增长变细形成胚孔，露出卵黄栓 (图1-p)；再经过1 h 36 min，胚孔完全闭合，形成弧形胚体，腹面包绕卵黄囊约3/5，到达胚孔封闭期 (图1-q)。
2.2.6 器官形成期
在受精21 h 43 min后，胚体发育形成6对肌节，卵黄囊被包裹约3/4，到达肌节出现期 (图1-r)；又经过6 h 54 min，胚体头部两侧出现椭圆状凸起泡，形成眼原基，到达眼基期 (图1-s)；随后经过31 min，眼原基中部下陷形成眼囊，包绕卵黄囊约4/5，进入眼囊期 (图1-t)；在受精35 h 34 min后，胚体尾部开始与卵黄囊分离，进入尾芽期 (图1-u)。又经过31 min，胚体尾部凸起形成尾泡，进入尾泡期，肌节增至16~18对 (图1-v)。受精37 h 51 min后，眼晶体形成 (图1-w)；又经过51 min后，尾部前方卵黄囊下陷，尾部完全脱离卵黄囊，可见尾鳍褶 (图1-x)；随着尾鳍和肌肉的发育，胚体开始扭动，可见心脏原基，胚胎卵黄囊内凹逐渐加剧 (图1-y)；在受精46 h 2 min后，胚体耳囊中形成清晰可见的耳石 (图1-z)，逐渐脱离卵黄囊，扭动加剧；又经过7 h后，胚体心脏开始跳动 (图1-aa)，达20~30次·min−1，胚体扭动频繁。
2.2.7 孵化期
在受精60 h 25 min后，长鳍光唇鱼受精卵进入孵化期，胚体长度明显增长，尾部环绕卵黄囊可延伸至眼囊与耳囊之间，胚体剧烈扭动翻滚 (图1-ab)，孵化期可持续约6 h；然后进入破膜期，一般胚体尾部先破膜露出，尾部扭动使头部脱膜 (图1-ac)，也有少部分胚胎先露出头部，出膜过程一般持续3~5 min；仔鱼出膜后即舒展身体，尾部微微上翘 (图1-ad)。
2.3 仔鱼发育
在水温22~24 ℃下，长鳍光唇鱼仔鱼发育经历21 d，21日龄 (Days after hatching, DAH) 进入稚鱼期。
2.3.1 卵黄囊仔鱼
初孵仔鱼体长 (8.33±0.23) mm，肌节34对，卵黄囊由前后两部分组成：前半部分卵黄囊呈圆球状，长径 (2.33±0.16) mm，短径 (1.67±0.13) mm；后半部分卵黄囊呈短棒状，长径 (3.09±0.16) mm，短径 (0.66±0.11) mm。仔鱼除头前部和尾部外，其余部分紧贴整个卵黄囊，卵黄囊约占体长的2/3。仔鱼胸鳍、背鳍、腹鳍以及肛门原基形成，尾椎微微上翘 (图1-ad)。
仔鱼出膜4 h后，卵黄囊逐渐被吸收，前半部分逐渐由圆球状变为橄榄球状，后半部分也延伸变细，开始出现血液流动 (图2-a)。
![]()
图 2 长鳍光唇鱼仔鱼发育注：a. 孵后4 h仔鱼；b. 孵后16 h仔鱼；c. 孵后24 h仔鱼；d. 孵后36 h仔鱼；e. 孵后48 h仔鱼；f. 3 DAH仔鱼； g. 4 DAH仔鱼； h. 5 DAH仔鱼；i. 6 DAH仔鱼；j. 7 DAH仔鱼；k. 8 DAH仔鱼；l. 9 DAH仔鱼；m. 10 DAH仔鱼； n. 11 DAH仔鱼；o. 12 DAH仔鱼； p. 13 DAH仔鱼；q. 14 DAH仔鱼；r. 15 DAH仔鱼；s. 21 DAH仔鱼。标尺=1 mm。Figure 2. Development of A. longipinnis larvaNote: a. Larv of 4 h after hatching; b. Larv of 16 h after hatching; c. Larv of 24 h after hatching; d. Larv of 36 h after hatching; e. Larv of 48 h after hatching; f. 3 DAH larv; g. 4 DAH larv; h. 5 DAH larv; i. 6 DAH larv; j. 7 DAH larv; k. 8 DAH larv; l. 9 DAH larv; m. 10 DAH larv; n. 11 DAH larv; o. 12 DAH larv; p. 13 DAH larv; q. 14 DAH larv; r. 15 DAH larv; s. 21 DAH larv. Bar=1 mm.出膜16 h后，仔鱼眼囊中开始出现黑色素沉积，各鳍原基进一步发育，鳍褶明显 (图2-b)。
出膜24 h后，鳔的雏形出现，眼囊中黑色素向中间汇集形成色素膜，眼睛发光灵动，在仔鱼腹部与卵黄囊之间开始出现散点状黑色素沉积 (图2-c)。
出膜36 h后，卵黄囊进一步被吸收，前后两部分分界不明显，呈流线型。在卵黄囊前方出现明显的心脏，血流明显，心跳速度80~85次·min−1。腹部黑色素沉积变多，头部出现黑色素沉积，透明度降低 (图2-d)。
出膜48 h后，仔鱼背部拱起，头部黑色素沉积成条带状，血管脉络增多，尾鳍褶开始分化 (图2-e)。
3 DAH仔鱼尾柄处和躯干上开始出现黑色素沉积，卵黄囊侧边出现了明显的黑色素沉积，表面缠绕了明显的血管组织，卵黄囊的吸收开始加快 (图2-f)。
4 DAH仔鱼黑色素进一步沉积，形成了5条横斑的雏形，肠道成形 (图2-g)。
5 DAH仔鱼出现了明显的鳔，增加了躯体和卵黄囊的间隔，卵黄囊进一步被吸收，肛门连通，多数可以平游 (图2-h)。
6 DAH仔鱼开始进食饵料，肠道中可见红色丰年虫，由内源性营养期到达混合营养期 (图2-i)。
7 DAH仔鱼卵黄囊被大幅度吸收，肠道内充满食物，鳔分化为前后二室，黑色素沉积明显增多 (图2-j)。
8 DAH仔鱼背鳍基本成形，尾鳍开始分叉，卵黄囊仅少量残余 (图2-k)。
2.3.2 晚期仔鱼
9 DAH仔鱼卵黄囊完全吸收，消化道充满食物，由混合营养期到达外源性营养期 (图2-l)。
10 DAH仔鱼胸鳍基本成形，背鳍和尾鳍鳍条分明，但腹鳍褶依然存在，尚未分化成为腹鳍 (图2-m)。
11~13 DAH仔鱼黑色素沉积逐渐增加，基本形成了以后的身体色彩分区，腹鳍褶向腹鳍逐步分化，臀鳍逐渐形成 (图2-n—2-p)。
14 DAH仔鱼腹鳍褶消失，分化为腹鳍 (图2-q)。
15 DAH仔鱼腹鳍鳍条明显，各鳍发育基本完毕 (图2-r)。
21 DAH仔鱼背鳍和尾鳍之间开始形成鳞片，进入稚鱼期 (图2-s)。
2.4 仔鱼生长
长鳍光唇鱼仔鱼全长和日龄的关系见图3。在2 DAH前仔鱼全长增长较快，而2 DAH后的卵黄囊仔鱼全长增长速度减慢，晚期仔鱼全长增长略有加快。各发育阶段的全长绝对生长率和全长特定生长率见表3，其中2 DAH前仔鱼阶段全长绝对生长率和全长特定生长率均最高，分别达1.10 mm·d−1和12.05%·d−1。
表 3 长鳍光唇鱼仔鱼各发育阶段的生长率Table 3. Growth rates of A. longipinnis larvae at different developmental stages发育阶段
Developmental stage经历时间
Experience time/d全长绝对生长率
AGR/(mm·d−1)全长特定生长率
SGR/(%·d−1)2 DAH前仔鱼 Prelarvae of 2 DAH 2 1.10 12.05 2 DAH后仔鱼 Larvae after 2 DAH 19 0.21 1.72 2~8 DAH仔鱼 Larvae of 2−8 DAH 6 0.16 1.49 卵黄囊期仔鱼 Yolk sac stage larvae 8 0.40 4.13 晚期仔鱼 Late larvae 13 0.23 1.82 仔鱼期 Larval stage 21 0.29 2.70 ![]()
图 3 长鳍光唇鱼仔鱼日龄和全长的关系Figure 3. Relationship between total length and days after hatching of A. longipinnis larvae3. 讨论
3.1 长鳍光唇鱼人工繁殖技术
本研究选用的亲本为经过驯养的3龄以上的长鳍光唇鱼，体质量均在0.3 kg以上，先进行强化培育再进行人工催产和受精，经过4次人工繁殖，平均孵化率约为80%。本次长鳍光唇鱼人工繁殖技术取得突破的关键在于：1) 亲本经过人工驯养，耐应激程度较高；2) 强化培育期间以硅藻配合鲤鱼饲料进行投喂，每天冲水刺激性腺发育，以及定期拉网锻炼；3) 合适的催产剂配方；4) 适宜的孵化条件 (水温、流速)。在软鳍新光唇鱼 (Neolissochilus benasi) 的研究中也提出了类似的观点：降低亲鱼应激程度、冲水刺激、驯养亲鱼和强化培育等措施都可以提高人工繁殖的成功率[18]。尤其长鳍光唇鱼具有应激易死和硅藻食性的特点，应该在人工繁殖中着重把握。在催产剂的配方上，本研究采用二次注射法，先对雌鱼注射1次LHRH-A3，然后采取LHRH-A3+DOM+HCG的组合，这与华泽祥等[12]在云南光唇鱼上所用催产剂配方 (LHRH-A2+DOM+HCG) 类似，可能是因为这种配方对光唇鱼属具有比较好的催产效果。
3.2 长鳍光唇鱼胚胎发育和仔鱼发育特点
鱼类胚胎发育是一个十分关键且脆弱的环节，受多种因素 (温度、光照、溶解氧、pH) 影响[19]，不同鱼类胚胎发育各具特点，这些特点也与各自生活习性息息相关。长鳍光唇鱼成熟卵呈金黄色，与光唇鱼属中的宽口光唇鱼[7]、光唇鱼[8-11]、云南光唇鱼[12-15]等颜色相同。长鳍光唇鱼受精卵为沉性卵、具弱黏性，卵径为 (2.42±0.06) mm，吸水膨胀后为(3.03±0.08) mm，卵周隙变大，受精卵密度变小，有利于随流水流动，是溪流性鱼类的特点之一[20]。长鳍光唇鱼受精卵卵膜较厚，呈淡黄色、半透明，富有弹性，有利于自然繁殖情况下保护其在河流底部石砾间滚动[21]。
长鳍光唇鱼受精卵胚胎发育过程历经胚盘形成期、卵裂期、囊胚期、原肠期、神经胚期、器官形成期和孵化期，共计7个阶段、28个时期，符合硬骨鱼类胚胎发育的基本过程[22]。孵化用时66~74 h，久于光唇鱼 (46 h) 和宽口光唇鱼 (49 h)，快于云南光唇鱼 (103 h) 和软鳍新光唇鱼 (120 h)。较长时间的胚胎发育期往往需要较长的流水江段，这与长鳍光唇鱼的卵膜特性也相对应[23]。在尾泡、眼晶体、耳石等器官的分化顺序上与其他同属鱼类基本一致，除光唇鱼外[11]。
仔鱼阶段是从孵化出膜到体表有鳞片生成，根据卵黄囊有无又分为卵黄囊仔鱼和晚期仔鱼[24]。鱼类成熟卵的卵径直接决定了卵黄囊的大小，长鳍光唇鱼成熟卵卵径为 (2.42±0.06) mm，与光唇鱼大小相当，大于同属的云南光唇鱼[12]、宽口光唇鱼[7]，以及同亚科的黑脊倒刺鲃 (Spinibarbus caldwelli)[25]、中华倒刺鲃 (S. sinensis)[26]、大鳞鲃 (Barbus capito)[27-28]等。而卵黄囊的大小与其被完全吸收，以及仔鱼开口摄食的时间呈正相关[29]。仔鱼完全以卵黄囊作为营养来源的阶段称为内源性营养阶段，在这期间仔鱼主要完成口、眼、消化道、血液循环系统和鳍的初步发育，建立了巡游模式，基本实现平游，具备了从内源性营养阶段转向外源性营养阶段的条件，可以开口摄食，但此时卵黄囊往往并没有被完全吸收，故而此阶段称为混合营养期[30]。不同鱼类内源性营养期时间不同，长鳍光唇鱼内源性营养期5 d，混合营养期3 d，卵黄囊仔鱼期共持续8 d，大于同属的光唇鱼[8-11]、云南光唇鱼[12-15]，和同亚科的黑脊倒刺鲃[25]、中华倒刺鲃[26]、大鳞鲃[27-28]等，以及黑尾近红鲌 (Ancherythroculter nigrocauda)[30]、岩原鲤 (Procypris rabaudi)[30]和太湖翘嘴红鲌 (Erythroculter ilishaeformis Bleeker)[31]等鲤科鱼类。比较不同鱼类之间的发育轨迹是解释物种早期发育阶段丰度、分布和种群生存能力差异的第一步[32]，较长的卵黄囊仔鱼期和混合营养期有利于长鳍光唇鱼建立初次摄食，能更好地适应环境。
3.3 长鳍光唇鱼仔鱼生长速率
长鳍光唇鱼仔鱼期全长增长速率并非完全和日龄呈正相关，与光唇鱼不同[8-11]。2 DAH前仔鱼阶段全长绝对生长率和全长特定生长率均最高，卵黄囊也被大幅吸收，由前圆球状后柱状的构型变为流线型。证明长鳍光唇鱼初孵48 h内，卵黄营养物质被大量吸收用来加速身体生长。而2 DAH后的卵黄囊仔鱼全长增长速度放缓，卵黄囊吸收速度下降，此时卵黄囊被用来支持消化系统发育，为进入外源性营养阶段做准备。卵黄囊吸收完毕后，进入外源性营养期，晚期仔鱼大量摄入饵料物质，全长增长再次加快。这样的仔鱼全长生长速率的变化趋势与扁吻鱼 (Aspiorhynchus laticeps)[16]、短须裂腹鱼 (Schizothorax wangchiachii)[33]基本一致。
4. 结论
本研究对3龄以上长鳍光唇鱼进行强化培育和人工催产，成功孵化3.80×104尾鱼苗，平均孵化率约为80%。长鳍光唇鱼受精卵呈金黄色，具弱黏性，卵膜淡黄色半透明具有弹性，胚胎发育过程历经胚盘形成期、卵裂期、囊胚期、原肠期、神经胚期、器官形成期和孵化期，共计7个阶段、28个时期，在水温22~24 ℃下经66~74 h孵化出膜，与同属鱼类发育特征相似。初孵仔鱼的卵黄囊由前后两部分组成，仔鱼经内源性营养期5 d，混合营养期3 d，卵黄囊仔鱼期共持续8 d。2 DAH前仔鱼阶段生长最快，随后全长增长速度放缓，卵黄囊吸收速度下降，仔鱼期共历时21 d，全长特定生长率为2.70%·d−1。
-
![]()
图 1 MGWR 模型各环境因子局部回归系数和显著性的空间分布
Figure 1. Spatial distribution of local regression coefficients and significance of each environmental factor in MGWR model
![]()
图 2 不同模型渔获潜能标准化预测及 2005—2019 年平均 CPUE 空间分布
Figure 2. Standardized catch potential prediction and average CPUE spatial distribution of different models from 2005 to 2019
表 1 解释变量间方差膨胀因子
Table 1 Variance inflation factor among explanatory variables
解释变量
Explanatory variable单位
Unit方差膨胀因子
VIF海表面温度
Sea surface temperature, SSTK 1.558 100 m 水深盐度
Sea water salinity at 100 m depth, S100‰ 1.090 净初级生产力
Net primary productivity, NPPmol·(m2·s)−1 1.858 55 m 东西向海水流速
Mean zonal (E-W) current velocity at
55 m depth, U55m·s−1 1.007 55 m 南北向海水流速
Mean meridional (N-S) current
velocity at 55 m depth, V55m·s−1 1.796 
下载: 导出CSV
表 2 GAM、GWR 和 MGWR 不同回归模型性能评价对比
Table 2 Comparison of statistical parameters of different linear regression models (GAM, GWR and MGWR)
参数
Parameter广义加性
模型
GAM地理加权
回归
GWR多尺度地理
加权回归
MGWR残差平方和 RSS 254 400.400 148.607 129.136 修正的赤池信息准则 AICc 1 477.598 1 572.437 1 287.873 拟合优度 R² 0.369 0.879 0.895 校正后的拟合优度
Adjusted R²0.273 0.846 0.871
下载: 导出CSV
表 3 MGWR 与 GWR 带宽对比结果
Table 3 Bandwidth comparison between classical MGWR and GWR models
变量
Variable多尺度地理
加权回归
MGWR占比
Proportion/
%地理加权
回归
GWR占比
Proportion/
%截距 Intercept 48 3.91 65 5.29 海表面温度 SST 54 4.40 65 5.29 100 m 水深盐度 S100 48 3.91 65 5.29 55 m 东西向海水
流速 U5548 3.91 65 5.29 55 m 南北向海水
流速 V551 227 99.84 65 5.29 净初级生产力 NPP 44 3.51 65 5.29
下载: 导出CSV
表 4 基于MGWR的各环境因子的局部系数统计描述
Table 4 Statistical description of MGWR local coefficient
变量
Variable均值
Mean标准差
SD最小值
Min.最大值
Max.正值比例
Positive ratio/%负值比例
Negative ratio/%显著性检验
P截距 Intercept −0.036 0.350 −0.552 0.939 38.975 61.025 1.00 海表面温度 SST 0.162 0.331 −0.555 1.551 66.802 33.198 <0.05 100 m 水深盐度 S100 0.574 0.952 −0.359 4.649 73.556 26.444 <0.05 55 m 东西向海水流速 U55 −0.180 0.504 −2.215 0.997 32.303 67.697 <0.05 55 m 南北向海水流速 V55 0.069 0.109 −0.142 0.553 80.797 19.203 <0.05 净初级生产力 NPP 0.172 0.351 −0.787 1.806 64.849 35.151 <0.05
下载: 导出CSV
-
[1] 陈新军, 郑波. 中西太平洋金枪鱼围网渔业鲣鱼资源的时空分布[J]. 海洋学研究, 2007(2): 13-22. doi: 10.3969/j.issn.1001-909X.2007.02.002 [2] WILLIAMS P, RUAIA T. Overview of tuna fisheries in the western and central Pacific Ocean, including economic conditions–2019[C]//16th WCPFC Scientific Committee Meeting (SC16). Pohnpei State, Federated States of Micronesia: WCPFC, 2020: 7.
[3] DUERI S, BOPP L, MAURY O. Projecting the impacts of climate change on skipjack tuna abundance and spatial distribution[J]. Global Change Biol, 2014, 20(3): 742-753. doi: 10.1111/gcb.12460
[4] 杨晓明, 戴小杰, 田思泉, 等. 中西太平洋金枪鱼围网渔场变动的预报分析[J]. 中国水产科学, 2016, 23(6): 1417-1425. [5] 王学昉, 许柳雄, 周成, 等. 中西太平洋金枪鱼围网鲣鱼自由鱼群捕获成功率与温跃层特性的关系[J]. 上海海洋大学学报, 2013, 22(5): 763-769. [6] 魏广恩. 北太平洋柔鱼渔场的时空分析与资源丰度的预测[D]. 上海: 上海海洋大学, 2018: 77 [7] 陈洋洋, 陈新军, 郭立新, 等. 基于BP神经网络的中西太平洋鲣鱼渔场预报模型构建与比较[J]. 广东海洋大学学报, 2017, 37(6): 65-73. [8] ROBINSON M, SEI-ICHI S. Ensemble modelling of skipjack tuna (Katsuwonus pelamis) habitats in the Western North Pacific using satellite remotely sensed data: a comparative analysis using Machine-Learning Models[J]. Remote Sens-Basel, 2020, 12(16): 2591. doi: 10.3390/rs12162591
[9] VENABLES W N, DICHMONT C M. GLMs, GAMs and GLMMs: an overview of theory for applications in fisheries research[J]. Fish Res, 2004, 70(2/3): 319-337.
[10] 杨彩莉, 杨晓明, 朱江峰. 不同类型厄尔尼诺事件中环境因子对中西太平洋金枪鱼围网鲣分布响应[J]. 南方水产科学, 2021, 17(3): 8-18. [11] DUQUE-LAZO J, GILS H V, GROEN T A, et al. Transfer ability of species distribution models: the case of Phytophthora cinnamomi in Southwest Spain and Southwest Australia[J]. Ecol Model, 2016, 320: 62-70. doi: 10.1016/j.ecolmodel.2015.09.019
[12] 姜珊, 杨晓明, 朱江峰. 中西太平洋金枪鱼围网鲣渔获率与海洋环境关系的时空异质性[J]. 中国水产科学, 2022, 29(5): 744-754. [13] FENG Y J, LIU Y, CHEN X J. Modeling monthly spatial distribution of Ommastrephes bartramii CPUE in the northwest Pacific and its spatially nonstationary relationships with the marine environment[J]. J Ocean U China, 2018, 17(3): 647-658. doi: 10.1007/s11802-018-3495-9
[14] 沈体雁, 于瀚辰, 周麟, 等. 北京市二手住宅价格影响机制: 基于多尺度地理加权回归模型(MGWR)的研究[J]. 经济地理, 2020, 40(3): 75-83. [15] FOTHERINGHAM A S, YANG W, KANG W. Multiscale geographically weighted regression (MGWR)[J]. Ann Am Assoc Geogr, 2017, 107(6): 1247-1265.
[16] YU H, FOTHERINGHAM A S , LI Z, et al. On the measurement of bias in geographically weighted regression models[J]. Spat Stat-Neth, 2020: 100453.
[17] 李斌, 邢汉发, 刘烨菲. 基于MGWR模型的城市景观与热环境关系分析[J]. 西安理工大学学报, 2021, 37(4): 517-525. [18] 周天军, 邹立维, 陈晓龙. 第六次国际耦合模式比较计划(CMIP6)评述[J]. 气候变化研究进展, 2019, 15(5): 445-456. [19] BUI D T, LOFMAN O, REVHAUG I, et al. Landslide susceptibility analysis in the Hoa Binh province of Vietnam using statistical index and logistic regression[J]. Nat Hazards, 2011, 59(3): 1413-1444. doi: 10.1007/s11069-011-9844-2
[20] HASTIE T J, TIBSHIRANI R J. Generalized Additive Models[M]. London: Chapman and Hall, 1990: 17.
[21] TSENG C T, SUN C L, YEH S Z, et al. Spatio-temporal distributions of tuna species and potential habitats in the Western and Central Pacific Ocean derived from multisatellite data[J]. Int J Remote Sens, 2010, 31(17/18): 4543-4558.
[22] SONG X L, MI N, MI W B, et al. Spatial non-stationary characteristics between grass yield and its influencing factors in the Ningxia temperate grasslands based on a mixed geographically weighted regression model[J]. J Geogr Sci, 2022, 32(6): 1076-1102. doi: 10.1007/s11442-022-1986-5
[23] 祝新明, 宋小宁, 冷佩, 等. 多尺度地理加权回归的地表温度降尺度研究[J]. 遥感学报, 2021, 25(8): 1749-1766. [24] JONES S, SILAS E G. Synopsis of biological data on skipjack Katsuwonus pelamis (Linnaeus) 1758 (Indian Ocean): FAO Fisheries Biology Synopsis No. 64 Species Synopsis No. 21[R]. Rome: FAO, 1963: 663-694.
[25] DRUON J N, CHASSOT E, MURUA H, et al. Skipjack tuna availability for purse seine fisheries is driven by suitable feeding habitat dynamics in the Atlantic and Indian Oceans[J]. Front Mar Sci, 2017, 4: 315. doi: 10.3389/fmars.2017.00315
[26] 周甦芳, 沈建华, 樊伟. ENSO现象对中西太平洋鲣鱼围网渔场的影响分析[J]. 海洋渔业, 2004(3): 167-172. [27] LEHODEY P M, BERTIGNAC M, HAMPTON J, et al. El Niño Southern Oscillation and tuna in the western Pacific[J]. Nature 1997, 389: 715-718.
[28] MEEHL G A. Characteristics of surface current flow inferred from a Global Ocean Current Data Set[J]. J Phys Oceanogr, 1982, 12(6): 538-555. doi: 10.1175/1520-0485(1982)012<0538:COSCFI>2.0.CO;2
[29] 郭爱, 余为, 陈新军, 等. 中国近海鲐鱼资源时空分布与海洋净初级生产力的关系研究[J]. 海洋学报, 2018, 40(8): 42-52. [30] 黄易德. 中西太平洋正鲣资源时空分布特性之研究[D]. 基隆: 台湾海洋大学, 2002: 83. [31] SUGIMOTO T, KIMURAA S, TADOKOROB K. Impact of El Niño eventsand climate regime shift on living resources in the western North Pacific[J]. Prog Oceanogr, 2001, 49: 113-127. doi: 10.1016/S0079-6611(01)00018-0
[32] 杨晓明, 戴小杰, 田思泉, 等. 中西太平洋鲣鱼围网渔业资源的热点分析和空间异质性[J]. 生态学报, 2014, 34(13): 3771-3778. [33] 杨晓明, 戴小杰, 王学昉, 等. 基于点格局的中西太平洋金枪鱼围网中鲣的空间格局特征[J]. 中国水产科学, 2017, 24(3): 633-639. [34] 杨胜龙, 周甦芳, 周为峰, 等. 基于Argo数据的中西太平洋鲣渔获量与水温、表层盐度关系的初步研究[J]. 大连水产学院学报, 2010, 25(1): 34-40. [35] 叶泰豪, 冯波, 颜云榕, 等. 中西太平洋鲣渔场与温盐垂直结构关系的研究[J]. 海洋湖沼通报, 2012(1): 49-55. [36] PICAUT J, IOUALALEN M, MENKES C, et al. Mechanism of the zonal displacements of the Pacific warm pool: implications for ENSO[J]. Science, 1996, 274(5292): 1486-1489. doi: 10.1126/science.274.5292.1486
[37] ZHANG F, ZHANG R H. Interannually varying salinity effects on ENSO in the tropical Pacific: a diagnostic analysis from Argo[J]. Ocean Dyn, 2015, 65(5): 691-705. doi: 10.1007/s10236-015-0829-7
[38] 唐未, 王学昉, 吴峰, 等. 基于最大熵模型模拟西印度洋剑鱼栖息地的时空分布[J]. 海洋学报, 2022, 44(10): 100-108. [39] 张小龙, 付东洋, 刘大召, 等. 基于EOF分析中西太平洋金枪鱼围网渔场的海洋环境[J]. 海洋学研究, 2019, 37(2): 81-94. [40] 唐浩, 许柳雄, 陈新军, 等. 基于GAM模型研究时空及环境因子对中西太平洋鲣鱼渔场的影响[J]. 海洋环境科学, 2013, 32(4): 518-522. -
期刊类型引用(10)
1. 朱文慧，谭桂芝，步营，李学鹏，励建荣，崔方超，檀茜倩，孟玉琼，马睿，郭晓华. 低温等离子体耦合微酸性电解水对三文鱼的保鲜作用. 中国食品学报. 2024(02): 228-238 . 
百度学术
2. 孙万青，岑剑伟，陈胜军，李春生，邓建朝，潘创，杨少玲，冯阳. 微酸性电解水对盐渍海蜇脱铝的效果研究. 中国食品学报. 2023(01): 160-170 . 百度学术
3. 韦丽娜，李来好，郝淑贤，黄卉，杨贤庆，相欢，赵永强，岑剑伟，魏涯. 渗透处理对冷冻干燥罗非鱼肉品质和肌原纤维蛋白的影响. 南方水产科学. 2023(02): 133-141 . 本站查看
4. 敖菲菲，方祥，梁钻好，陈海强，梁凤雪，黄成龙，余铭. 酸性电解水对罗非鱼片的杀菌工艺研究. 食品科技. 2021(03): 105-111 . 百度学术
5. 张宇，宋力，李达，王志敏，田瑞. 太阳能与辅助热源互补供暖仿真研究. 工程热物理学报. 2021(08): 2048-2059 . 百度学术
6. 孔金花，诸永志，葛庆丰，卞欢，闫征，刘芳，徐为民. 超声波协同微酸性电解水对小龙虾净化及品质的影响. 食品工业科技. 2021(21): 182-189 . 百度学术
7. 蓝蔚青，刘琳，孙晓红，赵亚楠，谢晶. 酸性电解水发生机理及在水产领域中的应用研究进展. 食品与发酵工业. 2020(04): 294-298 . 百度学术
8. 岑剑伟，于福田，杨贤庆，李来好，黄卉，赵永强，魏涯，郝淑贤，林织. 几种杀菌方法对罗非鱼片鱼肉品质的影响比较. 食品与发酵工业. 2020(07): 123-129 . 百度学术
9. 岑剑伟，于福田，杨贤庆，李来好，黄卉，魏涯，赵永强，林织. 微酸性电解水对罗非鱼片保鲜效果的研究. 食品与发酵工业. 2019(18): 209-214 . 百度学术
10. 石慧，杨少玲，吴燕燕，林婉玲，杨贤庆，黄卉. 卵形鲳鲹鱼片热风干燥条件优化及其品质特性研究. 食品与发酵工业. 2019(17): 129-135 . 百度学术
其他类型引用(7)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 17
粤公网安备 44010502001741号
