基于加权可用面积（WUA）时间序列的鱼类栖息地时间动态指标：气候变化情景下地中海调节河流的环境流量评估

《Ecological Indicators》：Indicators from fish habitat time series for climate change and environmental flows assessment

　　

基于加权可用面积（weighted usable area，WUA）的指标可量化鱼类栖息地时间动态。指标捕捉幅度、变异性、持续时间、频率和持续性。该框架支持基于栖息地的环境流量（environmental flow，e-flow）评估。关键生活史阶段（critical life-stage periods，CLPs）可增强对栖息地变化的检测。幼年barbel表现出最强的气候相关栖息地响应。

论文解读：研究人员针对气候变化下地中海调节河流鱼类物理栖息地的时间动态，构建了由WUA时间序列衍生的指标框架，并将其用于环境流量（e-flow）评估。研究背景是淡水生态系统生物多样性下降严重，河流栖息地受水文改变、破碎化、污染、过度开发与外来种等多重压力影响。自然流态范式强调流量大小、频率、持续时间、时机和变化率决定河流生态完整性；指示水文改变（Indicators of Hydrologic Alteration，IHA）可描述水文改变并支持e-flow评估，但不直接量化流量改变如何转化为物理栖息地可用性。IFIM/PHABSIM框架通过WUA连接流量与目标物种及生活史阶段的物理栖息地，但WUA时间序列很少被组织为紧凑、类似IHA、可跨物种、尺寸类和关键生活史阶段比较的指标套件。研究人员以西班牙Jucar河中游Madrigueras河段为案例，该河段受上游水库、农业需水和地表水-地下水交换影响，选择本地Eastern Iberian barbel（Luciobarbus guiraonis Steindachner, 1866）L1、L2、L3三个尺寸类和入侵bleak（Alburnus alburnus Linnaeus, 1758），评估CMIP6 SSP1–2.6、SSP3–7.0和SSP5–8.5下近期与中期栖息地指标变化。研究意义在于将栖息地评估由静态流量-栖息地曲线推进为对限制条件“多久、多频繁、多连续”出现的时间显式诊断。

关键技术方法不超过250字：研究人员采用集成建模链，包括偏差校正CMIP6气候投影、TETIS v9.1水文模拟、STIG-CROPROD水资源系统模型、SEFA水力学模拟和模糊规则系统（fuzzy rule-based systems，FRBS）栖息地建模。气候数据来自ISIMIP3b偏差校正CMIP6，选用五个全球气候模型（global climate models，GCMs）和ERA5-Land参考数据，空间分辨率约10 km。研究区为Jucar河Madrigueras段Los Frailes控制点附近200 m河段，代表地中海调节河流，含18个横断面。目标鱼为Eastern Iberian barbel三个尺寸类和bleak。关键生活史阶段（critical life-stage periods，CLPs）依据已发表繁殖和生态证据确定。水力学模拟50个流量（0.08–45 m³ s^-1），以水深、平均流速和cover建立WUA–Q关系。调节月流量经线性插值转为月WUA序列，计算六个指标：WUA_mean、WUA_min、WUA_CV、WUA_NM30、WUA_CM30、WUA_Ev30，阈值取30% WUA_max，并用Cliff's delta与5000次bootstrap评估效应。

研究结果：

3.1 WUA–Q relationships and monthly habitat distribution：通过水力学与FRBS栖息地模拟得出，WUA–Q关系因物种和尺寸类而异。bleak和barbel L1呈单峰响应，峰值低于3 m³ s^-1，分别为4652 m²和5075 m²，高于13 m³ s^-1后近零。barbel L2呈饱和响应，低流量陡增，10 m³ s^-1以上仍大于5400 m²。barbel L3呈延迟响应，3 m³ s^-1以下低于35 m²，45 m³ s^-1时1539 m²。CLP内月栖息地可用性差异超过一个数量级，L2常高于3200 m²，L3很少超过300 m²。

3.2 Temporal evolution of habitat indicators：研究人员在各自CLP内计算指标并比较历史期与未来期。栖息地赤字定义为月WUA低于30% WUA_max。

3.2.1 Eastern Iberian barbel：L1指标分布变化小，WUA_min变化最明显，历史51%年份最干月零可用栖息地，中期SSP1–2.6、SSP3–7.0、SSP5–8.5分别增至70%、90%、73%；WUA_NM30中期一致增加，赤字2个月年份增加9%；WUA_CM30无一致模式；中期97%年份仅1次赤字事件。L2历史栖息地高，47%年份处于最高WUA_mean四分位；未来最低四分位增加，SSP5–8.5中期达67%；WUA_min变化最尖锐，历史所有年份高于阈值1654 m²，未来所有组合出现零，中期分别20%、37%、73%；WUA_CV增加；历史无赤字，中期无赤字年份分别为80%、63%、27%，最高排放下30%年份有4个赤字月；WUA_Ev30中期1次事件分别为20%、37%、73%。L3的WUA_mean持续低，最大633 m²，最高四分位历史40%，中期降至37%、13%、3%，最低四分位由11%升至17%、23%、67%；WUA_min从未超过阈值462 m²，但零栖息地年份历史无，中期为3%、13%、47%；WUA_CV增加；阈值指标历史已受限，WUA_NM30四个月赤字历史54%，中期为50%、87%、97%；WUA_CM30类似；WUA_Ev30所有年份均为1次事件，无区分。

