结合GLDAS模型估计GRACE/GRACE-FO全球的RMS比率（信噪比）

在我们阅读文献的时候，经常出现如下类型公式：

上述公式出自 Chen J L, Wilson C R, Seo K W. Optimized smoothing of Gravity Recovery and Climate Experiment (GRACE) time‐variable gravity observations[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B6). 中的公式（6）在JGR: Solid Earth本篇文章网站显示所有引用这篇文章的学者论文

大多数文献中对于该公式的详细解释和具体操作是缺乏的，这里我想试着根据自己的理解大致实现这一公式。如有问题欢迎各位友友大胆批评指正！

1、数据准备

（1）陆地和海洋的掩膜数据（UTCSR网站：GRACE/GRACE-FO - Gravity Recovery and Climate Experiment (utexas.edu)）

（2）GLDAS_NOAH025_M.A.021数据从2002年4月至2009年12月，缺失月份和GRACE对应上

去除2004.01-2009.12的均值，代码中“gldas02_09.mat”即为该格网数据

（3）GRACE CSR RL06 2002.04-2009.12去平均，不经过任何滤波的质量数据（cs_mss02_09.mat），经gmt_cs2grid即可转为格网EWH

2、代码展示

clearvars -except gldas02_09;
addpath E:\Data\CSR\result\RL06;
load E:\Code\GLDAS\GLDAS2006m025\gldas02_09.mat
%%如何计算误差项，是除去季节项和趋势项之后剩余的部分么？
load cs_mss02_09.mat;cs_mass(:,:)=cs_mss02_09(17,:,:);%以2003年11月份数据为例
grid=gmt_cs2grid(cs_mass,0,0.25);
file_land='E:\Data\CSR\Mascon\data\CSR_GRACE_GRACE-FO_RL06_Mascons_v02_LandMask.nc';
file_ocean='E:\Data\CSR\Mascon\data\CSR_GRACE_GRACE-FO_RL06_Mascons_v02_OceanMask.nc';
% info = ncinfo('E:\Data\CSR\Mascon\data\CSR_GRACE_GRACE-FO_RL06_Mascons_v02_OceanMask.nc');
lon=ncread(file_land,'lon');lat=ncread(file_land,'lat');
L_val=ncread(file_land,'LO_val');L_msk=flipud(L_val');
O_val=ncread(file_ocean,'LO_val');O_msk=flipud(O_val');

grid_weight=zeros(720,1440);
for j=1:720
    for i=1:1440
        grid_weight(j,i)=cos(lat(j)*pi/180);
    end
end
land_weight=grid_weight.*L_msk;ocean_weight=grid_weight.*O_msk;
l_weight_sum=sum(sum(land_weight));o_weight_sum=sum(sum(ocean_weight));

temp1(:,:) = grid(:,:).*100;temp2(:,:) = gldas02_09(:,:,17)*0.1;
RMS_L=sqrt(sum(sum(((temp1-temp2).*L_msk).^2))./sum(sum(L_msk)));
RMS_O=sqrt(sum(sum(((temp1-temp2).*O_msk).^2))./sum(sum(O_msk)));
SNR=RMS_L./RMS_O;
SNR_log=10*log10(SNR);

3、结果与讨论

计算得到GRACE无约束解的RMS比率为1.5755；后期实验高斯300滤波的RMS比率为3.5005、Swenson滤波为2.4583。结果符合我们的认知，经过滤波处理后RMS比率有所提升，说明滤波处理能提高信噪比。但是上述代码仍有缺陷，没有完全与RMS比率公式匹配上，缺乏Err项的计算。对于Err误差项如何计算，我目前还未找到合适的方法，如果有明白其中细节的，欢迎私信交流（q：2629512206）

如果对你有帮助，多多转发多多集赞哦。Thanks！

参考文献：

[1] Chen J L, Wilson C R, Seo K W. Optimized smoothing of Gravity Recovery and Climate Experiment (GRACE) time‐variable gravity observations[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B6)

[2] 魏征强. GRACE时变重力数据球谐域滤波方法与质量分布反演应用研究[D]. 中国矿业大学, 2023. DOI:10.27623/d.cnki.gzkyu.2023.000516