全球冰川均衡调整(GIA)数据下载链接汇总

在处理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.jpghttps://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