从出流减少到生态退化：基于Kolmameter的河床淤填（colmation）阈值识别退化溪流

《Limnologica》：Thresholds in the dark – Colmation thresholds as a tool for identifying degraded streams

　　

Kolmameter可量化沉积物中间隙出流减少（interstitial outflow reduction）。Kolmameter出流减少为沉积物淤填（sediment colmation）提供实用代理指标。已识别间隙动物（hyporheic meiofauna）的淤填阈值。30%出流减少阈值可捕捉间隙群落更替（hyporheic community turnover）。当≥50%样本超过30%阈值时，河段被判定为退化（degraded）。

论文主体解读

细颗粒泥沙在河床累积是广泛但常被忽视的溪流生态压力。欧洲土地利用与水文形态变化增加泥沙负荷和细颗粒滞留；水框架指令（Water Framework Directive, WFD）旨在保护和改善地表水生态与化学状态并实现良好状态，但多数德国地表水仍未达良好，河床淤填（colmation，指<2 mm细颗粒泥沙在河床沉积物内和表层的入渗与累积）是缺口贡献之一。河床淤填可表现为表面细颗粒累积或内部孔隙入渗导致水力传导下降；虽属自然过程，人类活动可加剧并使其更持久。它降低沉积物渗透性、限制垂向水交换并破坏间隙带（hyporheic zone，河床下和侧方地表水与地下水混合的生态交错带），影响氧供应、营养循环和有机质处理。间隙动物（hyporheic meiofauna，<1 mm的无脊椎动物）利用间隙空间避难、摄食和繁殖；孔隙空间、氧穿透和间隙流减少使敏感生物栖息地不适，导致群落变化、丰度和分类丰富度下降、组成改变及自净和营养周转等功能受损。已有基于大型无脊椎动物的细泥沙评估工具和敏感响应指标，但间隙动物群落的明确阈值仍大多未被明确。生态阈值是环境梯度上小压力增加引起生物群落突然、非线性变化的点或带，对管理重要；阈值指示类群分析（Threshold Indicator Taxa Analysis, TITAN）可检测群落和类群响应。研究人员在一般退化（general degradation, GD）1—2级河流中，以出流减少为淤填代理识别生态阈值，检验丰度、分类丰富度和群落组成是否随出流减少一致变化，识别断点及敏感与耐受类群，为泥沙管理和WFD实施提供经验基准。

关键技术方法（不超过250字）：研究于2021年10月至2022年11月在德国Rhineland-Palatinate 10条河和奥地利Upper Austria 1条河开展；德国含5条类型5、5条类型7，奥地利对应类型5；每河设20 m河段、10个采样点，3次季节采样。关键技术包括：Kolmameter（一种测定河床沉积物出流减少的原位装置）量化出流减少（0.05 bar、5 s、3重复），压测管采集间隙水并现场测间隙溶解氧（interstitial dissolved oxygen, IDO）、pH、电导率（electrical conductivity, EC）和温度，73 μm筛分保存并鉴定间隙动物至较高分类阶元；以TITAN识别阈值（限于GD 1—2，1000次bootstrap，可靠性≥0.90），用置换多元方差分析（PERMANOVA）、多元离散度检验（PERMDISP）、Spearman相关、多元线性回归（MLR）和线性混合效应模型（LMM）分析，并以Youden's J统计量和层次bootstrap构建分类框架。

研究结果

3.1 分类单元对出流减少的响应（Taxa response to outflow reduction）
11条河324个采样点共记录32个分类单元；排除少于3个采样点出现的类群以及Collembola、Ephemeroptera、Plecoptera和Trichoptera后剩25个。GD 1—2级7条河205个采样点保留22个分类单元。TITAN识别出5个z^?指示类群（Cyclopoida、Ostracoda、Chironomidae、Macroturbellaria、Nauplii）和1个z⁺指示类群（Parastenocaridae）。bootstrap群落阈值为z^?组23.1%（95% CI：9.5—36），z⁺组33.4%（95% CI：29.6—39.3）；密度图显示约23%负峰和约33%正峰，表明敏感类群向耐受类群转变。按河流类型，类型5的fsumz^?阈值为24.1%，类型7为26.1%；类型5的fsumz⁺为33.1%，与聚合数据33.4%接近，类型间fsumz^?差异约2%且生态上可忽略，支持合并GD 1—2河流分析。因23.1%和33.4%不便于应用，进一步检验30%阈值；基于四次方根转化丰度的PERMANOVA显示30%阈值上下群落组成显著不同（Bray-Curtis相异性；p=0.011，R²=0.013），PERMDISP亦显著（p=0.023），提示群落离散度差异可能贡献该结果。

3.2 间隙指示类群的方向性丰度格局（Directional abundance patterns of hyporheic indicator taxa）
30%阈值上下总丰度差异显著（Mann-Whitney U检验：p=0.02），中位数分别为阈值下166个个体、阈值上122个个体。指示类群方向一致：Cyclopoida（p=0.008）、Macroturbellaria（p=0.027）和Ostracoda（p=0.004）在阈值下较高；Parastenocaridae（p=0.004）在阈值上较高；Chironomidae（p=0.065）和Nauplii（p=0.134）未显著。总丰度下降、多数指示类群一致响应及阈值组间显著差异，支持30%群落阈值的生态相关性及其在间隙群落监测管理中的潜在应用。

