一条河少了水，未必只是天旱了
2016年5月，加拿大育空地区的 Slims River 水位突然下降。这里的水很大一部分来自 Kaskawulsh 冰川，过去经 Slims River 流入北侧的 Kluane 湖，最终进入育空河水系。变化之后，原来供应这条路线的冰川融水大量转入另一侧的 Kaskawulsh River，再经 Alsek River 流向海洋。4
两条河没有移动着互相碰撞，水也没有翻过一座完整的山。发生变化的是冰川末端的通路：退缩、融化和侵蚀打开了通往较低湖面的新出口。原来隔开的排水路线由此接通，上游来水被重新分配。
这是一次有当年水位记录、地形测量和影像可对照的河流袭夺。它提供了一个直观入口，但普通山地的分水岭往往以另一种方式缓慢调整。要理解两者，先要知道分水岭究竟是什么。
边界由地形维持，也会被侵蚀改写
想象雨落在同一条山脊两侧。沿地表流走的那部分水，一边进入甲河，另一边进入乙河。把最终流向同一出口的地表区域圈起来，就是这里讨论的汇水区；两区之间的分界就是分水岭。这个定义针对地表排水，不保证地下水也服从同一边界。
分水岭并不一定是全山最高的山峰。关键是它两侧的坡向把水分开。河道一旦改变周围的地形，这条线也可能改变位置。典型的渐进过程可以拆成四步：1 2
- 一侧河道先降低。 出口基准面下降、岩石更易侵蚀或径流条件改变，都可能使下切增强。基准面是所研究河段向下游连接时受到的高程约束，不一定等于海平面。
- 坡脚降低，坡面开始调整。 河流不是在山顶直接“拉”分水岭；它降低邻近坡脚，改变坡面坡度、物料输运和侵蚀。若近分水岭一侧地表降低得更快，在其他条件相近时，边界可向另一侧退移。
- 原来的归属发生改变。 原先流向乙河的一小片坡地，可能转而流向甲河。若侵蚀切通一个原属乙河的上游河段，变化会从一片坡地扩展为整条支流的改接。
- 新的汇水反过来改变侵蚀。 面积增大通常意味着能收集更多径流，但这一步依赖降雨、入渗和蒸发条件；更大的流量又可能改变河床下切。失去汇水的一侧则面对另一组条件。
这条反馈没有要求水向上游倒流。向上游传播的是河床降低和坡面调整的影响，水本身仍沿当时可通行的下坡路线流动。
也不要把“侵蚀强的一侧获得面积”写成无条件定律。若两侧抬升不同，分水岭附近地表的净高程变化要同时算抬升与侵蚀；若岩层倾斜，一侧先切穿软层、再遇到硬层，迁移方向还可能反转。构造造成的水平搬运又是另一项过程。2 8
把出口降低，边界会不会跟着变
用一个受控例子把前面的过程走完。本库采用 Landlab 2.11.0 的流路累积、Fastscape河道下切和线性坡面扩散组件，地形每一步都更新，流路也重新计算。9
河道部分采用常见的简化关系：
$E$ 是下切速率，$A$ 是上游汇水面积，$S$ 是顺流方向的坡度大小，$K$ 汇集侵蚀效率等条件。本例用 $m=1/2$、$n=1$。它把汇水面积当作径流代理，因此默认空间上相同的有效产流条件；它不逐颗模拟砂石撞击河床。1
地表高程变化再加上均匀抬升和坡面输运，可概括为：
$D\nabla^2z$ 表示物料沿坡输运对地形的平滑作用。河道下切只在汇水面积达到0.04平方千米的节点启用；这个门槛、线性坡面关系以及不让淹水节点下切，都是模型选择，并非所有真实河流遵守的规则。
模型区域东西长5.9千米、南北长4千米，网格间距100米，两出口位于左右边界中央，其余边界封闭。对称山脊先在均匀条件下演化100万模型年，随后各组从完全相同的地形开始比较。这个准备过程用于取得共同起点，没有据此宣称地形已严格达到稳态。
接下来50万模型年中，一组保持两个出口高程不变，一组让左出口以每年0.5毫米下降，另一组改降右出口。还有一组同样降低左出口，但关闭河道下切，仍保留坡面输运与流路重算。均匀抬升为每年0.2毫米，$D=0.1$平方米/年，$K=5\times10^{-5}$每年。这些全是示意参数，没有拟合任何自然流域。
最终结果如下，比例的分母是有产流面积的内部网格，总计22.62平方千米：
| 条件 | 最终流向左出口的面积比例 | 改变归属的网格数 |
|---|---|---|
| 两出口不变 | 50.00% | 0 |
| 左出口下降250米 | 53.18% | 72 |
| 右出口下降250米 | 46.82% | 72 |
| 左出口下降，但关闭河道下切 | 50.00% | 0 |
| 左出口下降，时间步从500年减为250年 | 53.23% | 73 |
左侧增加的面积约0.72平方千米。这里可以确认的是：在指定模型里，出口变化经下切与坡面响应改变了排水分界。关闭下切后没有留下相同的最终净转移，但不是整个过程中都没有变化：左侧比例曾短暂达到51.28%，该组还出现过最多58个淹水节点，面积归属包含洼地溢流路由的处理。因而不能把这条对照说成“下切是任何边界变化的必要条件”，更不能把短时路由变化都视为山脊侵蚀。正常下切的左右两组全程没有淹水节点。
每个记录时点都检查全部内部网格能沿流路到达两个出口之一，两个出口累计面积之和等于总产流面积。时间步减半的末态差一个网格，即约0.044个百分点，支持这里的变化方向和量级没有依赖原时间步偶然出现。但没有做空间分辨率收敛研究，不能拿53.18%当自然系统的精确预测。图中的阶梯变化也不等于已经模拟出一次真实规模的支流袭夺。
这套模型还有一道明确的边界：它没有泥沙储存、河床覆盖和冰川。因而它不能独自解释“多下雨为什么有时反而迁移得慢”，也没有重现开头的2016年事件。
多雨的一侧为什么不一定一直占上风
Bonnet 的2009年实验把均匀抬升的硅粉湿料放在受控降雨下。先用均匀降雨形成近似对称地形，再施加横向降雨梯度。此后主分水岭向较干的一侧迁移，多雨侧伸长，干燥侧的一些流域缩小。3
值得注意的是，缩小的流域不只是被旁边“吃掉”：其中一些还分裂成两个流域。作者根据连续地形测量与沉积位置提出，径流减少后活动河道变窄，原洪泛地的一部分不再受同样的侵蚀，在持续抬升中变成新的分界。地形中分裂和抬升的变化是观察结果；雨幕遮挡使河道内部过程不能全程直接观察，具体的河道变窄机制包含解释。3
这个实验说明，气候不对称能够重组地形，也说明流域数量并不只会通过吞并而减少。它没有给出把实验小时数换成自然年数的可靠比例。作者在方法中明确提醒，材料、过程和尺度无法完整缩比，因此这类装置用于理解过程，不能直接当作缩小的真实山脉。
再看泥沙。水流携带的颗粒撞击裸露河床，可以成为磨蚀工具；但沉积物覆盖床面时，又会把基岩隔开。磨蚀模型因此要同时考虑颗粒撞击、单次撞击造成的磨损，以及床面有多少仍裸露。只把雨量增加写成 $K$ 增大，会隐去这个竞争。6
这也不是说“沙多就一定保护河床”。粒径、搬运能力、覆盖厚度和悬移质是否仍撞击床面都会改变结果。Lamb 等的2008年模型纳入悬移颗粒后，就改变了只考虑跳跃床载时对高搬运强度的部分预测；连“颗粒开始悬浮后磨蚀怎样变”都依赖具体机制，不能靠一个通用的抛物线替代。6
2026年 Yan 等的双侧沙箱研究正沿着这个问题推进：其摘要报告，在所研究的构造沉降与供沙条件下，增加降雨可通过供沙、覆盖和沉积反馈减慢分水岭迁移。这是值得追踪的新约束。本篇只取得摘要、网页预览和书目，未核全文方法与重复实验，因而不采用其中精确的百分比，也不把它提升为“降雨通常抑制迁移”的一般规律。7
2016年的事件，为什么能叫袭夺
渐进迁移改变边界位置；袭夺则强调原先属于一套排水系统的上游水流，接入另一套系统。两者可以相连，但不是同一个观测量。冰川事件的证据尤其有助于看清这个区别。4
Shugar 等对2016年事件使用了几种互补证据。Slims River 的水位记录在5月26日开始的四天内明显下降，原有日变化随后减弱；现场测量和无人机摄影测量显示冰川末端存在一条连接两侧湖泊的冰壁河道；卫星时间序列则把新通路与原有河谷、湖面变化放在同一地形中。摄影测量以地面控制点定位，作者对局部高程误差作了保守估计，约为1米。4
这样就能串起一条具体过程：冰川末端退缩、变薄，融水能够进入通往较低湖面的通路；水流和冰壁局部坍塌扩大通道；原供水湖降低约17米，原出水口被留在更高的位置，新路线得到保持。降低一条出口的高程门槛，与沿整条山脊缓慢侵蚀，是不同的触发方式。
单纯干旱可以让一条河水少，却不能单独解释新连接通道和水被转往另一水系；单纯水位传感器出错，也不能解释现场地形与影像中的一致变化。这些是同一事件的不同测量方式，不是多次独立袭夺实验。
时间和分母也需要分清。2016年9月3日测得 Slims River 流量约11立方米/秒；论文另引约130立方米/秒的改道前流量，却来自1970年的历史事件，不能直接拿两者计算“2016年下降了百分之多少”。影像里的窄水面同样不是流量计。
本次已经发生的是冰川融水排水路线重组。论文讨论 Kluane 湖未来可能改变排水状态、进一步转入另一水系，属于当时的预测，不能写成2016年已全部实现。作者将这次改道判断为难以逆转，也依赖冰川不重新推进封堵等条件。本篇不对2017年发表之后的区域全部变化另作更新判断。
一张地形图能告诉我们多少
很多古袭夺没有发生前后的仪器记录。研究者于是寻找地形中的痕迹：突兀的河道转向、与现有小河不相称的大谷地、废弃河谷，或分水岭两侧不协调的河道剖面。这些能定位值得调查的地方，但“看起来像被截走”仍不是完整历史。8
一个常用指标是 $\chi$。它把沿河距离按上游面积重新加权：
$x_b$ 是选定出口，$A_0$ 是固定参考面积。在均匀抬升、侵蚀效率相同、河道达到相应稳态等条件下，前面的下切关系会使高程与 $\chi$ 近似成直线。若相邻河道还共享相同出口高程，而且近山脊坡面段的差异足够小，比较两侧 $\chi$ 可以指出哪一种面积调整有助于两侧匹配。1
但 $\chi$ 不是直接测得的分水岭速度。Forte 和 Whipple 在2018年的受控模型中专门改变抬升分布或让倾斜硬岩层穿过侵蚀地表：分水岭可以趋于稳定而 $\chi$ 差异增大，也可以在岩层切换过程中出现 $\chi$ 对当前移动方向给错提示的阶段。2
原因并不神秘。河道高程不只取决于网络几何，还取决于抬升、岩性、降雨和出口高程。把这些差异省掉，再把所有剖面差异归给“正在争夺面积”，就把一个条件性的指标用成了事实。换一个不合适的出口基准，甚至可能改变指标对比的图像。
因此，应比较靠近分水岭的坡度、局部起伏、同一参考汇水面积处的河床高程，并核对岩性、抬升和出口条件。这些局部指标也受阈值坡面、尺度选择和地形误差影响，不能仅凭两根误差棒是否重叠就宣布移动或静止。2
| 想回答的问题 | 有帮助的证据 | 仍须防止的跳步 |
|---|---|---|
| 今天是否有不对称或调整趋势 | 地形指标、近分水岭侵蚀差异、岩性与基准面对照 | 不对称不自动等于正在迁移 |
| 曾经是否改接水系 | 古河谷连接、沉积物物源、河流阶地与年代互相吻合 | 一个急转弯不能唯一确定袭夺 |
| 何时发生、移动多快 | 有时间约束的地形或沉积记录，必要时与模型比较 | 现在的指标幅度不能直接换成历史速度 |
| 为什么发生 | 能改变候选机制的对照或独立约束 | 多种模型能拟合终态时，不能只挑最顺的一条故事 |
物源证据本身也要防止再搬运与混合：旧沉积物被后来河流挖出来，不等于它的矿物来自当时直接相连的上游。年代、搬运路线与沉积环境需要一起解释。这里给的是历史重建所需的证据结构，不宣称已经重建了某条古河的完整袭夺史。
哪些环节较稳，哪些还没有接上
较稳的是：排水分界可以变化；河道下切与坡面响应能够共同推动调整；通道一旦连通较低出口，也能造成快速重组。受控地形实验、数值过程模型和有现场记录的自然事件分别支持这些环节，它们承担的证据角色不同。
本文的整合判断是：流域边界既是侵蚀过程的结果，也是下一轮侵蚀的条件。 得到多少水会影响地形变化，地形变化又决定下一轮水从哪里来。这个关系可以产生自我加强，也可能被岩性、泥沙覆盖、坡面输运速率和新的基准面条件限制。
尚未接上的，是把这些环节放到一座真实山地后同时约束其强弱与时序。尤其需要区分：现在的地形是在响应当前气候和构造，还是仍保留较早条件的调整过程；新增供水带来的下切，是否被新增供沙抵消；一次袭夺后，哪些地形痕迹会很快被后续侵蚀抹掉。这些问题决定模型能否从“演示可能过程”走到“解释这个地方为什么这样变化”。
所以，理解分水岭移动的关键不只是找出哪条河更大、更陡，而是看清地形、出水通路与物料输运怎样依次改变彼此，再用能约束方向和时间的证据把过程接起来。