GRACE重力数据重建全球地下水：物理约束网络给出每年358立方公里的亏损量
核心概要
GRACE与GRACE-FO卫星只能测到整柱总水量，地下水必须靠扣除土壤水、雪和冠层储水才能分离，而这一步长期依赖的模型会低估最大的趋势。研究者用一套物理约束的分解LSTM，把每个格点的地下水信号拆成趋势、季节循环和残差，并把趋势限制在大气净供水之内、把流域质量守恒写进训练。重建在留出月份上的均方根误差为6.97厘米、Nash-Sutcliffe效率为0.69，在全部测试尺度上对监测井的趋势误差都低于一个独立水文模型和一个数据同化产品。
深度剖析
全球地下水储量变化现在有了一份连续、受物理约束的格点记录，覆盖2002年4月至2022年12月，分辨率为半度。 此前基于GRACE的地下水估计靠从总水量中减去陆面模型给出的土壤水、雪和冠层储水，而全球水文模型相对卫星观测会明显低估最大的下降与上升趋势，这份误差会被继承下来。 在留出的观测月份上，重建的均方根误差为6.97厘米，Nash-Sutcliffe效率为0.69；对独立监测井的趋势误差在单元、含水层、流域和大于20万平方公里流域四个尺度上都低于WaterGAP 2.2d和一个GRACE数据同化产品。
模型把每个格点的地下水异常写成趋势、季节循环与残差三项之和，使物理约束只作用在长期项上。 趋势项只由静态属性驱动，经双曲正切压缩后乘以以年为中心的时间轴，其斜率被限制在该格点大气净供水的1.5倍加1厘米每年之内，该约束的中位数约为8.33厘米每年。 训练时另加一项软约束，要求233个流域的预测平均趋势与观测流域趋势一致；在45轮训练的对照中，权重0.5把流域闭合误差从0.5285降到0.4172厘米每年，同时留出误差下降2.4%。
重建的亏损量被换算成海平面与极移两个行星尺度后果。 此前的地下水亏损估计分散在约126.5到550立方公里每年之间，彼此口径不同。 在显著下降的21632个格点上，总亏损为357.64±11.6立方公里每年，净储量变化为−36.25立方公里每年，相当于0.100±0.049毫米每年的海平面上升；相应的质量重分布给出极移6.92厘米每年、朝东经68.6度方向。
启示与展望
结果针对2002年4月至2022年12月的全球陆地半度格点，共52427个格点参与训练与评估；验证只使用GROW监测井汇编，且井数据只在学习下游进入。作者提出，随着GRACE Follow-On延续记录，可用同一物理约束框架桥接两代任务之间的空档。
给水度以假设值进入不确定度预算，取0.15或0.40，而非在卫星足迹尺度上实测；它在67.6%的格点上是主导误差项。在单格点尺度上，把推算的相对误差扫到约0.51时显著格点比例会骤降，这是显著性检验本身的性质；聚合到约6度的相关区后，显著比例在保守与独立两种假设下分别为61.4%和76.7%。极移速率对分析窗口敏感，缩短到2002至2019年时升至8.145厘米每年，方向只移动0.65度。