**OxyPOM：具有过程特异性温度敏感性的氧气与颗粒有机物动态生物地球化学模型**

《Ecological Modelling》：OxyPOM: A biogeochemical model for Oxygen and Particulate Organic Matter dynamics with process-specific temperature sensitivity

　　

溶解氧是水质和水生生态系统健康的基本指标。氧气浓度降低，尤其是缺氧（<4 mg L^?1）和无氧（<2 mg L^?1），会对好氧生物造成不利影响，并促进沉积物中温室气体、营养盐和毒素的释放。水系统中的氧气收支由生物地球化学（BGC）过程控制，在模型应用中这些过程充当氧气的源（如光合作用）或汇（如有机物矿化、呼吸、硝化及其他生物与非生物反应）。大气-水氧气交换（再曝气，re-aeration）根据水体温度依赖的饱和状态而充当源或汇。这些控制变量因其温度敏感性而表现出显著的季节性变化，导致主导过程在一年内发生转换。水温扮演着尤为关键的角色，因为长期变暖趋势和短期热浪均与观测到的脱氧速率密切相关。然而，在当前的模型应用中，温度变化对氧气相关过程的相对效应仍缺乏足够的约束。许多现有水质和生态系统模型将氧气动态描述为复杂程度不同的生物与非生物相互作用的副产物，仅有少数模型解析了直接支配氧气收支的过程。此外，其中若干模型与特定水动力模型耦合，缺乏开源实现和标准化接口，限制了其在不同水动力建模平台上的可及性、互操作性和可重用性。本文中，研究人员提出了生物地球化学模型OxyPOM（氧气与颗粒有机物模型），该模型在氧气相关过程以及低营养级的其他生物地球化学速率中纳入了各自的温度敏感性。研究人员的目的是将OxyPOM呈现为一个适用于量化驱动氧气和POM动态的生物地球化学与物理过程相对重要性的模型，并着重关注温度敏感性的作用。研究人员的开源实现采用水生生地化模型应用程序接口框架（FABM API），确保了可及性、与其他水生过程模型的互操作性、跨不同水动力模型的可重用性，以及从零维到三维域的可应用性。为突显特异性温度敏感性在生态系统动态中的影响，研究人员比较了两种温度敏感性公式情形：（1）所有过程共享单一敏感性参数；（2）每个过程拥有独立的参数。研究人员将OxyPOM耦合到通用海洋湍流模型（GOTM）中的一个理想化、一维、深度分辨的水动力装置中，以模拟代表具有季节性缺氧特征的典型河口情境下的氧气动态。

## 论文解读

**一、研究背景与现存问题**

溶解氧是衡量水生生态系统健康状况的核心指标。当溶解氧浓度降至缺氧（<4 mg L^?1）甚至无氧（<2 mg L^?1）水平时，好氧生物的生存受到严重威胁，同时沉积物中温室气体、营养盐和毒素的释放也会加剧。水生态系统中氧气收支由多种生物地球化学过程共同支配，包括光合作用（产氧源）、有机物矿化、呼吸作用、硝化作用（耗氧汇）以及大气-水界面的再曝气过程（根据饱和状态可为源或汇）。这些过程均表现出显著的季节性波动，而水温是驱动这种变化的关键因素。长期气候变暖与短期热浪事件均与观测到的水体脱氧现象密切相关。然而，现有大多数水质与生态系统模型在描述氧气动态时，要么将其视为复杂生物地球化学相互作用的副产物而未能充分解析，要么对不同过程的温度敏感性采用统一化处理，忽略了各过程在温度响应上的本质差异。此外，许多模型与特定水动力模型紧密耦合，缺乏开源实现和标准化接口（如FABM API），严重制约了模型的可及性、互操作性和跨平台可重用性。尤其值得注意的是，当前先进模型仍普遍忽略了若干已被证实具有温度敏感性的过程，包括病毒致病机制（pathogenesis）和浮游植物最适光强（optimal light intensity for phytoplankton）等。准确评估气候变暖背景下的氧气响应因此变得极为困难。

**二、研究内容与总体结论**

针对上述问题，研究人员开发了名为OxyPOM（Oxygen and Particulate Organic Matter）的生物地球化学模型。OxyPOM在氧气相关过程及低营养级生物地球化学速率中均赋予各自独立的温度敏感性参数（即过程特异性Q₁₀系数），涵盖光合作用、呼吸、矿化、硝化、再曝气、最适光强、捕食死亡率以及病毒介导的裂解与失活等关键过程。模型基于FABM API实现开源，确保与不同水动力模型的互操作性和从零维到三维域的应用灵活性。为揭示过程特异性温度敏感性的重要性，研究人员将OxyPOM耦合到理想化的一维GOTM（通用海洋湍流模型）水动力框架中，模拟典型河口季节性缺氧环境，并比较了两种参数化方案：一是所有过程共享统一的Q₁₀值（几何均值2.5）；二是各过程采用文献报道范围内的独立Q₁₀值（参考情形）。研究表明，温度敏感性的参数化方式对生态变量和氧气过程的季节动态具有显著影响，统一温度敏感性会减弱初级生产、有机物矿化与营养盐再生之间的耦合强度，导致夏季底层水体溶解氧偏高、碳相关库偏低、营养盐偏高的系统性偏差，并使绝对氧通量低估2至5倍。OxyPOM为研究变暖世界中的氧气动态和极端物理-海洋学事件（如海洋热浪）提供了灵活而有力的建模工具。

**三、关键技术方法概述**

