河流屏障移除后结构-功能生态响应的空间不对称与时间动态：全球meta分析及BACI设计启示

《Journal of Applied Ecology》：Spatially asymmetric and temporally dynamic ecological responses to river barrier removal: A global meta-analysis

　　

河流屏障移除（river barrier removal）在全球范围内日益被实施，以恢复河流连通性和生态系统完整性，然而生态响应仍高度多变，限制了基于证据的恢复规划与监测。为此，研究人员开展了一项全球meta分析，量化河流生态系统在屏障移除后的结构性和功能性响应（structural and functional responses），并明确考虑相对于原屏障的空间位置、时间动态和研究设计。研究人员使用多水平meta分析模型（multilevel meta-analytic models）整合了来自29项研究的803个观测，通过随时间变化的轨迹（BACI方法，before-after control-impact）以及移除后相对于对照河段的偏离（CI-after比较，control-impact after）评估响应。生态响应显示出清晰模式，但强烈依赖情境。结构指标在不同类群和河段位置间呈相反响应；例如，基于轨迹的效应量显示鱼类丰富度（fish richness）在下游下降（LRR = ?1.37），而大型无脊椎动物丰富度（macroinvertebrate richness）在原蓄水河段增加（LRR = 2.36）。功能响应较弱且一致性较低，可能既反映数据有限，也反映基于性状的指派指标相对于生态系统过程直接测量占主导。受影响河段与对照河段之间的生态差异随时间逐渐减小，模型预测轨迹在原蓄水河段约12个月、下游约60个月后接近对照条件，尽管不同河流系统的恢复轨迹仍高度可变。推断结果还强烈依赖研究设计，轨迹法效应量（BACI）比CI-after设计得到的较简单移除后效应量更一致地检测到可归因于屏障移除的生态变化。综合与应用：屏障移除启动生态恢复，而非立即实现恢复，其响应受空间背景、初始条件和移除后时间塑造。因此，监测应超出移除后第一年，并在可行时持续至少5年，尤其在下游河段，同时纳入原蓄水河段和下游河段，并优先采用BACI设计。结构指标，特别是鱼类和大型无脊椎动物，应通过生态系统过程的直接测量加以补充，如代谢（metabolism）、养分吸收（nutrient uptake）和有机质分解（organic-matter decomposition），以弥补当前证据基础中的主要空白，并为恢复成功提供更全面评估。

1 INTRODUCTION
河流生态系统在数百年间被坝、堰等人工屏障广泛改造，用于灌溉、水电、蓄水、航运和防洪，导致河流网络普遍碎片化。屏障被界定为改变流态、改造物理栖息地并形成上游蓄水河段的人工结构。屏障形成上下游相反环境：上游蓄水段流速降低、滞留时间增加、细沙沉积，呈静水特征；下游常缺沙且流态改变。这些物理环境作为生态过滤器塑造群落和过程。大量屏障已失效或不再发挥社会经济功能，移除成为恢复纵向连通性和生态系统完整性的策略，并在政策推动下快速增加。然而，生态响应高度可变，受原蓄水段与下游不同环境变化、时间动态、结构/功能指标及研究设计影响。研究人员提出四个假设：（i）移除在原蓄水段和下游均引发可测响应，下游短期干扰更强，原蓄水段发生静水-流水重组；（ii）结构和功能指标响应幅度/方向不同；（iii）受影响与对照河段差异随时间减小；（iv）前后-对照-影响（BACI, before-after control-impact）比移除后对照-影响（CI-after）提供更强、更一致推断。

2 MATERIALS AND METHODS
2.1 文献检索与研究选择
2025年4月在ISI Web of Science检索英文同行评议文章，无出版年份限制，并筛查参考文献。初获773篇，按标题、摘要和全文逐步筛选。纳入标准：比较受影响与对照河段评估屏障移除效应；报告至少一个结构或功能生态指标；提供均值及离散度/样本量以估计效应量。排除缺乏合适指标、设计不适或统计信息不足的研究。数据均为已发表数据，无需伦理批准。

