GRACE卫星观测到的与汶川Ms8.0地震有关的重力变化

利用GRACE卫星重力资料,计算了中国大陆及周边的卫星重力时变场和地表密度变化分布,获取了具有代表性的点位区域的每月重力变化时间序列.同时获得了WUSH、LHAS、KUNM、LUZH站相对于区域参考框架的GPS位移时间序列.卫星重力观测结果显示喜马拉雅弧形带的重力在2004年苏门答腊M_w9.3地震后快速下降, 2006~2008年尤为明显,西域地块西北边界带上震后重力下降也较为显著;而沿青藏高原北至东边界2007年出现明显的重力上升沿构造边界的弧形分布,且2008年南北地震带中南段重力上升变化显著.这些苏门答腊地震后的重力变化趋势到汶川地震发生后才开始改变.GPS位移结果显示四个台站均记录到苏门答腊大地震的同震信号,震后WUSH、LHAS、KUNM站水平位移向量出现明显的运动趋势改变,且一直持续到2008年汶川M_s8.0地震的发生.GRACE卫星揭示的青藏高原及周边地表质量的变化为解释汶川地震的动力机制提供了新的观测途径和资料.本文结合区域构造运动的特点和GPS位移,对GRACE观测的时变重力场特征及汶川地震的动力机制进行了初步解释和讨论.     

GRACE卫星重力场同震变化的经验正交函数分解:以日本M_W9.0地震为例

2011年3月11日日本东海发生MW9.0地震,造成日本岛整体东移、下沉并伴随巨大的质量重新分布。对于海底地震的震中区域,空间观测的GRACE卫星重力数据很好地弥补了GPS、InSAR等形变资料的缺失。利用GRACE卫星月重力场数据计算了地面0.5°×0.5°网格点上的重力变化时间序列,采用最小二乘拟合、经验正交函数(EOF)2种方法,提取了同震重力变化,结果显示震中两侧区域的重力变化呈两极分布,其中弧后区域重力下降,最大降幅约6μgal,海沟区域重力增加,最大增幅约3μgal。EOF方法避免了最小二乘拟合方法所需的地震发生时刻等先验信息,但卫星重力信号是由多种地球物理过程引起的重力变化的叠加,EOF结果的可靠性及其反映的真实物理来源往往随着事件的规模、观测时间的长短等而改变。文中第2,3,4主成分主要反映了非构造因素的影响,通过第1主成分空间一致性提取的同震重力变化与位错理论模型计算结果较为接近,因此较真实地反映了地震引起的重力变化特征。     

利用GRACE观测数据研究苏门答腊区域的黏滞性结构

本文尝试利用卫星重力观测资料和震后黏弹性松弛理论研究苏门答腊地区的区域流变结构,为更好地认识区域地球动力学环境提供依据．利用GRACE卫星重力资料,计算了2004年苏门答腊Mw9．3地震的同震及震后的重力变化．计算中使用平滑半径为500km的高斯滤波器．结果显示苏门答腊Mw9．3地震破裂东侧以陆地为主的上盘同震下沉量很大,造成约9x10＾-8mS＾－2的重力下降阶变,而西侧处于海水下的下盘重力同震上升较小约2x10＾－8mS＾－2,但其震后上升较快．流变结构对岩石的变形有很大的影响,是地球动力学数值模拟取得可靠科学结果的基础．本文尝试了利用卫星重力变化时间序列来反演苏门答腊地区的黏滞性结构;即基于GRACE时变重力场,利用自重力、黏弹性、平面分层模型模拟了该地震的同震和震后变形,并将获得的空间固定点的重力变化与GRACE重力场及点位时间序列进行比较,估计该地区的黏滞性系数在1．0x10＾18PaS的量级,且断层两侧的流变参数存在差异．最后结合苏门答腊区域的构造特点讨论了黏滞性系数的影响因素．     

重力卫星GRACE检测出2010年智利Mw8.8地震的同震重力变化

本研究通过对重力卫星GRACE观测数据的处理,采用去相关加300 km半径的高斯平滑,成功地提取了2010智利M_w8.8地震所产生的重力变化信号,最大变化幅度达到7 μGal,并且与位错理论计算结果具有较好的一致性.这是继GRACE检测出2004苏门答腊M9.3地震重力变化后的又一个卫星观测地震的例证,说明GRACE具有检测出M<9.0量级地震的能力,为利用GRACE研究地震以及其更广泛的应用提供了可靠的依据.     