研究采用的主要技术方法包括：OxyPOM模型开发与FABM API开源实现；基于van't Hoff规则的Q₁₀温度敏感性公式（τ_K(T) = Q_10,K^(T?T*)/10，参考温度T* = 20 °C）；一维理想化水动力模拟采用GOTM模型，设置20 m水深垂直分辨，气象强迫基于易北河河口（53°52′N, 8°42′W）气候学特征；模拟连续运行两年，仅采用第二年的输出数据进行对比分析，以避免初始条件带来的伪影；设置统一敏感性（所有Q₁₀ = 2.5）和过程特异性敏感性（参考情形，各Q₁₀见表2）两组数值实验；通过绝对差异（x_uni ? x_ref）、相对百分比差异D（D = 200 · (x?_uni ? x?_ref) / (x?_uni + x?_ref)）以及温度-相对差异二维密度图评估两种参数化的影响；所有数据分析与可视化在R语言环境中完成。

**四、研究结果**

**（1）生态变量的时空差异**

溶解氧、总浮游植物、DOC（溶解有机碳）、总POC（颗粒有机碳）以及营养盐（氨氮、硝酸盐、磷酸盐）的时空分布在不同温度敏感性方案下存在重要差异。统一敏感性情形下，7月至10月底层水体溶解氧偏高约100 mmol O₂ m^?3；参考情形中夏季浮游植物水华更为显著，7月至8月上层10 m总浮游植物偏低约10 mmol C m^?3；DOC在统一敏感性情形下各季节和深度均略低，7月至10月底层积累偏低10–40 mmol C m^?3；总POC在统一敏感性情形下偏低，底层差异可达100–200 mmol C m^?3。与碳相关变量相反，氨氮、硝酸盐和磷酸盐在统一敏感性情形下大部分时间均偏高。

**（2）病毒与氧气相关过程的温度敏感性响应**

垂直平均的病毒丰度在两种配置间差异较小，仅在7月至10月及4月的个别峰值处出现较大偏离，且与温度的关联较弱，仅在10–14 °C狭窄温度区间出现较大绝对差异。反硝化速率季节趋势紧随温度变化，夏季最高，但特异性敏感性情形全年大部分时间反硝化速率更低，仅在夏季意外偏高，相对差异范围为?100%至50%。矿化和硝化耗氧量均随温度变化，夏季速率最高，1月至4月统一敏感性情形下矿化与硝化均偏低，相对差异范围分别为?10%至100%和?5%至100%。光合产氧量与呼吸耗氧量并不严格紧随温度变化，两者均在夏季达到峰值，但相对差异与温度的关联较弱，表明在藻华期间，光照可获性、营养盐限制、捕食及病毒感染等非温度因素可能发挥主导作用。夏季所有产氧和耗氧过程在特异性敏感性情形下均更高（D > 0）。

**（3）氧通量差异**

夏季氧气通量的绝对量级在不同温度敏感性方案下差异巨大。统一敏感性情形的绝对氧通量仅为特异性敏感性情形的1/5至1/2。7月和8月统一敏感性情形中缺失了强脱氧事件的信号。两种情形在相对氧通量方面模式相近，但7月例外。光合作用是主要的氧来源，再曝气充当氧汇，呼吸是主导氧汇；自8月下旬起矿化和硝化的贡献增强，而此时再曝气转为氧来源。

**五、讨论与结论总结**

讨论部分指出，温度敏感性的参数化选择会通过直接和间接途径修饰整个生态系统的响应，不仅限于氧气浓度本身。统一温度敏感性可能削弱初级生产、有机物矿化与营养盐再生之间的耦合，这对设计温度操控实验具有重要参考意义。令人意外的是，将浮游植物生长的Q₁₀从参考情形中的1.7提升至统一情形的2.5后，夏季浮游植物浓度反而降低，这源于统一敏感性下藻华规模受抑进而通过下行效应减少了捕食、呼吸和矿化，降低氧消耗与营养盐释放。病毒致病过程对温度参数化仅表现出有限敏感性，与实验观察中病毒感染速率对温度的依赖性弱于宿主生长的现象一致。不同微生物途径（如反硝化菌与硝化菌）具有不同的热最适性，统一Q₁₀无法捕捉这些差异，导致温暖月份过程速率系统性低估2至5倍。此外，最适光强温度敏感性虽在模型中不常被纳入，却能对藻华发生时间与规模产生强烈影响。文献中报道的各过程Q₁₀值高度可变，构成模型应用中的主要不确定性来源之一。OxyPOM通过独立参数化的方式能够灵活重现从单调递增到钟形曲线等复杂的温度-生长响应模式，这种灵活性源于生长、呼吸和死亡等过程温度效应的独立表达与组合。模型当前的主要局限包括忽略硫、铁、锰等氧化还原过程、采用固定化学计量比假设、未显式包含沉积物-水界面过程，且结果来自理想化水动力装置，尚未经过针对真实系统的观测数据校准验证。然而，OxyPOM已涵盖解释远洋系统氧气动态所需的最相关过程，并提供了温度依赖性速率的灵活参数化方案。

研究结论可翻译如下：OxyPOM模型通过赋予每个与氧气相关的过程独立且特定的温度敏感性系数，克服了当前主流水质模型在该方面的不足。研究结果揭示，过程特异性温度敏感性对于准确模拟溶解氧、碳库和营养盐的季节动态以及氧通量量级具有决定性作用——若采用统一温度敏感性参数化，可能导致绝对氧通量被低估2至5倍。OxyPOM基于FABM API的开源架构确保了其在不同水动力模型和空间维度上的高可移植性，为研究气候变暖背景下水生生态系统的脱氧过程及极端事件影响提供了强有力的定量分析工具。