2.2 数据提取与指标分类
数据从正文、表格、作者或WebPlotDigitizer数字化提取。记录屏障身份、地理位置、出版年份及屏障高度和建造年份等特征。响应指标覆盖移除后对照河段（自由流动、上游未受屏障影响）与受影响河段（原蓄水段和下游）；如有移除前数据也提取。移除后时间以月记录。指标先分结构指标（描述群落组成、丰度或多样性）和功能指标（再分性状指标与生态系统功能指标）。结构指标按鱼类、大型无脊椎动物和生物膜分组；性状指标按运动策略、功能摄食群、底质偏好、流水偏好和耐受性状分组；生态系统功能指标包括凋落叶分解等过程。多个响应变量或多个采样时间均作为独立观测。

2.3 效应量与方差计算
效应量采用对数响应比（LRR, log response ratio）。CI-after效应量量化移除后受影响河段相对对照河段的偏离；BACI轨迹效应量在移除前数据可用时计算，量化受影响河段相对对照的时间变化。LRR为零分别表示无差异或变化等同，负/正值表示受影响河段较低/较高或相对下降/上升。方差按报告标准差和样本量计算并用于加权；标准差缺失时用数据集水平变异系数中位数近似，以避免排除报告不全研究。BACI方差假设站点内相关系数ρ=0.5。

2.4 统计分析
分析在R中完成，使用metafor包。CI-after和BACI分别分析。多水平meta回归按指标类型和效应量拟合，调节变量包括河段位置、目标变量和响应指标；屏障身份作为随机截距处理非独立性。用Akaike信息准则（AIC）比较候选固定效应模型，限制于至少3个观测的组合。最优模型用限制性最大似然（REML）重新拟合，预测效应量和95%置信区间使用CR2聚类稳健方差估计和Satterthwaite自由度。时间轨迹分析以LRR_CI-after为响应，将时间自移除后（月）作log₁₀转换并作为连续预测变量，结合结构/功能指标和河段位置选模。探索性分析将屏障高度分<5 m、5–15 m、>15 m三类，检验是否解释移除后响应。比较BACI与CI-after配对效应量差值和线性关系。敏感性分析排除标准差近似观测，漏斗图和回归检验评估小研究效应和发表偏倚。

3 RESULTS
3.1 数据集概况
最终数据集含23个河流系统的29项研究，803个效应量，其中结构指标480个、功能指标323个。44.2%来自原蓄水段，55.8%来自下游。68.5%效应量来自同时包含移除前后数据的研究，31.5%仅基于CI-after设计。下游比较多采用BACI设计（83.5%），原蓄水段有无移除前数据约各半。研究分布于北美、欧洲和亚洲七个国家，美国占一半以上。2012年后出版明显增加，75.9%研究发表于此后；早期侧重结构指标，功能指标自2021年以来增多。

3.2 结构响应
结构响应在类群、河段位置和研究设计间高度异质。BACI和CI-after模型均支持河段位置×目标生物交互，响应指标嵌套于目标生物。残差异质性仍高，但纳入调节因子改善模型。鱼类总体负响应，尤其下游；移除后效应量显示持续负偏离，原蓄水段更接近对照；轨迹分析显示下游负变化更强。多数观测在移除后3年内，反映早期重组。大型无脊椎动物总体正响应，原蓄水段轨迹显著增加，移除后效应量接近零，表明趋近对照；下游弱或中性。生物膜响应弱且不一致，移除后12–24个月无明显偏离，提示较快稳定，但研究少。总体结构响应可测但高度可变，取决于分类群和空间位置。

3.3 功能响应
功能指标响应较结构指标更弱、更不一致。BACI模型支持河段位置和功能类别的加性效应，CI-after模型支持二者交互。多数性状指标轨迹效应量接近零，少数显著：下游游泳者增加，原蓄水段耐受鱼类下降；移除后原蓄水段穴居者密度增加、滤食者比例降低。原蓄水段轨迹数据少。生态系统功能指标研究极少，凋落叶分解是唯一可定量分析的过程，移除后与对照无显著差异。功能响应证据有限，可能反映功能冗余，也可能反映经验覆盖不足。

