历史与路径依赖14 分钟

融水进入冰川底部,为什么有时使冰川加速,有时加速会减弱?

融水对冰川运动的影响,取决于它怎样改变底部承重,以及排水通路怎样随后调整。本文核对经典通道模型、2025年软床冰川研究与2026年南极钻孔观测,重算南极水压和定位记录,并用两个平行通道说明压力的快速响应与几何的缓慢变化。更多来水可以先提高滑动倾向,通路扩张后响应又减弱;但软床储水、沉积物形变和海洋边界提供了不同路径,不能仅凭一条季节速度曲线识别全部底部机制。

同样是水增加,冰川的反应可以随时间改变

春天开始融化,一些冰川会突然加快;到了夏天,来水仍很多,速度却可能回落。另一些冰川的夏季速度相对平稳,冬季暖事件反而能带来明显加速。

这些现象提醒我们,水的作用不只有“增加一点润滑”。水会改变冰与床面怎样承重,流动又会改变自己的通路:有的通道扩大,有的空隙闭合,有些水暂时储存,有些区域后来才接入排水网络。1

因此,需要同时问两件事:现在的来水怎样改变压力,过去的来水又把排水系统变成了什么样。 两个时刻的融水量一样,若之前形成的通路不同,速度响应仍可能不同。

不过,冰川底部并不是一种统一材料。下面可能是起伏岩石,也可能覆盖能变形的沉积物;出口可能通向空气、湖泊或海洋。理解一种冰川的过程之后,还要检查它能带到另一种冰川多远。

水压改变的,是接触处还需要承担多少载荷

冰很厚,自己的重量会压向底部。水进入底部连通空间后,可以承担其中一部分压力。一个基本变量是有效压力:

$$ N=P_i-P_w. $$

$P_i$为冰的上覆压力,$P_w$为相关底部空间的水压。$N$越低,水承担的压力越接近上覆冰压力;在合适的接触条件下,冰与床面更容易分离,或者沉积物颗粒之间的有效压紧减弱。1

这不意味着“压力差一小点,摩擦就只变一小点”。假如水压从上覆压力的93%变到97%,它只增加了四个百分点,但剩余压力差从7%降到3%,下降了一半以上。

2026年发表的南极Langhovde冰川研究,提供了一个具体例子。本篇按作者绘图中的几何和密度常数,重算2022年1月事件:事件前的钻孔水位中位数约为海拔30.87米,峰值约51.05米;相应有效压力由约0.341兆帕降到0.143兆帕。67

水位海拔需要先减去床面海拔,才得到支持水压的水柱高度,不能直接拿“海拔51米”去除以“漂浮水位海拔”。此外,淡水钻孔的水头与海水浮力几何使用的密度不同,两个“漂浮高度”不能随手互换。

$N$也不是一个可以包办全部摩擦的数字。床面形状、接触面积、颗粒重排、冰的变形和周围应力传递,都会参与决定运动。尤其是一个钻孔测得的压力,未必代表整片床面的平均状态。

通路会被挤小,也会被运动和流水撑大

冰会缓慢变形,因而冰下空隙不是固定管道。水压降低、有效压力升高时,周围冰通常更容易向通路内部蠕变,促使空隙闭合。

开口则有不同来源。

在起伏岩床上,冰滑过凸起,背面会留下空腔。运动不断打开空隙,蠕变又把它关小。若更高水压抑制闭合,空腔就能变大。

另一类通道主要靠流水的能量耗散融化冰壁而扩大,通常称为R通道。较多水流带来更多可用融蚀能量;通道扩大以后,在较低水压下也可能排走较大流量。这种几何调整需要时间。12

所以,突然来水与长期维持来水,不能用同一个固定管径来理解。突然增加时,原几何和原压差可能不足以输送来水,水压或储水量先响应;通路扩大后,同样的输送任务所需压力会改变。

