利用GPS垂直位移反演区域陆地水储量变化的TSVD-Tikhonov正则化方法

利用GPS垂直位移反演区域陆地水储量变化(TWSC)属于典型的病态问题,其关键是如何进行稳定求解并提高反演结果的精度和可靠性.本文引入TSVD-Tikhonov组合正则化方法对利用GPS垂直位移反演区域TWSC的病态问题进行求解,并以四川省TWSC反演为例进行分析与验证.首先,通过数值模拟对TSVD、Tikhonov和TSVD-Tikhonov正则化方法采用不同正则化参数选取策略(RMSE最小准则、GCV法和L-curve法)进行反演,结果显示基于TSVD-Tikhonov正则化反演的TWSC比单独使用TSVD或Tikhonov正则化反演结果的精度和可靠性更高,这三种正则化方法反演2005年1月至12月的TWSC差值的平均STD分别为14.97 mm、7.03 mm和5.04 mm.其次,利用中国地壳运动观测网络(CMONOC)的72个GPS测站的垂直位移数据,基于TSVD-Tikhonov正则化反演了四川省2010年12月至2021年2月的TWSC时间序列,结果表明GPS反演的TWSC与GRACE/GFO Mascon模型(JPL、CSR和GSFC)的空间分布特征及季节性变化符合较好,但其TWSC信号的振幅比GRACE/GFO Mascon模型更强.最后,采用广义三角帽方法(GTCH)融合不同类型的降水、蒸散发和径流数据,并根据水量平衡方程计算的dTWSC/dt序列(PER-dS/dt)对GPS反演的dTWSC/dt序列(GPS-dS/dt)和GRACE/GFO Mascon模型融合的dTWSC/dt序列(GRACE/GFO-dS/dt)进行验证,结果表明这三类dTWSC/dt序列的季节性变化符合较好,平滑后GPS-dS/dt和GRACE/GFO-dS/dt序列与PER-dS/dt序列的相关系数分别为0.78和0.87,但GPS相比GRACE/GFO对降水变化的响应更为敏感.本文研究证明了TSVD-Tikhonov组合正则化方法能够提高GPS垂直位移反演区域TWSC的精度和可靠性,同时也表明GPS观测数据对局部水质量负荷变化更为敏感,可作为GRACE/GFO反演区域TWSC的有益补充.