世界大型河流系统的自相似临界性

《Journal of Hydrology》：Self-similarity criticality of the world’s large river systems

　　

• •
  提出自相似有效比（SVR）来量化河网中尺度不变组织特性的持续程度。
• •
  建立一个框架，将河网的自相似临界区与能量响应联系起来。
• •
  SVR与黄金比例的吻合标志着分层能量分布的平衡状态。
• •
  提供适用于洪涝评估、侵蚀诊断、泥沙管理和水电规划的筛选尺度工具。

引言

河流系统作为分层网络，支配着全球水循环、物质传输和生态过程（Carraro和Altermatt，2022；Rodriguez-Iturbe等，2009；Widder等，2014；Wollheim等，2022），在塑造地球表面形态和维持人类文明方面发挥着关键作用（Perron等，2012；Wang等，2022；Wang等，2024）。其结构的自相似性——即子系统在不同尺度上复制整体特征的 fractal模式——已被确认为连接地貌演化与水文学功能的根本特性（Rodriguez-Iturbe等，1994；Rodriguez-Iturbe等，2009），提供了一种强大的尺度桥梁工具，用于解读河网拓扑、流域形态、水文连通性以及沿河级层次的能量耗散。河级作为统一指标，将分支结构、坡度变化、径流汇积和河功率重新分配联系起来（Horton，1945；Strahler，1952；Strahler，1957）。当分支拓扑、河道几何、集水面积和水文汇积在各级河级间遵循稳定的标度关系时，局部河段尺度或子流域尺度的观测便可在连贯的全流域组织框架内进行解读。Horton定律（Horton，1945）和Strahler河级排序系统（Strahler，1952；Strahler，1957）长期以来为描绘自然河网格局（Wang等，2019）和通过分叉比、河长比和流域面积比来量化这种自相似性提供了框架。

然而，新近证据揭示了一种分化现象：自相似模式可能在较高河级处突然崩塌，这一现象被称为"自相似临界性"（Chen等，2019）。自相似性的崩塌代表一个关键转折点：河网脱离自组织的分层模式，转而受特定流域约束主导，包括洪泛区蓄积、岩性变化、构造活动、冰冻圈变化以及人类干预。越过这一转折点，基于标度关系的水文与地貌预测将可靠性下降，并伴随更强的空间异质性。这种崩塌对传统模型提出了挑战，凸显了解析控制河网组织的非线性阈值的必要性，而这一空白阻碍了对流域响应气候变化和人为压力的预测。

自相似临界性假说认为，河流系统通过自组织向能量最优构型演化，平衡侵蚀力与物质传输效率（Rodríguez-Iturbe和Rinaldo，1997；Yang和Song，1979；?ukowski等，2022）。虽然最优河道网络和耗散结构原理等理论能够解释低河级的规律性（Abed-Elmdoust等，2016；Sarker等，2025；Sarker等，2019），但它们无法解释为何高河级会偏离。近期研究将归因于地理气候驱动因子：坡度阈值调控能量耗散率（Cohen等，2018；Hooshyar等，2017；Stronga和Mudda，2022），而降水—径流耦合调节径流分配（Feng和Gleason，2024；Messager等，2021）。有趣的是，黄金比例（φ = $\sqrt{\text{5}}/\text{2}\text{-}\text{1}/\text{2}\approx \text{0.618}$），在某些生物和物理系统中是一种通用的优化度量指标（Coldea等，2010；Maor，2003），被广泛认为是由分层堆积和尺度不变分支驱动的自组织系统的最优状态，并且已被初步与河流蜿蜒和分叉角联系起来（Stolum，1998），暗示了一种内在的能量最小化原理（Rodriguez-Iturbe等，1992）。然而，整合这些多尺度控制的定量框架仍然缺失，阻碍了流域韧性的预测建模。

当前研究仍存在关键局限性。第一，自相似性的空间范围通常仅通过目视粗略估计，缺乏标准化度量指标，导致关键阈值的全局模式未得到解决。特别是，缺乏一个统一的无量纲指标来量化不同背景流域中河网尺度不变自相似性的持续程度。第二，尽管坡度和径流被公认为主要的能量调节因子（Gernaat等，2017），但它们对临界性的协同效应仍未被量化，特别是在气候变化引起的水文转变条件下。第三，φ在结构优化中的作用存在争议，在不同流域中的实证验证十分有限。解决这些空白需要一个统一方法，结合高分辨率地理空间分析、能量—坡度耦合模型和非线性动力学理论。

本研究利用HydroATLAS数据集对66条全球大型河流（流域面积≥100,000 km²，补充表1）的自相似临界性进行了调查，旨在：（1）建立一个通用指标来量化河网中尺度不变组织特性的持续程度；（2）揭示自相似性崩塌的地理气候和能量阈值；（3）解析φ在能量优化网络构型中的调控作用。通过将结构临界性与水—地貌驱动力及能量响应联系起来，本研究将河网自相似性从一种结构性观察转化为基于物理的诊断框架，可广泛应用于环境变化下的洪涝风险评估、泥沙调控和水能开发（Best，2019）。

章节摘要

研究区域

我们分析了66个全球大型河流流域（流域面积≥100,000 km²），总计覆盖6.19×10⁷ km²（占全球河流汇水区面积的66%），入海径流量为6.29×10¹¹ m³·a^-1（占全球淡水入海通量的91%）。这些流域横跨六大洲，展现出多样的全球水气候和地貌特征（补充表1）。

自相似区与黄金比例吻合

通过设定Horton自相似假设的接受阈值（R²≥0.90，p＜0.05），所有河流在全流域尺度上均表现出统计上的普遍自相似性（补充图2）。进一步以河级为自变量分析分叉比（R_B），以划定河网的自相似区（图2）。按大洲汇总，中低河级河流（河级1–4）观察到R_B的收敛（均值范围为4.33–4.55）。

三种自相似临界区的机制解释

所提出的SVR量化了河网等级中尺度不变组织仍然有效的比例，而CS、CQ和CSP则刻画了该标度域上界处的水—地貌条件和能量状态。当以坡度主导的能量释放被以径流主导的能量传输所替代无法再维持河功率的平滑分层再分配时，自相似性发生崩塌，其

结论

本研究通过将结构临界性与水—地貌驱动力及能量响应联系起来，为66个全球大型河流流域建立了一个自相似临界性框架。我们证明，河网中的自相似结构模式在低至中河级（1–5）稳健收敛，但在高河级（6–10）越过关键阈值后突然崩塌，受地貌和水文耦合控制。所提出的自相似有效比（SVR）量化了

CRediT作者贡献声明

Zhenqi Li：方法论、调查、形式分析、数据整理。Yichu Wang：写作—审阅与编辑、写作—原稿撰写、可视化、验证、资源、方法论、资金获取、形式分析、数据整理、概念化。Jinren Ni：写作—审阅与编辑、指导。

利益冲突声明

作者声明不存在已知的竞争性经济利益或可能影响本文所报告工作的个人关系。

致谢

本研究由国家自然科学基金资助[批准号 52394233、52479057、U2243201]。

Zhenqi Li|Yichu Wang|Jinren Ni