在景观中何处采取措施以最小化河流系统的沉积和污染：一种多目标启发式方法

《Land》：Where to Act in the Landscape to Minimize Sedimentation and Contamination of the River System: A Multi-Objective Heuristic Approach

　　

造林可以减轻水、沉积物以及溶解或吸附的污染物向河流系统的输出，但确定有效的干预地点需要考虑多个与流动相关的标准及其空间相互作用。本文提出了一种多标准启发式方法，扩展了CAMF（基于元胞自动机的启发式最小化流动方法，Cellular Automata-based Heuristic for Minimizing Flow），该方法最初设计用于从栅格化景观中选择干预单元，以最小化目标地点的沉积物产量。研究人员将距离理想点（DIST2IP，Distance-to-Ideal-Point）算法集成到CAMF中，从而能够选择干预措施可以同时最小化两种或更多流动的单元。多标准CAMF被事后应用于日本福岛县受放射性污染的Niida河流域，以识别去污区内1000个单元，在这些单元上立即进行造林将最大程度地减少沉积物和残余¹³⁷Cs的输出。由DIST2IP选择的1000个最佳单元，占去污单元的4%，将减少22%的沉积物输出和6%的¹³⁷Cs输出。所选单元位于两个单标准优化所识别的单元的并集内，并且主要靠近水体，证实了在流动路径连接到河流系统之前阻断它们是最有效的。

**论文解读：基于多目标启发式优化的景观干预选址——以福岛Niida河流域沉积物与¹³⁷Cs污染缓解为例**

**研究背景与问题**
2011年福岛第一核电站事故导致超过2.7 PBq的放射性铯（¹³⁷Cs）沉降于景观表面，¹³⁷Cs作为长期暴露的关键放射性核素，与土壤颗粒紧密结合并随沉积物迁移，造成下游及近海区域持续污染。当地政府虽启动了去污计划（移除表土及植被），显著降低了表面污染，但此举加剧了土壤侵蚀，导致低污染沉积物进入河流系统，例如Niida流域在2016年悬浮沉积物负荷较治理前（2013年）增加了237%。造林作为一种快速恢复植被覆盖的措施，可有效减少侵蚀及沉积物和相关污染物的输出。然而，由于景观的空间异质性和空间交互作用（即某地点的干预会改变下游地点的状态和行为），干预措施的有效性高度依赖选址，且需在预算约束下最大化效益。已有单目标优化方法（如CAMF，基于元胞自动机的启发式最小化流动方法）仅能针对单一目标（如沉积物产量）进行选址，无法同时权衡沉积物与¹³⁷Cs输出。因此，研究人员开发了一种多目标启发式方法，以确定能同时最小化沉积物和¹³⁷Cs输出的最优造林地点。

**研究内容与结论**
本研究将距离理想点（DIST2IP）算法集成到CAMF中，提出多目标优化框架，并事后应用于日本福岛县受放射性污染的Niida河流域。通过模拟造林干预，识别去污区内1000个单元（约占去污单元的4%），评估其对沉积物输出和¹³⁷Cs输出减少的贡献。结果表明：DIST2IP选择的最优单元可减少22%的沉积物输出和6%的¹³⁷Cs输出；所选单元位于两个单目标优化（分别以沉积物和¹³⁷Cs为目标）所选单元的并集内，且主要靠近水体，证实阻断流路行前连接河流系统是最有效的策略。该研究发表在《Land》期刊，为受污染流域的多目标环境规划提供了实用妥协策略。

**关键方法**
研究采用以下关键技术方法（样本队列来自Niida流域，数据来源包括日本地理空间信息局的数字高程模型（DEM，Digital Elevation Model）和Worldcover 2020土地覆盖产品）：
1. **沉积物产生与传输模型**：使用修正通用土壤流失方程（RUSLE，Revised Universal Soil Loss Equation）计算年均土壤流失量；采用多流向（MFD，Multiple Flow Direction）模型和传输能力（TC，Transport Capacity）方法模拟沉积物传输与累积。
2. **¹³⁷Cs传输模型**：基于固体夹带系数和初始¹³⁷Cs沉降数据，估算沉积物中¹³⁷Cs浓度，并假设混合均匀，通过质量平衡计算累积活度。
3. **单目标优化（CAMF）**：迭代评估候选单元（去污区）的造林对沉积物产量或¹³⁷Cs产量的减少效果，按影响排序选择最优单元，显式考虑空间交互。
4. **多目标优化（DIST2IP-CAMF）**：在每次迭代中，将候选单元的归一化沉积物减少量（SYR，Sediment Yield Reduction）和¹³⁷Cs减少量（CYR，Cesium Yield Reduction）与理想点（理论最大值）的加权欧氏距离作为选择标准，权重相等。

**研究结果**
- **沉积物产量减少与¹³⁷Cs产量减少对比**：通过单目标（仅优化SYR或仅优化CYR）与DIST2IP多目标三种策略，分别造林100、200…1000个单元，比较累积SYR和CYR。基于SYR的选择始终获得最高SYR但最低CYR，反之亦然。DIST2IP提供中间结果，SYR和CYR均接近各自最大值，避免了单目标优化的严重失衡。例如，当选择100个单元时，DIST2IP的CYR仅比基于CYR的单目标优化低约0.28%，但SYR显著高于基于CYR的选择（约8-12%）。
- **空间分布**：DIST2IP所选单元（1000个）主要靠近河流，中位距离河流较近，且集中在中低海拔（约200-600 m a.s.l.），但高海拔（约600 m）也有显著集群。初始沉积物产量高的单元并非全部被选，说明沉积物路径和连通性同样重要。值得注意的是，DIST2IP所选单元完全包含于单目标优化所选单元的并集内，无独立于两者的单元，表明该案例中多目标妥协未引入新地点。
- **敏感性分析**：针对去污效率（¹³⁷Cs浓度降低40%、20%对比50%）和表土混合深度（2、3、4 cm对比5 cm）进行敏感性测试。结果表明，初始CYR和CYR变化较小，且所选100个单元集合完全一致（仅排序略有差异），说明CAMF-DIST2IP的输出空间格局对这些输入参数不敏感。当混合深度设为2 cm时，模拟初始CYR与观测值吻合。

**讨论与结论**
讨论部分指出：空间靶向干预需考虑场地内和场地外效应，去污实践虽降低了局部污染，但增加了下游侵蚀和污染输出。造林作为干预措施时，靠近河岸的单元因能拦截上游沉积物而效率最高，这与既有研究一致。多目标优化中，DIST2IP提供近最优平衡，且可通过调整权重反映不同管理优先级。方法学上，RUSLE和简单传输函数虽为简化，但用于识别模式优于精确值；CAMF可扩展至其他污染物（如氮、磷）或更多目标（如成本）。结论部分强调：在受¹³⁷Cs污染的流域（如Niida），基于多目标优化的靶向及时造林能同时缓解下游沉积物通量和¹³⁷Cs传输，增强长期生态安全。DIST2IP作为实用妥协策略，具有广泛适用性。