摘要：本文利用2003—2011年的GRACE RL05数据提取了苏门答腊地震（Mw9.3）引起的震后重力变化,发现断层两侧震后重力变化速率存在明显差异,断层下盘总体变化率为0.55μGal/yr,断层上盘为0.16μGal/yr.基于子断层叠加的编程思想,本文将Tanaka的黏弹球体位错理论配套计算程序（简称黏弹位错程序）加以改造,克服了其近场计算精度不足（甚至错误）的缺陷,可用来研究大地震引起的近场震后位移与重力变化.本文利用改造后的黏弹位错程序计算了2004年苏门答腊地震（Mw9.3）产生的同震重力变化,计算结果在空间分布和量级上均与利用弹性位错程序计算获得的结果一致,验证了我们对黏弹位错程序进行改造的正确性.最后,结合GRACE卫星观测数据,本文利用Tanaka的黏弹位错理论研究了苏门答腊地区的地幔黏性因子.结果表明,该地区地幔黏滞性具有显著的横向差异,当发震断层上下两盘的地幔黏滞性系数分别取8×1018 Pa·s和1×1018 Pa·s时,模拟的震后重力变化在总体空间分布和变化趋势上与GRACE卫星观测结果更接近.

关键词：球体地球位错理论 grace重力卫星 重力变化 黏滞性 

单位：中国地震局地震预测重点实验室(地震预测研究所) 北京100036 中国科学院计算地球动力学重点实验室 北京100049 大地测量与地球动力学国家重点实验室 武汉430077 东京大学地震研究所 东京113—0032