GRACE卫星重力在地震研究中的应用进展

地球重力场恢复与气候实验卫星(GRACE)在运行期间提供了大量的地球时变重力场观测数据,在大地测量、地球环境变化等领域有非常广泛的应用.在固体地球科学研究中,GRACE重力场数据被广泛应用于天然地震研究,由于地震过程中存在大范围的质量迁移,大型地震引起的重力变化可以被GRACE卫星探测到.同时GRACE记录的地震同震及震后长期的重力场变化对反演地震震源参数也起到了帮助作用.本文从GRACE卫星重力场在地震研究中的应用出发,在回顾了GRACE卫星重力在地震应用的主要进展的基础上,总结了以地震研究为应用目标的数据处理方法与流程,为地震同震及震后卫星重力数据处理提供了技术思路.本文结合前人在2004年M_W9.3苏门答腊地震、2010年M_W8.8智利地震、2011年M_W9.0日本东北地震的研究成果,针对地震重力变化的同震观测、震后观测、间接观测等领域,总结了GRACE卫星重力的主要应用进展,提出了其中涉及的GRACE数据处理与地震综合解释的主要研究问题.在总结技术研究进展的基础上,本文以2004年M_W9.3苏门答腊地震为研究对象,对GRACE卫星重力数据序列进行处理,得到该地震的同震重力变化特性,并以此为基础进行了地震同震重力变化的特征分析.在回顾和总结的基础上,本文结合GRACE-Follow on计划的优势,提出未来GRACE卫星应用于地震研究的发展展望.     

利用GRACE重力和重力梯度变化估计2012年苏门答腊地震断层参数

2012年4月11日苏门答腊北部附近海域发生M_W=8.6地震,国际上采用不同数据获得了该地震的断层滑动模型,但断层走滑性质存在较大差异.卫星重力GRACE（Gravity Recovery and Climate Experiment）观测覆盖震中区域,可以提供很好的断层参数估计约束.本文采用GRACE月重力场模型数据提取了此次地震同震重力和重力梯度变化,表明了北向分量的优越性.并与断层模型理论计算结果进行了比较,分析了T_xx（北-北向重力梯度）分量对断层参数的敏感性,以及基于粒子群算法以及Okubo位错模型反演了该地震断层参数.结果表明GRACE观测到的同震重力梯度变化在空间形态分布上与断层模型模拟结果符合较好,但振幅差异较大.重力和重力梯度的北向分量可以很好地压制条带误差,其中g_N（北向重力变化）和T_xx的误差最小,其次是T_xz（北-上向重力梯度）,误差最大和对条带误差抑制效果最差的是T_xy（北-东向重力梯度）.T_xx对断层的深度不敏感,对其余断层参数较敏感.基于GRACE反演得到的断层的走向角为113.63°,倾角为89.99°,滑移角为175.26°,平均滑移量为28.18 m,相应的矩震为8.71级,且此次地震的走滑性质为右旋走滑.     

利用震后黏弹性位错理论研究苏门答腊地震(M_w9.3)的震后重力变化

本文利用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卫星观测结果更接近.     

日本MW9.0地震震区及其周缘2002—2015年卫星重力变化时空特征

本文在考虑全球陆面数据同化系统陆地水储量变化影响后,利用2002年4月—2015年4月GRACE卫星RL05月重力场模型获取了2011年日本M_W9.0地震震中及其周边区域的重力场信息;然后给出了日本及其周边2003—2015年的年度累积重力变化和差分重力变化,并且利用经验正态函数方法深入分析了该地震过程中同震重力变化对区域重力场的贡献. 结果显示:日本M_W9.0地震前研究区域出现了幅值高达2×10-8 m/s2的异常重力变化,其同震效应的时间和空间特征均出现在第一模态,且同震重力变化和震后区域重力场变化特征显著,这充分表明该地震对区域重力场的影响显著.      

云南地区及邻域卫星重力变化及其与区域大震的关系

