GRACE研究海洋质量变化的不确定度分析_刘书豪

《GRACE研究海洋质量变化的不确定度分析_刘书豪》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

第51卷第3期华中科技大学学报（自然科学版）Vol.51No.32023年3月J.HuazhongUniv.ofSci.&Tech.(NaturalScienceEdition)Mar.2023DOI：10.13245/j.hust.230311GRACE研究海洋质量变化的不确定度分析刘书豪杨帆吴毅唐璐(华中科技大学物理学院，湖北武汉430074)摘要以全球海洋质量变化为案例，使用重力恢复与气候实验卫星(GRACE)及其后继卫星(GRACE-FO)重力场和相关数据，通过不同的后处理方法和模型修正的组合，计算了全球海洋质量变化．所使用的组合共计1296个，通过对这些组合的统计研究，计算了GRACE后处理的不确定度．结果表明：GRACE后处理不确定度在全球海洋上的数值结果约为0.64mm/a，并且存在一定的年周期性，这种量级的系统的误差值得引起GRACE科学应用的关注．根据这些结果评估了GRACE在研究全球海洋质量变化的可靠性，并给出相应的后处理方法建议．关键词海洋质量变化；GRACE；重力场；后处理；不确定度中图分类号P223；P229文献标志码A文章编号1671-4512(2023)03-0141-08UncertaintyanalysisofGRACEonglobalmeanoceanmasschangeLIUShuhaoYANGFanWUYiTANGLu(SchoolofPhysics，HuazhongUniversityofScienceandTechnology，Wuhan430074，China)AbstractTakingtheglobaloceanqualitychangeasanexample，gravityrecoveryandclimateexperiment(GRACE)andGRACEFollow-On(GRACE-FO)gravityfieldsandrelateddatawereutilizedtocalculatetheglobalmeanoceanmass(GMOM)changethroughthecombinationofdifferentpost-processingmethodsandmodelcorrection．Upto1296combinationsweretakeninthiswork．Throughthestatisticalstudyofthesecombinations，theuncertaintyofGRACEpost-processingiscalculated，andthenumericalresultsofGRACEpost-processinguncertaintyintheglobaloceanwasgivenforthefirsttime，whichisabout0.64mm/aandhasacertainannualperiodicity．SuchsystematicerrorofthismagnitudedeservestheattentioninGRACEapplication．Basedontheseresults，thereliabilityofGRACEinGMOMchangesresearchwasevaluated，andsomesuggestioninpost-processingisalsogiveninthisstudy．Keywordsoceanmassvariation；GRACE；gravityfiled；post-processing；uncertainty海平面的变化是衡量全球气候变化的一个重要为海洋领域的研究提供了一种全新的手段．2002指标，与人类活动息息相关．对海平面变化的研究年3月，重力恢复与气候实验(gravityrecoveryand也是地球系统变化研究中非常重要的一个课题，相climateexperiment，GRACE)重力卫星正式发射并关的研究在持续进行中，关于海平面变化的具体数投入使用，采集了届时直至2017年的全球时变重[3]值还存在较大争议，因此继续加强海平面变化的监力场数据；作为GRACE任务的后继，GRACE-[1-2]测乃至预测是值得重视的．FO(GFO)于2018年5月正式发射并投入使用．到目海洋质量变化是引起海平面变化的关键因素之前为止，GRACE和GFO重力卫星已经提供了长达一．传统的验潮站观测、卫星高程监测手段等并不约20a的全球时变重力场月解产品，这为全球海洋能直接得到海洋质量的信息，而在21世纪初出现质量变化的研究提供了关键的数据支撑．然而，在的高时空分辨率的卫星重力测量弥补了这一空白，将GRACE时变重力场产品应用于具体问题之前，收稿日期2021-12-21．作者简介刘书豪(1997-)，男，博士研究生；杨帆(通信作者)，副研究员，yfan_cge@hust.edu.cn．基金项目国家自然科学基金资助项目(41804016，41931074，42061134007)．

