种群生存力分析(Population viability analysis, PVA)是以物种生活史信息为参数，运用计算机模型来估计濒危物种种群未来大小和灭绝风险的一种方法，多用于模拟孤立小种群的生存动态过程，可以预测种群数量变化和灭绝概率(Brook et al, 2000a、b)。应用这一理论，发展了许多计算机应用程序，包括RAMAS(Ferson et al, 1989)、NEMESIS(Gilpin, 1993)、ALEX(Possingham et al, 1992)和Vortex (Lacy et al, 1993)等，其中，Vortex模型(漩涡模型)应用最为广泛。国内学者已经利用该模型对白鳍豚(Lipotes vexillifer) (张先锋等, 1994)、长江江豚(Neophocaena asiaeorientalis asiaeorientalis) (张先锋等, 1999)、普氏野马(Equus przewalskii) (裴鹏祖等, 2018)、川金丝猴(Rhinopithecus roxellana) (张宇等, 2018)、猕猴(Macaca mulatta)(楚原梦冉等, 2019)等种群动态进行模拟和研究。张先锋等(1999)采用漩涡模型(Vortex 7.3)模拟长江江豚的种群生存力，预测了长江江豚未来的种群动态。20年来，长江江豚种群结构发生了新变化，种群生态学研究也取得了新成果(杨健等, 2000; 于道平等, 2005; 张新桥, 2011; Kimura et al, 2012; Mei et al, 2014、2017; Zhou et al, 2018; Yuan et al, 2018)，有必要对长江江豚种群生存状况进行再评估和再认识。为此，本研究利用漩涡模型(Vortex 10.3.7.0)对长江江豚种群生存力及其致危因素进行分析，以期为未来种群管理和有效保护提供有益参考。

长时间以来，长江和海洋江豚的分类关系一直存在争议，曾被认为是同一个物种(3个亚种)，后被认定为2个物种：以台湾海峡为界，以南海域是印太江豚；以北海域包括长江流域为窄脊江豚(Neophocaena asiaeorientalis)。2种江豚形态上区分的显著特点是其脊背宽窄和疣粒数量多少。长江江豚是江豚的最濒危的一个亚种，仅生活在长江中下游干流和大型通江湖泊—鄱阳湖和洞庭湖中。但目前基于江豚种群基因组学最新研究结果显示，极度濒危的长江江豚种群应被作为一个独立的物种(Zhou et al, 2018)。长江江豚不仅是全世界唯一生存于淡水环境的鼠海豚科(Phocaenidae)动物，也是长江流域极可能唯一现存的的水生哺乳动物(2007年白鱀豚功能性灭绝)。自上世纪九十年代以来，长江江豚的数量呈现出持续性的衰减趋势，2018年7月正式发布《2017长江江豚生态科学考察报告》，这也是目前最新的长江江豚科考成果。江豚种群数量大幅下降的趋势得到遏制，估算长江江豚数量约为1012头，其中，鄱阳湖约为457头，是野外江豚分布密度最高的区域。长江江豚发生了一些新变化，种群趋于向受人为扰动较少的江段集中，过半数水域种群持续下降或没有发现种群，种群分布呈碎片化特征(先义杰, 2018)。因其受威胁的严峻程度，长江江豚已被按照国家一级保护动物进行保护。2018年，中国科学院深海科学与工程研究所海洋哺乳动物研究团队对长江江豚进行了全基因组测序，组装得到了2.22 Gb的基因组草图，对江豚基因组的研究奠定了基础(Yuan et al, 2018)。

本研究以1999年作为回顾分析的起点，2017年作为模拟验证时间节点和现状分析的起点，进行长江江豚种群生存力分析。为进一步提升精度和可信度，将迭代次数设定为10000，为前人同类研究的10倍。MVP定义为以95%的概率健康存活100年所需的最小种群数量。估算方法为：分析不同初始种群数量(Ni)在100年内的种群生存力，当种群存活概率达到95%时的Ni即为MVP(Lacy et al, 2009、2017)。

长江江豚与海江豚之间存在着显著而稳定的遗传分化，提示它们之间缺乏基因交流而出现了生殖隔离，长江江豚因此被认定为独立物种，长江江豚和海江豚最明显的差别是一个生活在淡水，另一个生活在咸水。通过大数据分析，发现长江江豚和海江豚在进化上显著不同，特别是基因上体现对水环境的适应性差异，表明了长江江豚和其他江豚之间，已经不可能出现基因交流和种群间混合，因此，不存在迁移扩散问题(Zhou et al, 2018)。