利用GRACE卫星重力场数据,计算并分析云南地区及邻域2003—2014年的卫星重力场变化特征,获得3个特征点上的月重力变化时间序列;结合该地区2003—2014年发生的M_S5以上地震分布,探讨区域卫星重力变化与区域大震的关系。结果表明:云南地区及邻域卫星重力场变化对区域M_S5以上地震有一定的反映,地震多发生在区域重力变化的高梯度带上,地震发生前后,区域出现与地震相对应的特征性卫星重力变化,可能是对壳幔物质调整的反映;区域大震前后,不同震区特征点的重力变化特征存在差异,说明分断层、分区域的重力场变化研究对地震预测有重大意义。     

芦山M_S7.0地震震前GRACE卫星时变重力场特征研究

基于GRACE卫星时变重力数据,利用小波多尺度分解方法,分析研究了芦山MS7.0地震前重力场的时空动态演化特征.结果表明,雅安地区震前卫星时变重力场具有明显的重力变化信息.卫星时变重力场在时间域表现为,累积重力场和差分重力场存在较强的规律性动态变化,以2003—2004年均值为基值的2009—2012年的变化特征与2005—2007年的变化特征具有较好的相似性,以2009年均值为基值的2009—2012年即芦山地震前4年,芦山震区的时变重力场呈现上升→减速上升→加速上升以及上升→轻微下降→快速上升的变化特征;卫星重力场在空间域表现为,震中位置由靠近重力场正变化的中心→接近重力场正变化边缘→接近重力正变化四象限分布中心的变化特征.因此,芦山地震的重力场动态变化信息和波场分离结果符合逆冲推覆兼走滑型地震引起的区域变形特征.     

Seismologic applications of GRACE time-variable gravity measurements

The Gravity Recovery and Climate Experiment (GRACE) has been measuring temporal and spatial variations of mass redistribution within the Earth system since 2002. As large earthquakes cause significant mass changes on and under the Earth’s surface, GRACE provides a new means from space to observe mass redistribution due to earthquake deformations. GRACE serves as a good complement to other earthquake measurements because of its extensive spatial coverage and being free from terrestrial restriction. During its over 10 years mission, GRACE has successfully detected seismic gravitational changes of several giant earthquakes, which include the 2004 Sumatra–Andaman earthquake, 2010 Maule (Chile) earthquake, and 2011 Tohoku-Oki (Japan) earthquake. In this review, we describe by examples how to process GRACE time-variable gravity data to retrieve seismic signals, and summarize the results of recent studies that apply GRACE observations to detect co- and post-seismic signals and constrain fault slip models and viscous lithospheric structures. We also discuss major problems and give an outlook in this field of GRACE application.     

Correction of Hydrological and Oceanic Effects from GRACE Data by Combination of the Steric Sea Level,Altimetry Data and GLDAS Model

In this study, a scheme to estimate oceanic and hydrological effects in the GRACE (Gravity Recovery and Climate Experiment) data is presented. The aim is to reveal tectonic signals for the case of the Sumatra earthquake on 26 December 2004. The variations of hydrological and oceanic effects are estimated with the aid of data set of GRACE, altimetry, World Ocean Atlas, and the GLDAS model for a period of January 2003 to December 2006. The time series of computed gravity changes over Sumatra region show some correlations to the deformation resulting from the earthquake occurred in December 2004. The maximum and minimum impacts of hydrological and oceanic effects on gravity changes are about 3 μGal in radial direction and–5 μGal in northward direction. The maximum and minimum amounts of gravitational gradient changes after the correction are 0.2 and–0.25 mE, which indicates the significant influences of hydrological and oceanic sources on the desired signal.     

GRACE detection of the medium- to far-field coseismic gravity changes caused by the 2004 MW9.3 Sumatra-Andaman earthquake

Large earthquakes cause observable changes in the Earth's gravity field, which have been detected by the Gravity Recovery and Climate Experiment (GRACE). Since most previous studies focus on the detection of near-field gravity effects, this study provides the results from the medium- to far-field gravity changes caused by the 2004 Sumatra-Andaman earthquake that are recorded within GRACE monthly solutions. Utilizing a spherical-earth dislocation model we documented that large-scale signals predominate in the global field of the coseismic gravity changes caused by the earthquake. After removing the near-field effects, the coseismic gravity changes show a negative anomaly feature with an average magnitude of -0.18×10^-8 m·s^-2 in the region ranging ~40° from the epicenter, which is considered as the "medium field" in this study. From the GRACE data released by Center for Space Research from August 2002 to December 2008, we retrieved the large-scale gravity changes smoothed with 3 000 km Gaussian filter. The results show that the coseismic gravity changes detected by GRACE in the medium field have an average of (-0.20€±0.06)×10^-8 m·s^-2, which agrees with the model prediction. The detection confirms that GRACE is sensitive to large-scale medium-field coseismic gravitational effects of mega earthquakes, and also validates the spherical-earth dislocation model in the medium field from the perspective of satellite gravimetry.     

