Page 60 - 《水产学报》2026年第3期
P. 60

3 期                                     水    产    学    报                                 50 卷

              繁殖过程 ，而后         2003  年三峡水库蓄水、2004       年       1    材料与方法
                       [3]
              葛洲坝下河势调整、2008           年葛洲坝下航道疏浚等
              一系列涉水工程又对中华鲟的自然繁殖产生了进                             1.1    基于中华鲟生活史的简化参数模型
              一步的影响      [4-5] ，虽然中央和地方采取了禁止商业                     基于中华鲟相关文献 ，本研究将中华鲟生
                                                                                      [2]
              捕捞、开展增殖放流、列为国家一级保护野生动                            活史过程简化为包含          16  个控制节点的参数模型。
              物、建立中华鲟保护区、季节性休渔、禁止科研
                                                               16  个控制节点为：繁殖洄游亲本数 → 繁殖洄游
              捕捞等一系列保护措施            [2, 6-7] ，中华鲟自然种群仍
                                                               群体的雌雄比 → 洄游过程中繁殖群体的生存率 →
              然持续萎缩 ，并于          2013  年出现自然繁殖中断 ，             繁殖洄游群体性腺发育成熟比例 → 产卵过程中的
                                                        [9]
                         [8]
              2017—2024  年已连续     8  年未监测到中华鲟自然繁               雌雄比阈值      (限制性指标) → 产卵窗口适配率 →
              殖，自然种群岌岌可危           [3, 10] 。                   产卵场繁殖参与率 → 产卵场容量              (限制性指标) →
                   中华鲟为何会陷入濒临灭绝的困境？其关键                         雌性亲鲟怀卵量 → 产卵受精-生存率 → 孵化-生存
              影响因素是什么？关于此问题，学界已有不少争                            率 → 仔稚幼鱼降河生存率 → 长江口适应调整期
              论。总体来看，有归因于大坝建设的 ，有归因                            生存率 → 海洋群体生存率的年龄分布 → 初次性
                                                [11]
              于过度捕捞的       [12] ，有归因于水文状况改变的            [13] ，  腺发育启动繁殖洄游的年龄分布 → 重复参与繁殖
              有归因于滞温效应的 ，有归因于产卵场底质退                            洄游的时间间隔分布，然后将该控制节点进行串
                                 [14]
                   [4]
              化的 ，还有其他等方面            [7, 15-16] 。对于葛洲坝救鱼       联建模   (图  1)。用  1980  年之前中华鲟在长江中的
              行动，有认为相关评估和建议是科学合理的                       [17] ，  已知参数率定其在海洋中的未知参数，并将所率
              有指出相关评估是存在严重误判的                 [16, 18] ，有提出    定的参数代入模型以应用于不同场景下中华鲟自
              增殖放流无效的 。似乎每个观点都指出了一部                            然种群演变的模型计算。
                             [11]
              分原因，但好像还都不足以充分解释当前困境，                                该简化参数模型是由各控制节点串联起来的
              更未能有效指导中华鲟拯救行动实践 。本研究                            循环迭代模型。在单循环计算中，除了                   2  个限制
                                                 [7]
              提出假设，在葛洲坝截流之后，中华鲟自然种群                            性指标是作为相应节点的判断项                (分段函数)，放
              中大多数个体可能已无法适应新的生境条件，其                            流指标作为相应节点的补充项               (加法)，其他指标
              自然生活史周期已无法支撑种群自持，其自然种                            均为相应节点的串联项           (乘法)。简化参数模型中
              群已开始陷入螺旋式衰退路径                (用一个形象生动            与葛洲坝截流相关的指标有繁殖洄游群体性腺发
              的比喻来描述，即陷入死亡螺旋)，后续的人类活                           育成熟比例、产卵窗口适配率、产卵场繁殖参与
              动干扰则进一步加速了其自然种群螺旋式衰退的                            率、产卵场容量，与三峡运行相关的指标有产卵
              速度。                                              窗口适配率，与放流相关的指标有人工放流规模、
                   为了检验中华鲟自然种群陷入死亡螺旋(陷入                        放流个体适应生存率，与洄游相关的指标有洄游
              螺旋式衰退路径) 的假设，并探讨中华鲟自然种群                          过程中繁殖群体的生存率、仔稚幼鱼降河生存率、
              改出死亡螺旋       (从螺旋式衰退路径中改出) 的必要                   长江口适应调整期生存率，与海洋生活史相关的
              条件，本研究基于中华鲟生活史和相关本底数据                            指标有海洋群体生存率的年龄分布、初次性腺发
              对中华鲟自然种群历史演变进行了简单的量化建                            育启动洄游的年龄分布、重复参与繁殖洄游的时
              模，通过模拟计算中华鲟自然种群在葛洲坝截流                            间间隔分布。将基于中华鲟生活史的简化参数模
              以来数十年的变化检验模型的可用性，然后模拟                            型表述为计算式：
              计算不同场景条件下中华鲟自然种群在过去和未                                                r (FM)
                                                                   N (mat)(F) = N (P) × [  ] ×r (sur−up) ×r (mat)
              来数十年的状况，进而检验中华鲟自然种群陷入                                              1+r (FM)
              死亡螺旋的假设和探讨不同拯救方案下中华鲟自                                                  1
                                                                   N (mat)(M) = N (P) × [  ] ×r (sur−up) ×r (mat)
              然种群改出死亡螺旋的可能性。                                                      1+r (FM)

                                                                     N (mat )(F)
                              
                                              N (mat )(F) ×r (adap ) ×r (partici ) ,
                                                                            ⩽ r (thres )(FM)( upper )
                              
                              
                                                                    N (mat )(M)
                              
                  N (poten )(B)(F) = 
                              
                                                                     N (mat )(F)
                              
                              N (mat )(M) ×r (thres )(FM)( upper ) ×r (adap ) ×r (partici ) ,  > r (thres )(FM)( upper )
                              
                              
                                                                     N (mat )(M)
              https://www.china-fishery.cn                           中国水产学会主办    sponsored by China Society of Fisheries
                                                            2
   55   56   57   58   59   60   61   62   63   64   65