3.3 出流减少、环境参数与间隙丰度的关系（Relationships between outflow reduction, environmental parameters, and hyporheic abundance）
Spearman相关显示出流减少与间隙温度正相关；含出流减少与采样轮交互项的MLR表明该关系随采样轮变化，仅2022年4月显著（estimate=0.289，p=0.007），2021年10月（estimate=?0.102，p=0.137）和2022年10月（estimate=0.091，p=0.448）不显著，说明温度相关具季节依赖性且春季最强。出流减少与分类单元数和总丰度负相关。阈值组间仅总丰度（p=0.02）和温度（p=1.4e^?5）显著；pH（p=0.4）、EC（p=0.2）、IDO（p=0.5）无显著差异，分类单元丰度仅边际趋势（p=0.055）。LMM显示所有环境变量在河流间均显著不同（p<0.001）。结果表明仅温度和总丰度与阈值组差异有意义，基本水化学与出流减少关系较弱。

3.4 管理分类框架（Classification framework for management）
候选决策阈值评估显示最佳cut-off约0.50（bootstrap均值=0.501；中位数=0.510），95% bootstrap区间0.30—0.81，平均Youden's J=0.344；Youden's J曲面呈近优区域而非单一锐点，0.45—0.50达到最大J的至少90%，定义为管理带。纳入采样不确定性后，层次bootstrap得到管理带上下界均值0.464和0.609，95% CI分别为0.3—0.781和0.40—0.90，支持约0.50的高分类性能中心区。因此，当河段中至少约50%样本超过30%出流减少阈值时，可视为退化（淤填）河段；3.4原文以“more than 50%”表述。

讨论部分总结：最小影响溪流中间隙群落对出流减少呈清晰阈值响应，与大型底栖参考方法一致；粗基质溪流在出流减少超过30%后群落强烈重组，TITAN敏感类群变化点23.1%、耐受类群33.4%。WFD未规定淤填阈值或是否随河流类型不同；本研究单一阈值适用于类型5和7，但可能不适用于更大或沙主导河流；因预算限制，30%阈值以同一数据集计算和验证，仍需额外河流独立验证；该阈值落在既有报告范围内，可捕捉敏感类群下降、耐受类群持续的功能重组起点。分类分辨率较粗可能掩盖种级差异；Parastenocaridae自然逐渐增加，负响应类群多较突然； harpacticoids钻得更深、小蠕虫形体型增强在淤填沉积物中移动，可能利于持续；Chironomidae随细泥沙下降但部分种类耐受较高泥沙，Macroturbellaria主要与低出流减少相关，Copepoda和Ostracoda随出流减少下降并与粗沉积物偏好一致。栖息地丧失可能是主要压力：细颗粒填充孔隙使出流减少增加、基质不可栖，并可能降低氧可用性，降低间隙带承载能力，解释总丰度下降及30%阈值以上敏感类群下降。温度、pH、EC和IDO影响次于出流减少；温度相关仅2022年4月明显，可能反映春季融雪和径流 mobilise 间隙沉积物并改变流路；pH、EC、IDO无阈值组差异且处于硅质和钙质溪流预期范围；河流间环境变量差异显著，说明其更多反映河流特征；理化参数仅在接近耐受限时影响分布，IDO足以支持多样类群；但生态状态、间隙群落和出流减少受多种互作或协变因素及未测栖息地特征影响，测试参数未涵盖全部驱动；结果仍表明出流减少是间隙群落变化的重要机制驱动，监测出流减少可作为早期检测泥沙压力的敏感指标，而标准水化学可能漏检早期生态响应。管理上，至少50%样本超过30%阈值提供监测管理基准；45%—50%管理带提示渐变而非锐断，应作为决策支持范围而非固定生态边界；bootstrap最优中心约0.50但CI宽0.30—0.81，0.50为最佳cut-off但处于不确定包络内；Youden's J透明可重复，未设成本比；该cut-off不是独立生态状态指标，而是管理导向指标而非WFD标准，应与其它WFD参数并用；即使GD 1—2站点也可能超过50%，因出流减少自然时空变异且局部淤填不一定改变总体分类；管理带内应谨慎并考虑补充采样或追查压力源；需考虑GD分类不确定性、河流类型和评估方法差异；GD 4和5各仅1条河，需更多严重退化溪流验证；细泥沙累积和出流减少自然随时间变化，建议全年重复采样；TITAN阈值基于GD 1—2减少自然背景变异，支持早期预警和适应性管理；结果支持间隙动物作为生物指示，因间隙带直接受出流减少影响且可能早于底栖群落响应；Kolmameter快速微创，可提高监测效率、标准化和客观性；出流减少不直接测淤填或细泥沙含量，但超过阈值提示敏感类群功能栖息地丧失，可作为优先恢复和评估管理效果的参考点。

研究结论翻译：受人类影响最小的粗基质溪流中间隙群落对河床淤填——此处以Kolmameter出流减少为代理——呈清晰阈值响应，出流减少超过30%后发生群落重组。敏感类群在该阈值以上下降，而更耐受类群持续存在，在德国西南部和奥地利粗沉积物溪流中形成一致的群落尺度格局。类群特异性响应表明沿出流减少梯度敏感性存在差异。总体而言，出流减少是间隙群落变化的主要驱动因素，而温度、pH、电导率和间隙氧等环境参数仅起次要作用。因此，30%阈值提供了稳健、具生物学依据的间隙栖息地泥沙诱导压力指标。当河流中至少半数样本超过该TITAN得出的出流减少阈值时，退化最易被检测到，这为早期生态损害检测提供了稳健基准。