大小和形状如何影响蓝藻菌落的垂直速度

《Water Research》：How size and shape affect the vertical velocity of cyanobacterial colonies

　　

蓝藻菌落常利用其浮力和大尺寸快速上浮形成密集的表面水华，这会降低水质、威胁生态系统和公众健康，并造成巨大的经济损失。然而，关于蓝藻菌落的形态如何控制其垂直速度仍知之甚少。研究人员对湖泊样品中微囊藻属（Microcystis）菌落进行了详细的单菌落水平三维形态表征，并在分层流中进行了受控的浮选实验。通过在垂直密度梯度中进行粒子追踪，研究人员独立量化了菌落形状和浮力密度在单个菌落水平上的贡献。研究结果表明，斯托克斯定律（Stokes' law）中的形状因子随菌落大小系统性地变化。因此，菌落的垂直速度不与菌落大小的平方成比例，而仅与1.13次幂相关，因为较大的菌落具有更不规则的形状，从而经历增强的阻力。因此，研究人员修正了常用的斯托克斯定律以考虑形状因子随尺寸的变化。有趣的是，在垂直迁移模型中实施该幂律关系显示微囊藻菌落迁移轨迹中存在广泛的混沌动力学。这些结果强调了蓝藻菌落形态可塑性的重要性，并可为有毒水华的预测模型和水动力控制策略提供信息。研究人员同时确定密度、形状因子和速度的方法广泛适用于淡水和海洋系统中具有复杂形状的悬浮聚集体。

**蓝藻菌落形态与垂直速度的关联：单菌落分辨率下的实验与模型研究**

**1. 研究背景与问题**

蓝藻水华在全球湖泊和水库中频发，尤其以微囊藻（*Microcystis*）属为典型代表，其形成的大规模表面水华会释放毒素、消耗溶解氧，严重威胁水质安全、生态系统平衡及公共健康，并造成巨大的经济损失。微囊藻菌落能够通过调节自身密度（依赖细胞内气体囊泡与碳水化合物压载物的平衡）实现垂直运动，从而获取光照或逃避高光强区域。这种垂直迁移能力是微囊藻在竞争中获得优势的关键生态性状。然而，菌落的垂直速度受其大小、形状和密度共同控制。经典的斯托克斯定律（Stokes' law）假设菌落为球形，并预测垂直速度与菌落等效直径的平方成正比（U ∝ D_ec²）。但大量观测表明，微囊藻菌落形态高度不规则，其形状随大小变化显著，且现有研究无法同时测量单个菌落的密度和形状因子，导致对垂直速度的预测存在较大偏差。Nakamura等（1993）和Li等（2016）虽指出形状因子随大小变化，但前者无法直接测量密度，后者则受限于手动操作且仅针对大菌落。因此，亟需发展一种能够在单菌落分辨率下独立测量密度、形状因子及垂直速度的方法，以揭示形态对垂直运动的真实影响，并改进现有水华预测模型和湖内控藻策略（如人工深混合）的设计。

**2. 研究方法概述**

本研究从荷兰的Braassemermeer湖（2024年8月6日采样）和't Joppe湖（2024年9月3日采样）采集表层水样，经35 μm网筛浓缩和低速离心分离获得富含微囊藻菌落的浮游组分（菌落丰度>98%）。主要技术方法包括：（1）**光学相干断层扫描（OCT）**：对压碎气体囊泡后的菌落进行三维扫描，获取菌落总体积（细胞+粘液层）的等效球直径（D_es）和等效圆直径（D_ec），并建立体积校正因子；（2）**共聚焦荧光显微镜与粘液染色（Alcian Blue 8GX）**：揭示细胞间排列和粘液层结构，并测量最近邻距离分布；（3）**气体分数测量**：利用毛细管压缩管计算每个细胞的气体囊泡体积分数（φ_g,c）和细胞直径；（4）**浮选速度测量**：在Percoll线性密度梯度中追踪单个菌落的上升轨迹，结合力平衡方程（含Richardson数修正）同时拟合出菌落的形状因子（ψ）和密度（ρ_a）。此外，对压碎气体囊泡后的菌落进行密度测量以验证方法稳健性。

**3. 研究结果**

