在处理GRACE数据时，由于GRACE监测的信号包含地表的质量迁移信号和固体地球物理信号，因此研究地表的质量迁移时需要扣除固体地球的信号，目前最主要的时冰川均衡调整（冰后回弹）的信号。具体的关于冰川均衡调整的信息可参考我的另外一个专栏文章：
GIA研究进展 - 哔哩哔哩 (bilibili.com)
下面我将目前最常使用的一些GIA模型的下载链接进行汇总。
（1）ICE系列
具体的数据可以参考Peltier的个人网站：Home Page: W.R. Peltier
在网站上可下载到历史的冰雪覆盖厚度数据。CE-6G_C results:
Files are available for the following times (Ka)
26.0, 25.0, 24.0, 23.0, 22.0, 21.0,
20.5, 20.0, 19.5, 19.0, 18.5, 18.0, 17.5, 17.0, 16.5, 16.0,
15.5, 15.0, 14.5, 14.0, 13.5, 13.0, 12.5, 12.0, 11.5, 11.0,
10.5, 10.0, _9.5, _9.0, _8.5, _8.0, _7.5, _7.0, _6.5, _6.0,
_5.5, _5.0, _4.5, _4.0, _3.5, _3.0, _2.5, _2.0, _1.5, _1.0, _0.5, _0.0
其中的每一个文件都是包含1°×1°空间格网的冰雪厚度数据。A file containing the ice-thickness for all times from 122Ka to the present is IceThick.nc.gz.下面我们读取这个数据并进行展示：
Global Attributes:
Model = 'ICE-6G_D'
Acknowledgement = 'Please cite both PAPER1 and PAPER2'
PAPER1 = 'Argus, D.F., Peltier, W.R., Drummond, R. and Moore, A.W.(2014). The Antarctica component of postglacial rebound model ICE-6G_C (VM5a) based upon GPS positioning, exposure age dating of ice thicknesses, and relative sea level histories. Geophys. J. Int., 198(1), 537-563, doi:093/gji/ggu140'
PAPER2 = 'Peltier, W.R., Argus, D.F. and Drummond, R. (2015). Space geodesy constrains ice-age terminal deglaciation: the ICE-6G_C (VM5a) model. J. Geophys. Res.-Solid Earth, 120, 450-487, doi:10.10002/2014JB011176'
Dimensions:
Time = 122
Lat = 180
Lon = 360
Variables:
stgit
Size: 360x180x122
Dimensions: Lon,Lat,Time
Datatype: single
Attributes:
name = 'IceT'
units = 'm'
long_name = 'Ice thickness'
Time
Size: 122x1
Dimensions: Time
Datatype: single
Attributes:
name = 'Time'
units = 'Ka'
long_name = 'Kyrs BP'
Lat
Size: 180x1
Dimensions: Lat
Datatype: single
Attributes:
axis = 'Y'
units = 'degrees_north'
standard_name = 'latitude'
Lon
Size: 360x1
Dimensions: Lon
Datatype: single
Attributes:
standard_name = 'longitude'
units = 'degrees_east'
modulo = 360
topology = 'circular'
axis = 'X'file = 'IceT.I6F_C.131QB_VM5a_1deg.nc';
ncdisp(file)
lon = ncread(file,'Lon');
lat = ncread(file,'Lat');
Time = ncread(file,'Time');
stgit = ncread(file,'stgit');
[x,y] = meshgrid(lon,lat);
for i = 1:12
subplot(4,3,i)
O.lon = x;
O.lat = y;
O.rg = stgit(:,:,i+10)';
wzq_plot(O)
end每个年代的冰雪的覆盖空间分布
下面是可以下载使用的数据：
ICE-6G_D results:
Rate of radial displacement (UP) on a 0.2x0.2 grid: drad.12mgrid.nc
Horizontal Velocity (East & North) on a 0.2x0.2 grid: Hvel.12mgrid.nc
图片引用自：GIA & Trends | Get Data – GRACE Tellus
（2）Caron模型
具体的模型信息可参考这篇文章：https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2017GL076644
这里给出了数据下载的网站：Glacial Isostatic Adjustment Simulation | Solid Earth | VESL | JPL | NASA
进入如下的界面
点击下载按钮，即可下载不同阶数的Stokes Coefficients Rate
（3）Paulson07 GIA 模型
https://commons.wikimedia.org/wiki/File:PGR_Paulson07_big.jpg
https://commons.wikimedia.org/wiki/File:PGR_Paulson07_big.jpg
The ICE-5G (VM2_L90), Geruo13 (a/b/c) and Paulson07 models are very similar in predicting the spatial distribution of velocity because they all use the same ice model: ICE-5G
（4）Geruo13 GIA模型
关于Geruo13 和Paulson07 GIA模型，有需要可以直接给论文作者发邮件获取。此处我提供Paulson07_mass和Geruo13_uplift
链接：https://pan.baidu.com/s/1VwAp9mkwxPzywywqHkkEuQ?pwd=ucas
提取码：ucas欢迎点赞收藏！
参考文献
Richard Peltier, W., Argus, D. F., & Drummond, R. (2018). Comment on “An assessment of the ICE‐6G_C (VM5a) glacial isostatic adjustment model” by Purcell et al. Journal of Geophysical Research: Solid Earth, 123(2), 2019-2028.
Caron, L., E.R. Ivins, E. Larour, S. Adhikari, J. Nilsson and G. Blewitt (2018), GIA model statistics for GRACE hydrology, cryosphere and ocean science, Geophys. Res. Lett., 45, doi: 10.1002/2017GL076644.Paulson A, Zhong S J, Wahr J .2007. Inference of mantle viscosity from GRACE and relative sea level data[J]. Geophysical Journal International, 171 (2) : 497–508. DOI:10.1111/gji.2007.171.issue-2
Geruo A, Wahr J, Zhong S J .2013. Computations of the viscoelastic response of a 3-D compressible Earth to surface loading:An application to Glacial Isostatic Adjustment in Antarctica and Canada[J]. Geophysical Journal International, 192 (2) : 557–572. DOI:10.1093/gji/ggs030
