卫星测高数据、GRACE数据时间序列中移除GIA信号

我们知道GIA，也称PGR，是地球物理中极其重要的一个概念，其会影响全球的地球物理观测信号。以下文章介绍了基本的概念。下图展示了全球GIA的信号分布（GIA & Trends | Get Data – GRACE Tellus）

GIA研究进展 - 哔哩哔哩 (bilibili.com)

下面我们介绍如何在卫星测高数据和GRACE数据中移除GIA。

（1）卫星测高数据

A_tt = aviso_msla.tt;
for ii=1:length(A_tt)
    cs_tmp         = aviso_msla.rg(:,:,ii).*aviso_mask;
    % remove PGR from aviso grid
    aviso_msla_gia.rg(:,:,ii) = cs_tmp - (A_tt(ii)-mean(A_tt)).*aviso_gia.*aviso_mask;
    disp(ii)
end

我们可以看到改正前后，对于海平面变化的影响相对较小。

（2）GRACE数据

A_tt = LLZ_CSR.tt;
for ii=1:length(A_tt)
    cs_tmp_grace         = LLZ_CSR.rg(:,:,ii);
    % remove PGR from aviso grid
    grace_gia.rg(:,:,ii) = cs_tmp_grace - (A_tt(ii)-mean(A_tt)).*O1.rg;
    disp(ii)
end

我们可以看到，经过GIA改正的结果与未改正还是存在较大差异。

♥欢迎点赞收藏♥