软床上还会有另一套开闭机制:水可能侵蚀沉积物,形成宽而浅的沟道,或在颗粒孔隙和多条互通路径中移动;沉积物自身又会变形。不能仅凭“存在水道”,就认定它全部由冰壁融化维持。

一个模型怎样把两种开口放在同一条式子里

Schoof在2010年的模型中,用一条面积收支式表达空腔与融蚀通道的共同部分:23

$$ \frac{\mathrm dS}{\mathrm dt} =c_1Q\Psi+u_bh-c_2N^3S. $$

$S$为通道截面积,$Q$为流量,$\Psi$为沿流向的水力势梯度。右边三项依次表示流水耗散带来的融蚀开口、冰滑过床面凸起带来的开口,以及蠕变闭合。$u_b$是滑动速度,$h$概括凸起尺度;系数包含潜热、冰的变形性质等。

模型还用一个湍流通量关系连接流量、截面与梯度。在正向流动时,它可写为:

$$ Q=c_3S^{5/4}\Psi^{1/2}. $$

固定水力梯度,并让通道充分调整到稳态,可以算出$N$随$Q$怎样变化。在小流量一侧,滑动开口相对重要,增大稳态流量对应较低的有效压力;在大流量一侧,融蚀开口变重要,较大流量反而可以对应较高有效压力,也就是较低水压。

本篇按原文参数重算了这个关系。其斜率改变处约为每秒0.218立方米,特定数值取决于这组参数;它不是自然界所有冰川共用的来水门槛。

按Schoof2010单通道关系重算,固定水力梯度为512帕/米,其他条件见计算记录。曲线比较充分调整后的稳态,不能用它直接代替突然增加来水时的瞬时变化。
按Schoof2010单通道关系重算,固定水力梯度为512帕/米,其他条件见计算记录。曲线比较充分调整后的稳态,不能用它直接代替突然增加来水时的瞬时变化。

这个图解释的是一种条件性反转。它还没有告诉我们一个真实网络里哪条通道会变大、附近空腔会不会接通,更没有证明任何已经观察到的夏季减速都由同一机制造成。

压力先变,通道随后变:一个实际算过的最小网络

为把时间过程讲完整,本篇用同类通量和开闭关系,建立了两个并排的通道。它们共享上游压力,通向同一个出口;总流量由输入决定,较大通道可以分到更多水。

初始两个截面只相差约2%。总输入先维持每秒0.4立方米,随后增加到每秒1.2立方米。这个模型忽略通道自身的储水,因而水压立即调整到能让输入流出的值;它演示的是压力与几何调整,不是积水体积的复算。

在一个对照中,滑动造成的开口速度固定。在另一个对照中,加入一条明确的近似反馈:有效压力降低时,滑动代理增大,滑动开口也随之增大。代理定义为初始有效压力除以当前有效压力,代表特定恒定驱动假设下的相对响应,并非真实冰面速度预报。

模型确实出现了先增强、后减弱:四天内提高输入时,含反馈版本的相对滑动代理峰值约1.319,随后回到接近原水平;到第300天约为0.992。最终输入仍是原来的三倍。

把增加过程放慢到40天,同一模型里的代理峰值只有约1.027。两者最终输入速率相同,不过在整个计算期间的累计输入量不同;这里没有进行等水量对照。

本库两个平行通道的示意计算。上图为输入,中图为基于压力的滑动代理,下图为四日增加、含反馈版本的两个截面。代理不是现场测量;通道数量、位置和大尺度几何预设,未计算完整冰流或软床侵蚀。
本库两个平行通道的示意计算。上图为输入,中图为基于压力的滑动代理,下图为四日增加、含反馈版本的两个截面。代理不是现场测量;通道数量、位置和大尺度几何预设,未计算完整冰流或软床侵蚀。

图中还出现了通道之间的分化。一条略大的通道分到更多水,融蚀开口增强,后来继续扩大;另一条分到的水相对减少。这种反馈并不是从一开始就一定占优:本篇的稳定性计算显示,低输入时,两个通道的微小差异会衰减;较高输入时,差异能够放大。