长江江豚为混交制，雌性可能先于雄性性成熟，雌性首次生育为4~6龄，雄性为4.5~7龄(郝玉江等, 2006)。Shirakihara等(1993)发现，2头21龄和23龄的个体可以怀孕和哺乳。长江江豚无论野外还是人工环境下均无产双胞胎记录，因此，在模拟中设定每胎最大产仔数为1。长江江豚新生婴豚的性比为1∶1，雌性首次生育年龄为4龄，雄性为5龄，最高繁殖年龄约为20龄。该物种的种群增长应为密度制约型，即参加繁殖的雌性数占全部成年雌性数的比例[P(N)]随种群大小(N)的变化而变化。Fowler(1981)认为，描述这种变化强度的参数B取值2能更好地模拟密度制约型种群的增长，描述雌性交配率降低的参数A取值也为2。假设N接近环境容纳量K时，繁殖雌性的比例P(K)为25%；当N接近0时，P(0)为70%(张先锋等, 1999)。

张先锋等(1992)研究指出，长江江豚繁殖周期一般为2年，年生殖率为20%。2020年1月公开报道显示，湖北石首天鹅洲迁地保护区内，江豚种群5年内净增长108%。因此，推测对仍处于种群数量下降趋势的长江江豚来说，目前可能的繁殖率只相当于其理论值的一半左右，即成年雌豚生殖率可能为20%~25%。目前，关于长江江豚流产率和幼豚夭折率及流产或幼豚夭折等对其繁殖周期的影响还不全面，进一步深入研究有助于提升相关参数估算的精度。

张先锋等(1999)研究发现，长江江豚0~1和1~2年龄组的死亡率为20%，其他年龄组的死亡率为15%。20世纪90年代以来，长江江豚种群下降速率每年约为6.3%，并推测幼体死亡率较高可能是决定其种群动态的重要因素(梅志刚等, 2011)。在模拟中设定0~1年龄组的死亡率为20%~30%，1~2年龄组的死亡率为20%，其他年龄组为15%。

张先锋等(1999)设置2种灾害类型，发生频率为10%，以死亡率均值为存活率下降率(95%)，用成年雌性死亡率表示繁殖率下降率(95%)。2016年，习近平总书记在重庆召开推动长江经济带发展座谈会，提出要把修复长江生态环境摆在压倒性位置，“共抓大保护、不搞大开发”。农业农村部发布《关于长江流域重点水域禁捕范围和时间的通告》，通告指出，长江干流和重要支流除水生生物自然保护区和水产种质资源保护区以外的天然水域，最晚自2021年1月1日0时起实行暂定为期10年的常年禁捕，期间禁止天然渔业资源的生产性捕捞。人类不利于长江江豚保护的活动减少，但像2008年的大雪灾害、2019年的鄱阳湖大干旱极低水位以及意外事故等灾害仍然存在。综合病害等生理风险和意外等生态风险2种灾害类型，同时，长江江豚种群存在碎片化分布趋势，影响程度更大。为此，以1999年作为回顾分析的起点，2017年为模拟验证时间节点，设置灾害2种类型，发生频率为10%，以死亡率均值为存活率下降率(95%)，成年雌性死亡率表示繁殖率下降率(95%)。以2017年作为现状分析的起点，设置灾害2种类型，发生频率均为10%，以死亡率均值为存活率下降率(90%)，成年雌性死亡率表示繁殖率下降率(90%)。

张先锋等(1999)设置长江江豚的初始头数为2000头，环境容纳量为5000头，每5年环境容纳量将减少10%。为此，以1999年作为回顾分析的起点，2017年为模拟验证时间节点，也设置初始头数为2000头，环境容纳量为5000头。环境容纳量每年减少2%~3%，连续减少15年。以2017年作为现状分析的起点，设置初始头数为1012头，环境容纳量为4000头，每年增加1%，连续变动10年。

对于小种群而言，近交衰退是影响其长期存活的重要因素之一。Vortex软件通过模拟有无近亲繁殖时的种群动态来说明其对种群的影响。Ralls等(1988)研究了40个哺乳动物种群致死等价系数，认为每个二倍体有3.14个致死基因当量。因此，设置3.14作为致死等价系数。

