Kemp's ridley 海龟体型波动与 DWH 溢油事件后群体结构的关联研究：一项基于 Rancho Nuevo 庇护所 2018-2023 年巢居雌龟直线背甲长度分析

《Agronomy》：Towards Stress-Resilient Canola via Genetic Engineering Approaches Ali Ijaz Ahmed, Aldrin Y. Cantila and Sheng Chen

本研究聚焦濒危的 Kemp's ridley 海龟种群在 2010 年“深水地平线”（Deepwater Horizon, DWH）溢油事件后的恢复状态。为解决巢数监测数据中断后群体结构信息缺失的难题，研究人员在墨西哥 Rancho Nuevo 庇护所开展了为期六年（2018-2023）的巢居雌龟直线背甲长度（SCL）监测。研究发现，雌龟平均 SCL 存在显著的年际波动，而非持续下降，其中 2020 和 2022 年的低值可能与 DWH 事件后成熟的初次筑巢个体（neophytes）有关。该成果为理解大型环境灾难对长寿海洋脊椎动物的长期亚致死效应提供了关键表型证据，并强调了在传统丰度数据缺失时，形态测量学可作为监测种群结构与补充动态的互补性重要工具。

　　

在墨西哥塔毛利帕斯州（Tamaulipas）绵延的海岸线上，有一处名为 Rancho Nuevo 的圣地，这里是世界上体型最小的海龟——坎普氏雷德利海龟（Kemp's ridley sea turtle, Lepidochelys kempii）最重要的筑巢家园。这种海龟的命运堪称一部跌宕起伏的现代保护史诗。它们曾因过度捕猎和渔业副渔获物而在 20 世纪 80 年代濒临灭绝，随后在严格的巢穴保护和强制性使用海龟排除装置（Turtle Excluder Devices, TEDs）等措施下，种群数量在 2010 年前呈现出指数级增长的复苏迹象。然而，2010 年发生的“深水地平线”（Deepwater Horizon, DWH）原油泄漏事件，给这片海域的生态系统带来了沉重一击，也使得海龟的复苏进程骤然停滞。此后，种群增长陷入停滞和波动，传统的监测指标——年度巢穴计数——虽然基础，但似乎已不足以揭示种群内部可能发生的深刻变化。

科学家们意识到，要理解一个种群的健康状况，不仅要看“有多少”，还要看它们是“谁”。对于海龟这类长寿的脊椎动物，体型大小是一个绝佳的生命史与环境经历的“记录仪”。雌龟的体型，特别是其直线背甲长度（Straight Carapace Length, SCL），与繁殖力（通常体型越大，产卵量越高）、性成熟年龄等关键种群参数紧密相关。更重要的是，筑巢雌龟群体体型分布的变化，可以揭示种群结构的动态，比如每年新加入的初次筑巢者（neophytes，通常体型较小）与有经验的回归筑巢者（remigrants，通常体型较大）的比例变化。这就像通过观察一支军队中新兵和老兵的比例，来判断其兵员补充和战斗力的变化。那么，在经历了 DWH 这样的重大环境灾难后，那些在泄漏后环境中成长的年轻海龟，其体型和成熟年龄是否受到了影响？种群的体型结构是否发生了系统性改变？这些问题对于评估 DWH 事件的长期生态影响、理解当前种群恢复停滞的原因至关重要。

为此，研究人员在 2018 年至 2023 年的筑巢季节，在 Rancho Nuevo 庇护所的核心筑巢区开展了一项系统性研究。他们选择在“arribada”（一种成百上千只雌龟同步筑巢的现象）高峰期，对正在产卵的雌龟进行采样测量。通过使用精度达 ±0.1 厘米的不锈钢游标卡尺，严格遵循前人确立的测量协议，科研人员收集了共计 191 只筑巢雌龟的 SCL 数据。为了确保研究结果的统计稳健性，他们不仅计算了描述性统计量，还深入检验了数据的分布特性（如偏度）和方差-均值关系（应用泰勒幂定律，Taylor's Power Law），以验证使用参数检验（如方差分析，ANOVA）的合理性。在确认数据满足近似正态分布和方差齐性（homoscedasticity）假设后，他们才进行后续的差异分析和趋势检验。

研究结果揭示了几个关键发现：

• •
  验证统计假设：对历年 SCL 数据的频率分布分析显示，所有年份的偏度系数（G₁）绝对值均 ≤ 0.5，表明分布近似对称。同时，泰勒幂定律分析得到的指数（b = -1.08）不显著，且决定系数（R²= 0.021）可忽略不计，证实了方差独立于均值，即数据具有方差齐性。这为后续使用参数检验和直接比较算术平均值奠定了坚实的统计基础。
• •
  平均 SCL 的年际变化：单因素方差分析（one-way ANOVA）显示，不同年份间的平均 SCL 存在高度显著差异（p < 0.001）。其中，2020 年（59.01 ± 1.79 厘米）和 2022 年（60.68 ± 1.47 厘米）的平均 SCL 显著低于其他年份（2018、2019、2021 和 2023 年，其均值在 61.75 至 62.77 厘米之间）。这种模式在展示数据分布和统计分组的箱线图与小提琴图中得到清晰体现，直观地展示了 2020 和 2022 年体型分布的左移。
• •
  检验持续线性趋势：简单线性回归分析表明，从 2018 年到 2023 年，平均 SCL 并未呈现显著的线性下降趋势（p = 0.78）。这一结果否定了“平均 SCL 在此期间持续下降”的先验假设，支持了“波动而非趋势”的解释。

结论与讨论综合来看，本研究通过严谨的统计分析，证实了在 2018 至 2023 年间，Rancho Nuevo 筑巢的 Kemp's ridley 海龟雌龟体型存在显著的年际波动，但并未形成持续的、单向的体型变小趋势。这种波动模式最合理的生物学解释是每年筑巢群体中初次筑巢者（体型较小）与回归筑巢者（体型较大）比例的变化。尤为重要的是，2020 年和 2022 年观测到的较低平均 SCL，在时间上与 DWH 溢油事件后首批成熟并加入筑巢群体的海龟（大约在 2008-2014 年间孵化）相吻合。这表明，DWH 事件可能通过多种亚致死机制（如慢性生理应激、觅食地退化、种内竞争加剧等）影响了相关群体在幼年和亚成年阶段的生长发育轨迹，导致它们以较小的体型达到性成熟。

这项研究的意义深远。首先，它将形态测量学（morphometrics）提升为一个在传统丰度监测数据（如巢数统计）中断或受限时，监测种群结构和补充动态的、可靠且必不可少的工具。其次，研究结果将 DWH 溢油事件对海龟种群的影响，从一个潜在的死亡率事件，延伸到了可量化的、长期的表型（体型）效应层面，为评估大型环境灾难对长寿物种的生态遗产提供了具体证据。最后，研究呼吁在未来加强综合性的监测体系，包括建立个体标记重捕（mark-recapture）项目以区分初次与再次筑巢者，开展将幼体发育环境与成体表型联系起来的整合研究，并迫切需要恢复跨国、跨机构的协作数据共享框架。尽管面临监测挑战和环境压力，但持续、科学的保护努力和像本研究所展示的这类严谨的科学研究，仍是引导这种极度濒危的海龟物种在充满变化和不确定性的墨西哥湾中继续生存和恢复的希望所在。