干扰持续时间主导死亡风险，而累积能量赤字塑造地中海贻贝的繁殖表现

《Marine Environmental Research》：Disturbance duration governs mortality risk while cumulative energetic deficit shapes reproduction in Mediterranean mussels

　　

摘要专业翻译 海洋贻贝是基础物种和生态系统工程师，日益暴露于海洋热浪、低氧事件和食物限制之下，然而表面相似的事件可在不同种群中产生截然不同的后果。本研究检验这种差异是否可由干扰持续时间、累积能量赤字与内部能量状态之间的相互作用来解释。研究人员利用地中海本地种紫贻贝（Mytilus galloprovincialis）和入侵种Brachidontes pharaonis的动态能量预算（Dynamic Energy Budget，DEB）模型，在四种环境强迫情景下，将繁殖表现分解为达到性成熟（maturity）的概率和成熟条件下的繁殖投入（reproductive investment）两个组分。在相等的累积能量赤字下，摄食抑制的不间断持续时间（uninterrupted duration）主导死亡风险和总体繁殖表现，而累积赤字决定存活个体的繁殖投入量。存活率在一个狭窄的持续时间范围内呈非线性下降，该范围的位置由未受干扰条件下维持的储备余量（reserve margin）排序。脆弱性因此取决于事件持续时间与能量状态的耦合，而后者随局部食物条件、季节和生命阶段而变化。本地种更接近其饥饿边界（starvation boundary），始终比维持更大储备缓冲的入侵种更脆弱。这些结果将储备余量确定为补充失败（recruitment failure）和死亡风险的一个可测量的候选预测因子。更广泛地，它们揭示了为何仅凭累积暴露指标可能遗漏驱动存活损失的时间连续性。这一个体层面的机制为以下现象提供了可检验的解释：表面相似的沿海干扰为何能产生不均衡的补充失败，并促成贻贝种群中空间异质的大规模死亡。

论文解读文章

一、研究背景与问题

沿海海洋系统的干扰格局通常以事件的频率、强度和持续时间来描述，学界已有大量理论组织生态系统对它们的响应。持续型扰动（press）与脉冲型扰动（pulse）的区分构成了这一工作的核心框架：持续型扰动是持续存在的压力，生物体生活于其中；脉冲型扰动是事件，有开始有结束，生物体承受其所失去的。然而，这一文献主要在群落或结构化种群层面运作，将干扰事件表示为人口统计率的瞬时变化。给定的累积干扰量是以少数长事件还是许多短事件传递，很少被解析；即便被涉及，考察的也是频率与强度之间的权衡，而非单个事件的持续时间。但一个压力事件的持续时间，相对于受影响生物体的生理时间尺度，从能量角度预期是决定性的。野外观察也指向同一方向：2003年和2022年地中海热浪相关的大规模死亡事件中，死亡率与异常热条件持续的时间相关，而不仅与其强度相关；加利福尼亚洋流热浪的响应同样跨类群依赖于事件持续时间。

贻贝是温带和地中海沿岸的基础物种，构建生物栖息地，过滤并耦合水柱与底栖环境，支撑依赖其床体形成的生物多样性。其种群的持续取决于两个人口统计枢纽：个体能否存活至繁殖以更新种群，以及存活者投入配子的数量以设定繁殖能力。本地种与入侵种在这些枢纽对日益加剧的干扰的响应上存在差异时，干扰格局本身便成为群落重组和变暖下入侵的动因。预测此类转变需要机制性的、性状层面的解释，说明哪些种群最脆弱及其原因，而非局限于单一地点或终点的描述。

生物体通过动用储备来应对能量供给的短缺；能否在中断中存活以及之后能否繁殖，取决于储备在生理极限被跨越前能支撑赤字多久。动态能量预算理论提供了框架，描述同化能量如何储存于储备中，并在固定分配规则下在体躯维持与生长一侧、成熟与繁殖另一侧之间分配。环境胁迫进入此类模型时体现为生物能量参数的修改，对贻贝而言，亚致死胁迫主要作用于摄食通量和体躯维持成本。一个由能量预算构建的繁殖响应不必是单一类型的，因为生命历程中的繁殖是两个条件事件的产物：达到成熟，以及成熟后投入配子。干扰可作用于任一组分——可减少成熟个体的繁殖能量分配，产生分级损失；也可通过将储备压低至足够低、足够久，在首次繁殖前触发死亡，从而阻止个体达到成熟。仅记录最终繁殖输出的终点指标混淆了二者。区分它们并识别干扰格局的哪个特征支配各组分，正是本研究的目的。

二、研究设计与方法