3.4 移除后生态恢复的时间轨迹
两类指标和两个空间背景均呈恢复轨迹，受影响河段随时间趋近对照。下游初始负响应更强，收敛更慢；时间自移除后是强预测因子。模型预测在原蓄水段约12个月、下游约60个月接近零。功能指标初始偏差较小。尽管如此，残差异质性仍大，说明恢复轨迹在河流系统间变异很大。探索性分析未发现屏障高度类别<5 m、5–15 m、>15 m间移除后响应显著差异，不支持屏障高度作为一致预测因子。

3.5 研究设计比较
CI-after与BACI配对效应量正相关（n=457，r=0.54），方向大体一致，但CI-after效应量普遍低于BACI（均值差?0.239，斜率0.42，R²=0.29），倾向于低估正负轨迹效应。结构指标差异大于功能指标。仅66.7%配对比较方向一致，约三分之一符号相反。漏斗图无清晰不对称，回归检验无一致小研究效应证据，表明发表偏倚不太可能实质性影响结果。

4 DISCUSSION
4.1 概述与综合
本全球meta分析表明，河流屏障移除启动可测生态变化和恢复，但响应既非均一也非即时。三大模式：响应随相对于屏障的位置变化，但方向和幅度取决于指标；结构指标响应强于功能指标，同时暴露功能过程和生态系统过程证据缺口；响应呈清晰时间轨迹，受影响河段逐步趋近参照条件，但恢复速率和幅度在河流系统间差异大。结果强调移除是生态重组过程，而非瞬时恢复，受局部栖息地遗留、干扰强度、生物性状和地理背景塑造。

4.2 空间结构化生态响应
结果支持第一假设。所有指标和空间背景均检测到响应，很少中性，表明移除持续改变生态系统结构和功能。移除恢复纵向连通性和流水水文地貌，同时干扰蓄水段和下游既有条件。下游鱼类的短期负响应应理解为更广泛生态响应的一部分，而非恢复失败。响应表现为改善、延迟恢复或暂时下降，取决于生态遗留、泥沙动态、栖息地重构、扩散限制和局地环境。

4.3 结构与功能响应差异
结果支持第二假设。结构指标响应更强、更一致，但分类群间不同。鱼类多呈负或中性短期响应，下游尤甚，反映栖息地不稳定、泥沙移动和位移。大型无脊椎动物多正响应，原蓄水段尤明显，与快速再定殖、短生活史和流水条件环境过滤有关。功能指标较弱，可能因功能冗余；少数显著响应与恢复流水条件的水文地貌变化一致。当前证据基础不平衡：功能指标多为性状指标，生态系统过程直接测量少，凋落叶分解是唯一充分代表过程，限制更广泛推断。

4.4 时间轨迹与研究设计的作用
结果支持第三和第四假设。移除后早期响应强，常为负，随后随水文地貌稳定而趋近对照；下游恢复更慢。BACI比CI-after更清晰，因为考虑移除前基线。CI-after虽捕捉方向性恢复信号，但常低估效应量，且约三分之一配对方向不一致，原因在于未考虑受影响河段移除前可能与对照不同。屏障高度未显示一致差异，可能被河谷坡度等混淆；需更全面报告屏障高度、年龄、泥沙储存、移除策略和地形背景。

4.5 对河流恢复评估与监测的启示
建议将屏障移除作为动态过程评估。监测应覆盖原蓄水段和下游，并超出移除后第一年；在可行时至少持续5年，尤其下游，但12和60个月为模型预测而非通用阈值。优先采用BACI设计；若无法获得移除前监测，CI-after结果应谨慎解释。结构指标中鱼类和大型无脊椎动物互补，适合评估变化；功能评估应补充直接过程测量，如生态系统代谢、养分吸收和有机质分解。改善屏障特征报告有助于识别恢复成功的物理驱动因子。在大型恢复政策推动下，合适的时间尺度、空间覆盖、指标和实验设计对适应性管理和最大化生态效益至关重要。