高分辨率集成水文建模揭示河流网络扩张与收缩的滞回动力学

《Water Resources Research》：Hysteretic Dynamics of Stream Network Expansion and Contraction Revealed Through High-Resolution Integrated Hydrological Modeling

　　

随着景观在湿润与干燥状态之间循环，河流网络随之扩张与收缩，驱动物质输出与生物地球化学转化，对水质和生态系统功能具有重要影响。由于直接测量活动河网的范围（L）具有挑战性，径流（Q）常被用作替代指标。然而，Q-L关系仍存在不确定性，尤其是在直接观测耗时费力的大型流域中。研究人员利用能够以高空间和时间分辨率解析淹没模式的基于过程的集成水文模型，研究了面积为75 km2的流域及其子流域的河网范围动态。在无需流域特定校准的情况下，使用ATS集成水文模型进行的模拟与观测的日径流数据以及观测的河流持续性状态（常流vs.非常年性）吻合良好。模拟结果表明，仅在小型流域中存在唯一的、同步的Q-L关系。具体而言，同步的Q-L数据因滞回效应而呈现显著离散，且该离散程度随流域面积增大而增加。然而，当转而考虑Q与L的异步事件尺度峰值时，出现了一种无滞回效应的稳健关系。该事件尺度的异步关系在模拟的峰值范围内可较好地近似为幂律关系，并在不同大小的嵌套子流域间保持一致，且不受前期湿度、季节和风暴细节的影响。在基流条件下，Q-L标度关系同样可用幂律近似，但指数更大，表明径流产生机制不同。这些结果为理解更大、更复杂流域中的河流扩张与收缩提供了更稳健的框架。

论文解读：高分辨率集成水文建模揭示河流网络扩张与收缩的滞回动力学

一、研究背景与问题

全球范围内，大量河流（估计占60%–80%）并非持续流动，而是在降雨事件、季节变化及长期降水、蒸散发和地下水水位趋势的驱动下，在湿润与干燥状态之间转换。这种动态性使得河流网络成为高度动态的系统，其扩张与收缩直接影响河流向大气的净气体通量、景观与河流系统之间的连通性，以及流域内水和物质的停留时间，对流域生物地球化学、淡水生态学、农业区养分管理及水政策与法规具有根本性意义。

径流（Q）与活动河网长度（L）之间的关系是理解河流网络动力学的关键。早期研究将L视为流域水文响应的首要控制因素，但随时间推移，L日益被视作动态诊断变量。最常用的Q-L关系为幂律形式，其指数代表网络弹性系数，野外实测的弹性系数通常在0.1–0.4范围内。然而，Q-L关系仍存在重大不确定性：部分建模研究发现其偏离幂律关系；Gamma函数被提出作为幂律的替代形式；且Q-L关系可能并非一一对应，同一径流可对应多个活动河网长度。事件尺度的滞回现象已被注意到，但在以往小流域研究中较为微弱。最大的不确定性在于Q-L关系对时间和空间尺度的依赖性——传统调查方法劳动密集且局限于小流域（一般<5 km2），遥感与无人机受冠层覆盖限制，传感器仅能指示局部湿润状态而无法获取总河网长度。因此，亟需一种能够在更大时空尺度上显式模拟Q-L关系的方法。

二、研究内容与意义

研究人员利用ATS集成水文模型，对位于田纳西州阿巴拉契亚山谷与山脊褶皱冲断带的东福克杨树溪（EFPC）流域（75 km2）及其子流域开展了高分辨率模拟研究。该模型能够显式解析窄小的季节性河道，无需流域特定校准即可模拟河网扩张与收缩动态。研究揭示了事件尺度滞回现象随流域面积增大而加剧，并提出了基于异步事件峰值的稳健Q-L幂律关系，为理解更大、更复杂流域中的河流扩张与收缩提供了更稳健的框架。

三、关键技术方法