1•142•华中科技大学学报（自然科学版）第51卷要经过一系列数据处理过程，这一过程称作后了某些后处理步骤对结果的影响，如文献[13]分析处理．了不同的滤波方法对估算全球平均海洋质量变化结比如，GRACE(-FO)Level-2产品以球谐系数果的影响，文献[14]则讨论不同信号泄漏恢复手段的形式给出了地球重力场的分布，目前最新版本对各流域结果的影响．这些研究局限于分析某一特Release06(RL06)的空间分辨率约为300km，但定步骤所用方法对结果的影响大小，而没有对整个GRACE(-FO)解算的重力场不包含一阶项，因此须后处理流程做全面的分析．用其他手段的观测值作为补充[4]；另外，GRACE本研究将以全球海洋质量变化作为案例，使用(-FO)对C20和C30的解算精度较低，因此须要使用不同的后处理手段，得到相应的全球海洋质量变化卫星雷达测距(SLR)数据将C20和C30项替换独立于结果，分析了不同后处理手段对结果的影响，进而[5-6]给出GRACE的后处理不确定度的量化结果．这种GRACE(-FO)的精度更高的观测值．量化结果将明确回答GRACE是否在海洋应用上具此外，由于给出的产品，尤其在高阶系数中，[7]有可靠性，为海洋学研究提供重要的参考．包含相关性噪声，其在空间分布上表现为南北条带误差，因此在使用GRACELevel-2重力数据解释1数据与方法地表质量变化之前，须要对其做一些相应的滤波处理，同时由于滤波所造成的信号泄漏也需要一定的1.1GRACE卫星重力数据反演全球海洋质量变化手段去修复．1.1.1　重力场数据最后，在使用这些数据解释海洋质量变化之GRACELevel-2产品数据主要由美国德克萨斯前，还须要做一些相应的模型改正，如GIA模型改大学空间研究中心(CenterforSpaceResearch，[8][9][10]正、地震模型改正、海洋潮汐模型和CSR)、喷气推进实验室(JetPropulsionLaboratory，[11]GMAM模型等．JPL)和德国地学中心(Helmholtz-CentrePotsdam-显然，上述不同的后处理方法对结果造成的影GermanResearchCentreforGeosciences，GFZ)三家响并不一致；因此，GRACE(-FO)重力场的不确定机构提供．除了上述三家机构提供的GRACE时变度除了由观测手段、测量误差带来的信号不确定度重力场数据外，也使用了国际时变重力场联合服务之外，还有上述后处理方法引起的不确定度，二者(InternationalCombinationServiceforTime-Variable对于所要研究问题最终结果的影响都是不可忽略GravityFields，COST-G)产品以作比较．在应用的．对于前者而言，由于缺乏其他观测数据的验GRACE重力场数据于研究区域之前，还须要对证，对GRACE以及GFO的噪声水平的估计是比较GRACE数据做相应的后处理，以及必要的模型修[12]困难的；然而，对于后者来说，由后处理引起的正，图1给出了大致的后处理流程，其中具体的步结果差异，则是可以量化的．已有的一些研究分析骤将在后几个小节介绍．图1GRACELevel-2重力场数据简化处理流程在重力场低阶项数据的替换步骤中，使用官方提出的DDK滤波方法也可以很好地减弱这样的提供的TN-13文档的一阶项以及TN-14的C20项和误差．C30项作为GRACE低阶项的替换．去相关滤波并不能很好地处理高纬度地区的[17]1.1.2　滤波方法“条带”，还须进一步引入高斯滤波器、各向异性[18][19]Swenson等发现了GRACE产品在高阶系数的高斯滤波器、扇形滤波器等低通滤波器来进一奇次项之间和偶次项之间存在相关误差，其在空间步抑制噪声．域则表现为“南北条带误差”，并提出了相应的去使用上述不同种类和参数的去相关滤波和低通[7]相关滤波方法，该方法将特定次不同阶的系数作滤波组合来抑制高阶系数的相关误差并比较其[15-16]多项式拟合，并从原始系数中扣除．Kusche等差异．