这让“排水能力改变”有了一个具体过程,而不是预先指定一句“夏天自动排得更好”。但这个最小网络还省略了节点储水、沉积物侵蚀、通道迁移及完整应力传递。它能证明所列机制在明确条件下可发生,不能直接估计Langhovde的通道大小或年速度。

原2010年论文也有一个相似的边界:用于展示冰流变化的速度是从水压推算的输出,未完整反馈到原排水模型的滑动开口速度。本篇另列的反馈对照只是补上其中一个环节,并不等于建立了完整的冰川动力学模型。

南极新观测:水、压力与运动的联系有多直接

Sugiyama等2026年报告了2021/22夏季在Langhovde冰川的钻孔和定位观测。2022年1月初的温暖与降雨期间,现场看到表面积水进入裂缝;钻孔水压随后明显升高,冰面也加速。接地段还出现向上位移,为底部接触发生改变提供额外线索。6

这些证据比单看气温与卫星速度相关更直接:它把表面来水、底部状态和运动放在了同一事件中。不过,从表面裂缝到床面的完整流水路线并未被逐段追踪;抬升也需要与垂向应变、床面坡度等解释比较。原文用接地与漂浮位置的差别及运动变化来约束这些替代解释。

本篇从公开定位坐标重新计算了局部速度,使用有实际观测覆盖的时间窗口,不填补长缺测。7 下图比较12小时与24小时窗口;它采用本篇明确的局部线性斜率方法,不冒充原作者未随包提供的平滑函数。

2022年1月公开观测的重分析,阴影表示1月2日至6日事件期间。速度由已有定位坐标计算,虚点线为24小时窗口、实线为12小时窗口;曲线中断保留观测覆盖不足。下图按作者公开换算得到钻孔水位,原表缺值没有被一律填满。
2022年1月公开观测的重分析,阴影表示1月2日至6日事件期间。速度由已有定位坐标计算,虚点线为24小时窗口、实线为12小时窗口;曲线中断保留观测覆盖不足。下图按作者公开换算得到钻孔水位,原表缺值没有被一律填满。

以1月1日均速为本篇固定参考,三种时间窗口得到以下事件峰值增幅:

计算窗口 接地段GNSS1 漂浮段GNSS2
6小时 约47% 约37%
12小时 约27% 约19%
24小时 约16% 约13%

加速的日尺度表现仍在,但较短窗口会保留更多短时变化与定位噪声,较长窗口会平滑峰值。这些差别不是三个独立试验,也不能挑最高一行当成更真实的加速幅度。

这也是为什么不能把图上的尖峰当成一个精确到分钟的机械响应时间。原文某些缺测段使用插值重建;其中第一次事件的部分速度还跨越缺测,而压力传感器尚未提供那次事件的同步记录。本篇对水压—运动关系的重点放在第二次有数据的事件。89

这里记录的是冰面运动,包含底部运动与冰内变形及应力传递。漂浮段能够随上游变化加速,并不意味着该位置也有相同的冰—床滑动。类似地,事件后水压下降,也不能仅凭下降曲线就断言已直接观察到R通道扩张。

潮汐可以传来压力,海水不必一起到达测点

Langhovde靠近海洋,出口边界也会变化。作者观察到钻孔水位中有与海潮相联系的周期成分。接地段钻孔所见水的盐度低于仪器检出限,这与压力通过连通水系传递、测点仍主要是淡水的情形相容。6

压力变化的传播和盐水颗粒的输运是两件事。前者不要求海水在同一时间、以同样速度到达钻孔。反过来,一个钻孔没有检出盐分,也不能替整个接地区域所有时间都排除海水进入。

本篇复算了作者补充图使用的潮位比较。需要保留一项处理:公开脚本先将潮汐时间平移2.5小时,之后才进行对应;漂浮区域钻孔BH2203与潮位的相关系数约0.965,进一步采用作者的四分之一小时相对错位后约0.970,均与其图中数值吻合。78