目前，没有对长江江豚的捕获和补充，故该参数设置为0。

通过Vortex模拟计算可知，以1999年作为回顾分析的起点，8种不同情景下具体结果见表 1。长江江豚种群未来100年的内禀增长率r为–0.0622~ –0.0231，周限增长率λ为0.9397~0.9701，净生殖率R_o为0.5747~0.8211，平均世代时间雌性为8.55~8.90年，雄性为9.41~9.73年，平均现存种群数量为8.20± 6.08~184.16±139.75头，平均种群数量为2.03±4.36~ 183.92±139.81头，遗传多样性为0.6990±0.1455~ 0.9585±0.0301。灭绝概率为0.0013~0.7828，平均首次灭绝时间为84.4~95.4年，具体见表 2。以2017年长江江豚数量约为1012头作为模拟验证来看，情景2和情景4与实际情况不吻合，为小概率情景。情景3模拟结果与实际调查情况契合度最好，情景1次之，情景5也较理想。情景5和情景7种群存活概率达到95%以上。

表 1(Tab.1)

表 1 长江江豚种群生存力回顾分析的不同模拟情景 Tab.1 Different simulation scenarios of the retrospective analysis of the viability of Yangtze finless porpoise populations 项目Items 情景1 Scenario 1 情景2 Scenario 2 情景3 Scenario 3 情景4 Scenario 4 情景5 Scenario 5 情景6 Scenario 6 情景7 Scenario 7 情景8 Scenario 8 环境容纳量变动年限Changing years of environmental capacity 15 15 15 15 15 15 15 15 环境容纳量年变动大小Annual variation of environmental capacity (%) –2 –2 –3 –3 –3 –3 –3 –3 0~1龄组死亡率 Mortality in the 0~1 age group (%) 30 30 30 30 25 25 20 20 成年雌豚生殖率Breeding rate (%) 25 20 25 20 25 20 25 20 初始头数Initial population size, Ni 2000 2000 2000 2000 2000 2000 2000 2000 初始年份Initial year 1999 1999 1999 1999 1999 1999 1999 1999

表 1 长江江豚种群生存力回顾分析的不同模拟情景 Tab.1 Different simulation scenarios of the retrospective analysis of the viability of Yangtze finless porpoise populations

表 2(Tab.2)

表 2 长江江豚种群生存力回顾分析结果 Tab.2 Results of the retrospective analysis on the survival ability of Yangtze finless porpoise populations 项目Items 情景1 Scenario 1 情景2 Scenario 2 情景3 Scenario 3 情景4 Scenario 4 情景5 Scenario 5 情景6 Scenario 6 情景7 Scenario 7 情景8 Scenario 8 r –0.0381 –0.0622 –0.0381 –0.0622 –0.0304 –0.0549 –0.0231 –0.0479 λ 0.9626 0.9397 0.9626 0.9397 0.9701 0.9466 0.9772 0.9532 Ro 0.7184 0.5747 0.7184 0.5747 0.7697 0.6158 0.8211 0.6568 平均世代雌性 Average generation female 8.68 8.90 8.68 8.90 8.61 8.84 8.55 8.77 平均世代雄性 Average generation male 9.53 9.73 9.53 9.73 9.47 9.67 9.41 9.61 N-extant±SD 40.30 ±35.44 8.20 ±6.08 40.57 ±34.86 8.22 ±6.51 86.84 ±70.77 11.74 ±9.73 184.16 ±139.75 17.72 ±15.71 N-all±SD 37.64 ±35.60 2.03 ±4.36 37.60 ±35.09 2.10 ±4.57 85.75 ±70.96 5.42 ±8.53 183.92 ±139.81 12.23 ±15.19 PE 0.0683 0.7828 0.0754 0.7746 0.0128 0.5629 0.0013 0.3232 GeneDiv±SD 0.8801 ±0.0829 0.6990 ±0.1455 0.8822 ±0.0803 0.6990 ±0.1411 0.9290 ±0.0533 0.7533 ±0.1282 0.9585 ±0.0301 0.8028 ±0.1132 平均TE Mean TE 93.3 84.4 93.5 84.5 94.4 87.9 95.4 90.5 2017年预测头数±SD Predicted value of 2017±SD 1021.36 ±260.19 662.07 ±178.28 1015.90 ±262.61 661.71 ±181.40 1167.91 ±304.61 756.99 ±206.54 1329.71 ±343.37 854.52 ±230.54 注：r: 内禀增长率; λ: 周限增长率; Ro: 净生殖率; N-extant: 平均现存种群数量; N-all: 平均种群数量; PE: 灭绝概率; GeneDiv: 遗传多样性; TE: 灭绝时间(年)。表 5同 Note: r: Annual internal growth rate; λ: Finite rate of increase; Ro: Net reproduction rate; N-extant: Average number of existing populations; N-all: Average population number; PE: Extinction probability; GeneDiv: Genetic diversity; TE: Extinction time (years). Same in the Tab. 5

