IF50.5｜格陵兰冰盖夏季融水储存的时空特征

GNSS监测基岩垂直位移揭示冰盖周边夏季融水缓冲储存，7月达峰后逐渐释放。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

格陵兰冰盖（GrIS）是当前全球海平面上升的最大单一贡献者，主要因北极放大效应导致加速融化。融水径流是2000年后冰盖质量损失的最大来源（约55%）。融水在流向海洋的途中，可暂时储存于冰面湖、粒雪层、冰内空腔、冰床界面或地下水中，这种缓冲储存（BWS）影响冰盖演化，例如通过改变冰下水压影响冰流速度。然而，由于冰盖内部水文过程复杂且观测稀疏，对BWS的时空变化、储存时长及其对径流的影响仍缺乏定量约束。

科学问题

卡在哪里： 现有观测手段各有局限：冰雷达可探测冰下水的存在但无法定量其体积和压力；卫星重力（GRACE）能监测总水量变化但空间分辨率粗，难以区域化评估；区域气候模型（如RACMO2.3p2）对融水径流的模拟存在系统偏差，尤其在异常温暖年份可能低估径流或高估粒雪持水能力。因此，缺乏一种能高时空分辨率地约束冰盖内部融水储存和释放过程的方法。

本文要回答： 利用格陵兰GNSS网络（GNET）连续监测的基岩弹性垂直位移，量化格陵兰冰盖BWS的时空演化，评估区域气候模型径流估算的系统误差，并确定不同区域的水储存时间。

技术路线

作者主要用了：GNSS垂直位移、GRACE卫星重力、区域气候模型RACMO2.3p2、解析模型、蒙特卡洛模拟。

作者利用22个GNET基岩GNSS站的连续垂直位移观测，扣除大气、非潮汐海洋、陆地水、SMB等负荷信号后，得到残余垂直位移，其季节性下沉和回升反映了BWS的积累与释放。为量化BWS，构建了一个解析模型，假设BWS随时间指数衰减，并引入年际径流校正因子，通过非线性最小二乘拟合GNSS时间序列，同时估计水储存时间Tst和径流缩放因子。该设计将GNSS对近场质量变化的高敏感性（约200 km半径内贡献80%）与区域气候模型的径流输出相结合，从而分离出BWS信号并评估模型误差。

图 1　Fig. 1。图1展示融水季节前后冰盖及基岩垂直位置的示意。a为融水季节前，水量最小，基岩实际位置与背景模型一致；b为融水季节早期，液态水快速积累但向海洋排放少，基岩实际位置仅轻微抬升（棕色箭头），而背景模型（不考虑水积累）显示快速抬升（白色箭头），两者分离。该图直观说明了BWS对基岩位移的影响机制。

核心亮点

亮点 1｜BWS季节循环及空间差异

所有GNET站显示一致的季节模式：2月至4月缓慢下沉，5月加速，7月达峰，之后逐渐回升至次年2月。该信号解释为BWS在融水季节的积累和随后的释放。空间上，东南部站点在7月后快速回升，表明BWS快速流失；而西部和东北部站点回升较慢，持续至9-10月。这种差异反映了不同区域水文体制的不同，可能与地形、积累率和消融带宽度有关。

图 2　Fig. 2。图2为22个GNET站残余垂直位移的年平均周期。所有站均显示2-4月缓慢下沉，5月加速，7月达峰，之后回升。红色阴影为1σ不确定性。图中还叠加了冰流速图，显示站点多位于高流速区域附近。该图揭示了BWS信号的普遍性和空间一致性。

亮点 2｜BWS定量及年际变化

通过解析模型拟合，得到各站BWS引起的垂直位移时间序列。校准后的BWS在2012年极端融化年达到最大，如KAGA站位移达14 mm。2010年也是高BWS年份。校准后的BWS年际变化比未校准的更大，表明区域气候模型可能低估了温暖年份的径流。

图 3　Fig. 3。图3比较了GNSS垂直位移（红色）与GRACE卫星重力反演的质量异常（黑色）的月平均序列。两者均显示5-7/8月质量增加，随后下降至次年2月。GRACE尺度反转以便比较。两者高度一致，验证了GNSS信号的真实性。

亮点 3｜径流缩放因子的温度依赖性

研究发现，温暖年份（如2010、2012）的径流缩放因子显著大于普通或寒冷年份，尤其在北部和东北部站点。缩放因子与夏季温度异常呈正相关（R=0.42）。对于最暖年份，径流需上调约20%。这表明区域气候模型在温暖年份可能低估融化或高估粒雪持水能力，导致实际径流高于模拟值。

图 4　Fig. 4。图4展示了基于校准BWS估计的解析模型得到的垂直位移时间序列。各站时间序列显示明显的年际变化，2012年峰值突出，尤其在南部和西南部站点。该图直观呈现了BWS的年际动态。

亮点 4｜水储存时间的区域差异

估算的水储存时间Tst在大多数站点为3-13周，平均约8周。东北部和西部平均约9周，而东南部平均仅31天（约4.5周）。最长的储存时间出现在西北端的MARG站，达129天。东南部较短的储存时间可能与该地区高积累率、陡峭地形和狭窄消融带导致的快速排水有关。

图 5　Fig. 5。图5显示温暖年（a）、普通年（b）和寒冷年（c）的径流缩放因子。温暖年缩放因子普遍大于1，寒冷年小于1，误差棒基于100次蒙特卡洛模拟的四分位距。该图支持缩放因子与温度相关的结论。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： GNSS垂直位移、GRACE卫星重力、区域气候模型RACMO2.3p2、解析模型、蒙特卡洛模拟
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事冰盖质量平衡、水文过程或海平面上升研究的科研人员，以及关注GNSS在地球物理应用中的方法开发者。

带走一句

GNSS基岩位移可高分辨率地约束冰盖内部水文过程，为改进区域气候模型的径流估算提供独立校准手段，尤其对温暖年份的模型偏差敏感。