研究人员采用了以下主要技术方法：使用ATS（高级陆地模拟器）集成地表-地下水文模型，该模型基于3D Richards方程求解变饱和地下水流，并以2D扩散波近似求解地表径流，通过通量和压力连续性实现地表-地下耦合；采用流线对齐混合多面体网格技术，沿河道中心线布置细长四边形单元，其余区域采用粗化三角剖分，从而高效高分辨率地模拟河道走廊动态；利用Watershed Workflow工作流工具包构建网格并映射地表与地下数据集；气象强迫数据来自Daymet（1 km分辨率），模拟时段涵盖2012–2020水文年，采用三步自旋-up流程（20,000天稳态运行、15年循环稳态、8年瞬态模拟），以小时分辨率进行模拟。模型评估采用修正Kling-Gupta效率（mKGE）指标，并与两个流量站的日径流观测数据及51个位置的野外流态分类进行对比。

四、研究结果

4.1 模型评估与置信度建立

模拟与EFK5.4位置观测径流的mKGE值为0.84（2013–2019水文年），与EFK16.2位置的mKGE值为0.72，表明模型性能优异。模型正确识别了80%的常流点位和71%的非常年性点位，考虑到这些点位多被有意选在常流与非常年性过渡区域附近，模型表现尤为突出。

4.2 全分布式集成水文模型揭示河网持续性的时空动态

模拟结果显示，NHDPlus-HR绘制的142 km流线中约15 km（一般位于源头区）在模拟期间从未流动，表明活动河网长度不受网格单元布置约束，而主要由地貌控制和水文条件决定。河网在12月至4月期间普遍更为扩张，在6–10月干旱期基流条件下收缩至约47 km，冬季风暴事件期间可达约128 km的扩张上限。连续活动河网平均占总活动河网的81%（范围75%–87%）。一阶河流主导整体扩张-收缩动态，二阶间歇性河段主要影响季节性趋势，三阶以上河段基本持续活跃。

4.3 事件尺度滞回及其对Q-L关系的影响

高分辨率模拟揭示了全流域清晰的顺时针事件尺度滞回现象，滞回指数中位数为0.46，四分位距为0.12。滞回源于L与Q对降水强迫的响应时间不同——河网扩张主要通过沿短坡面的浅层地下水流在活动河网远端发生，而径流是整个流域的集成响应，响应更为缓慢。滞回导致同步Q-L关系呈现显著离散，幂律拟合的预测价值被削弱。通过分离事件峰值与基流时段，研究人员发现事件驱动的异步Q-L峰值关系呈现更紧密的幂律标度（弹性系数0.19），而基流时段的同步关系也符合幂律但指数更高（0.37），两者在统计上显著不同，且均不随季节变化。

4.4 流域面积对Q-L标度、滞回及拟合质量的影响

对EFPC全流域及三个嵌套子流域（EFK5.4、Bear Creek、Mill Branch）的分析表明，较小流域呈现更清晰的幂律标度和较弱的滞回，较大流域则因河网激活与出口径流之间的时间滞后加剧滞回，导致同步Q-L数据离散度增大。相比之下，基于事件尺度的峰值-峰值Q-L关系在各流域间显著降低了离散度，弹性系数在三个较大子流域间范围为0.18–0.22，而最小子流域（MB）为0.37。滞回指数分布随流域面积增大而系统性增大，表明滞回强度与空间尺度密切相关。

五、讨论与结论

研究结果表明，同步的Q-L唯一关系仅存在于小型流域。在HUC12尺度（40–160 km2）下，径流与河网范围对降水响应的不同滞后时间产生了事件尺度滞回，导致同步Q-L数据广泛离散。研究将事件驱动的扩张与收缩同基流驱动的河网激活区分开来，识别出异步事件峰值Q-L之间稳健的幂律关系，该关系在不同大小的嵌套子流域间保持一致，不受前期湿度、季节和风暴细节影响。基流条件下Q-L标度指数更大，表明径流产生机制不同。研究结论指出，事件峰值与基流的区分性标度规律为理解不同水文条件下的河网动态提供了更清晰的框架，而同步Q-L数据中的滞回效应可能掩盖结构性信号并限制可解释性。尽管研究指出幂律关系在有限范围内可能被过度解读，但低离散度的异步峰值Q-L关系本身是重要发现，不同弹性系数值仍可靠地指示了不同的河网激活动力学。