研究人员使用标准DEB模型，包含四个状态变量：结构体积（V）、储备能量（E）、成熟度（E_H）和繁殖缓冲（E_R）。同化遵循II型功能响应，速率经五参数Arrhenius函数进行温度校正。存活率S(t)为逐小时存活概率的乘积，饥饿风险（starvation hazard）采用分段线性激活规则。繁殖响应分解为E(R) = P_M·Q_M，其中P_M为达到成熟年龄的存活概率，Q_M为成熟后经条件存活加权的繁殖投入。模型参数化针对本地种紫贻贝和入侵种Brachidontes pharaonis，四种环境强迫情景源自三个地中海沿岸环境：的里雅斯特（亚得里亚海北部）、巴勒莫开放与巴勒莫沿岸（共享温度序列但食物供给不同，平均食物密度分别为0.19与0.49）以及西西里西部的马尔萨拉泻湖。本地种在巴勒莫开放强迫下未能在未受干扰对照中达到繁殖，故该组合被排除，共形成六种物种-强迫组合。干扰设计为规则的摄食中断序列，每年包含1至52个事件，每个事件将同化通量p_A降低σ_F比例；当事件持续时间超过事件间起始间隔时，连续事件融合为持续型扰动。研究区分了孤立脉冲域与持续型扰动域，主要分析在孤立脉冲域进行。计算了施加赤字D_imp（外生的，衡量格局设计所取走的能量）与实现赤字D_real（内生的，衡量受扰个体实际未能同化的能量）。统计上，各响应组分对log₂事件持续时间和log施加赤字进行回归，报告各预测因子的偏方差解释量而非显著性检验，因为该集合是确定性的。存活转变宽度量化为D₉₀?D₁₀，临界持续时间定义为D50。生存转变在三种替代风险公式（陡峭Logistic、平滑指数、宽分段线性）及饥饿参数的三乘三变化下重新计算。发育阶段分析采用三个初始发育状态（定居者、中间、近成熟）在共同日历时刻重新初始化并与两个季节起始（冬季、夏季）交叉，以分离阶段与季节效应。

三、研究结果

3.1 持续时间与累积能量赤字

在相等的累积能量赤字下，持续时间和施加赤字支配响应的不同组分。回归各组分于log持续时间和log施加赤字，持续时间主导达到成熟的概率（本地种偏方差解释量0.45，入侵种0.19）和总体表现（0.48和0.26），施加赤字贡献甚微（低于0.03）。对于条件投入——给定达到成熟后的繁殖投入——排序反转：施加赤字携带较大的偏方差（本地种0.28），持续时间较小（0.05）。在狭窄赤字带内，达到成熟的概率仍随持续时间急剧下降——中赤字带内从6小时时的0.91降至72小时时的0.38——而在受控持续时间下，条件投入按赤字带分离。当实现赤字替代施加赤字时结果不变。生物学上，短暂中断可存活且主要损失部分繁殖投入，而相同能量短缺作为一次长不间断时段传递，才使个体在繁殖前越过不可逆点：贻贝是否死亡由赤字无缓解持续多久决定，存活者繁殖多少由总能量损失多少决定。

3.2 繁殖响应的分解

所有未在时间范围内达到成熟的轨迹在成熟前已达到数值灭绝标准：480个单元中，401个成熟，79个在成熟前死亡，无一个在时间范围结束时存活但未成熟。因此1?P_M是成熟前死亡率而非延迟成熟。两个响应组分在形式上不同：未达到成熟概率随持续时间非线性增加而非在单一阈值处跃变，以10%至90%最大死亡率的区间表示，转变宽度本地种为23小时，入侵种为46小时。成熟个体中的条件投入在共同支持范围内随持续时间下降；分段拟合优于线性拟合，早期断点位于约12至18小时，远早于存活转变。将总繁殖损失L = 1?P_MQ_M分解为存活效应、繁殖效应与交互项，存活项在转变之外占主导。30小时最大摄食胁迫下，本地种存活效应为+0.93，条件繁殖效应为+0.15，交互项为?0.14。实际含义是：一旦事件长至足以致死，大部分繁殖损失源于个体从未成熟而非存活者繁殖减少。

3.3 储备余量与临界持续时间（D50）

在六个物种-强迫组合中，临界中断持续时间D50按对照储备余量排序：所有储备余量较小的组合在更短或相等的持续时间跨越转变。本地种在其强迫情景下维持约0.08至0.11的储备密度，入侵种为0.12至0.35。由于六个组合非独立，此排序作为描述性排序而非验证过的预测因子呈现。生物学上，未受扰时贻贝携带的备用储备设定其可吸收的中断时长。推断上，更接近能量边缘的种群更早达到致死阈值，使静息储备余量成为可在干扰前测量以标记最暴露种群的性状。

3.4 对风险公式的敏感性

基线饥饿风险在e_thr处有不可微点，因此非线性存活转变原则上可能由该规则强加。在陡峭Logistic、平滑指数和覆盖两倍储备范围的宽分段线性规则下重新计算，转变在每种情形下均持续存在。宽规则提高基线死亡率，因为其风险达到低余量组合的静息储备密度，但随持续时间的转变仍然存在，且组合按储备余量的排序在所有公式下保持或加强。在各自对照上归一化后，四种公式显示相似的转变位置和大致重叠的宽度，归一化D50在14至21小时内。稳健结果是储备消耗生成的非线性转变及其按储备余量的排序；精确D50取决于风险形式。存活崩溃的突然性是储备消耗本身的属性，而非死亡建模方式的伪迹。

3.5 对发育阶段的依赖性（与季节分离）

