19 - grace数据处理 - 补充 - 地下水储量计算过程分解 - 冰后回弹(GIA)改正
![](https://img-blog.csdnimg.cn/direct/74a6fee40c8e4f50959e646fe2da7bea.jpeg)
0 引言
由水量平衡方程可以将地下水储量的计算过程分解为3个部分,第一部分计算陆地水储量变化
、第二部分计算地表水储量变化
、第三部分计算冰后回弹改正
、第四部分计算地下水储量变化
。本篇简单介绍下第三部分的内容,主要是冰后回弹改正量的计算 。有需要对冰后回弹的内容进一步了解的可以看文章。
1 gia数据处理过程
前文有介绍过GIA的数据形式和数据处理策略,对于球谐系数的GIA数据要进行球谐系数截断、高斯滤波、球谐展开到等效水厚等处理过程,可以得到对应区域的冰后回弹速率;对于网格形式的GIA数据,要先进行球谐分析得到球谐系数再进行高斯滤波、球谐展开等得到区域冰后回弹变化速率。本文使用ICE-6G球谐系数形式的GIA数据计算青藏区域的冰后回弹效应。主要程序如下: