《Journal of Freshwater Ecology》：From field observations to conditional simulation: Cladophora glomerata in Middle Atlas streams, Morocco
编辑推荐：
本研究考察了Cladophora glomerata覆盖度的空间和季节变化、其环境关联性以及在残差重抽样下的稳定性。数据集包含来自上乌姆埃尔雷比亚河网20个固定点位的240次月度观测。分析结合了含点位和月份效应的回归、按点位聚类的推断、主成分分析（PCA）以及
本研究考察了Cladophora glomerata覆盖度的空间和季节变化、其环境关联性以及在残差重抽样下的稳定性。数据集包含来自上乌姆埃尔雷比亚河网20个固定点位的240次月度观测。分析结合了含点位和月份效应的回归、按点位聚类的推断、主成分分析（PCA）以及AICc模型比较。所选模型用于校准包含10,000次实现的蒙特卡洛模拟。观测协变量以及拟合的点位和月份效应保持固定，而完整的年度残差序列在各河流内的点位间进行重抽样。平均覆盖度为50.5%（范围5.0%至71.7%）。点位和月份共同解释了转换后覆盖度方差的86.2%。温度具有最强的调整后关联性，但在约19.5°C处的最大值无法与季节性物候区分开来。硝酸盐和正磷酸盐未显示出受支持的独立单峰关系，流速也未显示中间最适值。所有12个观测到的月平均值均落在其条件模拟区间内。在96.7%至100%的实现中，环境系数保留了其观测到的符号。这些结果描述了在观测设计下残差重抽样的稳定性。它们并未确立生理阈值，也未验证对其他点位或年份的预测。
论文解读
一、研究背景与目的
丝状绿藻在河流和浅水水域的底栖生境中具有重要生态影响，其藻垫可改变近底水流、截留细颗粒物、提供食物和栖息地，并导致溶解氧和pH值的日变化。然而，过度生长也会损害栖息地质量、娱乐功能和水资源利用。Cladophora glomerata是一种常见于光照硬质基底上的丝状绿藻，但其分布广泛，难以用单一环境变量解释。现有研究虽已描述了摩洛哥中阿特拉斯山脉及上乌姆埃尔雷比亚流域的大型植物和藻类格局，但未能分离不同点位间的稳定差异、共同的年度周期以及C. glomerata覆盖度与环境条件之间的剩余关联。因此，研究人员旨在解决三个问题：第一，在一个完整年度周期内，C. glomerata覆盖度在不同点位和月份间如何变化？第二，在调整点位和月份效应后，哪些环境关联仍然存在，并且硝酸盐、正磷酸盐或流速是否在观测范围内呈现最大值？第三，残差重抽样能否保留月度覆盖度模式和环境关联的方向？
二、关键技术方法
研究人员采用的主要技术方法包括：基于20个固定点位、每月一次、持续12个月（2022年3月至2023年2月）的系统采样设计，共获得240次观测；使用含点位和月份固定效应的线性回归模型，并采用按点位聚类的CR1协方差估计进行推断；应用主成分分析（PCA）描述九种标准化环境变量的联合结构；利用校正后赤池信息准则（AICc）比较全部512种候选模型组合；最终基于最低AICc模型校准蒙特卡洛模拟，通过在同一条河流内的点位间重抽样完整年度残差序列，运行10,000次实现评估条件稳定性。
三、研究结果
3.1 空间和季节变化
年平均C. glomerata覆盖度为50.5%（标准差17.2%，范围5.0%至71.7%）。四个子流域的平均覆盖度分别为：Amengous 34.7%、Fellat 42.7%、Ichbouka 58.0%、Oum Er-Rbia干流57.7%。仅含点位的模型解释了36.5%的方差，仅含月份的模型解释了49.7%的方差，而点位和月份共同解释了86.2%（调整R2=0.842）。
3.2 原始相关性
未经调整的相关分析显示，与覆盖度相关性最强的是水温（r=0.922）、溶解氧（r=-0.836）和电导率（r=0.550）。流速、硝酸盐和正磷酸盐的相关性不显著，r值分别为0.080、0.025和0.111。
3.3 调整点位和月份后的关系
在调整点位和月份效应后，水温的二次项关系仍显著（F?,??=55.44，Holm校正后p<0.001，偏R2=0.813），曲线呈凹形并在19.50°C达到最大值。流速的二次项也显著（F?,??=15.76，Holm校正后p<0.001，偏R2=0.386），但曲线呈凸形，无中间最适值。硝酸盐和正磷酸盐均未显示受支持的独立单峰关系。溶解氧、亚硝酸盐和电导率在Holm校正后仍显著，而pH和铵则不显著。
3.4 多变量环境结构
前三个主成分分别解释了环境方差的34.8%、26.5%和16.9%，累计78.2%。包含七个主成分项的模型在点位和月份之外解释了额外变异（F?,??=49.65，p<0.001，偏R2=0.679）。PC1、PC12以及PC1×PC2交互作用在聚类推断下显著。
3.5 AICc模型比较
最低AICc模型包含水温、流速、亚硝酸盐和电导率，加上点位和月份效应后解释了97.8%的方差，其Akaike权重为0.140。累积权重显示，温度和亚硝酸盐的支持度最高（均为1.000），其次是流速（0.842）和电导率（0.834），正磷酸盐中等（0.420），硝酸盐最低（0.097）。
3.6 模拟结果
在10,000次实现中，模拟的年平均覆盖度为50.580%（条件95%区间50.501%-50.652%），观测值为50.544%。所有12个月的观测均值均落在其条件95%模拟区间内。温度和亚硝酸盐的系数在100%的实现中保留了观测符号，流速线性项为96.7%，流速二次项为99.7%，电导率为99.96%。
四、讨论与结论
讨论部分指出，覆盖度呈现出明显的年度模式，点位和月份效应解释了大部分方差，这与着生藻类现存量为生产与损失累积结果的解释一致。温度保持最强的调整后关联，但19.5°C的最大值不能视为生理最适温度，因为月份效应反映了日历平均模式而无法将温度与辐照度、水文、藻垫年龄等季节性过程实验分离。多变量分析和模型比较表明，没有任何单一模型明显占优，环境关联是在考虑了共享的季节模式和点位差异后估计的。研究未发现硝酸盐和正磷酸盐的区域性最适浓度，这不排除营养盐限制的可能性，但表明瞬时浓度不能可靠预测各河流中的藻类覆盖度。溶解氧的负相关不代表藻类偏好低氧水体，而是光合作用和呼吸作用的共同反映。流速效应取决于剪切应力、群落抵抗力和扰动历史，月度平均流速无法捕捉这些机制。研究局限性包括仅覆盖一个年度周期、测量的是视觉覆盖度而非生物量或生长速率、缺乏辐照度、浊度、底质粒径、放牧、河道宽度、流量和营养盐负荷等数据，以及未保留凭证标本或分子确认。
研究结论部分翻译如下：Cladophora glomerata覆盖度在不同点位和年度内变化强烈，并与多个环境梯度相关联。温度具有最强的调整后关联，但在约19.5°C处的拟合最大值无法与年度物候分离。硝酸盐和正磷酸盐未显示出受支持的独立单峰关系，流速也未显示中间最适值。一项以20个点位月度观测为校准、包含10,000次实现的模拟，在残差重抽样下保留了季节模式和大多数系数的符号。这些发现描述了在观测设计内的稳定性，并未确立生理阈值或对研究网络和年份之外的预测性能。
