Analysis of ontogenetic phenotype changes of Enoploteuthis chunii beak
-
摘要: 角质颚作为头足类重要的硬组织，记录了其丰富的个体生长信息。在西北太平洋海域，富山武装乌贼 (Enoploteuthis chunii) 是众多海洋哺乳动物和鱼类的食饵，掌握角质颚的生长变化规律不仅有助于了解其个体发育期摄食生态的变化，还有助于进一步研究其在海洋生态系统中的地位和作用。根据2019年3月和9月于西太平洋海域采集的91尾样本，结合胴长数据从角质颚长度及色素沉积分析了富山武装乌贼个体发育期角质颚的表型变化规律。结果表明，随着胴长增加，角质颚各部长度与色素沉积比例均呈正相关，角质颚各部长度在胴长31 和41 mm处跃变，上颚脊突、侧壁及下颚侧壁生长较快；上颚头盖、下颚头盖和下颚侧壁在胴长31 mm处色素沉积发生跃变，上颚侧壁在胴长41 mm处色素沉积跃变明显。ANOVA分析发现，富山武装乌贼在不同个体发育期，角质颚表型变化差异显著 (P<0.05)。研究认为，通过对角质颚生长过程中形态大小、色素沉积等表型变化进行分析，可为富山武装乌贼个体发育期摄食转变的研究提供基础。Abstract: Cephalopod beak is an important hard tissue that records lots of growth information. In the Northwest Pacific Ocean, Enoploteuthis chunii is the prey of many marine mammals and fishes. Understanding phenotype changes in the growth of beak helps to know E. chunii ontogenetic trophic ecology variation and to investigate its ecological role in the system. In this study, the ontogenetic changes in beak phenotype were investigated by analyzing the variation of beak length and pigmentation of 91 squids, which were collected in the Northwest Pacific Ocean in March and September in 2019. AIC analysis shows that beak length and pigmentation of the squids had positive relationship with increasing mantle length. The beak length saltated in 31 and 41 mm. The upper crest length (UCL), lateral wall length (ULWL) and lower lateral wall length (LLWL) had high linear regression slope, which indicates that they grew faster than other beak parts. The proportion of pigmentation in the upper hood length (UHL), lower hood length (LHL) and LLWL saltated in 31 mm, but that in ULWL saltated in 41 mm. ANOVA analysis shows significant phenotypic changes of E. chunii beak (P<0.05). The study of ontogenetic phenotype changes such as size and pigmentation provides references for research on ontogenetic feeding shift.
-
Keywords:
- Enoploteuthis chunii /
- Beak /
- Pigmentation /
- Growth characteristics /
- Western Pacific
-
珠江是中国南方第一大河，地处热带、亚热带季风气候区，年径流量高于3.3×1011 m3，水资源总量居全国第二，流域内海拔跨度大于3 000 m，地貌多样，包括高原、山地、丘陵和三角洲平原等。异质多样化的生境，孕育了珠江水系丰富的生物资源，鱼类物种多样性在全球大河流域生态系统中排名第四，特有种类数排名第二[1]。截至2015年7月，珠江水系记录鱼类达650种(亚种，包括河口鱼类)[2]，特有鱼类约126种，国家级水产种质资源保护品种约25种[3]，是中国重要的淡水渔业生产基地和水生生物资源基因库。
珠江水系鱼类生态学研究基础薄弱，陆奎贤[4]于1981—1983年首次开展了珠江水系渔业资源全面性调查，李捷等[5]分析了2005—2008年西江河段鱼类多样性，帅方敏等[3]分析了2015年珠江水系鱼类空间格局。2015年以来珠江水系渔业管理制度和水利水电工程建设已发生了重大调整，例如珠江水系规划设立自然保护区约30个，国家级水产种质资源保护区20个；自2017年开始延长禁渔期，每年3月1日—6月30日在珠江流域实行全面禁渔[6]；浔江上游的大藤峡水利枢纽工程于2015年9月全面开工。2015年以来珠江水系鱼类群落时空格局尚处于空白状态。
西江作为珠江水系的干流和主体部分，发源于云南曲靖市马雄山，由南盘江、红水河、黔江、浔江和西江等河段组成，全长2 214 km。本研究选取了浔江和西江河段作为珠江水系鱼类群落研究的代表江段。2016—2018年笔者研究团队在西江水域开展了鱼类种群和水环境调查，通过分析鱼类群落时空格局和环境影响，补充更新了西江鱼类群落动态信息，服务珠江水系生态保护。
1. 材料与方法
2016年春季(3月)和夏季(6月)，2017年冬季(2月)、夏季(7月和8月)和秋季(10月)，2018年冬季(1月)、春季(4月)、夏季(7月)和秋季(10月)，在西江桂平、濛江镇、梧州、郁南、德庆和肇庆江段开展了渔获物采样调查(图1)，通过在各站点收集20船次的渔获物，采集渔具为三重刺网和地笼，每种渔具分别收集10船次，作业时间为8~10 h；在采样站点的干流中泓线区域，利用YSI同步监测表层水体的水温、盐度、溶解氧含量、pH、氧化还原电位、溶解性固体总量、电导率等，用透明度盘测定水体透明度，采集表层水体在实验室测定叶绿素a含量；在电子地图上测定采样点至河口距离，从全国水雨情网站(http://xxfb.mwr.cn/)下载采样点的径流量数据。
![]()
图 1 西江鱼类种群采样和水环境监测站点Figure 1. Sites for fish sampling and environment variables monitoring in West River在PRIMER 6.0中计算每次采样渔获鱼类的物种数(S)、香浓威纳指数(H')和均匀度指数(J')。其中
$H' = - \sum {{P_i}\ln {P_i}} $ ，Pi为渔获鱼类中种i的个体数比例；$ J' = H'/\log S $ ，S为物种数。利用方差分析(ANOVA)和非参数检验方法比较分析了不同季节和站点的S、H'和J'的差异，分析过程在R软件中完成。群落间的物种组成表现为物种丢失或物种获得，若仅含1个过程，则物种数少的群落是另一个群落的嵌套子集，群落变化为群落嵌套；若包含2个过程，则群落间存在物种周转，群落变化呈现为物种周转[7-8]。本研究分别分析了不同季节和不同站点间的群落变化过程。通过相对重要性指数(IRI)确定鱼类物种的重要性，IRI=(N%+W%)×F%，其中N%为某物种个体数占渔获鱼类个体总数量的百分比例，W%为某物种质量(生物量)占渔获鱼类总质量的百分比例，F%为物种出现的站点数占调查总站数的百分比例。本研究中将IRI>1 000的物种定为优势种。
为去除稀有种影响，取出现频次≥3的鱼类N%、W%数据，在PRIMER 6中利用置换多因素方差分析(PERMANOVA)比较了不同季节和站点中鱼类群落结构的变化，并利用百分比相似性分析(SIMPER)解析了导致群落结构差异的重要物种。对出现频次≥3的鱼类N%、W%数据，在Hellinger转换后进行DCA排序，最长梯度长度均<3，在环境相关性分析中适合选择冗余分析(Redundancy analysis，RDA)[9]。为分析重要鱼类物种的环境梯度要求，取最大值>10%的鱼类N%、W%数据，进行RDA分析。对环境因子数据进行了log (x+1)转换，为避免环境因子间的共线性，在RDA分析中去除具有最高膨胀因子(VIF=14.88)的电导率环境因子，并使用双终止原则前向选择[10]筛选出对物种分布具有重要影响的环境变量。
2. 结果
2.1 物种组成多样性
2016—2018年在西江桂平至肇庆江段，共采集鱼类12目29科99种(附录，详见http://dx.doi.org/10.12131/20190142的资源附件)，其中鲤形目3科52种，鲤科46种，鳅科4种，条鳅科2种；鲈形目7科16种；鲇形目7科15种；鲱形目2科2种；胡瓜鱼目1科2种；合鳃鱼目2科3种；颌针鱼目1科1种；鳗鲡目1科2种；脂鲤目2科2种；鳉形目1科1种；鲻形目1科2种；鲀形目1科1种；杂交鲟1尾，未定种鱼类2尾。
每次采样鱼类物种数S为10~45 (图2)，时空差异不显著(非参数检验，P>0.05)。H'为0.90~2.93 (图2)，时空差异不显著(方差分析，P>0.05)。J'为0.30~0.89 (图2)，时空差异不显著(非参数检验，P>0.05)。春季共采集鱼类51种，夏季共采集77种，秋季共采集60种，冬季共采集44种，相邻季节间群落变化包含了物种获得和物种丢失两个过程，季节间鱼类群落呈群落周转格局。在桂平江段共采集鱼类61种，在濛江镇共采集51种，在梧州共采集52种，在郁南共采集39种，在德庆共采集40种，在肇庆共采集42种，相邻江段间群落变化包含了物种获得和物种丢失两个过程，西江桂平至肇庆江段鱼类群落呈群落周转格局。
![]()
图 2 西江平均鱼类物种数、香浓威纳指数和均匀度指数空间分布Figure 2. Spatial distribution of average fish species number, Shannon-Wiener index and evenness index in West River2.2 鱼类优势种时空分布
鱼类IRI表明，西江主要优势种为鲮(Cirrhinus molitorella)、广东鲂(Megalobrama terminalis)、赤眼鳟(Squaliobarbus curriculus)和鲤(Cyprinus carpio)等(图3)。优势种组成具有一定的季节波动性，春季主要优势种为广东鲂、鲤、鲢(Hypophthalmichthys molitrix)和赤眼鳟等，夏季主要优势种为鲮、赤眼鳟和广东鲂等，秋季主要优势种为鲮、赤眼鳟和草鱼(Ctenopharyngodon idella)等，冬季主要优势种为鲮、广东鲂和鲤等。鲮IRI在夏、秋和冬季均为物种间最高值：鲮秋季IRI最高，优势地位最高，春季其IRI降低，低于广东鲂和赤眼鳟。广东鲂在一年四季均为优势种，且在冬季具有最高的IRI。赤眼鳟IRI表现为夏季>冬季>秋季>春季。鲤在冬季具有较高的IRI，春季和夏季次之，但秋季IRI<1 000。草鱼、䱗(Hemiculter leucisculus)、鲢和瓦氏黄颡鱼 (Pelteobagrus vachelli) 的IRI仅在个别季节高于1 000，总体优势度偏低。
![]()
图 3 西江鱼类优势种相对重要性指数季节动态变化Figure 3. Seasonal index of relative importance variation of dominant fish species in West River不同站点的鱼类优势种组成具有一定的差异(图4)。桂平江段主要优势种为鲮、鲤和尼罗罗非鱼；濛江镇的主要优势种为赤眼鳟、鲤和鲮；梧州江段的主要优势种为鲮、广东鲂和赤眼鳟；郁南江段的主要优势种为鲮和广东鲂；德庆江段的主要优势种为赤眼鳟、鲮和瓦氏黄颡鱼；肇庆江段的主要优势种为广东鲂、鲢和鲤。
![]()
图 4 西江鱼类优势种相对重要性指数空间动态变化Figure 4. Spatial index of relative importance variation of dominant fish species in West River优势种的N%和W%体现了时空分布差异(表1)。鲮N%在夏、秋、冬季为种间最高，鲮W%在秋季为种间最高，鲮N%和W%均在秋季达到峰值，夏季次之，春季最低，其中夏、秋、冬季的郁南以及夏、秋季的梧州分布有较高资源量(N%和W%)的鲮。广东鲂W%为冬季>春季>夏季>秋季，其N%季节变化同样是冬季最高，秋季最低。赤眼鳟W%在夏季最高，尤其是德庆江段，春季最低；其N%在冬季最高，春季最低。鲤W%在冬季的濛江镇和肇庆江段最高，其次是春季的濛江镇和肇庆江段以及夏季的桂平江段，秋季W%属季节间最低；其N%在冬季最高，秋季最低。
表 1 西江优势种鱼类的数量百分比、质量百分比、出现站位百分比和相对重要性指数Table 1. N%, W%, F% and relative importance index of dominant fish species in West River物种
Species季节
Season数量百分比
N%质量百分比
W%出现站位百分比
F%相对重要性指数
IRI鲮 C. molitorella 春季 7.54 4.78 83.33 1 026.92 夏季 16.73 18.00 100.00 3 472.68 秋季 25.88 25.53 100.00 5 140.18 冬季 22.76 15.46 100.00 3 822.01 广东鲂 M. terminalis 春季 10.52 20.18 83.33 2 558.13 夏季 10.74 15.80 75.00 1 990.54 秋季 9.32 9.83 62.50 1 197.11 冬季 12.29 24.83 80.00 2 969.70 赤眼鳟 S. curriculus 春季 8.14 6.10 100.00 1 424.15 夏季 10.03 20.40 91.67 2 789.60 秋季 11.14 10.13 100.00 2 126.82 冬季 13.02 9.04 100.00 2 206.35 鲤 C. carpio 春季 7.20 14.00 83.33 1 766.70 夏季 3.62 8.89 91.67 1 146.91 秋季 2.99 7.60 87.50 926.81 冬季 8.96 25.58 100.00 3 453.70 草鱼 C. idella 春季 0.90 7.83 50.00 436.81 夏季 2.05 7.69 75.00 730.31 秋季 4.87 16.04 100.00 2 090.89 冬季 1.46 2.92 80.00 350.15 鲢 H. molitrix 春季 8.43 24.10 50.00 1 626.66 夏季 5.52 16.62 66.67 1 476.20 秋季 1.80 6.43 87.50 720.34 冬季 0.88 7.59 80.00 677.74 䱗 H. leucisculus 春季 12.64 3.52 66.67 1 077.56 夏季 9.82 1.40 58.33 654.81 秋季 12.19 2.04 75.00 1 067.71 冬季 6.22 1.67 80.00 631.27 瓦氏黄颡鱼 P. vachelli 春季 13.35 7.37 66.67 1 381.01 夏季 7.64 5.79 66.67 895.24 秋季 1.99 0.71 87.50 236.14 冬季 7.95 2.08 80.00 802.65 2.3 西江鱼类群落结构分析
PERMANOVA分析表明不同站点间鱼类N%存在显著性差异(表2，P<0.05)，其中桂平与梧州和肇庆江段存在显著差异(PERMANOVA比对分析，P<0.05)。桂平—梧州间鱼类N%的差异主要体现在鲮、银(Squalidus argentatus)、瓦氏黄颡鱼、广东鲂、赤眼鳟、尼罗罗非鱼、银飘鱼(Pseudolaubuca sinensis)、卷口鱼(Ptychidio jordani)、粗唇鮠 (Leiocassis crassilabris)和鳊(Parabramis pekinensis，表3，SIMPER分析)。桂平—肇庆间鱼类N%的差异主要体现广东鲂、银、鲮、鲤、鲢、黄颡鱼(Pelteobagrus fulvidraco)、花鰶 (Clupanodon thrissa)、尼罗罗非鱼、赤眼鳟、卷口鱼、䱗、银飘鱼和粗唇鮠 (表3，SIMPER分析)。
表 2 不同季节和站点间的鱼类数量百分比、质量百分比组成的PERMANOVA分析Table 2. PERMANOVA analysis of fish number percentage and weight percentage variation among seasons and sites方程
Function变异源
Source自由度
df平方和
SS均方
MS伪F
Pseudo-FP (蒙特卡罗检验)
P (MC)数量百分比 (季节+站点+季节×站点)
N% (season+site+season×site)季节 3 7 768.8 2 589.6 1.331 5 0.176 站点 5 22 074 4 414.8 2.270 0 0.002 季节×站点 11 23 407 2 127.9 1.094 1 0.321 残差 11 21 393 1 944.8 总和 30 76 267 质量百分比 (季节+站点+季节×站点)
W% (season+site+season×site)季节 3 10 173 3 391.1 2.219 7 0.013 站点 5 18 012 3 602.5 2.358 0 0.004 季节×站点 11 20 862 1 896.5 1.241 4 0.167 残差 11 16 805 1 527.7 总和 30 69 358 表 3 西江鱼类数量百分比空间结构差异源的SIMPER分析Table 3. SIMPER analysis of fish number percentage difference for spatial structure in West River物种
Species桂平—梧州 GP−WZ 物种
Species桂平—肇庆 GP−ZQ 平均相异性
Av. Diss贡献比
Contrib/%平均相异性
Av.Diss贡献比
Contrib/%鲮 C. molitorella 10.09 13.48 广东鲂 M. terminalis 8.47 10.35 银鮈 S. argentatus 5.73 7.65 银鮈 S. argentatus 5.72 6.99 瓦氏黄颡鱼 P. vachelli 5.44 7.28 鲮 C. molitorella 5.63 6.89 广东鲂 M. terminalis 5.37 7.17 鲤 C. carpio 4.40 5.38 赤眼鳟 S. curriculus 4.24 5.66 鲢 H. molitrix 4.27 5.22 尼罗罗非鱼 O. niloticus 4.08 5.45 黄颡鱼 P. fulvidraco 4.17 5.10 银飘鱼 P. sinensis 3.92 5.25 花鰶 C. thrissa 4.09 5.00 卷口鱼 P. jordani 3.01 4.02 尼罗罗非鱼 O. niloticus 4.05 4.96 粗唇鮠 L. crassilabris 2.66 3.56 赤眼鳟 S.curriculus 3.46 4.22 鳊 P. pekinensis 2.49 3.33 卷口鱼 P. jordani 3.09 3.77 鲫 C. auratus 2.26 3.01 䱗 H. leucisculus 2.98 3.64 斑鳠 H. guttatus 2.03 2.72 银飘鱼 P. sinensis 2.97 3.63 䱗 H. leucisculus 2.01 2.69 粗唇鮠 L. crassilabris 2.67 3.27 PERMANOVA分析表明不同季节和不同站点间鱼类W%均存在显著性差异(表2，P<0.05)。PERMANOVA比对分析表明夏—冬季间鱼类W%存在显著差异(P<0.05)，主要源于鲤、鲮、广东鲂、赤眼鳟、鲢、草鱼和尼罗罗非鱼W%的季节差异(表4)。桂平—梧州鱼类W%存在显著差异(PERMANOVA比对分析，P<0.05)，主要源于鲮、广东鲂、鲤、鲢、尼罗罗非鱼、草鱼、鳡(Elopichthys bambusa)和卷口鱼W%的空间差异(表4)。
表 4 西江鱼类质量百分比时空结构差异源的SIMPER分析Table 4. SIMPER analysis of fish weight percentage difference for spatial-temporal structure in West River夏季—冬季 Summer−Winter 桂平—梧州 GP−WZ 物种
Species平均相异性
Av. Diss贡献比
Contrib/%物种
Species平均相异性
Av. Diss贡献比
Contrib/%鲤 C. carpio 11.56 17.79 鲮 C. molitorella 10.26 14.49 鲮 C. molitorella 9.82 15.12 广东鲂 M. terminalis 8.84 12.48 广东鲂 M. terminalis 8.15 12.55 鲤 C. carpio 6.12 8.64 赤眼鳟 S. curriculus 6.54 10.06 鲢 H. molitrix 5.61 7.93 鲢 H. molitrix 5.78 8.89 尼罗罗非鱼 O. niloticus 5.23 7.38 草鱼 C. idella 2.88 4.44 草鱼 C. idella 4.62 6.52 尼罗罗非鱼 O. niloticus 2.71 4.18 鳡 E. bambusa 3.65 5.16 瓦氏黄颡鱼 P. vachelli 2.16 3.32 卷口鱼 P. jordani 2.80 3.95 花鰶 C. thrissa 1.82 2.80 赤眼鳟 S. curriculus 2.51 3.55 银鲴 X. argentea 1.35 2.07 鳙 A. nobilis 2.47 3.49 2.4 群落环境影响分析
西江桂平至肇庆江段距河口190~520 km，水体盐度0.08~0.18，pH 7.03~8.28，水温17.5~29.7 ℃，溶解氧5.37~9.12 mg·L−1，氧化还原电位10.4~399.4 mv，透明度20~150 cm，溶解性固体总质量浓度0.11~0.22 g·L−1，叶绿素a质量浓度0.28~6.50 mg·L−1，径流量2 033~22 585 m3·s−1。水温、透明度和径流量季节变化显著(P<0.05，Kruskal-Wallis非参数检验，图5-A—C)。Kruskal-Wallis多重比较检验表明，夏季和秋季水温显著高于春季和冬季(P<0.05，图5-A)，夏季透明度显著低于春季和冬季(P<0.05，图5-B)，夏季径流量显著高于春季和冬季(P<0.05，图5-C)。水温、透明度和径流量的空间变化见图5-E—F，t检验表明长洲坝水利枢纽上游(桂平和濛江镇)与下游(梧州、郁南、德庆和肇庆)间的径流量存在显著差异(P<0.05)。而且透明度和径流量数据存在显著的幂函数关系(图6)，两者的log对数存在显著的线性相关关系(P<0.001)。水体透明度的波动变化，受控于水文调度和河流径流量的变化。
![]()
图 5 季节和站点间水温、透明度和径流量的波动变化Figure 5. Variation of water temperature, transparency and runoff among seasons and sites![]()
图 6 透明度和径流量的幂函数关系Figure 6. Power function relationship between transparency and average runoff鱼类N%-环境因子的RDA分析(图7-A)中，所有环境因子共解释物种数据变异的49.92%，其中第一轴解释了13.00%，第二轴解释了11.71%，主要影响环境因子为距河口距离、水温、径流量、透明度和溶解性固体总量，通过前进式筛选的环境因子为距河口距离和透明度(P<0.05)。RDA三序图上物种在环境变量轴上的投影，反映了物种分布的环境梯度；物种和原点的连线与环境变量轴之间的夹角，反映了物种多度和环境梯度间的相关性[11]。重要物种N%-环境因子的RDA分析表明，瓦氏黄颡鱼、卷口鱼、粗唇鮠、尼罗罗非鱼、银飘鱼、银和南方拟䱗 (Pseudohemiculter dispar) N%与距河口距离正相关，与径流量负相关，主要分布在研究河段上游、低径流量、高透明度的低温水域(图7-B)；鲤和黄颡鱼N%与TDS负相关，分布在TDS含量低的水域；花鰶、鲢和广东鲂N%与距河口距离负相关，与径流量正相关，分布在高径流量的下游江段；鳊、黄尾鲴(Xenocypris davidi)、草鱼和舌虾虎鱼(Glossogobius giuris)分布在研究河段偏下游的高径流量的高温水域；鲮N%与TDS及水温正相关，分布在TDS含量高的高温水域；䱗N%与溶解性固体总量及距河口距离正相关，分布在研究河段上游的TDS含量较高的水域；赤眼鳟和翘嘴鲌(Culter alburnus)位于三序图的原点位置，对环境条件具有较弱的选择倾向。
![]()
图 7 西江鱼类数据-环境变量RDA分析图鱼类N% (A、B)和鱼类W% (C、D) (△代表采样点，+代表物种)；黑色箭头代表环境变量: Chl a. 叶绿素a含量；Dis.. 采样点至河口的距离；DO. 溶解氧含量；ORP. 氧化还原电位；Runoff. 采样前一个月内的平均径流量；Sal.. 盐度；Tra.. 透明度；TDS. 溶解性固体总量；Tem.. 水温；红色字母为物种学名缩写，Cluth. 花鰶；Cteid. 草鱼；Eloba. 鳡；Squcu. 赤眼鳟；Psesi. 银飘；Parpe. 鳊；Hemle. 䱗；Psedi. 南方拟䱗；Megte. 广东鲂；Culal. 翘嘴鲌；Culre. 海南鲌；Xenar. 银鲴；Xenda. 黄尾鲴；Arino. 鳙；Hypmo. 鲢；Squar. 银鮈；Ostsa. 纹唇鱼；Cirmo. 鲮；Cirmr. 麦鲮；Ptyjo. 卷口鱼；Cypca. 鲤；Pelfu. 黄颡鱼；Pelva. 瓦氏黄颡鱼；Leicr. 粗唇鮠；Latma. 花鲈；Oreni. 尼罗罗非鱼；Glogi. 舌虾虎鱼；Rhigi. 子陵吻虾虎鱼Figure 7. RDA plot of fish data and environmental variables in West RiverFish N% (A, B) and fish W% (C, D) (△. Sampling site,+. Species); black arrows represent environmental variables: Chl a. Chl a contents; Dis.. Distance between sampling site and the Pearl River estuary; DO. Dissolved oxygen; ORP. Oxidation-reduction potential; Runoff. Mean daily runoff during the month before the sampling day; Sal.. Salinity; Tra.. Water transparency; TDS. Total dissolved solids; Tem.. Water temperature; red letters indicate species, Cluth. C. thrissa; Cteid. C. idella; Eloba. E. bambusa; Squcu. S. curriculus; Psesi. P. sinensis; Parpe. P. pekinensis; Hemle. H. leucisculus; Psedi. P. dispar; Megte. M. terminalis; Culal. C. alburnus; Culre. C. recurviceps; Xenar. X. argentea; Xenda. X. davidi; Arino. A. nobilis; Hypmo. H. molitrix; Squar. S. argentatus; Ostsa. Osteochilus salsburyi; Cirmo. C. molitorella; Cirmr. C. mrigala; Ptyjo. P. jordani; Cypca. C. carpio; Pelfu. P. fulvidraco; Pelva. P. vachelli; Leicr. L. crassilabris; Latma. Lateolabrax maculatus; Oreni. O. niloticus; Glogi. G. giuris; Rhigi. Rhinogobius giurinus鱼类W%-环境因子的RDA分析(图7-C)中，所有环境因子共解释物种数据变异的49.61%，其中第一轴解释了15.30%，第二轴解释了13.01%，主要环境影响因子为水温、透明度、径流量和距河口距离，通过前进式筛选的环境因子为水温和距河口距离(P<0.05)。重要物种W%-环境因子的RDA分析表明，鲮和草鱼W%与水温和径流量正相关，它们分布在高水温、高径流量、低透明度、叶绿素含量高的水域；广东鲂和鲢W%与距河口距离负相关，倾向于分布在下游江段的高溶氧水域(图7-D)；鲤W%与水温和径流量负相关，倾向于分布在高透明度的低温水域；瓦氏黄颡鱼和赤眼鳟W%与距河口距离正相关，倾向于分布在研究河段的偏上游水域；䱗、麦鲮(Cirrhinus mrigala)和花鰶分布在中等偏高的水温和径流量的环境梯度中；花鲈、鳡、鳙和海南鲌分布在低径流量、低温的下游水域；粗唇鮠、卷口鱼和尼罗罗非鱼分布在中游江段的低温水域。
3. 讨论
3.1 物种多样性呈降低趋势
西江作为珠江水系的干流部分，具有郁江、柳江、北流河、桂江、贺江、罗定江和新兴江等支流汇入，年净流量占珠江总径流量的70%以上，西江水生态系统在很大程度上代表了珠江水生态系统。
在西江桂平至肇庆江段共采集鱼类12目29科99种，鲤科鱼类物种数占46.5%，采集到濒危鱼类花鳗鲡(Anguilla marmorata) 1尾和易危物种花鰶和长臀鮠 (Cranoglanis bouderius)多尾，未采集到鲟形目、赤魟 (Hemitrygon akajei)、鲥(Tenualosa reevesii)和舌鳎(Cynoglossidae)等名贵珍稀物种，确证了其在西江的濒危状态。20世纪80年代西江分布鱼类136种[4-5]，李捷等[5]于2005—2008年在西江梧州至肇庆江段采集鱼类13目29科96种，本研究在上述同一江段采集鱼类11目24科76种，没有采到银鱼科种类和鲽形目种类，西江鱼类物种多样性呈降低趋势。部分种类未被采集到，也反映了渔获物采样调查的随机性，今后应增加调查站点和调查频次。
研究江段外来物种入侵严重，本研究共采集到外来物种5目8科8种，分别为条纹鲮脂鲤(Prochilodus lineatus)、短盖巨脂鲤(Piaractus brachypomus)、麦鲮(C. mrigala)、革胡子鲇(Clarias gariepinus)、斑点叉尾鮰(Ictalurus punctatus)、豹纹翼甲鲶(Hypostomus plecostomus)、食蚊鱼(Gambusia affinis)和尼罗罗非鱼(O. niloticus)。在全球范围内，伴随河流生境受干扰程度加剧，例如水坝拦截加强了原分离生态系统的联系，促进了外来物种的成功定殖[12]，生物入侵已成为影响河流生物群落和生态系统的重要因素[13]。入侵种通过食物竞争、捕食关系和栖息地竞争等，导致土著种群生态位萎缩甚至种群灭绝[14]。故西江土著鱼类的保护和管理，一定程度上应从控制外来生物入侵着手。
3.2 鱼类群落的时空结构
西江桂平至肇庆江段鱼类物种组成和群落结构呈现出显著的时空演替。物种组成时空分布变化明显，为季节和空间的群落周转格局，同张东等[15]研究的河流中、下游间的分类群变化趋势一致。优势种组成和群落结构均存在显著季节差异，鱼类优势度的季节更替以鲮最为显著，夏、秋、冬季鲮为第一优势种，而在春季鲮优势度较低；夏—冬季间群落结构的季节差异，主要源于鲤、鲮、广东鲂、赤眼鳟、鲢、草鱼和尼罗罗非鱼W%的差异。
河流生态系统是一个具连续性和异质性特点的线性结构[16-17]，沿着这一线性结构理化环境梯度的改变，鱼类区系组成和群落结构往往发生空间演替[18]。本研究中圆筒体形的鲤科种类银、蛇和卷口鱼是近上游桂平和濛江镇江段的特有优势种，侧扁体形的鲱形目种类花鰶是下游肇庆江段的特有优势种。与本研究类似，欧洲河流上、下游间的鱼类区系和体形特征也发生了有规律的变化：从上游至下游“鲑科—鲤科”的变化[18]，以及从上游圆筒体形至中下游的侧扁体形的变化[8]。桂平江段与梧州、肇庆江段群落结构的显著差异，主要源于鲮、广东鲂、银等物种的空间分布差异。流域内鱼类群落空间变化与环境变化密切相关，连续性变化的环境条件下，群落多呈嵌套式空间格局，但跳跃式变化的环境条件，会导致群落呈周转格局[8]。因此，本研究中桂平江段较高的物种多样性，以及桂平—肇庆江段间鱼类分类群的周转格局，一定程度上由环境条件的非连续性跳跃变化引发。
3.3 影响鱼类群落的主要环境因子
物种组成和群落结构的时空动态主要由环境因子(人类干扰、气候和水文因子等)在时空尺度上的异质性引起[3]。本研究中RDA分析表明影响鱼类分布的主要因子为距河口距离、水温、透明度和径流量。
距河口距离通过影响鲮、广东鲂、银、鲤、鲢、尼罗罗非鱼等鱼类在不同江段间的分布，导致了桂平江段与梧州和肇庆江段间的群落结构差异。桂平与梧州江段间的长洲水利枢纽被推测为重要影响因素。水坝通过改变河流水文状况和河流连通性等影响鱼类群落[12]。2007年8月长洲水利枢纽开始下闸蓄水[19]，限制了洄游性鱼类广东鲂的生殖洄游，其上溯洄游路线已受到严重干扰，2012年以来，每年的6—7月在梧州江段形成大规模的广东鲂渔汛，而长洲坝上游的桂平和濛江镇江段的广东鲂的分布受到限制。长洲坝也改变了西江水文状况，将部分河段由动态环境转变为静态环境，桂平江段形成了低径流量的缓流库区，由于鲤具有中等水深和缓流静水区的倾向性[20]，吸引了鲤聚集于桂平江段并形成优势种。笔者发现自2004年底长洲水利枢纽初步截流以来，其下游鲤种群资源量显著低于上游江段。长洲水利枢纽阻隔引起的河流水文状况的改变和栖息地的片段化，是导致桂平与梧州、肇庆江段鱼类群落结构差异的重要因素。
水温是决定河流生态环境的重要因子，影响其他环境因子，另外作用于鱼类的分布、生长、繁殖和迁移分布等[21]。河流水温的季节性波动显著影响了水生物种的分布(河流连续性理论[17])，如高温促使鱼类生理需氧量增加同时又降低水体中溶氧饱和度，导致鱼类分布受限或死亡[16]。多数鱼类存在特定的适宜温度范围，适宜温度范围和河流水温共同决定了鱼类的时空分布[22]。鲮的季节性分布与其适温范围密切相关，鲮适宜水温为22~26 ℃[23]，本研究中夏、秋季平均水温分别为(27.86±1.62) ℃和(27.57±1.65) ℃，适宜于鲮繁殖产卵和广泛分布[24]；冬、春季平均水温分别为(18.25±0.87) ℃和(17.92±0.18) ℃，不利于鲮的广泛分布，但冬季郁南江段分布有大量鲮，初步推测郁南江段是鲮的越冬场。冬季鲤优势度较高，源于其较强的耐低温能力，其适宜温度为0~30 ℃，甚至可耐受0~4 ℃达数月之久[25]。RDA分析表明鲤、广东鲂、鲮和赤眼鳟具有不同的水温分布梯度(图7)，鲮分布在较高水温，赤眼鳟和广东鲂次之，鲤在低水温环境中具有较高的N%、W% (图7-B、D)。因此，水温的季节性波动显著影响了鱼类物种的N%、W%和IRI，导致了群落结构的季节性差异。
径流量和流速也是影响鱼类分布和繁殖产卵的重要环境因子。本研究中长洲坝上、下游江段间径流量的显著差异，以及对应站点间鱼类群落结构的显著差异，表明长洲坝已改变了西江水文环境，并作用于了鱼类分布。水坝拦截引起的水位、流速和径流量等水文状况的改变，是水坝影响鱼类群落的最广泛性作用机制[12]。另外，四大家鱼的繁殖产卵对径流量具有一定的水文需求，易雨君[26]模拟的四大家鱼栖息地适合度模型显示，日流量涨幅在800 m3·s−1时具有较高的家鱼产卵场适合度；珠江鳙大量繁殖的基础水文需求为持续2 d以上高于5×103 m3·s−1的径流量[27]。本研究中草鱼、鲢和鳙在夏、秋季具有较高的W%，同时夏、秋季具有持续整月的高于5×103 m3·s−1的径流量(图5-B)，日流量涨幅高于800 m3·s−1的月均天数分别为8.9 d和5 d，该水文条件满足了家鱼产卵需求，适于草鱼、鲢和鳙的分布和繁殖。
透明度与水深、径流量、水体理化因子以及水生植物生物量密切相关[28]，并作用于鱼类的分布、摄食和繁殖活动[3]。与多数江河鱼类群落和湖泊鱼类群落研究一致[3,28-29]，透明度对西江鱼类分布具有重要影响。透明度对鱼类分布的影响，与径流量关系密切，涨水引起营养盐浓度增大，但也会导致水体透明度降低[30]。本研究中水体透明度的波动变化，受控于水文调度和河流径流量的变化，两者协同影响鱼类分布(图6)。
综上，2016—2018年在西江桂平至肇庆江段共采集鱼类12目29科99种，与前期调查比较，鱼类多样性降低，外来种增多。鱼类群落物种组成变化为季节和空间的周转格局；物种优势度时空变化显著，优势种鲮、广东鲂、赤眼鳟和鲤的优势度均具有明显的季节和空间变化(图3、图4)；群落结构时空差异显著，主要表现为夏—冬季间的群落结构差异，以及桂平和梧州、肇庆江段间的群落结构差异；影响物种分布的主要环境因子为距河口距离、水温、透明度和径流量。西江鱼类群落空间断裂显著，一定程度上是长洲水利枢纽的拦截阻隔所致，应加强水坝过鱼通道管理等工作，以提高西江上、下游的连通性。
-
![]()
图 1 角质颚形态参数示意图
A. 头盖长 (HL)；B. 脊突长 (CL)；C. 喙长 (RL)；D. 侧壁长 (LWL)；E. 翼长 (WL)。
Figure 1. Scheme of morphometric measurements for beak
A. Hood length (HL); B. Crest length (CL); C. Rostrum length (RL); D. Lateral wall length (LWL); E. Wing length (WL).
![]()
图 3 角质颚各项形态参数与胴长关系
Figure 3. Relationship between each beak morphometric parameter and mantle length
![]()
图 4 角质颚色素沉积面积占比与胴长的关系
Figure 4. Relationship between proportion of beak pigmentation area and mantle length
![]()
图 5 角质颚各项形态参数随胴长变化图
Figure 5. Variation of beak morphometrics variation with mantle length
![]()
图 6 角质颚各部分色素面积占比随胴长变化图
Figure 6. Variation of proportion of beak pigmentation area with mantle length
表 1 富山武装乌贼采样信息
Table 1 Sampling information of E. chunii
采样海域
Sampling
sea area采样日期
Sampling
date样品数量
Sample
size胴长范围
Mantle length
range/mm148°00'E, 33°00'N 2019年3月21日 25 32~62 148°00'E, 31°99'N 2019年3月23日 4 43~49 148°10'E, 38°76'N 2019年9月5日 62 18~45 
下载: 导出CSV
表 2 富山武装乌贼角质颚形态参数值
Table 2 Morphometric parameters of beak of E. chunii mm
形态指标
Morphologic
index最大值
Maximum最小值
Minimum均值
Mean标准差
Standard
deviation上头盖长 UHL 6.16 3.34 4.48 0.86 上脊突长 UCL 8.68 4.57 6.14 1.18 上喙长 URL 2.18 0.87 1.45 0.31 上侧壁长 ULWL 7.09 3.56 4.95 0.98 上翼长 UWL 3.36 1.55 2.38 0.50 下头盖长 LHL 3.25 1.87 1.47 0.32 下脊突长 LCL 3.43 1.69 2.32 0.44 下喙长 LRL 2.60 1.01 1.60 0.38 下侧壁长 LLWL 6.45 3.33 4.56 0.92 下翼长 LWL 4.12 1.90 2.60 0.52
下载: 导出CSV
表 3 角质颚表型变化信息
Table 3 Information of beak phenotype change
胴长组
Mantle length
group/mm上颚 Upper beak 下颚 Lower beak 喙部
Beak侧壁色素沉积面积
Pigmentation area of
lateral wall翼部色素沉积面积
Pigmentation area
of wing侧壁脊
Ridge of
lateral wall侧壁色素沉积面积
Pigmentation area of
lateral wall18~31 较长 24% 无沉积 不明显 42% 32~41 尖且长 36% 较小 略明显 64% 42~62 尖且长 51% 较大 明显 65%
下载: 导出CSV
-
[1] WATANABE H, KUBODERA T, MOKU M, et al. Diel vertical migration of squid in the warm core ring and cold water masses in the transition region of the western North Pacific[J]. Mar Ecol Prog Ser, 2006, 315: 187-197. doi: 10.3354/meps315187
[2] HONG S Y. A short note on Enoploteuthis (Paraenoploteuthis) chunii (Cepphlopoda: Enopoteuthidae) from the Korean water[J]. Anim Syst Evol Divers, 1992, 8(2): 331-335.
[3] AHN S, KIM H W, KIM S. Diet of long-beaked common dolphin (Delphinus capensis) in the East Sea, Korea[J]. Anim Cells Syst, 2014, 18(5): 340-350. doi: 10.1080/19768354.2014.955823
[4] SHIMOSE T, KIM H W, KIM S, et al. Ontogenetic diet shift of age-0 year Pacific bluefin tuna Thunnus orientalis[J]. J Fish Biol, 2013, 82(1): 263-276. doi: 10.1111/j.1095-8649.2012.03483.x
[5] KEAR A J. Morphology and function of the mandibular muscles in some coleoid cephalopods[J]. J Mar Biol Assoc UK, 1994, 74(4): 801-822. doi: 10.1017/S0025315400090068
[6] UYENO T A, KIER W M. Functional morphology of the cephalopod buccal mass: a novel joint type[J]. J Morphol, 2005, 264(2): 211-222. doi: 10.1002/jmor.10330
[7] FRANCO-SANTOS R M, VIDAL E A G. Beak development of early squid paralarvae (Cephalopoda: Teuthoidea) may reflect an adaptation to a specialized feeding mode[J]. Hydrobiologia, 2014, 725(1): 85-103. doi: 10.1007/s10750-013-1715-2
[8] 刘必林, 林静远, 方舟, 等. 头足类角质颚稳定同位素研究进展[J]. 海洋渔业, 2018, 40(2): 242-248. doi: 10.3969/j.issn.1004-2490.2018.02.014 [9] 方舟, 陈新军, 陆化杰, 等. 头足类角质颚研究进展Ⅰ——形态、结构与生长[J]. 海洋渔业, 2014, 36(1): 78-89. doi: 10.3969/j.issn.1004-2490.2014.01.013 [10] 方舟, 陈新军, 陆化杰, 等. 头足类角质颚研究进展Ⅱ——运动、分类与生态[J]. 海洋渔业, 2014, 36(2): 183-192. doi: 10.3969/j.issn.1004-2490.2014.02.013 [11] 刘必林, 陈新军, 方舟, 等. 利用角质颚研究头足类的年龄与生长[J]. 上海海洋大学学报, 2014, 23(6): 930-936. [12] 刘必林, 刘娜, 李建华, 等. 智利外海茎柔鱼角质颚微结构及其年龄与生长研究[J]. 南方水产科学, 2020, 16(1): 62-68. doi: 10.12131/20190116 [13] 林静远, 刘必林, 王立权. 基于角质颚微结构的剑尖枪乌贼的日龄与生长[J]. 华东理工大学学报(自然科学版), 2019, 45(5): 775-782. [14] 方舟, 陈新军, 瞿俊跃, 等. 北太平洋柔鱼角质颚形态及生长年间差异[J]. 上海海洋大学学报, 2020, 29(1): 109-120. doi: 10.12024/jsou.20181202485 [15] 陆化杰, 陈子越, 宁欣, 等. 中国南海西沙群岛海域鸢乌贼角质颚色素沉积变化[J]. 生态学杂志, 2020, 39(5): 1600-1608. [16] 陈炫妤, 陆化杰, 王洪浩, 等. 西北印度洋鸢乌贼角质颚色素沉积特性分析[J]. 动物学杂志, 2020, 55(4): 468-476. [17] 胡贯宇, 陈新军, 方舟. 个体生长对秘鲁外海茎柔鱼角质颚形态变化的影响[J]. 水产学报, 2016, 40(1): 36-44. [18] 刘必林, 陈新军. 印度洋西北海域鸢乌贼角质颚长度分析[J]. 渔业科学进展, 2010, 31(1): 8-14. doi: 10.3969/j.issn.1000-7075.2010.01.002 [19] 刘必林, 陈新军, 方舟, 等. 基于角质颚长度的头足类种类判别[J]. 海洋与湖沼, 2015, 46(6): 1365-1372. [20] 韩青鹏, 陆化杰, 陈新军, 等. 南海北部海域中国枪乌贼角质颚的形态学分析[J]. 南方水产科学, 2017, 13(4): 122-130. doi: 10.3969/j.issn.2095-0780.2017.04.015 [21] 李建华, 陈新军, 陈芃, 等. 赤道海域茎柔鱼角质颚形态及生长特征分析[J]. 海洋渔业, 2018, 40(4): 385-393. doi: 10.3969/j.issn.1004-2490.2018.04.001 [22] 方舟, 金岳, 胡飞飞, 等. 南海北部秋季杜氏枪乌贼角质颚形态和生长特征[J]. 南方水产科学, 2018, 14(4): 17-26. doi: 10.3969/j.issn.2095-0780.2018.04.003 [23] 何静茹, 陆化杰, 陈炫妤, 等. 冬春季西北印度洋鸢乌贼角质颚外部形态及生长特性[J]. 应用生态学报, 2020, 31(8): 2775-2784. [24] 刘必林, 林静远, 陈新军, 等. 西北太平洋柔鱼角质颚微结构及其生长纹周期性研究[J]. 海洋与湖沼, 2016, 47(4): 821-827. doi: 10.11693/hyhz20160300044 [25] 方舟, 陈新军, 李建华. 西南大西洋公海阿根廷滑柔鱼角质颚色素变化分析[J]. 水产学报, 2013, 37(2): 222-229. [26] 胡贯宇, 陈新军, 方舟. 秘鲁外海茎柔鱼角质颚色素沉积及影响因素的初步研究[J]. 海洋湖沼通报, 2017(2): 72-80. [27] ARKHIPKIN A I. Diversity in growth and longevity in short-lived animals: squid of the suborder Oegopsina[J]. Mar Freshw Res, 2004, 55: 341-355. doi: 10.1071/MF03202
[28] YOUNG R E. Aspects of the natural history of pelagic cephalopods of the Hawaiian mesopelagic-boundary region[J]. Pac Sci, 1995, 49(2): 143-155.
[29] UCHIKAWA K, SAKAI M, WAKABAYASHI T, et al. The relationship between paralarval feeding and morphological changes in the proboscis and beaks of the neon flying squid Ommastrephes bartramii[J]. Fish Sci, 2009, 75: 317-323. doi: 10.1007/s12562-008-0036-2
[30] CAMARILLO C S, SALINAS Z C, LAVANIEGOS B E, et al. Food in early life stages of Dosidicus gigas (Cephalopoda: Ommastrephidae) from the Gulf of California, Mexico[J]. J Mar Biol Assoc UK, 2013, 93(7): 1903-1910. doi: 10.1017/S0025315413000398
[31] FIELD J C, ELLIGER C, BALTZ K, et al. Foraging ecology and movement patterns of jumbo squid (Dosidicus gigas) in the California current system[J]. Deap-Sea Res II, 2013, 95(15): 37-51.
[32] 方舟, 陈新军, 陆化杰, 等. 北太平洋两个柔鱼群体角质颚形态及生长特征[J]. 生态学报, 2014, 34(19): 5405-5415. [33] CASTRO J J, HERNÁNDEZ-GARCÍA V. Ontogenetic changes in mouth structures, foraging behavior and habitat use of Scomber japonicus and Illex coindetii[J]. Sci Mar, 1995, 59(3/4): 347-355.
[34] 贾涛, 李纲, 陈新军, 刘必林, 等. 9−10月秘鲁外海茎柔鱼摄食习性的初步研究[J]. 上海海洋大学学报, 2010, 19(5): 663-667. [35] 刘必林, 陈新军. 头足类角质颚的研究进展[J]. 水产学报, 2009, 33(1): 157-164. [36] MISEREZ A, SCHNEBERK T, SUN C J, et al. The transition from stiff to compliant materials in squid beaks[J]. Science, 2008, 319(5871): 1816-1819. doi: 10.1126/science.1154117
[37] 林静远, 刘必林, 金宵. 茎柔鱼角质颚的机械强度特性[J]. 上海海洋大学学报, 2020, 29(3): 385-391. doi: 10.12024/jsou.20190502690 [38] MING C W, KWANG T S, SHIANG L H, et al. Food partitioning among three sympatric odontocetes (Grampus griseus, Lagenodelphis hosei and Stenella attenuate)[J]. Mar Mammal Sci, 2012, 28(2): E143-E157. doi: 10.1111/j.1748-7692.2011.00501.x
-
期刊类型引用(12)
1. 朱书礼，陈蔚涛，武智，夏雨果，杨计平，李跃飞，李捷. 基于环境DNA技术的珠江中下游鱼类多样性初步研究. 南方水产科学. 2024(01): 120-129 . 
本站查看
2. 李朝晖，单楠，王鹏伟，包萨茹，庞博，徐德琳，陈星瑶. 呼伦湖流域鱼类群落结构特征及与物理生境特征关系研究. 生态与农村环境学报. 2024(02): 233-244 . 百度学术
3. 夏雨果，陈蔚涛，李新辉，朱书礼，李捷，李跃飞. 广东省内陆江河鱼类多样性. 南方水产科学. 2024(04): 34-45 . 本站查看
4. 何佳洋，吴志强，黄亮亮，陈中兵，王冬杰，孙扬言，文浩，宋琦亮. 广西内陆水域外来鱼类入侵风险分析. 生物安全学报. 2023(02): 168-177 . 百度学术
5. 武智，李跃飞，朱书礼，杨计平，夏雨果，李捷. 基于渔业声学调查的珠江东塔产卵场鱼类栖息地适宜性研究. 南方水产科学. 2023(03): 11-18 . 本站查看
6. 王珍，蔡金妤，吉红香，张瀚，冯小丽. 广东省主要干流河道自然岸线保有率及分布情况研究. 广东水利水电. 2023(08): 99-103+114 . 百度学术
7. 李育森，赵会宏，王庆，王大鹏，罗永巨，周磊，黄世慧，余梵冬. 大藤峡大坝蓄水对浔江鱼类群落的影响. 生态学报. 2023(17): 7140-7156 . 百度学术
8. 谷先坤，刘燕山，唐晟凯，李大命，张彤晴，曾庆飞，毛志刚，陈辉辉，谷孝鸿. 网围拆除后东太湖鱼类功能群时空变化及影响因子. 湖泊科学. 2022(02): 569-581 . 百度学术
9. 朱书礼，陈蔚涛，李新辉，李捷，李跃飞. 柳江鱼类群落结构及多样性研究. 水生生物学报. 2022(03): 375-390 . 百度学术
10. 冯杰，吴志强，黄亮亮，高明慧，何佳洋，王一菲，李明锶. 西江来宾江段鱼类群落结构季节变化特征. 四川动物. 2022(06): 681-694 . 百度学术
11. 刘燕山，谷先坤，唐晟凯，李大命，曾庆飞，陆建明，毛志刚，殷稼雯，陈辉辉，沈冬冬，何浩然，张彤晴. 东太湖鱼类群落结构特征及其与环境因子的关系. 生态学报. 2021(02): 769-780 . 百度学术
12. 夏雨果，李跃飞，朱书礼，李捷，李新辉. 珠江流域草鱼和鲢单位捕捞努力量渔获量时空分布特征及温度影响. 南方水产科学. 2021(01): 10-16 . 本站查看
其他类型引用(16)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 28
粤公网安备 44010502001741号