**3.1 微尺度与介尺度特征：体积与投影尺寸的错配**
OCT扫描显示，大菌落常呈现分支结构和内部隧道，其等效圆直径（D_ec）系统性高估了等效球直径（D_es），两者遵循幂律关系：D_es/L_d = (D_ec/L_d)^f_d/3，其中L_d = 77 ± 20 μm，f_d = 2.36 ± 0.10（N=66，R²=0.89）。粘液染色和共聚焦显微镜证实粘液填充细胞间隙并形成薄外鞘，细胞间最近邻距离分布宽广（平均NND/d₁ = 1.35，SD = 0.55）。气体分数测量显示，两个湖泊的每个细胞气体体积相似（约2.0–2.3 μm³），但细胞直径略有差异（Braassemermeer湖平均4.6 μm，'t Joppe湖平均4.1 μm）。这些结果揭示了投影几何导致的体积低估偏差，并提供了菌落微观结构的关键参数。

**3.2 浮选速度受菌落密度和尺寸依赖形状因子的影响**
通过粒子追踪实验，研究人员发现形状因子（ψ）与等效圆直径（D_ec）高度正相关（Pearson r = 0.87，p < 0.001）。小菌落（D_ec ? 100 μm）的ψ ≈ 1，而大菌落（D_ec = 971 μm）ψ可达10.44。ψ与D_ec遵循幂律：ψ = (D_ec/L_s)^f_s，其中两湖合并数据的最佳拟合参数为f_s = 0.87 ± 0.03，L_s = 67 ± 3 μm。菌落密度（ρ_a）对尺寸的依赖较弱（Pearson r = 0.16，p = 0.03），且两湖之间无显著差异（Mann–Whitney U检验，p = 0.71）。气体囊泡压碎后密度无显著尺寸依赖，且分布方差与完整菌落相似，验证了密度测量方法的可靠性。基于这些结果，研究人员修正了斯托克斯定律，得到垂直速度与D_ec的1.13次幂成正比（U = [L_s²(ρ_w - ρ_a)g / (18η_w)] · (D_ec/L_s)^1.13），该修正优于经典平方律，并能更好地拟合已发表的多组数据（R² = 0.19 vs. 0.10）。

**3.3 生态意义：垂直迁移模式的混沌动力学**
将修正的斯托克斯定律应用于标准垂直迁移模型（Feng等，2025），发现大菌落（D_ec > 710 μm）才表现出周期性日迁移，而中等菌落（D_ec = 350 μm）在经典定律下为周期性迁移，在修正定律下则变为混沌迁移。小菌落（D_ec = 50 μm）在两种模型中均呈现混沌模式。此外，迁移模式对形状分形指数（f_s）、细胞体积分数（φ）和细胞气体分数（φ_g,c）等形态参数高度敏感，表明形态可塑性是影响水华动态的关键因素。

**4. 讨论与结论**

研究结果表明，微囊藻菌落的形状因子随尺寸系统性增加，导致垂直速度的尺寸依赖指数远低于经典平方律（1.13 vs. 2），这修正了以往模型中普遍采用的球形假设。该发现对水华预测具有重要意义：用修正模型预测的迁移速度更慢，使得更多菌落表现出混沌迁移模式，从而影响其在光照和水柱中的分布。此外，这一结果可指导人工深混合系统的设计，避免因过高估计浮选速度而导致过度混合和能源浪费。研究同时指出，形态参数（如形状分形指数、细胞体积分数）受生长条件和形态种组成影响，未来需进一步比较不同微囊藻形态种及其他蓝藻菌落的形状因子。方法学上，本研究所建立的同时测量密度、形状因子和速度的技术，可推广至海洋雪、藻类-微塑料异质聚集体等复杂形状悬浮颗粒的垂直输运研究。

**研究结论**（翻译自原文第4节）：本研究通过提供单菌落分辨率的形状因子和浮力密度精确测量，推进了对微囊藻浮选速度的现有认识。研究人员证明菌落形状随其大小系统性变化，因此菌落的垂直速度不与菌落大小平方成比例，而是以远低于经典斯托克斯定律预测的指数增长。此外，研究人员显示菌落的水动力阻力可从其最大线性尺寸估计，该参数易于通过显微镜获取。将新的尺寸依赖形状因子应用于垂直迁移模型，预测大菌落迁移速度更慢，从而扩大了呈现复杂迁移模式的菌落尺寸范围，并将大菌落限制为每天一个迁移周期。研究结果表明，考虑菌落大小和形状可提高对蓝藻菌落垂直迁移的理解，该迁移在蓝藻水华形成中起关键作用。