公开脚本口径的潮位比较。BH2201采用1月13–25日,BH2203采用1月25日至2月6日;图中已包含作者对潮汐时间的2.5小时平移。点数是时间网格数,不是彼此独立的重复。
公开脚本口径的潮位比较。BH2201采用1月13–25日,BH2203采用1月25日至2月6日;图中已包含作者对潮汐时间的2.5小时平移。点数是时间网格数,不是彼此独立的重复。

直接使用未平移的时间标签,BH2203相关约0.508。这个对照不否定它与潮汐的联系,说明的是相位与时间处理会影响“相关有多高”。在时钟、对齐和沿途响应还未完全分开时,本篇不把2.5小时或额外15分钟解释为已经测准的海水流动时间,也不用大量连续时间点给出过度精确的独立样本显著性。

软床研究为什么不能全套进“夏天形成大通道”

Hart等2025年的研究比较了四个温暖海洋性环境中的软床冰川,提出从较集中通道到多通道分布式排水的连续变化。4

其中一种图景是:夏季来水分散在许多互通路径与储水空间里,冬季又分次释放;另一种则有较集中的排水、较突出的春季响应和不同的夏季变化。沉积物的透水性、可侵蚀性、湖泊水位及冰川尺度,都可能影响它们。

这里有一个特别容易误读的观测:有些加速事件发生时,埋在沉积物里的探针水压反而下降。它不必与“高界面水压促滑动”矛盾。运动引起的卸载和颗粒剪胀,本身可能改变局部孔压;一个传感器测到的是循环中的一部分,而不是整个床面同一时刻的统一状态。

不过,这项研究的证据并非四套完全同步、同样完备的测量。两个冰川没有对应的倾斜探针数据,部分过程是由速度代理推断;一个地点的年度速度使用其他年份建立的关系重建;挪威地点的卫星比较使用了附近另一条冰川。不同图还存在跨年份拼接。4

因此,这些季节模式值得作为机制假设和分类入口,尚不能把它们当成一套已经在陌生冰川上独立验证的识别规则。卫星12日平均与一个GNSS点的日变化,提供的信息也不同。分类中的5℃或6℃是按地点设定的季节、事件门槛,不是“冰到这个温度才融化”。

研究给出的约15%滑动、70%变形、其余无明显变形,也是按代理识别的时间份额。它不能翻译为“70%的总位移来自软床变形”。

一个算术例子就能看出区别:若互不重叠的滑动阶段占15%时间、速度为10单位,而变形阶段占70%时间、速度为0.2单位,那么前者贡献约91%的位移。这个例子不是冰川估计,只说明时间比例不能替代速度积分。真实冰川的不同运动分量还可能同时发生,识别难度更高。

到一条新的冰川上,需要哪些证据才能继续判断

观察到速度随融水增加而升高,可以提出水文机制,但还不能选择唯一过程。更有区分力的是把测量与具体推论接起来。

想解释的变化 有帮助的证据 不能直接替代的判断
来水是否改变底部承重 水压、冰厚与来水时间,多个测点的接触线索 只由气温推算全床压力
排水系统是否随季节调整 压力、流量、储水、示踪或地球物理信息的共同变化 只看夏季速度降低就宣布通道化
是否发生软床形变 局部应变或倾斜、孔压、速度及空间位置 用某一阶段持续久推断位移贡献大
海洋边界是否参与 潮位、多个钻孔与接地状态,明确时间对齐 把压力相位差等同盐水行程
短时事件是否造成长期净加速 足够长的速度与通量记录,包含事件后的恢复 用一个峰值代表全年冰量损失

本题能回答的是,融水改变运动存在明确的物理环节与现场证据;排水几何的适应也能够让持续较大来水对应较弱的后续加速。具体冰川走哪条路线,则取决于床面材料、连通和储水状态、出口条件以及来水过程。