表 2 长江江豚种群生存力回顾分析结果 Tab.2 Results of the retrospective analysis on the survival ability of Yangtze finless porpoise populations

情景5的模拟结果与实际调查情况契合度较理想，种群灭绝概率为0.0128。因此，情景5可以作为阶段性的科学研究及精准管护的目标。为此，求取情景5条件下的最小可存活种群(Minimum viable population, MVP)，初始头数为1000开始模拟，按照需要每次增减5头，不难发现，情景5时，MVP为1305头，具体结果见表 3。同理，采取渐进逼近的方法发现，情景5条件下，满足MVP为1305头的最小环境容量为2300头。

表 3(Tab.3)

表 3 特定情景下长江江豚种群MVP模拟分析 Tab.3 Simulation analysis of MVP of Yangtze finless porpoise population under specific scenarios 初始头数Initial population size, Ni 1325 1320 1315 1310 1305 PE 0.0460 0.0484 0.0476 0.0540 0.0549

表 3 特定情景下长江江豚种群MVP模拟分析 Tab.3 Simulation analysis of MVP of Yangtze finless porpoise population under specific scenarios

首先，在回顾分析中情景3基础上，以2017年作为现状分析的起点，设置初始头数为1012头，环境容纳量为5000头和4000头，发现灭绝概率分别为0.7012和0.6950。在情景5的基础上，设定环境容纳量为4000头，环境容纳量分别每年变动–1%、0%和1%，连续变动10年，发现灭绝概率分别为0.4635、0.4592和0.4588。随后，环境容纳量设定为4000头，每年增加1%，连续变动10年，以2017年作为起点进行不同模拟情景的生存动态分析，具体见表 4。长江江豚种群100年的内禀增长率r为–0.0733~ –0.0207，周限增长率λ为0.9293~0.9795，净增值率R_o为0.5235~0.8413，平均世代雌性为8.36~8.83年，雄性为9.23~9.66年，平均现存种群数量5.49±3.58~ 110.81±108.10头，平均种群数量为0.11±0.80~ 107.91±108.10头，遗传多样性为0.5801±0.1671~ 0.9169±0.0659。灭绝概率为0.0264~0.9860，平均首次灭绝时间为68.4~91.9年，具体见表 5。

表 4(Tab.4)

表 4 长江江豚种群生存力现状分析的不同模拟情景 Tab.4 Different simulation scenarios for the analysis of the current viability of Yangtze finless porpoise population 项目Items 情景1 Scenario 1 情景2 Scenario 2 情景3 Scenario 3 情景4 Scenario 4 情景5 Scenario 5 情景6 Scenario 6 情景7 Scenario 7 情景8 Scenario 8 0~1龄组死亡率 Mortality in the 0~1 age group (%) 25 25 25 25 30 30 30 30 成年雌豚生殖率 Breeding rate of adult female (%) 20 25 27.5 30 20 25 27.5 30 初始头数Initial population size, Ni 1012 1012 1012 1012 1012 1012 1012 1012 初始年份Initial year 2017 2017 2017 2017 2017 2017 2017 2017

表 4 长江江豚种群生存力现状分析的不同模拟情景 Tab.4 Different simulation scenarios for the analysis of the current viability of Yangtze finless porpoise population

表 5(Tab.5)