2第3期刘书豪，等：GRACE研究海洋质量变化的不确定度分析•143•1.1.3　信号泄漏的处理地震可以在短时间内引起大范围的质量迁移，不同滤波手段的引入虽然很好地抑制了高阶噪研究海洋质量变化时这一信号须要被扣除．现在一声，但不可避免地对高阶信号部分也有所削弱，在般认为，GRACE可以观测到8级以上的地震信息．空间上则表现为信号的偏差(bias)和泄漏(leakage)，[9]Tang等在利用GRACE数据分析全球海洋质量变前者指研究区域内的信号往外流失，后者则是指外化之前，利用最小二乘的方法对GRACE观测信号部区域的信号流至内部．不同区域信号的偏差和泄拟合出某一点的地震信息并扣除．参照文献[9]的漏影响程度不一．对于偏差的修复，一般可以根据方法，分析了GRACE期间强于里氏8.5级的地震，不同手段预估出一个尺度因子来修正信号；对于泄扣除了相应的拟合信号，具体为漏而言，文献[17]提出了一种简单的处理手法(这里■c0(t

3•144•华中科技大学学报（自然科学版）第51卷[11]值约为-0.3mm/a．水高的贡献，其中GSFC给出的结果在2012年之前海平面高度的变化可以细化为海洋质量的变化没有C20的值，这一段时间内使用CSR机构的数据和比容海平面的变化，卫星测高系统可以得到海平代替．面的高度变化，而“BOA_Argo”网格资料所提供的海水温盐数据则可以得到其中比容海平面变化的部分．因此，将卫星测高得到的海平面变化扣除比容海平面，可以得到海洋质量变化，这一结果理论上应该与GRACE所观测到的海洋质量变化保持一致．2方法验证2.1　联合卫星测高和海水温盐数据得到的全球海洋质量变化图2给出了海平面高度变化和比容海平面变化和联合卫星测高数据和海水温盐数据得到的海洋质量变化，表3列出了其数值结果．由于CMEMS提图3　使用COST-G数据并做不同的低阶项替换后，前三阶供的卫星测高月解不包含2020年5月份以后的数对全求海洋质量变化的贡献据，当计算测高和Argo数据时，给出的是2005—表4C20和C30对海洋质量变化等效水高趋势的影响mm/a2019年间的结果．数据来源C项贡献C项贡献/10-22030COST-G0.176-2.65CSR0.211-2.77GFZ0.116-2.60JPL0.141-2.49GSFC0.120-2.622.3　滤波方法和信号泄漏恢复策略的影响图4和表5给出了不同滤波和信号泄漏恢复方法在海洋区域上质量变化的结果．表5中滤波1~滤波4分别为300km各向同性高斯滤波、500km图22005—2019年间，使用测高卫星和海洋温盐数据估计的全球海平面和海洋质量变化趋势表3联合卫星测高和Argo数据得到的海平面变化海平面和海水质量线性趋势/周年振幅/周年相位/(°)变化(mm/a)mm卫星测高海平面3.855.89144.50Argo比容海平面0.894.126.35测高-Argo2.969.37161.642.2　重力场低阶项的影响图3给出了2005—2020年期间，COST-G给出的重力场在不同低阶项替换的情况下前三阶对全球海洋质量变化的贡献．表4分别给出了不同机构以图4　选取不同滤波和信号泄漏恢复策略得到的全球海洋及GSFCSLR估计的C20和C30项对于海洋质量等效质量变化趋势

4第3期刘书豪，等：GRACE研究海洋质量变化的不确定度分析•145•各向同性高斯滤波、(300，800)km各向异性高斯滤波和(300，300)km扇形滤波，泄漏恢复a~c分别是正向建模法、使用扣除300km缓冲带的海洋区域和使用文献[17]提出的信号泄漏恢复方法并乘以相应的尺度因子．图4则给出了其中部分组合策略的时间序列．图52005—2020年间，使用GAA模型零阶项计算的GMAM信号表5滤波与信号泄漏恢复策略对海洋质量等效水高GSM海洋信号扣除这一部分信号．变化的影响mm/a计算海洋质量变化时，陆海交界处的地震信号策略泄漏恢复a泄漏恢复b泄漏恢复c也是不可忽视的系统误差，这一部分信号可以由滤波11.8912.181.782GRACE观测信号拟合而得．以2012年苏门答腊地滤波21.832.021.79滤波31.842.151.793震为例，图6给出了地震范围内一点(2.5°N，93.5°滤波41.862.171.78E)的GRACE观测信号以及相应的拟合信号，实际注：上标1表示图4中的组合策略1；上标2表示图4中的应用中，将表2所示的地震影响范围内的GRACE组合策略2；上标3表示图4中的组合策略3．信号扣除拟合而得地震信号．表6给出了2005—2.4　模型修正的影响2020年间地震信号对海洋质量变化的趋势影响．对于全球海洋而言，从表1可以看出：GIA的影响是尤为重要的一点．然而，各GIA模型的误差来自于GIA模型本身，因此所有案例均使用文献[24]给出的GIA模型．GMAM信号如图5所示，其代表了解算GSM模型时将0阶项置零时引入的误差，这一部分误差图62005—2020年间，苏门答腊区域一点(2.5°N，93.5°E)的可以从GAA产品导出．实际分析时，将GRACE-重力场信号及拟合出的2012年地震信号表62005—2020年期间地震对海洋质量变化的影响mm/a地震全球海洋扣除300km缓冲带的海洋区域Nias20050.0090.006Maule20100.0090.005Tohoku-Oki20110.0220.025Sumatra20120.0290.018鄂霍茨克海20130.0130.009地震总效应0.0950.0673不确定度分析取如表7所示的方法和参数，其中，模型修正则是“有”或者“无”．因此，最终形成的组合数量为第2节所述的每一步策略对结果都有不可忽略4×3×4×3×3×3=1296种，以上多组后处理的影响．由于后处理流程的方法前后相互影响，每流程将会相互比较得到结果．最终的结果由图7~一步骤对最终结果的影响并不是线性的．除了取决图10和表8给出，其中为了与测高-Argo结果比较，于自身的参数外，还取决于其他步骤策略以及对应表8里的内容均是2005—2019年间的相应结果．参数的选择．对GRACE数据处理的每一步骤，选图7为这1296个案例给出的全球海洋质量在表7使用的后处理方法和参数数据选择低阶项替换高斯滤波信号泄漏恢复模型修正去相关滤波COST-G仅替换1阶项300km高斯滤波正向建模GIA：ICE-6G_D(固定)不使用CSR替换1阶项、C20500km高斯滤波次滤波法GAD模型P3M5GFZ替换1阶项、C20和C30(300，800)km各向异性陆海缓冲区法GMAM模型滑动窗口拟合JPL—(300，300)km扇形—地震模型固定窗口拟合

5•146•华中科技大学学报（自然科学版）第51卷研究时间段内的结果，其中：黑色实线是测高-Ar‐go得到的海洋质量变化结果；红、绿实线是所有组合得到全球海洋质量变化趋势的上、下包络线；中间虚线则是触及到包络线的组合，其余组合的结果则不在图中显示；触及到包络线的部分结果的时间序列以虚线给出，其他部分的结果则不在图里给出．同时，为了规避一些少数极端的结果，图8给图9GRACE后处理不确定度结果出了类似于图7，但根据最终计算的趋势大小，扣除了前后各5%的样本的结果，其中GRACE结果仅包含上下包络线．从图7、图10和表8可以看出：GRACE在海洋上的结果与测高-Argo的结果总体基本一致，但是在细节上存在一些不符合的现象，如相对于符合较好的周年振幅而言，GRACE与Argo-测高的结果在线性趋势上并十分相符．一方面从数据本身所代表的物理信息而言，Argo数据本身不包含2000m水深以下的信息，同时，Ar‐go浮标也并没有对全球范围内的海洋区域都覆盖图10　不同方法计算全球海洋质量变化趋势频次图[25]到；并且通常在海洋区域，扣除了已有的模型后(90%样本)和进行一些改正之后，才认为得到的信号是由海洋趋势，尤其在GFO期间，后处理不确定度表现得质量变化引起的，但是这些模型和相应的改正本身比GRACE期间更为明显，这可能与2016年以后存在较大的不确定度．[26]GRACE和GRACE-FO期间的单加速度计模式有关．由表8可见：通过计算90%样本的最大、最小趋势的差，可以得到后处理得到全球海洋质量增长趋势不确定度约为0.64mm/a．另外，经统计计算，所计算的样本中，线性趋势的标准差约为0.24mm/a．表82005—2019年间，GRACE后处理最大、最小趋势(90%样本)图7GRACE后处理不确定度海平面和海水质量线性趋势周年振幅/周年相位/(°)变化(mm/a)mm测高-Argo2.969.37161.64GRACE最大趋势2.3010.78175.82GRACE最小趋势1.667.91173.04GRACE平均趋势1.949.54174.53GRACE中位趋势1.939.77174.62不确定度0.641.482.78图8　剔除最大、最小的5%线性趋势后得到海洋质量除了线性趋势之外，图9还显示出后处理不确变化趋势定度具有年周期性．从拟合结果可见：其振幅为图9给出了图8中上包络线与下包络线的差值1.48mm，并且在每年的3月份左右，GRACE的后及其拟合的线性趋势、年周期结果，其中：黑色实处理不确定度较小，而9月份左右则达到最大峰线是图8中上下包络线的差值；红色虚线则是其拟值．这表示在每年3月份左右GRACE在海洋上的合的线性趋势和年周期信号，后面称这个差值为结果相对更为可信．“后处理不确定度”．从图9的结果可以看出：在研除了上述对时间序列的分析之外，图10还给究的这段时间内，后处理不确定度随时间有增大的出了上述1296种的结果所得到的线性趋势在扣除

6第3期刘书豪，等：GRACE研究海洋质量变化的不确定度分析•147•最大、最小的5%样本后的频次图．通常而言：期GRACE后处理不确定度的研究对评估GRACE反望这样的统计结果呈现高斯分布是因为这样能够将演地表水质量结果有着重要的意义，首次在全球海后处理不确定度解释为随机误差．但是，从图10洋区域做了一定的分析，后续将进行更为广泛细致可以看出：对于不同组合得到的线性趋势呈双峰分的研究．布，这说明这里面有一些主导的关键物理机制．经参考文献统计分析发现：这样的双峰差异主要来自于后处理的信号泄漏恢复步骤中，是否使用陆海缓冲带的方[1]赵宗慈，罗勇，黄建斌．全球变暖和海平面上升[J]．法，因此图10中的结果进一步以使用缓冲带和非气候变化研究进展，2019，15(6)：700-703．[2]张通，俞永强，效存德，等．IPCCAR6解读：全球和缓冲带的信号泄漏恢复方法加以区分．从表1可以区域海平面变化的监测和预估[J]．气候变化研究进看出：使用陆海缓冲带后得到的GIA趋势比整个海展，2022，18(1)：12-18．洋区域计算的GIA区域更大．因此可以认为，这样[3]TAPLEYBD，WATKINSMM，FLECHTNERF，et的双峰结果很大一方面来自于GIA模型在不同定义al．ContributionsofGRACEtounderstandingclimate的海洋区域的结果差异．实际上，即使考虑了GIAchange[J]．NatureClimateChange，2019，9(5)：模型在全球海洋区域和300km缓冲带的海洋区域358-369．的差异，这样的双峰现象虽然得以改善，然而仍然[4]SWENSONS，CHAMBERSD，WAHRJ．Estimating是存在的，从这一角度来看，现有的GRACE后处geocentervariationsfromacombinationofGRACEand理流程存在系统误差，如何改善或提出更有效可行oceanmodeloutput：ESTIMATINGGEOCENTERVARIATIONS[J]．JournalofGeophysicalResearch：的后处理方法，是值得关注的一个话题．SolidEarth，2008，113(B8)：B08414．[5]CHENGM，TAPLEYBD．Variationsintheearth'sob‐4结论latenessduringthepast28years：variationsintheearth'soblateness[J]．JournalofGeophysicalResearch：Solid以全球海洋质量变化为案例，使用GRACE和Earth，2004，109(B9)：B09402．GRACE-FO时变重力场数据，通过不同的后处理[6]LOOMISBD，RACHLINKE，WIESEDN，etal．方法组合得到了GRACE后处理引起的海洋质量增ReplacingGRACE/GRACE-FOC30withsatellitelaser长趋势和年周期振幅的最大、最小值和不确定度．ranging：impactsonantarcticicesheetmasschange[J]．GRACE在海洋上的线性后处理不确定度约为0.64GeophysicalResearchLetters，2020，47(3)：mm/a，标准差约为0.24mm/a，并且其周期性显GL085488．示，GRACE在每年3月的可信度相对优于9月份的[7]SWENSONS，WAHRJ．Post-processingremovalof可信度．总体来看：GRACE在海洋应用上的后处correlatederrorsinGRACEdata[J]．GeophysicalRe‐searchLetters，2006，33(8)：GL025285．理不确定度在一个较稳定的范围，因此认为[8]CARONL，IVINSER，LAROURE，etal．GIAmod‐GRACE的重力场数据是可靠的、有意义的．elstatisticsforgracehydrology，cryosphere，andocean从结果也可以看出：GRACE在海洋上的结果science[J]．GeophysicalResearchLetters，2018，45与测高-Argo结果在2015年之前较为相符，后续时(5)：2203-2212．间则符合得并不理想，这一方面由于Argo数据和[9]TANGL，LIJ，CHENJL，etal．SeismicimpactofGRACE数据各自在测量原理上不能完美地表示整largeearthquakesonestimatingglobalmeanoceanmass个海洋的信号或将海洋信号与其他信号(如固体地changefromGRACE[J]．RemoteSensing，2020，12球信号等)分离开，另一方面这与GIA模型的不确(6)：935-950．定性也有很大的关系．[10]UEBBINGB，KUSCHEJ，RIETBROEKR，etal．Processingchoicesaffectoceanmassestimatesfrom通过以上研究可以看出：GRACE后处理不确GRACE[J]．JournalofGeophysicalResearch：Oceans，定度是不可忽略的，并且呈现一定的规律，说明这2019，124(2)：1029-1044．些处理方法存在系统性误差，因此对后处理研究还[11]CHENJL，TAPLEYB，SeoK，etal．Improvedquan‐须要进一步深入，才能得到更加准确的GRACE反tificationofglobalmeanoceanmasschangeusing演结果．除了海洋的研究之外，GRACE在其他小GRACEsatellitegravimetrymeasurements[J]．Geophys‐尺度流域等的研究可能存在更大的后处理误差，这icalResearchLetters，2019，46(23)：13984-13991．将会对一系列的科学应用带来更大的挑战．[12]CHENJL，TAPLEYB，TAMISIEAME，etal．Error

7•148•华中科技大学学报（自然科学版）第51卷assessmentofGRACEandGRACEFollow-onmassageerrorinGRACE-observedlong-termicemasschange[J]．JournalofGeophysicalResearch：Solidchange：acasestudyinwestantarctica[J]．JournalofEarth，2021，126(9)：JB022124．Geodesy，2015，89(9)：925-940．[13]郗慧，张子占，陆洋，等．利用GRACE监测全球海[21]LANDERERFW，SWENSONSC．Accuracyof水质量变化时滤波处理的影响分析[J]．大地测量与地scaledGRACEterrestrialwaterstorageestimates：accu‐球动力学，2016，36(5)：380-385．racyofGRACE-TWS[J]．WaterResourcesResearch，[14]VISHWAKARMAB，DEVARAJUB，SNEEUWN．2012，48(4)：WR011453Whatisthespatialresolutionofgracesatelliteproducts[22]GERUOA，WAHRJ，ZHONGS．Computationsofforhydrology?[J]．RemoteSensing，2018，10(6)：theviscoelasticresponseofa3-DcompressibleEarthto852-869．surfaceloading：anapplicationtoglacialisostaticadjust‐[15]KUSCHEJ．Approximatedecorrelationandnon-isotro‐mentinAntarcticaandCanada[J]，GeophysicalJournalpicsmoothingoftime-variableGRACE-typegravityfieldInternational，2013，192(2)：557-572．models[J]．JournalofGeodesy，2007，81(11)：[23]PELTIERWR，ARGUSDF，DRUMMONDR．733-749．Spacegeodesyconstrainsiceageterminaldeglaciation：[16]KUSCHEJ，SCHMIDTR，PETROVICS，etal．theglobalICE-6G_C(VM5a)model：globalglacialiso‐DecorrelatedGRACEtime-variablegravitysolutionsbystaticadjustment[J]．JournalofGeophysicalResearch：GFZ，andtheirvalidationusingahydrologicalmodel[J]．SolidEarth，2015，120(1)：450-487．JournalofGeodesy，2009，83(10)：903-913．[24]PELTIERWR，ARGUSDF，DRUMMONDR．Com‐[17]WAHRJ，MOLENAARM，BRYANF．Timevariabili‐menton“anassessmentoftheICE-6G_C(VM5a)gla‐tyoftheEarth'sgravityfield：Hydrologicalandoceaniccialisostaticadjustmentmodel”byPurcelletal．：TheeffectsandtheirpossibledetectionusingGRACE[J]．ICE-6G_C(VM5a)GIAmodel[J]．JournalofGeophysi‐JournalofGeophysicalResearch：SolidEarth，1998，calResearch：SolidEarth，2018，123(2)：2019-2028．103(B12)：30205-30229．[25]LIH，XUFH，ZHOUW，etal．Developmentofa[18]HANSC，SHUMCK，JEKELIC，KUOCY，etal．globalgriddedArgodatasetwithBarnessuccessivecor‐Non-isotropicfilteringofGRACEtemporalgravityforrections：anewglobalgriddedargodataset[J]．Journalgeophysicalsignalenhancement[J]．GeophysicalJournalofGeophysicalResearch：Oceans，2017，122(2)：International，2005，163(1)：18-25．866-889．[19]ZHANGZZ，CHAOBF，LUY，etal．Aneffective[26]BEHZADPOURS，MAYER-GÜRRT，KRAUSSS．filteringforGRACEtime-variablegravity：fanfilter[J]．GRACEfollow-onaccelerometerdatarecovery[J]．Jour‐GeophysicalResearchLetters，2009，36(17)：L17311．nalofGeophysicalResearch：SolidEarth，2021，125[20]CHENJL，WILSONCR，LIJ，etal．Reducingleak‐(5)：JB021297．