冰川记住过去来水的一种方式,就是把它写进了当前的排水通路和接触状态。 这份历史可以帮助解释响应为什么变化,但需要用水压、流量与运动的证据辨认,不能只从一条速度曲线把整个底部结构反推出来。

可回查:南极逐小时水位换算、定位速度重分析、窗口敏感性、潮位比较、单通道稳态、两通道稳定性、两通道时序、时间与位移例子、计算汇总及边界。

来源与依据

  1. S1 领域机制综述

    Modelling water flow under glaciers and ice sheets

    Gwenn E. Flowers · 2015

    Proc R Soc A471,20140907;公开发表PDF,PMC4991255 · 相关全文已核

    基础守恒/势能式3.1–3.5、排水要素分类和channels/cavities段落

    区分流量与水头、几何容量与实际储水;多种硬床/软床要素,不把R通道当全部排水系统。

  2. S2 原始通道模型

    Ice-sheet acceleration driven by melt supply variability

    Christian Schoof · 2010

    Nature468,803–806;2010-12-08在线、12-09卷期 · 相关全文已核

    主文式1/2、Fig1稳态/稳定性、网络与变化输入、Methods Summary参数

    本库重算单通道稳态关系及Qc;两个平行通道为缩减示意,不称重现原20000导管网络。单通道固定梯度稳态与暂态压力关系分开。

  3. S3 模型推导与限定

    Schoof2010 — Supplementary Information

    Christian Schoof · 2010

    · 相关全文已核

    单导管动力学、网络式44–49与节点储水、28–29页式58–62速度推算和反馈限制

    原基准忽略导管储水;uslide为被动输出,ub未反馈,不能称完全耦合。大指数与近零N的过大速度原作者也承认不现实。本文只补标量开口反馈并保持N>0。

  4. S4 原始多地点观测与机制解释

    Seasonal velocity patterns provide insights for the soft-bed subglacial hydrology continuum

    Jane K. Hart; Nathaniel R. Baurley; Amy Bonnie; Benjamin A. Robson; Graeme Bragg; Kirk Martinez · 2025

    Communications Earth & Environment6,223;2025-03-21;PMC11928315 XML · 相关全文已核

    Results/Discussion;Table1–4;Methods测量/季节门槛/阶段划分

    四地点非四套同步全观测:Skal跨年温度关系重建,Fjalls/Breid无tilt代理,Briks卫星用邻近Nigards;部分夏季跨年。77/83%等是分类吻合比例,不是Pearson;15/70%是时间非位移。季节5/6℃门槛不等于冰熔点。

  5. S5 旧观测数据入口与说明

    GLACSWEB — Skálafellsjökull public dataset documentation

    Jane K. Hart and collaborators · 2021

    · 相关全文已核

    Probe/flow README、目录更新时间及旧数据范围

    取得2009/10探针与较早输入/输出表;目录支持2021相邻研究,不是完整2025四冰川数据。本文未重跑旧表,不把未取得的最新四地数据当已复现。

  6. S6 原始钻孔与冰面观测

    Acceleration of an Antarctic outlet glacier driven by surface meltwater input to the base

    Shin Sugiyama; Ken Kondo; Masahiro Minowa; Akira Watanabe · 2026

    Nature Communications;2026-05-06;PMC13357834正式XML · 相关全文已核

    Results水压/速度、Discussion融水/抬升与海洋联系、Methods定位及水头测量

    重点第二次2022年1月事件,第一事件没有同步压力且部分GNSS缺测;裂缝全路径未直接追踪,抬升比较替代解释。近海区水压相关不能自动给出全南极年度净冰量损失;盐度未检出限于该观测。

  7. S7 原始论文公开处理数据包

    Field data from Langhovde Glacier and MATLAB code for plotting — V1

    Shin Sugiyama; Ken Kondo; Masahiro Minowa; Akira Watanabe · 2026

    Mendeley Data V1;公开CC BY4.0;2021/22野外观测 · 相关全文已核

    公开目录、数据日期、许可与19个已选数据/代码文件

    2026-05-01发布,原论文参考文献写2021不当数据发表日期。18个dat/m文件按公开SHA验证,root README另得;图间同名文件SHA同一,存一份。不是原始卫星载波相位和所有影像。

  8. S8 原始补充图表

    Sugiyama2026 — Supplementary Information

    Shin Sugiyama and coauthors · 2026

    · 相关全文已核

    FigS4潮汐与相关/相移、FigS6定位平滑与缺测;S1机制示意和S2温度范围

    BH3约.965/.970可重算;相移不自行解释为海水流动时间。缺测插值不是实测数据;未用环境生物图像承重冰流。

  9. S9 公开审稿与作者回复

    Sugiyama2026 — Transparent Peer Review File

    Anonymous referees; Shin Sugiyama and coauthors · 2026

    · 相关全文已核

    缺测插值说明、海拔/水柱混淆、第一事件无水压、检测限及外推

    审稿赞同不当独立数据;回正式稿与源表。未找到2.5小时与平滑函数定义的完整说明,不自行编造原因。

  10. S10 公开数据/绘图代码

    BH2201 public pressure series

    Shin Sugiyama and coauthors · 2026

    Mendeley V1;文件SHA及路径见计算汇总 · 相关全文已核

    逐分钟时间/压力;图代码偏置9.2dbar与1.0197m/dbar

    原表已处理且含NaN,保留缺值;水头与有效压力据绘图几何换算,未核原始仪器校准。

  11. S11 公开数据/绘图代码

    GNSS1 public positioning output

    Shin Sugiyama and coauthors · 2026

    Mendeley V1;文件SHA及路径见计算汇总 · 相关全文已核

    接地段坐标、时间与数据间隔

    WGS84 ECEF转局部EN,再用明确6/12/24h线性窗口;不跨长缺测,未复现作者local_regression。

  12. S12 公开数据/绘图代码

    GNSS2 public positioning output

    Shin Sugiyama and coauthors · 2026

    Mendeley V1;文件SHA及路径见计算汇总 · 相关全文已核

    漂浮段坐标、时间与数据间隔

    漂浮段水平速度不等于此位置冰床滑动;平滑敏感性完整保留,不挑更大峰值当真实响应。

  13. S13 公开数据/绘图代码

    Syowa tide public series

    Shin Sugiyama and coauthors · 2026

    Mendeley V1;文件SHA及路径见计算汇总 · 相关全文已核

    30秒时标/海面高度;代码减全段均值和-2.5h

    同时保留未平移与作者平移后相关;序列自相关显著,不把网格点数当独立重复或简单p值。

  14. S14 公开数据/绘图代码

    BH2203 public pressure series

    Shin Sugiyama and coauthors · 2026

    Mendeley V1;文件SHA及路径见计算汇总 · 相关全文已核

    列13压力、时标、.089MPa偏置与101.97m/MPa

    25Jan–6Feb原代码时间网格比较;海水边界处淡水钻孔高度与海平面不必相同。

  15. S15 公开数据/绘图代码

    Figure3 plotting code

    Shin Sugiyama and coauthors · 2026

    Mendeley V1;文件SHA及路径见计算汇总 · 相关全文已核

    几何、压力换算、GNSS缺测/样条、h=12/24与外部local_regression

    主文写1h Gaussian而外部函数未在包中取得,不能从参数h自行确定完整原滤波核;本库另给明确算法。固定压力偏置不证明之前未有时变校正。

  16. S16 公开数据/绘图代码

    FigureS4 plotting code

    Shin Sugiyama and coauthors · 2026

    Mendeley V1;文件SHA及路径见计算汇总 · 相关全文已核

    潮汐整体-2.5h、BH3另加.25h相对错位、插值范围

    实际重算BH1 r.497944与BH3 .964901/.970365;原始时标BH3 .508188。不把所用相移当独立物理传播测量。

最后更新 2026-09-12

来源