表 5 长江江豚种群生存力现状分析结果 Tab.5 Analysis results of current viability of finless porpoise population in the Yangtze River 项目Items 情景1 Scenario 1 情景2 Scenario 2 情景3 Scenario 3 情景4 Scenario 4 情景5 Scenario 5 情景6 Scenario 6 情景7 Scenario 7 情景8 Scenario 8 r –0.0660 –0.0416 –0.0308 –0.0207 –0.0733 –0.0493 –0.0386 –0.0287 λ 0.9361 0.9593 0.9697 0.9795 0.9293 0.9519 0.9621 0.9717 Ro 0.5609 0.7011 0.7712 0.8413 0.5235 0.6543 0.7198 0.7852 平均世代雌性 Average generation female 8.76 8.54 8.44 8.36 8.83 8.61 8.51 8.43 平均世代雄性 Average generation male 9.60 9.40 9.31 9.23 9.66 9.46 9.37 9.30 N-extant±SD 6.32 ±4.47 18.18 ±17.56 42.13 ±42.94 110.81 ±108.10 5.49 ±3.58 11.67 ±10.17 22.16 ±22.07 51.10 ±52.51 N-all±SD 0.34 ±1.58 10.08 ±15.65 35.83 ±42.24 107.91 ±108.10 0.11 ±0.80 3.75 ±7.66 14.50 ±20.52 45.58 ±51.98 PE 0.9589 0.4588 0.1523 0.0264 0.9860 0.6985 0.3551 0.1098 GeneDiv±SD 0.6113 ±0.1620 0.7745 ±0.1303 0.8532 ±0.1026 0.9169 ±0.0659 0.5801 ±0.1671 0.7191 ±0.1483 0.7976 ±0.1225 0.8677 ±0.0981 平均TE Mean TE 73.4 87.1 89.8 91.9 68.4 83.5 88.0 90.2

表 5 长江江豚种群生存力现状分析结果 Tab.5 Analysis results of current viability of finless porpoise population in the Yangtze River

在现状分析情景6的基础，设定性比雄性化、雌性化偏离10%，即雄性所占比例分别为0.55、0.50和0.45，模拟计算得出，长江江豚种群100年的内禀增长率r分别为–0.0608、–0.0493和–0.0386，周限增长率λ分别为0.9410、0.9519和0.9621，净增值率R_o分别为0.5889、0.6543和0.7198，平均世代雌性分别为8.71、8.61和8.51年，雄性分别为9.56、9.46和9.37年，平均现存种群数量分别为7.30±5.69、11.67±10.17和22.42±22.58头，平均种群数量分别为0.76±2.66、3.75±7.66和14.42±20.39头，遗传多样性分别为0.6490±0.1597、0.7191±0.1483和0.7950± 0.1237。灭绝概率为0.9157、0.6985和0.3664，平均首次灭绝时间为76.5、83.5和88.0年。

本研究以1999年作为回顾分析的起点，2017年为模拟验证时间节点，在综合分析历史文献资料以及最新研究成果的基础上，围绕环境容纳量、0~1龄组死亡率、成年雌豚生殖率设定了回顾分析8种不同情景，长江江豚种群100年灭绝概率为0.0013~0.7828，平均灭绝时间为84.4~95.4年。以2017年长江江豚数量约为1012头作为模拟验证来看，情景2和4与实际情况不吻合，情景3模拟结果与实际调查情况吻合度最好，情景1次之，情景5也较理想。情景5和7的种群存活概率达到95%以上。提示上述情景3中模拟参数与历史生境及长江江豚种群特性的契合度较好，在该情景下进行深入系统分析具有较强的理论和实际价值。通过Vortex模拟计算可知，环境容纳量为5000头和4000头，发现灭绝概率分别为0.7012和0.6950。进一步设定环境容纳量为4000头，环境容纳量分别每年变动–1%、0%和1%，连续变动10年，发现灭绝概率分别为0.4635、0.4592和0.4588，提示环境容纳量小范围变动对长江江豚种群生存力影响不大。从回顾分析结果来看，长江江豚的灭绝概率最可能为0.0754，且不高于0.5629。平均灭绝时间最可能为93.5年，且不少于87.9年。0~1龄组死亡率最可能为30%左右，成年雌豚生殖率最可能为25%左右。杨光等(1998)研究指出，长江江豚0~2龄组死亡率为0.308，繁殖率可高达0.7。张先锋等(1999)研究发现，长江江豚0~1龄组的死亡率为20%，年生殖率为20%。因此，模拟计算确认的长江江豚0~1龄组死亡率最可能为30%左右，成年雌豚生殖率最可能为25%左右，在合理范围之内。通过进一步分析，不难发现，成年雌豚生殖率和0~1龄组死亡率与种群生存力关系密切，模拟计算出，成年雌豚生殖率为25%时，长江江豚0~1龄组死亡率由30%分别下降为25%和20%时，灭绝概率由0.0754下降为0.0128和0.0013。0~1龄组死亡率为30%时，成年雌豚生殖率由25%下降为20%时，灭绝概率由0.0754增高为0.7746。因此，长江江豚种群所面临灭绝“漩涡”，直接原因应该是“一低一高”，即成年雌豚生殖率低，而0~1龄组死亡率高，但除了长江江豚自身生理生态以及进化演化适应方面的原因，长江江豚的生境干扰、破坏，特别是适合母子豚栖息、抚幼以及幼豚生长发育生境的数量减少和质量下降，导致其种群被动碎片化分布，个体水平可能有利于减小生存压力。但从种群水平来说，繁殖、遗传方面等风险加剧，“一低一高”现象持续恶化。另一方面，饵料鱼类资源下降，特别是渔汛等减少，捕食成本大幅提高，而长江江豚作为淡水水生哺乳动物，能量需求相对较大，特别是求偶、怀孕和哺乳等特定生理阶段，饵料相对不足，必将严重影响繁殖活动，并对幼仔生存带来巨大挑战，最终可能让长江江豚陷入种群灭绝“漩涡”。在回顾分析中，求取了情景5，即最可能达到情景下的MVP，为1305头，且满足该条件的最小环境容纳量为2300头。提示在进行长江江豚保护时，长江干流、鄱阳湖以及洞庭湖等水域在加强重点水域重点保护的同时，必须强化长江中下游水域大范围的“互联互通”，单一特定水域难以承载长江江豚的持续生存繁衍。

