01 - GRACE数据处理理论介绍回顾

GRACE数据处理理论介绍回顾

  好像在每一篇介绍GRACE数据处理的硕博论文中，都能看到类似的理论😀

  大地水准面被认为是一个重力等位面，那么重力场的变化就和大地水准面变化是有关联的，二者可以相互表征。地球表层物质的迁移与重新分布会引起局部重力场的变化，这一变化又作用于GRACE双星系统，使星间距离产生偏差从而得到了测量数据。利用GRACE月重力场球谐系数的变化量来表示大地水准面的变化($\Delta N$)的公式如下:
$$\Delta N(\theta ,\lambda )=a\sum _{l=1}^{\infty }\sum _{m=0}^l{\overline{P}}_{lm}(cos\theta )(\Delta C_{lm}cos(m\lambda )+\Delta S_{lm}sin(m\lambda ))$$

  a为地球的平均半径6371km；$\theta$ 和 $\lambda$ 为地心余维和地心经度；l和m为球谐系数的阶次；$\Delta C_{lm}$和$\Delta S_{lm}$ 是完全规格化球谐系数；${\overline{P}}_{lm}$ 为勒让德函数。

  引起大地水准面(重力场)变化是由一定体积的物质密度变化($\Delta \rho (r,\theta ,\lambda )$)引起的，根据上边的式子可以推出物质密变化与球谐系数的关系式：

$$\left\{\begin{array}{r}\Delta C_{lm}\\ \Delta S_{lm}\end{array}\right\}=\frac{3}{4\pi a{\rho }_{ave}(2l+1)}\int \Delta \rho (r,\theta ,\lambda ){\overline{P}}_{lm}(cos\theta )\times (\frac{r}{a}{)}^{l+2}\left\{\begin{array}{r}cos(m\lambda )\\ sin(m\lambda )\end{array}\right\}sin\theta d\theta d\lambda dr$$
${\rho }_{ave}$ 为地球平均密度 $\approx 5517km/m^3$ ；$\Delta \rho (r,\theta ,\lambda )$ 表示物质密度的变化量；

  假说物质密度变化主要发生在地球表层，那表层的各圈层的质量交互就是引起该密度变化的主要因素，大气、海洋、冰川及水等。地球表层物质交互主要集中在约10-15km厚度内，那么通过积分可以得到这一厚度内的面密度变化：
$$\Delta \sigma (\theta ,\lambda )=\int \Delta \rho (r,\theta ,\lambda )dr$$

  由于变化区域相较于地球半径a太小了，$(r/a{)}^{l+2}\approx 1$ 所以上边公式就可以进一步化简为：
$${\left\{\begin{array}{r}\Delta C_{lm}\\ \Delta S_{lm}\end{array}\right\}}_{surf}=\frac{3}{4\pi a{\rho }_{ave}(2l+1)}\int \Delta \sigma (\theta ,\lambda ){\overline{P}}_{lm}(cos\theta )\left\{\begin{array}{r}cos(m\lambda )\\ sin(m\lambda )\end{array}\right\}sin\theta d\theta d\lambda$$

  引起大地水准面变化主要有两个因素，其一为上文所谈到的地球表面薄层质量变化，其二为地表负荷变化，地球表层物体的质量长期对地球产生负荷作用，导致固体地球产生形变，故计算大地水准面变化也需要考虑地表负荷变化的影响，即：
$${\left\{\begin{array}{r}\Delta C_{lm}\\ \Delta S_{lm}\end{array}\right\}}_{solid}=\frac{3k_l}{4\pi a{\rho }_{ave}(2l+1)}\int \Delta \sigma (\theta ,\lambda ){\overline{P}}_{lm}(cos\theta )\left\{\begin{array}{r}cos(m\lambda )\\ sin(m\lambda )\end{array}\right\}sin\theta d\theta d\lambda$$
$k_l$ 为负荷勒夫数

因此总的大地水准面变化为：

$$\left\{\begin{array}{r}\Delta {\hat{C}}_{lm}\\ \Delta {\hat{S}}_{lm}\end{array}\right\}={\left\{\begin{array}{r}\Delta C_{lm}\\ \Delta S_{lm}\end{array}\right\}}_{solid}+{\left\{\begin{array}{r}\Delta C_{lm}\\ \Delta S_{lm}\end{array}\right\}}_{surf}$$
$\Delta {\hat{C}}_{lm}$ 和 $\Delta {\hat{S}}_{lm}$ 地表质量异常球谐展开的无量纲的规格化球谐系数

接着对面密度($\Delta \sigma (\theta ,\lambda )$)做球谐展开就能得到GRACE重力场模型计算地球表层密度变化的关键公式了：
$$\Delta \sigma (\theta ,\lambda )=\frac{a{\rho }_{ave}}{3}\sum _{l=0}^{\infty }\sum _{m=0}^l{\overline{P}}_{lm}(cos\theta )\frac{2l+1}{1+k_l}(\Delta C_{lm}cos(m\lambda )+\Delta S_{lm}sin(m\lambda ))$$

再通过$\Delta H=\Delta \sigma /{\rho }_w$ 就将面密度变化转化为了等效水高变化，${\rho }_w$ 表示水的密度，那么公式就可以转化到最终形态：
$$\Delta H(\theta ,\lambda )=\frac{a{\rho }_{ave}}{3{\rho }_w}\sum _{l=0}^{\infty }\sum _{m=0}^l{\overline{P}}_{lm}(cos\theta )\frac{2l+1}{1+k_l}(\Delta C_{lm}cos(m\lambda )+\Delta S_{lm}sin(m\lambda ))$$