滇西地区GPS时间序列中陆地水载荷形变干扰的GRACE分辨与剔除

利用“中国大陆构造环境监测网络”在云南西部地区的13个连续GPS观测站和法国空间大地测量研究组Space Geodesy Research Group）的GRACE时变重力场资料,定量分析了该区域陆地水载荷所产生的非构造形变的量值和变化特点,探讨了利用GRACE分辨和剔除GPS观测中陆地水负荷所引起的非构造形变干扰的依据和模型.结果表明：滇西地区GPS坐标变化时间序列的垂向分量中,普遍包含有明显的年周期非构造形变波动,高值可达12mm,其中约42%源于陆地水迁徙变化所引起的负荷形变;通过主成份分析方法所获取的区域GPS共模误差与GRACE陆地水载荷形变序列的相关性高达0.87,若以GRACE扣除陆地水负荷形变,则滇西地区GPS网共模误差可消除约64%,且物理机制明确.然而,由于目前的GRACE只能有效分辨大约400km范围内陆地水载荷的整体变化,所以对于各GPS站点更加局部化的陆地水负荷非构造形变干扰,尚无法进行有效分辨.     

Large-scale aseismic creep in the areas of the strong earthquakes revealed from the GRACE data on the time variations of the Earth’s gravity field

The refinement of the accuracy and resolution of the monthly global gravity field models from the GRACE satellite mission, together with the accumulation of more than a decade-long series of these models, enabled us to reveal the processes that occur in the regions of large (M_w≥8) earthquakes that have not been studied previously. The previous research into the time variations of the gravity field in the regions of the giant earthquakes, such as the seismic catastrophes in Sumatra (2004) and Chile (2010), and the Tohoku mega earthquake in Japan (2011), covered the coseismic gravity jump followed by the long postseismic changes reaching almost the same amplitude. The coseismic gravity jumps resulting from the lower-magnitude events are almost unnoticeable. However, we have established a long steady growth of gravity anomalies after a number of such earthquakes. For instance, in the regions of the subduction earthquakes, the growth of the positive gravity anomaly above the oceanic trench was revealed after two events with magnitudes M_w=8.5 in the Sumatra region (the Nias earthquake of March 2005 and the Bengkulu event of September 2007 near the southern termination of Sumatra Island), after the earthquake with M_w=8.5 on Hokkaido in September 2007, a doublet Simushir earthquake with the magnitudes M_w = 8.3 and 8.1 in the Kuriles in November 2006 and January 2007, and after the earthquake off the Samoa Island in September 2009 (M_w=8.1). The steady changes in the gravity field have also been recorded after the earthquake in the Sichuan region (May 2008, M_w = 8.0) and after the doublet event with magnitudes 8.6 and 8.2, which occurred in the Wharton Basin of the Indian Ocean on April 11, 2012. The detailed analysis of the growth of the positive anomaly in gravity after the Simushir earthquake of November 2006 is presented. The growth started a few months after the event synchronously with the seismic activation on the downdip extension of the coseismically ruptured fault plane zone. The data demonstrating the increasing depth of the aftershocks since March 2007 and the approximately simultaneous change in the direction and average velocity of the horizontal surface displacements at the sites of the regional GPS network indicate that this earthquake induced postseismic displacements in a huge area extending to depths below 100 km. The total displacement since the beginning of the growth of the gravity anomaly up to July 2012 is estimated at 3.0 m in the upper part of the plate’s contact and 1.5 m in the lower part up to a depth of 100 km. With allowance for the size of the region captured by the deformations, the released total energy is equivalent to the earthquake with the magnitude M_w = 8.5. In our opinion, the growth of the gravity anomaly in these regions indicates a large-scale aseismic creep over the areas much more extensive than the source zone of the earthquake. These processes have not been previously revealed by the ground-based techniques. Hence, the time series of the GRACE gravity models are an important source of the new data about the locations and evolution of the locked segments of the subduction zones and their seismic potential.