本研究中，从0~1龄组死亡率和成年雌豚生殖率的角度设定了8种不同模拟情景，环境容纳量设定为4000头，每年增加1%，连续变动10年，以2017年作为现状分析的起点，从上述回顾分析的结果来看，现状分析中的情景6和情景7是最可能情景，即长江江豚种群灭绝概率可能为0.3551~0.6985，平均灭绝时间可能为83.5~88.0年。而情景1和情景5是小概率情景，情景2和情景3是可能性较低情景，情景4和情景8是相对理想情景。长江江豚新生幼豚性比(雄性所占比例)对模拟结果具有较大影响，性比雄性化、雌性化偏离10%，灭绝概率由0.6985分别增加为0.9157和下降为0.3664。提示长江江豚新生幼豚性比研究对于长江江豚保护具有特别意义，通过比较相关数据，现状分析中情景6新生幼豚性比雌性化10%，与情景7模拟结果相当，即在情景6时，新生幼豚性比雌性化10%与成年雌豚生殖率提高10%，效果是相当的，提示人为诱导甚至控制长江江豚新生幼豚性比可能是进行长江江豚抢救性保护的措施之一。当种群面临灭绝漩涡风险时，重点保护、修复母子豚栖息区等特定自然生境，缓解甚至逆转长江江豚成年雌豚生殖率低，而0~1龄组死亡率高的现状，有利于长江江豚的保护，当然这需要对长江江豚生理生态学、行为生态学、繁殖生态学等进行深入而系统的研究。需要特别指出的是，我国实施伏季休渔制度已有20多年，在缓解捕捞对渔业资源带来的巨大压力和资源养护方面发挥了积极作用(胡芷君等, 2019)，长江十年禁捕等战略性政策的落地，必将有效增加长江江豚饵料资源，对其保护意义重大。本文模型使用的参数是由系统整理前人长江江豚研究结果获得的，特别是繁殖率和最高繁殖年龄的设定是整体权衡的结果，单一参数可能存在较大误差，但参数组合误差相对较小，且本研究首次运用迭代1万次的方法进行相关模拟和估计，迭代次数为前人研究的10倍，一定程度上减小了系统误差，能够为长江江豚种群生存力分析及相关濒危机制探讨，特别是长江江豚种群的保护提供一定借鉴和管护依据。同时，有报道使用灰色预测模型，分析影响资源丰度变化因子，并建立多种资源丰度预报模型，以此反映该种类的资源丰度动态，应用效果较好，提示多种模型协同应用可能是未来种群生存力分析的发展方向之一(闫永斌等, 2019)。