利用GRACE估算全球海质量变化

在研究海平面变化时，通常可以使用GRACE反演得到全球海洋质量变化。而由于GRACE反演得到的质量迁移包含众多的要素，因此诸如由于GIA产生的地表垂直位移会造成全球海洋的变化，因此需要从GRACE数据中扣除GIA的影响。同时需要考虑回加GAD模型。本专栏使用的经过反演的全球海洋质量变化时间序列的基本处理方法如下：

上面我们看到了这一处理方式：# Ocean Kernel: Total ocean with 300km coastal buffer zone。这是为了避免陆地信号的泄露对海洋质量估计的影响。一般我们使用的海洋掩膜如下所示：

然而为了减小陆地信号泄露的影响，我们需要对海岸线数据进行300 km和500 km的缓冲区处理，下面是缓冲之后的效果：

运行结果

接着在GRACE GSM数据中回家GAD的球谐系数，反演得到全球的质量变化之后扣除GIA的影响，根据不同的GIA模型，可以得到不同的结果。我们此处分析了CSR、GFZ、JPL和ITSG18四个数据类型得到的全球海洋质量变化的结果：

带有误差区间阴影的时间序列。

最后，附上全球海洋300 km和500 km缓冲区的掩膜数据的下载链接：

欢迎交流学习！

参考文献：