3.2.2 Bleak：除WUA_min外，各指标分布变化小。WUA_mean最低四分位历史29%，中期17%–30%；最高四分位历史34%，中期20%–23%，无排放梯度。WUA_min历史40%年份零，中期63%、77%、77%，近期40%、44%、52%；高于阈值1396 m²的年份不超过4%。WUA_CV中期略降。阈值指标变化小，WUA_NM30历史1、2、3个月赤字分别40%、31%、23%，中期2个月最常见，为60%、47%、50%。WUA_Ev30历史97%年份单事件，中期93%–100%。

3.3 Changes in habitat indicators across climate-change scenarios：通过Cliff's delta比较历史与未来，正效应表示更限制性栖息地。barbel L2的36个比较均至少小效应，28个显著，含7个大效应和15个中效应；L3有10个显著，3个大效应限于幅度指标，18个可忽略；bleak和L1分别仅2个和1个显著，30个和32个可忽略。L2响应最强，WUA_mean由中期SSP1–2.6的δ=0.331增至SSP5–8.5的δ=0.819，近期SSP5–8.5已δ=0.536；WUA_min由δ=0.466增至δ=0.757；三个阈值指标相同，δ=0.192至0.733。L3最大效应在幅度指标，WUA_mean中期SSP3–7.0 δ=0.477、SSP5–8.5 δ=0.777；WUA_min中期SSP5–8.5 δ=0.651；WUA_Ev30δ=0。Bleak仅WUA_min在中期SSP3–7.0和SSP5–8.5显著，δ=0.410和0.309；WUA_Ev30|δ|≤0.051。

讨论总结：研究人员认为，物理栖息地是随时间变化的信号，而非单一静态关系；WUA–flow曲线仍能表示栖息地幅度，但传播调节流量序列最有信息之处在于静态适宜性相似的尺寸类发生分化。L2历史可用性最高且未来退化最强，L3已接近赤字阈值，进一步退化主要在幅度，L1均值和变异性无方向变化而最小值明显移动。CLP限制使年均会合并的维度得以分离。不对称气候响应不在于是否响应，而在于哪个维度响应及传播多远；L2是唯一全指标传播的类，其WUA–Q单调且不归零，最小值零对应无流，变化是夏季窗口从历史无赤字转为大多数年份赤字，暴露是下界。bleak和L1因WUA–Q单峰，低流量和高流量两端均可低于阈值，最小值零可因干涸或高流量，指标含义不同。管理层面上，指标可直接用于现有e-flow框架，无需额外野外调查，并有助于识别本地类最长期限制；L2结果表明历史基线可能不再代表未来，最小月可用性应优先监测。局限包括30%阈值影响分辨率、不同CLP差异描述的是相关期暴露而非共同时间质量、曲线来自单一河段、鱼调查数据不连续难以验证种群响应。

结论翻译：本研究表明，月度WUA时间序列可被凝练为可解释指标，用以描述地中海河段鱼类栖息地时间动态。通过将WUA–Q关系转换为在关键生活史阶段计算的六个互补指标，并依据已发表繁殖和生态证据，该框架将水文改变指标的逻辑扩展到栖息地领域，提供了对限制物理栖息地条件的幅度、变异性、持续时间、频率和持续性的紧凑描述。应用于气候变化情景后，指标揭示了物种和尺寸类间的显著差异。暴露集中在L2 barbel类，其CLP历史任何年份无栖息地赤字，而在最严重中期情景下多数年份出现赤字，所有36个比较均非可忽略。L3 barbel类在整个繁殖窗口已经受限，阈值指标接近上限，进一步退化主要通过幅度检测。L1 barbel类和入侵bleak中，响应限于CLP内最小栖息地可用性月份，在这两类中反映的是偏离低流量最优而非栖息地收缩。检测到的转变是从无栖息地赤字到L2 CLP中经常出现赤字，而非从偶发限制到持续限制；解析这一转变需要按物种和尺寸类定义关键生活史阶段。若干指标也随三种情景的辐射强迫变化，L2最一致，幅度指标中L3亦如此，表明即使在调节河段中水库运行和需水过滤气候信号，这些指标仍对排放水平保持敏感性。这些结果强调，在气候非稳态下评估e-flows时，应解析限制条件的持续时间和持续性；指标行为取决于尺寸类相对赤字阈值的位置及其适宜性关系形状，解释前应检查二者。该框架提供了可转移的一线诊断，用于评估气候驱动的物理栖息地风险，适用于调节、自然或天然化流态，并将通过生物监测数据验证得到强化。