Fish community structure in mainstream of lower reaches of Jinsha River (Shuifu-Yibin section)
-
摘要:
金沙江下游干流 (水富—宜宾段) 是多种长江珍稀特有鱼类的幼鱼庇护所,也是大型水利工程建设较为密集的地区之一,其鱼类生物多样性受到严重威胁。为评估其鱼类资源现状,于2021—2023年在安边镇、三块石及柏溪镇3个江段开展了10次调查。结果显示,共采集鱼类8目20科50属69种,其中鲤形目最多 (48种);长江上游特有鱼类17种,国家重点保护野生动物7种,外来物种8种。群落优势种以铜鱼 (Coreius heterodon)、瓦氏黄颡鱼 (Tachysurus vachelli)、粗唇鮠 (Leiocassis crassilabris) 为主。生物多样性指数中仅Margalef丰富度指数呈现出时间和空间的显著性差异 (p<0.05)。ABC曲线结果表明,鱼类群落整体处于中度干扰状态,但大个体优势种有增加趋势。聚类分析结果显示,在56.71%的相似性水平上,2021和2022年的鱼类种类较为相似,归为一组;而2023年的种类与历史记录较为相似,归为另一组;与2021—2022年相比,2023年金沙江下游干流 (水富—宜宾段) 的鱼类种类有所增加。建议后续持续开展鱼类增殖放流、控制外来鱼类入侵以及加强科研监测等措施,进一步推动金沙江及长江鱼类种群数量的恢复。
Abstract:The lower reaches of the Jinsha River (Shuifu-Yibin section), a nursery habitat for various rare and endemic fish species of the Yangtze River, are one of the areas with dense construction of large hydrop projects. However, the fish biodiversity in this region is under severe threat. To assess the status of its fish resources, we had conducted ten surveys from 2021 to 2023 in three river sections: Anbian Town, Sankuashi and Baixi Town. The results show that a total of 69 fish species belonging to 8 orders, 20 families and 50 genera were collected, with the Cyiformes being the most abundant (48 species). There were 17 endemic fish species from the upper Yangtze River, 7 species of nationally protected wildlife, and 8 invasive species. The dominant species in the community included Coreius heterodon, Tachysurus vachelli, and Leiocassis crassilabris. Among the biodiversity indexes, only the Margalef richness index showed significant temporal and spatial differences (p<0.05). The ABC curve indicates that the fish community was in a state of moderate disturbance, but with a trend of increasing dominance of larger individuals. Cluster analysis shows that at a similarity level of 56.71%, the species from 2021 and 2022 were grouped together due to their similarities, while the ones in 2023 were more similar to historical records so they were grouped together. Compared with the results from 2021 to 2022, there had been an increase in fish species in the mainstream of the Jinsha River (Shuifu-Yibin section) in 2023. It is recommended that measures such as continuous fish stocking, controlling the invasion of non-native fish species, and strengthening scientific monitoring be implemented to further promote the recovery of fish populations in the Jinsha River and the Yangtze River.
-
柔鱼 (Ommastrephes bartramii) 属大洋性头足类,广泛分布于三大洋温带和亚热带海域[1],是海洋生态系统的重要组成部分,目前其商业性开发利用主要集中在北太平洋海域[2]。日本从1974年开始商业捕捞[3],中国于1994年开始进行商业性开发[4]。作为重要的大洋性经济柔鱼类,柔鱼资源的可持续开发利用越来越受到国内外研究者的重视。目前国内外学者对柔鱼的种群结构[5-7]、年龄生长[8-11]、摄食[12]等生物学特性进行了研究,发现其个体生长在不同种群、性别及季节存在明显差异。柔鱼生长迅速,生命周期通常为1年[13],其生长速度受生物和非生物影响。
胴长与体质量关系 (Length-mass relationship, LMR) 能够通过胴长来估计体质量,反之亦然,可以反映物种的生长模式 (等速还是异速生长),也可以反映物种的种内、种间关系[14]。了解胴长与体质量关系对头足类的种群划分和资源评估至关重要[15]。已有许多研究表明,头足类的生长会受到生物和非生物 (季节、地理和环境条件等因子) 的影响,因而其关系式中的参数存在一定的异质性。如Vasconcelos等[16]对阿尔加维海岸的乌贼 (Sepia officinalis) 生物学特性进行了研究,发现其胴长和体质量关系存在性别和地理位置上的差异。
在针对参数异质性研究中,线性混合模型 (Linear mixed model, LMM) 是一种应用极其广泛的统计学模型[17]。其线性预测同时包含固定效应和随机效应,并具有能研究多样型数据的优点[18],已广泛应用于多个领域。Ma等[19]使用线性混合模型对中国北部沿岸小黄鱼 (Larimichthy spolyactis) 的体长、体质量关系进行了研究,发现其体长、体质量在不同年份和海域存在显著性差异。瞿俊跃等[20]在对东海的蓝点马鲛 (Scomberomorus niphonius) 研究中使用了线性混合模型,模型中考虑了不同性别和月份的随机效应。以往研究发现不同群体的柔鱼生长特征亦有所不同,但性别和月间所产生的影响程度仍不清楚。因此,本文根据近些年我国鱿钓船在北太平洋生产调查期间所获得的渔获样本,对不同柔鱼群体的生长特征进行研究,并运用线性混合模型分析其性别与时间的差异,以期为北太平洋柔鱼的种群生长差异及资源评估与管理提供数据支撑和参考依据。
1. 材料与方法
1.1 数据来源
选取2016—2018年的柔鱼样本进行分析,样本主要采集于北太平洋海域 (172°08'E—176°25'W、38°42'N—39°43'N,150°42'E—167°00'E、38°00' N—43°13'N,图1),采集时间为5—11月 (表1),采样时间涵盖了不同群体柔鱼的主要生活史过程[21]。采集的柔鱼样本经冷冻保存运回实验室进行生物学分析。
表 1 不同月份柔鱼样本基本信息Table 1. Basic information of O. bartramii samples in different months群体
Population年份
Year月份
Month样本数量
Number of sample胴长范围
Range of mantle length/mm东部群体 Eastern population 2016—2017 5 394 163~500 6 456 212~516 西部群体 Western population 2016—2018 7 124 187~290 8 283 195~357 9 331 199~386 10 165 190~377 11 105 215~421 1.2 研究方法
1.2.1 群体划分
根据Yatsu等[5]的划分方法,170°E以东范围内为东部群体 (秋生群),170°E以西范围为西部群体 (冬春生群)。同时结合鱿钓船捕捞的时间和个体大小,以170°E为界,共获得东部群体样本850尾,西部群体样本1 008尾。
1.2.2 基础生物学测量
对解冻后的样本进行生物学测定,包括胴长 (Mantle length, ML)、体质量 (Body mass, BM)、性别、性腺成熟度等。使用皮尺测定样本的胴长,精确至1 mm;使用电子天平测定体质量,精确至0.1 g。根据Lipinski和 Underhil[22]将其性腺成熟度划分为Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ期,其中Ⅰ和Ⅱ为未成熟,Ⅲ期以上为成熟个体。
1.2.3 数据处理与分析
1) 采用t检验,分析不同群体间柔鱼的差异;2) 采用频率分析法分析样本胴长及体质量的组成,组间距分别为30 mm和300 g,绘制频率分布图,确定优势胴长组和体质量组,计算优势组的比例;3) 采用方差分析 (ANOVA)[23] 和多重比较 (LSD法)[24] 对不同月份柔鱼的差异进行检验,并对不同性别柔鱼胴长的差异进行检验,以分析不同因子对柔鱼生长的影响;4) 利用线性混合效应模型[25],建立胴长-体质量关系:
$$ m=a{L}^{b} $$ (1) 式中m为体质量 (g),L为胴长 (mm),a是条件因子,b是异速生长参数 (代表柔鱼胴长方向,体高方向和体宽方向3个维度的生长速度[26])。
1.3 异质性分析
本文使用胴长和体质量关系模型及9个线性混合模型来描述东西部群体柔鱼的胴长和体质量关系 (表2),其中9个LMM模型将月份和性别的影响作为参数a和b的随机效应,以探讨时间及生物学对柔鱼生长的影响。
表 2 不同群体柔鱼胴长-体质量关系模型的拟合效果指标Table 2. Fitting effect index (AIC) of mantle length-body mass relationship models for different O. bartramii populations种群
Population模型缩写
Model abbreviation模型
Model对数转化
Log-transformed赤池信息量准则
AIC均方根误差
RMSE东部群体
Eastern populationLMR m = a×Lb ln(m)=ln(a)+b×ln(L) −466.402 0.183 265 RSS m=a×L(b+RSS) ln(m)=ln(a)+(b+RSS)×ln(L) −472.338 0.181 617 RSI m=a×exp(RSI)×Lb ln(m)=[ln(a)+RSI]+b×ln(L) −472.27 0.181 641 RSI&S m=a×exp(RSI)×L(b+RSS) ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L) −469.057 0.181 428 RMS m=a×L(b+RMS) ln(m)=ln(a)+(b+RMS)×ln(L) −489.107 0.179 784 RMI m=a×exp(RMI)×Lb ln(m)=[ln(a)+RMI]+b×ln(L) −492.772 0.179367 RMI&S m=a×exp(RMI)×L(b+RSI) ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L) −595.137 0.166 468 RS&MS m=a×L(b+RSS+RMS) ln(m)=ln(a)+(b+RSS + RMS)×ln(L) −494.452 0.178 558 RS&MI m=a×exp(RSI)×exp(RMI)×Lb ln(m)=[ln(a)+RSI + RMI]+b×ln(L) −497.694 0.178 211 RS&MI&S m=a×exp(RSI)×exp(RMS)×L(b+RSI+RSI) ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L) −596.138 0.165 394 西部群体
Western populationLWR m = a×Lb ln(m)=ln(a)+b×ln(L) −1 795.63 0.099 003 RSS m=a×L(b+RSS) ln(m)=ln(a)+(b+RSS)×ln(L) −1 797.15 0.098 614 RSI m=a×exp(RSI)×Lb ln(m)=[ln(a)+RSI]+b×ln(L) −1 797.31 0.098 603 RSI&S m=a×exp(RSI)×L(b+RSS) ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L) −1 796.72 0.098 512 RMS m=a×L(b+RMS) ln(m)=ln(a)+(b+RMS)×ln(L) −1 807.04 0.097 626 RMI m=a×exp(RMI)×Lb ln(m)=[ln(a)+RMI]+b×ln(L) −1 807.41 0.097 615 RMI&S m=a×exp(RMI)×L(b+RSI) ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L) −1 807.59 0.097 215 RS&MS m=a×L(b+RSS+RMS) ln(m)=ln(a)+(b+RSS + RMS)×ln(L) −1 808.14 0.097 347 RS&MI m=a×exp(RSI)×exp(RMI)×Lb ln(m)=[ln(a)+RSI + RMI]+b×ln(L) −1 808.56 0.097 333 RS&MI&S m=a×exp(RSI)×exp(RMS)×L(b+RSI+RSI) ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L) −1 808.09 0.096 976 注:第一列为柔鱼的不同群体,第二列为第三、四列模型的缩写 RSS、RSI、RSI&S、RMS、RMI、RMI&S、RS&MS、RS&MI 和 RS&MI&S 分别为斜率 b 的性别随机效应、截距 ln(a) 的性别随机效应、截距 ln(a) 与斜率 b的性别随机效应、斜率 b 的月份随机效应、截距 ln(a) 的月份随机效应、截距 ln(a) 与斜率 b 的月份随机效应、斜率 b的性别与月份随机效应、截距 ln(a) 的性别与月份随机效应,以及截距 ln(a) 与斜率 b 的性别与月份随机效应Note: The first column shows different O. bartramii populations. The second column shows the abbreviations of models detailed in the third and fourth columns. RSS, RSI, RSI&S, RMS, RMI, RMI&S, RS&MS, RS&MI和RS&MI&S indicate random effects on slope (b) in sexual difference, random effects on intercept [ln (a)] in sexual difference, random effects on intercept [ln (a)] and slope (b) in sexual difference, random effects on slope (b) in months, random effects on intercept [ln (a)] in months, random effects on intercept [ln (a)] and slope (b) in months, random effects on slope (b) in sexual difference and months, random effects on intercept [ln (a)] fin sexual difference and months, and random effects on intercept [ln (a)] and slope (b) in sexual difference and months, respectively. 线性混合模型是一类非常重要的统计模型,既含有固定效应也含有随机效应。一般来说,如果抽取某一个因子的若干个特定水平,目的是对这些水平的效应进行估计或比较,这个因子通常会当成固定效应;如果因子的若干水平是来自总体的随机抽样,目的是要通过该样本推断总体,那么该因子通常被当作随机效应,其模型可表示为[18]:
$$ y=X\tau +Zu+e $$ (2) 式中y为所有观测值构成的向量,τ为所有固定效应构成的向量,X为固定效应的关联矩阵,u是所有随机效应构成的向量,Z为随机效应的关联矩阵,e为随机残差向量。
本研究以性别和月份作为随机效应进行模型拟合,基于随机效应的不同组合,共建立9个线性混合模型 (表2),其随机效应分别来自性别 (RSS、RSI和RSI&S)、月份 (RMS、RMI和RMI&S)、性别和月份 (RS&MS、RSS&MI和RS&MI&S)。
使用F检验评估各因子对柔鱼胴长-体质量关系的显著性[18]。利用赤池信息量准则 (Akaike Information criterion, AIC) 来比较10个模型的拟合程度[27],同时利用均方根误差 (Root mean square error, RMSE) 来判断最适模型;其值越接近0,表示模型的预测值与样本的真实值差异越小,即模型的拟合效果越好。所有模型的构建过程均采用R软件的 “lme4”程序包来完成[28]。
2. 结果
2.1 胴长与体质量组成
对东、西部群体比较发现,2个群体间胴长、体质量存在极显著性差异 (P<0.01),因此分别对两个群体的生物学特征进行分析。发现东部群体柔鱼雌雄比较高 (7.02∶1);西部群体雌雄比较低 (1.03∶1)。
经t检验,东、西部群体雌性与雄性柔鱼胴长、体质量均存在极显著性差异 (P<0.01)。结果显示,东部群体雌性个体胴长、体质量为163~516 mm、222~4 167 g,优势胴长、体质量为240~270 mm和330~420 mm、300~600 g和900~2 100 g,占总数的80.10%、83.70%,平均胴长、体质量为351.04 mm、1 366.10 g;雄性个体胴长、体质量范围为212~280 mm、292~665 g,优势胴长、体质量为210~270 mm、600~900 g,占总数的98.11%、97.17%,平均胴长、体质量为246.00 mm、395.20 g (图2)。
西部群体雌性个体胴长、体质量范围为185~421 mm、170~2 328 g,优势胴长、体质量为210~330 mm、300~900 g,占总数的85.35%、67.185%,平均胴长、体质量为268.55 mm、610.51 g;雄性个体胴长、体质量范围为190~395 mm、169~1640 g,优势胴长、体质量为210~330 mm、300~900 g,占总数的91.93%、79.84%,平均胴长、体质量为258.77 mm、538.14 g (图3)。比较两个群体的渔获物组成发现,东部群体雌雄数量相差较大,而西部群体则相近。东部群体柔鱼个体明显大于西部群体。
方差分析表明,东部群体和西部群体月份间胴长、体质量均存在极显著性差异 (P<0.01)。结果显示,5月东部群体柔鱼雌性个体较大,而雄性个体则相反。西部群体柔鱼雌雄个体表现出相同的趋势,在7—10月间随着月份的增加,雌雄个体逐渐增大,10月的个体最大;11月则略微下降。总的来看,各月份雌性个体均大于雄性个体 (图4)。
图 4 柔鱼不同月份胴长体质量分布图a. 东部群体胴长分布图;b. 东部群体体质量分布图;c. 西部群体胴长分布图;d. 西部群体体质量分布图Figure 4. Distribution of mantle length and body mass of O. bartramii in different monthsa. Distribution of mantle length for eastern population; b. Distribution of body mass for eastern population; c. Distribution of mantle length for western population t; d. Distribution of body mass for western population2.2 胴长和体质量参数分析
由于不同月份东、西部群体柔鱼胴长和体质量,以及雌、雄个体的胴长和体质量均存在极显著性差异,因此,本研究对不同性别柔鱼胴长和体质量关系 (图5),以及不同月份柔鱼胴长和体质量关系分别进行讨论 (表3)。
图 5 不同群体柔鱼体质量-胴长关系a. 东部群体雌性体质量-胴长关系;b. 东部群体雄性体质量-胴长关系;c. 西部群体雌性体质量-胴长关系;d. 西部群体雄性性体质量-胴长关系Figure 5. Mantle length-body mass relationship for different O. bartramii populationsa. Relationship between BM and ML of female for eastern population; b. Relationship between BM and ML of male for eastern population; c. Relationship between BM and ML of female for western population; d. Relationship between BM and ML of male for western population表 3 不同月份东西部群体柔鱼胴长与体质量关系Table 3. Mantle length-body mass relationship for different O. bartramii populations in different months群体
Population月份
Month胴长与体质量关系
Mantle length-mass relationship样本数量
Number of samples判定系数R2
Coefficient of determinationP 东部群体
Eastern population5 m =1.9×10−3L2.295 7 394 0.670 6 0.001 6 m =1×10−5L3.160 1 456 0.973 8 0.001 西部群体
Western population7 m = 9×10−6L3.205 3 124 0.889 5 0.001 8 m = 2×10−5L3.090 9 283 0.934 5 0.001 9 m = 2×10−5L3.059 7 331 0.952 0 0.001 10 m = 1×10−5L3.178 3 165 0.972 8 0.001 11 m = 5×10−5L2.902 8 105 0.947 7 0.001 胴长体质量关系模型的结果显示,东、西部群体柔鱼胴长体质量关系具有明显的性别和时间差异。从不同性别来看,东部群体雌性的a值比雄性小,b值则与之相反;西部群体亦然 (表4)。从不同月份来看,东部群体5月的a值比6月大,b值则与之相反。西部群体a值在11月最大,8月和9月次之,7月最小;b值月份分布由高到低依次为7月>10月>8月>9月>11月 (表3)。东部群体中a和b值的性别差异范围分别为0.069 68和0.367 26;而a和b值的月份差异分别为0.306 10和0.822 80,说明东部群体的时间差异大于性别差异。西部群体中a和b值的性别差异范围分别为0.027 45和0.156 62;而a和b值的月份差异分别为0.056 84和0.320 31,说明西部群体的时间差异大于性别差异。
表 4 不同性别东西部群体柔鱼胴长与体质量关系Table 4. Mantle length-body mass relationship of O. bartramii of different sex for different populations群体
Population性别
Sex胴长与体质量关系
Mantle length-mass relationship样本数量
Number of samples判定系数R2
Coefficient of determinationP 东部群体
Eastern population雌性 m=3×10−5L2.981 8 744 0.892 9 0.001 雄性 m=2×10−4L2.621 0 106 0.794 3 0.001 西部群体
Western population雌性 m=1×10−5L3.159 9 512 0.968 1 0.001 雄性 m =2×10−5L3.103 5 496 0.941 5 0.001 2.3 模型选择
本研究通过建立1个广义线性模型 (Generalized linear model, GLM) 和9个线性混合模型 (Linear mixed model, LMM) 来描述胴长和体质量关系 (表2),其中9个LMM模型分别讨论了不同月份影响下a与b的随机效应和不同性别影响下a与b的随机效应。结果表明,东部群体中,同时具有性别和月份对两个参数的随机效应的线性混合模型对柔鱼胴长-体质量关系的拟合效果最好,AIC为−596.138,RMSE为0.165 394;西部群体亦是如此,其中,AIC为−1 808.09,RMSE为0.096 976 (表2)。综合发现,性别和月份随机效应模型均优于一般的LMR模型,证明东、西部群体柔鱼均存在月份和性别上的差异。
3. 讨论
本研究采集的柔鱼样本中,东部群体5—6月雌性个体胴长范围为163~516 mm,雄性为212~280 mm。西部群体7—11月雌、雄性个体胴长范围分别为185~421 mm和190~395 mm。总体来看,东部群体柔鱼个体明显大于西部群体,这与黄洪亮等[29]对1996—2001年柔鱼渔获物胴长的研究结果一致。Ichii等[30]认为,东、西部群体柔鱼大小的差异可能是由过渡区叶绿素锋区和海表面温度的季节性变化引起,东部群体孵化时位于生产力较高的海域,西部群体则位于低产区。Katugin[31]对柔鱼样本进行种群基因差异性分析,发现东西部群体柔鱼在等位基因频率分布和遗传变异水平上存在显著性差异。因此,东、西部群体柔鱼大小的差异还可能与内在基因表达有关。
根据LMM (RS&MI&S) 的结果,北太平洋东部群体柔鱼条件因子a (0.000 021~0.000 46) 和西部群体条件因子a (0.000 013~0.000 017) 的估计值显著大于1996—2001年黄洪亮等[29]研究得到的东部群体a值 (0.000 006) 和西部群体a值 (0.000 004)。Hile[32]指出鱿鱼的种间和种内生长模式具有明显的差异性,其b介于2.5~4.0,本研究中东部群体柔鱼的b估计值 (2.731 4~3.050 51) 和西部群体的b估计值 (3.093 3~3.13536) 恰好在此范围内,且东部群体的值小于西部群体。东部群体柔鱼的b固定值 (2.763 4) 小于3,说明东部群体柔鱼呈负异速增长趋势,随着胴长的增加,柔鱼的肥满度降低,身体往细长方向发展;而西部群体柔鱼的b固定值 (3.117 19) 大于3,说明西部群体柔鱼呈正异速增长趋势,随着胴长的增加,柔鱼的肥满度升高,胴宽也增大。导致这一差异的原因可能是由于饵料与环境因素的不同。有研究表明,由于黑潮与亲潮的交汇区含有丰富的饵料,可能会造成太平洋西部海域柔鱼的摄食强度高于东部海域[33-34]。另外,海水温度变化对头足类的生长也有一定影响[35],东、西部群体的季节性差异会导致其生活水域温度的差异,因此,水温的不同可能会造成其生长模式的不同,从而导致不同群体b值的差异。
根据LMM (RS&MI&S) 的结果,本研究中,东、西部群体各月份的胴长与体质量均有极显著的差异性。东部群体柔鱼5月b (2.259 7) 小于3,呈现负异速增长;而6月b (3.1601) 大于3,呈现正异速增长。由于东部群体雄性样本较少,因此b的变化主要取决于雌性的生长变化。东部群体中,雌性个体5月到达亚北极边界海域,6—7月洄游至高生产力的亚北极锋区的南部海域[11],从而使得6月雌性柔鱼的摄食强度更高。这一结果可解释本研究中东部群体b的变化规律。西部群体中7—10月b均大于3,而11月b (2.9028) 最低,小于3,反映出11月的柔鱼较其他月份更偏向细长。研究表明,西部群体柔鱼8—10月主要分布在40°N—46°N亲潮前锋区及其周围海域的索饵场,其饵料丰富;10—11月后柔鱼开始向南进行产卵洄游,在产卵的过程中会导致柔鱼的肥满度降低[21]。这一现象也可解释西部群体柔鱼b的变化。
东、西部群体柔鱼雌雄的胴长与体质量均有极显著差异,存在显著的雌雄二态性。本研究中东部群体柔鱼雌雄比高达7.02∶1,而西部群体仅为1.03∶1,群体间性比差异显著 (P<0.01),这与东、西部群体洄游路线不同有关[30,36]。西部群体的雌、雄个体均生活在同一海域,有着相似的洄游路线;而东部群体中雌、雄个体的洄游路径有很大差异,雌性个体洄游方式和西部群体类似,先向北洄游,待个体成熟后再往南洄游至产卵场;雄性个体在整个生长过程中一直栖息于北太平洋副热带海域,直到7月开始向南产卵洄游。本研究发现,东、西部群体柔鱼中的雌性异速生长参数b均大于雄性,表明雌性的肥满度高于雄性,这一结果与乌贼 (Sepia officinalis)[16]、茎柔鱼 (Dosidicus gigas)[37] 相似。Domingues等[38-39]认为不同性别间乌贼的生长差异与生殖器官的发育和质量有关。乌贼在接近成熟时,雌性个体在生殖上投入的精力比在体细胞生长上投入的精力多,从而导致雌性的生殖器官变得相对较重 (相当于体质量),因此雌性的形状丰满度比雄性高。由此可推测,不同性别柔鱼生长差异可能与生殖器官的发育有关。
本文中任意一个线性混合模型的拟合程度均优于一般的LMR模型,其中,包含参数a和b的性别和月份随机效应LMM (SS&MI&S) 拟合效果最好,更能体现出生物学上的异质性。简单的LWR模型,常构建单因子回归模型,对多因子复合影响的模型能力较差[20]。而混合效应模型能把性别和时间的随机效应在单个模型中体现出来,从而更有效、更合理地估计时间和生长的差异,进一步证实了此模型在异质性研究中的优势。北太平洋柔鱼存在着不同的群体,在对其进行资源评估时需分开讨论。个体生长特征的异质性可能直接导致其种群动力学的变化,从而影响资源评估的准确性[40]。因此,本文根据不同年份、不同群体的柔鱼数据构建胴长体质量关系模型,并考虑到月份与性别的差异,从而更好地评估其生长特征,进而推断其资源量。由于本文仅考虑了性别和月份的随机效应,存在一定不足;未来的研究应充分考虑季节、水域和环境因子等对不同群体柔鱼生长的影响,更全面地研究不同群体柔鱼个体生长特征的变化。
-
表 1 金沙江下游干流 (水富—宜宾江段) 渔获物调查采样点信息
Table 1 Information of fish sampling sites in mainstream of lower reaches of Jinsha River (Shuifu-Yibin section)
序号
No.调查样点
Sampling site经度
Longitude/(°E)纬度
Latitude/(°N)1 安边镇上下游1 km范围内
1 km upstream and downstream of Anbian Town104.430 8 28.625 9 2 三块石上下游1 km范围内
1 km upstream and downstream of Sankuaishi104.489 7 28.660 9 3 柏溪镇上下游1 km范围内
1 km upstream and downstream of Baixi Town104.539 3 28.687 9 表 2 金沙江下游干流 (水富—宜宾段) 不同年份和样点鱼类优势种的相对重要性指数
Table 2 Index of relative importance of main fish species of different time and sample in mainstream of lower reaches of Jinsha River (Shuifu–Yibin section)
物种Species 年份 Year 样点 Sample 2021 2022 2023 安边镇
Anbian Town三块石
Sankuaishi柏溪镇
Baixi Town长江鲟 Acipenser dabryanus 807.2 — 871.6 — 2 890.5 576.5 无须高墨头鱼 Ageneiogarra imberba — — 907.9 — — — 鲤 Cyprinus carpio — 622.9 — — — — 䱗 Hemiculter leucisculus — 1 312.4 704.5 — 1 129.8 — 寡鳞飘鱼 Pseudolaubuca engraulis 612.2 — — — 156.4 — 团头鲂 Megalobrama amblycephala — — — — 542.2 — 宜昌鳅鮀 Gobiobotia filifer — — 1 042.6 1 238.9 — — 圆口铜鱼 Coreius guichenoti 546.1 — — — — — 铜鱼 Coreius heterodon 1 508.1 2 693.3 6 658.5 5 627.7 1 709.1 4 938.0 蛇鮈 Saurogobio dabryi 946.6 — — — 564.8 — 瓦氏黄颡鱼 Tachysurus vachelli 8 160.3 9 632.8 3 725.6 5 042.5 4 520.4 8 086.0 粗唇鮠 Leiocassis crassilabris 2 078.8 1 534.6 — 693.5 — 1 572.4 梭鲈 Sander lucioperca — 1 767.4 — — 1 717.4 — 表 3 金沙江下游干流 (水富—宜宾段) 渔获物组成情况
Table 3 Composition of fishery catches in mainstream of lower reaches of Jinsha River (Shuifu-Yibin section)
物种
Species尾数
Number质量
Mass/g尾数百分比
Number
percentage/%质量百分比
Mass
percentage/%相对重要性
指数
IRI群落地位
Community
status斑点叉尾 鮰 Ictalurus punctatus 17 1 236.0 0.38 0.26 63.8 III 长江鲟 Acipenser dabryanus 39 47 760.3 0.88 9.89 2 152.6 I 杂交鲟 Acipenser sp. 1 1 710.0 0.02 0.35 12.5 III 胭脂鱼 Myxocyprinus asiaticus 9 1 939.8 0.20 0.40 80.5 III 长薄鳅 Leptobotia elongata 15 825.7 0.34 0.17 50.8 III 红唇薄鳅 Leptobotia rubrilabris 1 41.7 0.02 0.01 1.0 IV 中华沙鳅 Sinibotia superciliaris 81 829.0 1.82 0.17 199.1 II 宽体沙鳅 Sinibotia reevesae 86 1 329.0 1.93 0.28 220.7 II 泥鳅 Misgurnus anguillicaudatus 2 8.7 0.04 0.00 1.6 IV 大鳞副泥鳅 Paramisgurnus dabryanus 2 32.4 0.04 0.01 1.7 IV 戴氏山鳅Claea dabryi 1 7.2 0.02 0.00 0.8 IV 犁头鳅 Lepturichthys fimbriatus 30 306.2 0.67 0.06 73.7 III 短身金沙鳅 Jinshaia abbreviata 4 45.5 0.09 0.01 3.3 IV 中华金沙鳅 Jinshaia sinensis 2 17.4 0.04 0.00 3.2 IV 中华倒刺鲃 Spinibarbus sinensis 8 4 530.8 0.18 0.94 37.3 III 白甲鱼 Onychostoma simum 7 4 726.8 0.16 0.98 75.7 III 鲮 Cirrhinus molitorella * 1 1 014.0 0.02 0.21 7.7 IV 无须高墨头鱼Ageneiogarra imberba 41 21 591.0 0.92 4.47 718.7 I 四川裂腹鱼 Schizothorax kozlovi 2 84.8 0.04 0.02 2.1 IV 岩原鲤Procypris rabaudi 25 8 373.1 0.56 1.73 535.5 I 鲤 Cyprinus carpio 14 8 994.6 0.31 1.86 507.8 I 鲫 Carassius auratus 16 2 085.59 0.36 0.43 184.6 II 宽鳍鱲 Zacco platypus 19 208.3 0.43 0.04 47.0 III 马口鱼 Opsariichthys bidens 3 94.3 0.07 0.02 5.8 IV 鲢 Hypophthalmichthys molitrix 1 1414.0 0.02 0.29 10.5 III 贝氏䱗 Hemiculter bleekeri 35 1 026.1 0.79 0.21 66.6 III 䱗 Hemiculter leucisculus 244 5 243.74 5.48 1.09 1 313.2 I 张氏䱗 Hemiculter tchangi 2 144.5 0.04 0.03 5.0 IV 似鳊 Pseudobrama simoni 2 73.5 0.04 0.02 4.0 IV 四川华鳊 Sinibrama taeniatus 1 42.2 0.02 0.01 1.0 IV 寡鳞飘鱼 Pseudolaubuca engraulis 82 3 024.4 1.84 0.63 164.5 II 银飘鱼 Pseudolaubuca sinensis 4 226.8 0.09 0.05 4.6 IV 翘嘴鲌 Culter alburnus 6 3 284.3 0.13 0.68 162.9 II 团头鲂 Megalobrama amblycephala* 30 14 604.3 0.67 3.02 862.6 I 厚颌鲂 Megalobrama pellegrini 3 3 186.5 0.07 0.66 48.5 III 黄尾鲴 Xenocypris davidi 8 4 617.5 0.18 0.96 75.7 III 圆吻鲴 Distoechodon tumirostris 7 9 606.5 0.16 1.99 71.5 III 高体鳑鲏 Rhodeus ocellatus 2 16.4 0.04 0.00 3.2 IV 宜昌鳅鮀 Gobiobotia filifer 211 3 657.03 4.74 0.76 1 099.3 I 鳅鮀 Gobiobotia pappenheimi 1 9.0 0.02 0.00 0.8 IV 异鳔鳅鮀 Xenophysogobio boulengeri 13 153.7 0.29 0.03 32.4 III 银鮈 Squalidus argentatus 19 233.7 0.43 0.05 31.7 III 唇䱻 Hemibarbus labeo 36 3 257.5 0.81 0.67 247.2 II 圆筒吻鮈 Rhinogobio cylindricus 5 333.6 0.11 0.07 6.0 IV 长鳍吻鮈 Rhinogobio ventralis 1 140.0 0.02 0.03 1.7 IV 吻鮈 Rhinogobio typus 26 4 598.5 0.58 0.95 409.6 II 圆口铜鱼 Coreius guichenoti 35 10 147.6 0.79 2.10 673.6 I 铜鱼 Coreius heterodon 451 176 264.7 10.13 36.49 12 431.6 I 蛇鮈 Saurogobio dabryi 97 3 125.3 2.18 0.65 565.2 I 斑点蛇鮈 Saurogobio punctatus 9 570.2 0.20 0.12 21.3 III 裸腹片唇鮈 Platysmacheilus nudiventris 5 40.4 0.11 0.01 4.0 IV 长须黄颡鱼 Tachysurus eupogon 1 40.0 0.02 0.01 1.0 IV 光泽黄颡鱼 Tachysurus nitidus 5 225.1 0.11 0.05 10.6 III 黄颡鱼 Tachysurus fulvidraco 28 884.46 0.63 0.18 54.1 III 瓦氏黄颡鱼 Tachysurus vachelli 2199 7 5316.7 49.39 15.59 19 495.4 I 粗唇鮠 Leiocassis crassilabris 365 16 843.7 8.20 3.49 3 505.6 I 长吻鮠 Leiocassis longirostris 11 796.5 0.25 0.16 68.7 III 切尾拟鲿 Pseudobagrus truncatus 3 86.6 0.07 0.02 5.7 IV 大鳍鳠 Hemibagrus macropterus 2 240.7 0.04 0.05 6.3 IV 白缘䱀 Liobagrus marginatus 10 110.83 0.22 0.02 24.8 III 中华纹胸鮡 Glyptothorax sinensis 8 26.7 0.18 0.01 12.3 III 鲇 Silurus asotus 4 3745.8 0.09 0.78 115.4 II 大口鲇 Silurus meridionalis 6 349.41 0.13 0.07 13.8 III 子陵吻虾虎 Rhinogobius giurinus 9 34.8 0.20 0.01 27.9 III 梭鲈 Sander lucioperca 25 26 282 0.56 5.44 1 200.4 I 齐氏罗非鱼 Coptodon zillii 1 63.2 0.02 0.01 1.2 IV 蓝鳃太阳鱼 Lepomis macrochirus 2 159.2 0.04 0.03 2.6 IV 鳜 Siniperca chuatsi 7 726.6 0.16 0.15 20.5 III 斑鳜 Siniperca scherzeri 4 305.8 0.09 0.06 5.1 IV 注:*. 外来种;I. 优势种;II. 重要种;III. 常见种;IV. 偶见种。 Note: *. Exotic specie;I. Dominant species; II. Important species; III. Common species; IV. Occasional species. 表 4 金沙江下游干流 (水富—宜宾段) 不同时间和样点鱼类物种多样性F、G和G-F指数
Table 4 F-index, G-index, G-F-index of different time and samples in mainstream of lower reaches of Jinsha River (Shuifu-Yibin section)
指标
Index年份 Year 样点 Sample 2021 2022 2023 安边镇
Anbian Town三块石
Sankuaishi柏溪镇
Baixi Town科 Family 12 13 19 17 14 16 属 Genus 24 27 48 43 30 35 种 Species 32 32 60 51 37 48 F指数 F-index 5.262 6.232 9.256 8.509 5.686 6.226 G指数 G-index 3.103 3.249 3.808 3.696 3.669 3.644 G-F指数 G-F index 0.410 0.479 0.589 0.566 0.355 0.415 表 5 2000—2023年金沙江下游干流 (水富—宜宾段) 优势种组成变化
Table 5 Changes of dominant species composition in mainstream of lower reaches of Jinsha River (Shuifu-Yibin section) during 2000–2003
物种
Species年份 Year 2000a 2007b 2008b 2008—2011c 2011d 2021e 2022e 2023e 长江鲟 Acipenser dabryanus — — — — — + — + 长薄鳅 Leptobotia elongata — — — — + — — — 中华沙鳅 Sinibotia superciliaris — — — — + — — — 无须高墨头鱼 Ageneiogarra imberba — — — — — — — + 鲤 Cyprinus carpio — — — + — — + — 鲫 Carassius auratus — — — — + — — — 鲢 Hypophthalmichthys molitrix — — — — + — — — 䱗 Hemiculter leucisculus — — — + — — + + 寡鳞飘鱼 Pseudolaubuca engraulis + — — — — + — — 宜昌鳅鮀 Gobiobotia filifer — + + — — — — + 异鳔鳅鮀 Xenophysogobio boulengeri — + + — — — — — 圆口铜鱼 Coreius guichenoti + +* +* +* +* + - — 铜鱼 Coreius heterodon — — — — + + + +* 蛇鮈 Saurogobio dabryi — + + + — + — — 瓦氏黄颡鱼 Tachysurus vachelli +* + + +* + +* +* + 粗唇鮠 Leiocassis crassilabris — — — — — + + — 鲇 Silurus asotus — — — + + — — — 大口鲇 Silurus meridionalis — — — + + — — — 梭鲈Sander lucioperca — — — — — — + — 注:a. 参考文献刘清等[20] (2005);b. 参考文献危起伟等[40] (2012);c. 参考文献高少波等[17] (2013);d. 参考文献李雷等[21] (2013);e.本次调查;+. 出现;—. 未出现;*. 绝对优势物种。 Note: a. Referenced from Liu, et al [20] (2005); b. Referenced from Wei, et al[40] (2012); c. Referenced from Gao, et al[17] (2013); d. Referenced from Li, et al[21] (2013); e. This survey; +. Detected species, —. Undected species; *. Absolute advantage species. -
[1] NORTHCOTE T G. Fish in the structure and function of freshwater ecosystems: a "top-down" view[J]. Can J Fish Aquat Sci, 1988, 45: 361-379. doi: 10.1139/f88-044
[2] 周向峰. 长江上游江安段鱼类资源调查研究[D]. 重庆: 西南大学, 2023: 2. [3] 夏雨果, 陈蔚涛, 李新辉, 等. 广东省内陆江河鱼类多样性[J]. 南方水产科学, 2024, 20(4): 34-45. doi: 10.12131/20230251 [4] 陆孝平, 富曾慈. 中国主要江河水系要览[M]. 北京: 中国水利水电出版社, 2010: 63-96. [5] 王忠原, 朱挺兵, 胡飞飞, 等. 十年禁渔初期金沙江干流鱼类群落结构及物种多样性[J]. 中国水产科学, 2024, 31(3): 286-300. [6] 周扬, 王宁, 温潇雄, 等. 长江上游梯级水电开发对区段水温累积影响及生态效应研究[J]. 水力发电学报, 2024, 43(11): 1-16. doi: 10.11660/slfdxb.20241101 [7] 孙赫英, 隋晓云, 何德奎, 等. 金沙江流域鱼类的系统保护规划研究[J]. 水生生物学报, 2019, 43(S1): 110-118. doi: 10.7541/2019.174 [8] 吴江, 吴明森. 金沙江的鱼类区系[J]. 四川动物, 1990, 9(3): 23-26. [9] 张春光, 杨君兴, 赵亚辉, 等. 金沙江流域鱼类[M]. 北京: 科学出版社, 2019: 20-23. [10] 吴江, 吴明森. 关于金沙江石鼓到宜宾段鱼类资源的概况及其利用问题[J]. 西南师范学院学报(自然科学版), 1985(1): 80-87. [11] CHEN D Q, XIONG F, WANG K, et al. Status of research on Yangtze fish biology and fisheries[J]. Environ Bio Fish, 2009, 85: 337-357. doi: 10.1007/s10641-009-9517-0
[12] 熊飞, 郭祺, 张伟, 等. 金沙江下游向家坝库区鱼类群落结构空间格局[J]. 水生态学杂志, 2024, 45(4): 82-91. [13] 胡兴坤, 邵科, 阙延福, 等. 金沙江下游宜宾段产漂流性卵鱼类早期资源动态变化[J]. 水生态学杂志, 2023, 44(2): 73-80. [14] 高少波, 唐会元, 乔晔, 等. 金沙江下游干流鱼类资源现状研究[J]. 水生态学杂志, 2013, 34(1): 44-49. [15] 孟宝, 张继飞, 叶华, 等. 长江上游珍稀特有鱼类国家级自然保护区鱼类产卵场功能现状分析及保护启示[J]. 长江流域资源与环境, 2019, 28(11): 2772-2785. [16] 陈大庆, 常剑波, 顾洪宾. 金沙江一期工程对保护区生态环境的影响与对策[J]. 长江科学院院报, 2005, 22(2): 21-24. doi: 10.3969/j.issn.1001-5485.2005.02.006 [17] 高少波, 唐会元, 陈胜, 等. 金沙江一期工程对保护区圆口铜鱼早期资源补充的影响[J]. 水生态学杂志, 2015, 36(2): 6-10. [18] 曹文宣. 长江上游水电梯级开发的水域生态修复问题[J]. 长江技术经济, 2019, 3(2): 5-10. [19] 龚廷登, 杨伟阶, 何滔, 等. 金沙江水富段水生生物的群落结构及水质监测[J]. 淡水渔业, 2014, 44(4): 25-34. doi: 10.3969/j.issn.1000-6907.2014.04.005 [20] 刘清, 苗志国, 谢从新, 等. 长江宜宾江段渔业资源调查[J]. 水产科学, 2005, 24(7): 47-49. doi: 10.3969/j.issn.1003-1111.2005.07.015 [21] 李雷, 危起伟, 吴金明, 等. 长江宜宾江段渔业资源现状调查[J]. 长江流域资源与环境, 2013, 22(11): 1449-1457. [22] 林俊强, 彭期冬, 任杰, 等. 赤水河与金沙江下游河段鱼类生境条件的相似性分析[J]. 淡水渔业, 2014, 44(6): 93-99. doi: 10.3969/j.issn.1000-6907.2014.06.016 [23] 李浩林, 赵亚辉, 张洁, 等. 金沙江下游与长江上游珍稀、特有鱼类国家级自然保护区鱼类物种多样性比[J]. 淡水渔业, 2014, 44(6): 104-108. doi: 10.3969/j.issn.1000-6907.2014.06.018 [24] 孙立元, 危起伟, 张辉, 等. 基于水声学的长江上游向家坝至宜宾江段鱼类空间分布特征[J]. 淡水渔业, 2014, 44(1): 53-58. doi: 10.3969/j.issn.1000-6907.2014.01.010 [25] 李祥艳 , 田辉伍, 蒲艳, 等. 长江上游宜宾江段鱼类早期资源现状研究[J]. 渔业科学进展, 2022, 43(4): 93-104. [26] 张觉民, 何志辉. 内陆水域渔业自然资源调查手册[M]. 北京: 中国农业出版社, 1991: 242-289. [27] 杨志, 乔晔, 张轶超, 等. 长江中上游圆口铜鱼的种群死亡特征及其物种保护[J]. 水生态学杂志, 2009, 2(2): 50-55. [28] 丁瑞华. 四川鱼类志[M]. 成都: 四川科学技术出版社, 1994: 1-641. [29] 褚新洛, 陈银瑞. 云南鱼类志上册[M]. 北京: 科学出版社, 1989: 1-377. [30] 褚新洛, 陈银瑞. 云南鱼类志下册[M]. 北京: 科学出版社, 1990: 1-313. [31] 危起伟, 吴金明. 长江上游珍稀特有鱼类国家级自然保护区鱼类图集[M]. 北京: 科学出版社, 2015: 1-156. [32] 马思琦. 鱼类分类和功能多样性研究:以济南地区鱼类为例[D]. 大连: 大连海洋大学, 2020: 12. [33] 蒋志刚, 纪力强. 鸟兽物种多样性测度的G-F指数方法[J]. 生物多样性, 1999, 7(3): 220-225. doi: 10.3321/j.issn:1005-0094.1999.03.010 [34] WARWICK R M. A new method for detecting pollution effects on marine macrobenthic communities[J]. Mar Biol, 1986, 92(4): 557-562. doi: 10.1007/BF00392515
[35] YEMANE D, FIELD J G, LESLIE R W. Exploring the effects of fishing on fish assemblages using Abundance Biomass Comparison (ABC) curves[J]. ICES J Mar Sci, 2005, 62(3): 374-379. doi: 10.1016/j.icesjms.2005.01.009
[36] 李圣法. 以数量生物量比较曲线评价东海鱼类群落的状况[J]. 中国水产科学, 2008, 15(1): 136-144. doi: 10.3321/j.issn:1005-8737.2008.01.018 [37] CLARKE K R, GORLEY R N. PRIMER V6: user manual/tutorial PRIMER-E[M]. Plymouth, UK: PRIMER-E, 2006: 57-68.
[38] 郭朋军, 江新琴, 俞存根, 等. 舟山沿岸渔场春秋季鱼类群落结构特征分析[J]. 渔业科学进展, 2020, 41(4): 1-11. [39] 唐会元, 朱其广, 金瑶, 等. 长江支流大宁河巫溪段鱼类群落结构的年际变化及保护对策[J]. 水生态学杂志, 2023, 44(6): 63-71. [40] 危起伟. 长江上游珍稀特有鱼类国家级自然保护区科学考察报告[M]. 北京: 科学出版社, 2012: 81-227. [41] 刘飞, 林鹏程, 黎明政, 等. 长江流域鱼类资源现状与保护对策[J]. 水生生物学报, 2019, 43(S1): 144-156. doi: 10.7541/2019.177 [42] 高明慧, 何力劲, 翁士创, 等. 武江乐昌段鱼类组成及物种多样性的年际变化[J/OL]. 水生态学杂志, 2024[2024-08-06]. [43] 段辛斌, 刘绍平, 熊飞, 等. 长江上游干流春季禁渔前后三年渔获物结构和生物多样性分析[J]. 长江流域资源与环境, 2008, 17(6): 878-885. doi: 10.3969/j.issn.1004-8227.2008.06.010 -
期刊类型引用(4)
1. 韩霈武,王岩,方舟,陈新军. 北太平洋柔鱼不同群体耳石日增量对海洋环境的响应研究. 海洋学报. 2022(01): 101-112 . 百度学术
2. 韩霈武,王岩,方舟,陈新军. 基于耳石信息的北太平洋东西部群体柔鱼日龄和生长特征. 生态学杂志. 2022(02): 393-403 . 百度学术
3. 薛薇,李楠,方舟. 剑尖枪乌贼个体生长年间差异及其影响因素. 上海海洋大学学报. 2022(03): 658-669 . 百度学术
4. 朱文斌,朱海晨,王雅丽,张亚洲,卢占晖,崔国辰. 基于线性混合效应模型的日本鳀幼鱼叉长-体重关系异质性. 应用生态学报. 2021(12): 4532-4538 . 百度学术
其他类型引用(6)
计量
- 文章访问数: 19
- HTML全文浏览量: 0
- PDF下载量: 0
- 被引次数: 10