三个发育初始状态在共同日历时刻重新初始化并与两个季节起始交叉，D50在四个共享物种-强迫组合中的三个随发育阶段增加，第四个持平而非下降。效应一致但温和——新定居与近成熟状态之间相差一到两个持续时间水平——远小于与季节时序混淆的朴素比较所暗示的。D50因此是储备余量与个体距成熟距离的联合属性，但发育组分是调制而非主导项。事件发生的时节比其打击的个体发育状态更重要。耐受性主要由贻贝持有多少储备决定，其次才由生命阶段决定。

3.6 摄食与维持胁迫

摄食通量的降低驱动存活转变；维持成本的升高产生繁殖的分级降低，但在所探索范围内很少产生成熟前死亡。由于两个轴不在共同能量尺度上，此为模型探索。

3.7 种间比较与长持续时间下的单元选择

限于两个共享强迫情景，本地种在死亡率发生的每个持续时间携带较大的成熟前死亡率，仅在最长持续时间两者均接近完全损失时才收敛于入侵种。条件投入降低在物种间混合。长持续时间下仍达到成熟的格局集合缩小——168小时后每物种仅剩一或两个——因此存活格局上的条件均值是选择性子集，极端持续时间下任何表面恢复反映的是该子集的组成变化而非生理恢复。

3.8 结构响应

结构长度随持续时间逐渐下降，比繁殖下降更缓：在成熟个体中，本地种在最长共同支持持续时间损失约一半条件繁殖投入但仅约五分之一壳长。

四、讨论与结论

研究证实了三个假设：H1成立——相等施加赤字下，一次长中断杀死而若干短中断不杀死，赤字大小设定存活者的投入；H2成立——存活非线性崩塌于窄持续时间带，而存活者投入逐渐下降；H3对储备余量成立，它排序六种组合的临界持续时间；对发育阶段，依赖性存在但弱于事件季节时序。核心结果是角色的分离：赤字无间断运行多久决定个体能否活到繁殖，赤字总移除多少能量决定存活者到达后的投入量。这是同一干扰格局的两个不同属性，作用于生命历程的两个不同部分，单独任一都不能描述响应。

存活响应是加速转变而非固定阈值，其位置由储备消耗生成，且在连续、可微和宽泛的风险公式下均持续，故非亚模型函数形式的伪迹。储备余量在未受扰条件下维持，排序各组合的D50，且一旦季节时序受控，D50随接近成熟适度增加。脆弱性是储备余量与距成熟距离的联合属性，储备余量为主导项。

对干扰指标的含义：累积暴露指标如度加热周（Degree Heating Week）将给定累积暴露视为等价而不论其时间结构。对存活组分此等价失效——重要的是赤字的最大不间断持续时间而非积分暴露；对条件投入组分，累积测度更为相关。预期繁殖影响的指标需同时保留事件的连续性（针对存活）与累积量级（针对条件投入）。结构响应强化此点：壳长下降小于条件繁殖投入，故基于大小的监测终点可能低估亚致死繁殖影响。

如果脆弱性由中断持续时间与个体储备状态的耦合联合设定，干扰结果不能仅从物理事件推断。相同中断对接近饥饿边界的个体可能致死，对携带较大储备余量的个体则可存活，而该余量随近期食物条件、季节和发育阶段变化。这为大规模死亡事件中观察到的部分空间和时间异质性提供了可检验的个体层面能量假说。本地种风险更高：凡持续时间致死之处，紫贻贝携带较大成熟前死亡率。模型给出的解读是能量性的：本地种在未受扰条件下更接近其饥饿边界，较短不间断中断即足以将其推越；入侵种维持较大储备缓冲，需要更长的中断。若沿海干扰持续延长，此不对称性应逐步有利于入侵种。

结论：沿生命历程分解个体能量预算的繁殖响应，将非线性存活转变（由能量赤字的不间断持续时间支配）与条件繁殖投入的分级降低（由累积赤字支配）分离。存活转变的位置由生物体维持的储备余量排序，并受其距成熟距离的调制。四项实际后果：第一，在固定窗口内累加压力的暴露指标不能分离一次长中断与移除相同能量的若干短中断，监测计划应记录事件的最长不间断间隔及其总量级；第二，基于大小的终点低估损害，追踪壳长将错过已发生的繁殖失败；第三，由于转变按未受扰条件下测量的储备余量排序，干扰前的储备状态测量提供事件到来前对种群暴露排序的途径，食物贫乏地点或低储备季节的种群最需关注；第四，两物种在此余量上不同，延长干扰格局预期将平衡推向入侵种。这些后果映射到海洋战略框架指令的监测义务：D1C3标准涵盖种群人口统计特征包括体型结构及繁殖力和存活率，D1C6涵盖基础物种构建的栖息地类型状况。结果显示D1C3的大小组分和D1C6的扩展组分可在繁殖输出已下降时保持接近参考条件，故贻贝床可能在结构上被评估为良好环境状态而其更新能力已受损。两项低成本的补充可弥合此差距：在压力层记录温度、氧或食物异常的最长不间断持续时间及其累积量级，以及在已按指令1和6采样的哨兵床中将基于大小的指标与储备状况或配子输出配对。