观潮 · 模式预报与气象 AI 观察
海洋施铁肥的碳账：三维拉格朗日斑块模型算到 42 天 690 吨
> 一个跟着水团走的三维拉格朗日斑块模型，把"给海面施铁肥能固多少碳"这件事从"假设斑块内部均匀"推进到"斑块一边被拉伸稀释、一边发生生化反应"。模型把水团近似成椭球，分辨四种改变它尺寸与形状的物理机制，用雷诺分解同时解出示踪物浓度的空间均值与二阶矩，并新引入一个分压差示踪量来算海气碳通量。对南大洋施铁实验 SOIREE 的 42 天模拟给出：前 13 天吸收 32 吨碳，实测是 31 吨，而一个针对 SOIREE 调优但不算斑块物理的模型给出 130 吨；42 天累计 690 吨，其中只有约 5% 发生在前两周。作者写明：水华晚期的直接观测缺失，而模型显示大部分生物量增长正发生在这段。
缺的是哪一块
给表层海洋补限制性营养盐会刺激浮游植物水华，进而吸收碳，这被当作一种潜在的二氧化碳移除手段提过很多次。但真要算清这笔碳账，难在两点。
观测侧：局地尺度的水华动力学要求从拉格朗日视角追踪那一团水随时间的生化演化，而可用的实测既稀疏又不同步。
模型侧：现有的多数施肥实验碳吸收模型，假设正在演化的水团内部是均匀的。近几年的工作指出，物理稀释、水平平流以及由此产生的斑块内空间异质性，对积分生产力有重要的定量后果——这一条被均匀假设抹掉了。
方法/系统怎么搭的
斑块近似成一个三维椭球。 宽、长、厚三个方向分别跟随演化的水团；参考系跟着斑块平移并水平旋转（倾斜忽略不计）。作者识别出四种改变三维水团尺寸和形状的物理机制：水平应变（沿某个水平方向拉伸的倾向）、辐散（斑块变平并径向扩张）、水平扩散、垂直扩散。水平与垂直扩散率分开处理，因为两者量级差一个数量级，而且水平搅拌与扩散的空间尺度依赖性是被动态解析的。
解三维平流扩散方程，得到体积增长。 解出来可以看到：水平扩散对体积膨胀的作用受"周长面积比 × 厚度"调制，垂直扩散的贡献则由"面积厚度比"调制。体积的相对增长率被定义为卷入率。
示踪物主控方程 + 雷诺分解。 斑块内示踪物浓度变化的三个驱动是：卷入浓度不同的周围水体、斑块内部的物理混合、示踪物之间的生化相互作用（比如吸收碳）。为了显式处理空间异质性，作者用雷诺分解，解出浓度的一阶矩（空间均值）和二阶矩（空间方差与协方差）。
新引入的分压差示踪量。 这是本文相对前人框架的增量：斑块内部与上方大气之间的二氧化碳分压差 Δp。它由三项驱动——海气交换、生物吸收碳、海表温度变化带来的溶解度强迫。浮游植物群落动力学沿用前人的框架。
两个积分指标。 拉格朗日碳通量 LCF，定义为 τ 天窗口内穿过不断演化的斑块表面吸收的碳质量异常（吨），积分量是斑块内部平均通量与水平周围水体通量之差；拉格朗日初级生产 LPP，在斑块体积上积分生物项相对周围水体的异常。
结果数字
对 SOIREE 的 42 天复算
SOIREE 是第一次在南大洋极地海域做的施铁实验，用惰性示踪剂 SF₆ 做拉格朗日设计：施铁后实地追踪了大约两周，之后靠卫星又跟了 40 天。
量 · 模型 · 观测 / 他家模型
第 8 天溶解铁浓度 · 约为水平周围的 23 倍、下方水体的 6 倍 · —
铁均值回到背景 · 约 25 天 · —
水华达峰 · 约 21 天，平均浓度接近周围水体的 16 倍 · —
前 13 天 LCF · 32 吨 · 实测 31 吨；SWAMCO 130 吨
前 13 天 LPP · 约 420 吨 · 沿 SF₆ 控制体测到的藻类碳增量 400 吨；他家评估 1200–1400 吨
42 天 LCF · 690 吨，其中只有约 5% 发生在前两周 · —
42 天 LPP · 1640 吨 · —
铁的空间方差约 4 天后就稳定下来，浮游植物的方差则要几周才达峰。当铁均值回落到背景水平，铁与浮游植物的相关性转负，反映的是浮游植物对局地铁亚斑块的耗竭。
和 SWAMCO 的对比值得单独说：SWAMCO 是针对 SOIREE 调优的生物地球化学模型，但不算水团物理。两个模型都低估了前 13 天斑块内 Δp 平均异常相对实测数据的水平——这是作者直说的差距，不是某一边单独的问题。他家对初级生产增量的评估高得多（1200–1400 吨），作者的判断是那些研究大概用了更大的积分体积。
约 3000 次集合试验的三条结论
稀释并不只是坏事。 在不分辨二阶矩（也就是回到均匀假设）时，LCF 随稀释增强单调下降——任何周围水体的侵入都会降低从大气吸收的二氧化碳。但分辨异质性之后，空间异质性提升了施肥后的 LPP。作者的解读是：稀释造出的示踪物异质性反过来能增强浮游植物生态系统的生物量增长。这条对实验设计有直接含义——海洋里那些"稳定、连贯"的区域，未必是最大化群落生产力的最佳投放点。
初始投放面积存在最优值。 均匀的斑块在营养盐集中在小于 70 km² 时更高产；把异质性算进去之后，给更大的斑块施肥（70 到 150 km²）带来更强的生态响应。这个最优值由物理驱动因子的尺度律与依赖于尺寸的浓度梯度共同决定。
垂直过程的方向性。 固定垂直扩散率、改变初始水平辐散的量级与符号，在所有水平稀释条件下 LCF 约提升 40%（图 5a、5c）；结论章把这一结果表述为"混合层以下的上升流增强以及相伴的水平辐散会提升水华期的净碳吸收"。LPP 的反应不同——它对体积效应敏感，体积缩小抵消了浮游植物浓度的上升，净变化可忽略。至于垂直扩散率的量级本身，作者对碳通量和生物量的影响描述是"更复杂"，没有给出单调结论。
原文图
图 1：(a) t 时刻一个三维椭球状拉格朗日斑块的快照，宽、长、厚分别为 W(t)、L(t)、T(t)；某一示踪物在斑块内的分布记为 u_i(x,t)，斑块四周、上方、下方的分布分别记为 h_i(x,t)、b_i(x,t)、a_i(x,t)。(b) 水平与垂直扩散、应变、辐散对一个拉格朗日斑块的作用：两个方向的扩散都通过稀释内容物增大斑块体积，而应变与辐散分别沿水平方向和垂直方向拉伸斑块。
图 2：SOIREE 施肥实验的模拟。实线是模型变量随时间的演化：(A) 斑块的宽与长，(C) 铁与浮游植物的平均浓度，(D) 两者的浓度方差，(B) 铁与浮游植物的空间相关。C 图用了左右两个纵轴，因为铁与浮游植物浓度的量级差得多。圆点是实测，叉号是卫星遥感观测。浮游植物生物量按铁当量表示。
图 3：模拟的 pCO₂ 异常（Δp）及其在 SOIREE 期间与铁、浮游植物的关系。灰色小点是实测 Δp，下方那团点对应施肥斑块内部的观测，上方那团是斑块周围水体的观测。实线是模型变量随时间的演化：(A) 模型 Δp 平均浓度对照 SWAMCO 的预测，(B) Δp 浓度方差，(C) 铁与 Δp 的空间相关，(D) 浮游植物与 Δp 的空间相关。A 图中的点划线是斑块外的模拟 Δp，对应上方那团观测点。
图 4：42 天后拉格朗日碳通量（LCF）与拉格朗日初级生产（LPP）随斑块稀释因子和初始注入面积的变化。红色点线是分辨了示踪物浓度二阶矩的配置，蓝色点线是忽略二阶矩（即均匀假设）的配置，黄色星号是 SOIREE 位置。
图 5：42 天后 LCF 与 LPP 随稀释因子的变化，叠加不同垂直动力条件，四幅面板都是分辨异质性的配置。黑线是图 4 中的红曲线，深红与青色曲线分别对应改变初始水平辐散 δ₀、垂直扩散率 κ₀ᵛ 的模拟输出，黄色星号是 SOIREE 位置。
作者承认的限制
- 水华晚期没有直接观测。 SOIREE 实地追踪只做了约两周，而模型结果显示大部分生物量增长发生在这之后。作者说这一段的直接观测（可能要靠拉格朗日自主平台）会是非常有价值的约束。
- 施肥水团的识别与追踪本身是根本挑战。 作者明说，水团总体积的不确定性是文献里报道的积分量彼此不一致的一个主要来源。
- 亚中尺度上的垂直过程只是推测。 原文的用词是"我们推测"（We speculate）这类垂直过程在亚中尺度上作用更大，因为那里垂直运动越来越重要。
- 真实中尺度与亚中尺度流动作用在不同斑块尺寸上的系统分析留待将来，本文没有做。
- 两个模型都低估了前 13 天的 Δp 平均异常，这是模型与实测之间尚未弥合的差距。
一句话判断
这笔碳账能不能算准，卡点不在生物地球化学的参数化，而在于有没有把水团本身的物理演化——被拉伸、被稀释、内部变得不均匀——一起算进去；本文给出的 32 吨对 31 吨，说明算进去之后量级就对上了。