GRACE重力卫星探测南极冰盖质量平衡及其不确定性

Ice sheet mass balance in Antarctica measured by GRACE and its uncertainty

2002年GRACE重力卫星的成功发射为南极冰盖质量平衡的研究提供了重力探测的新纪元.本文利用美国德克萨斯大学CSR公布的2003年1月到2013年12月期间的RL05版本GRACE月重力场数据,采用最优平均核函数法和组合滤波法两种GRACE后处理方法反演了南极冰盖质量的时空变化.结果表明:在2003—2013年期间南极冰盖物质平衡呈明显的负增长状态,质量变化趋势为-163±50Gt/a(GW13)、-129±41Gt/a(IJ05)、-81±27Gt/a(W12a),加速度为-8±10Gt/a2,质量消融的主要区域分布在西南极阿蒙森海岸和南极半岛的北部.另外本文还重点探讨了可能影响到估算结果的各项误差及不确定性,分析结果显示影响南极冰盖质量平衡估算结果的最大误差源为GIA改正.通过假设检验和信息准则对时间序列分析中拟合参数的合理选取进行了探讨和分析,在联合周年项、半年项和S2、K2、K1潮汐混频项进行拟合分析时发现K1项对拟合结果的加速度影响比其他周期项稍大,尽管考虑该项的合理性因当前GRACE数据时间序列长度有限而无法确切证实,但K1项的影响值得后续关注.对比两种GRACE后处理方法的结果发现:当采用的数据时间跨度一致,误差改正方法相同,两种相异的后处理方法,其估算结果也具有较好的一致性.

As a critical component of the cryosphere,the Antarctic Ice Sheet（AIS）has strongconnection with the sea level change and global climate change.Accurate quantification of the current spatial and temporal mass changes of AIS is very important to improve our understanding and prediction of its response and contribution to global change.The Gravity Recovery and Climate Experiment（GRACE）mission has provided new and useful observations to detect AIS mass balance since its launch in March 2002.There are significant differences among the GRACE estimates of the total mass change.The big difference is due in part to considerable uncertainty in the accuracy of glacial isostatic adjustment（GIA）signals,and also due to use of different time spans,different versions of GRACE products and different GRACE post-processing methods.Using 124 monthly GRACE gravity field solutions of Release 5（RL05）produced at the Center for Space Research（CSR）of the University of Texas,Austin,spanning the interval from January of 2003 through December of 2013,the mass balance of AIS is estimated by two postprocessing ways：the optimizing averaging kernel method（also named VW）and the combined filter method（the first step is called P5M11 decorrelation filter to remove correlated noise by fitting and subtracting a fifth-order polynomial to even and odd coefficient pairs at spherical harmonic orders eleven and above,the second involves smoothing with a 250 km Gaussian filter）.A detailed error analysis is provided including consideration of leakage-in,leakage-out,and errors in modeling mass variations of the atmosphere,ocean and GIA.In addition,a statistical model selection criterion is employed in computation of trends from mass variation time series,and the impact of K1 tidal alias is analyzed.The results reveal that during 2003—2013,the total mass of the ice sheet decreased significantly at change rates of-163±50,-129±41and-81±27Gt/a for three GIA models：GW13,IJ05,W12 a.There was a distinct region with mass loss in the Amundsen Sea Embayment of West Antarctic ice sheet and the Northern Antarctic Peninsula,while an increasing mass gain was concentrated in the Dronning Maud Land and the Enderby Land of East Antarctic ice sheet.Furthermore,we use hypotheses testing and information criteria evaluation to select the best trend model fitting together with sinusoidal functions of annual（365.0-d）and semi-annual（181.0-d）signals and the S2（161.0-d）,K1（2725.4-d）and K2（1362.7-d）tidal aliases.We found that K1 tidal alias has a potential to falsify the acceleration estimates.Although it is not good enough to confirm the K1 tidal alias based on an eleven-year time-series,the impact of K1 tidal alias deserves further notice.By comparing the quantities of total mass balance computed by the two different processing methods and three different GIA models in the Antarctica,we find that the differences are less than 15Gt/a between two processing methods,but the largest difference is about 80 Gt/a between different GIA models.The analysis of the uncertainty of GRACE′s estimation of AIS mass balance indicates that the largest source of error is the GIA correction.Our results indicate that during January 2003 to December 2013 the contribution of AIS to sea level rise was about＋0.34±0.11mm/a.Significant mass loss increases were limited to the basin that contains Pine Island Glacier along the Amundsen Sea coast of West Antarctica.During the analyzed time period,the total mass acceleration was-8±10Gt/a2,equivalent to＋0.02±0.03mm/a2 sea level rise.Results of analysis point to the conclusion that when using agiven GRACE data set with same error correction,the differences of total mass changes are not highly dependent on which post-processing strategies to be used but on the different GIA models.Therefore,a more accurate GIA model is the key for determining Antarctic ice mass change from GRACE in the present and future.