龙源期刊网 http://www.qikan.com.cn
2011年日本9。0级地震前后GPS基线时间序列分析
作者：张风霜 占伟 孙东颖 来源：《地震研究》2012年第02期
摘要：利用1999年初至2011年9月中国地壳运动观测网络基准站及中国大陆周边包含日本境内的IGS基准站的资料，分析了2011年3月11日日本
级地震前日本境内及周边
IGS基准站基线时间序列的变化情况。讨论了该地震对中国大陆GPS基准站的同震影响，并对地震后中国大陆GPS基线时序结果的变化情况进行了总结分析。通过分析认为，日本大震后我国东部地区部分GPS基线开始转折拉伸，因而一定程度上缓解了这些基线的趋势性缩短，但是东部地区仍然有部分GPS基线的变化趋势基本不受日本地震的影响，近几年以来的趋势性异常至今仍然存在。
关键词：GPS基准站；基线时间序列；构造形变场；日本 中图分类号： 0引言
2011年3月11日，日本本州岛北部以东海域发生
级地震，此次地震及其引发的海啸
级地震
文献标识码：A文章编号：1000-0666（2012）02-0190-11
造成了极大的破坏和严重的人员伤亡。据日本官方统计，实际死亡人数超过2万7千人，本州岛东北部地区大量沿海村镇被地震引发的海啸卷入海中。地震对福岛核电站4个核反应堆造成破坏，导致大气和海水遭受放射性污染。
由日本境内GPS观测数据计算得到的同震位移场结果显示，本州岛北部地区存在显著的东向位移，水平位移量最大达5～6 m（尹继尧，朱元清，2011）。日本
级地震对我国的
GPS基准站同震影响有多大？日本境内及周边的IGS基准站基线时间序列是否包含了此次大震的孕震形变信息？本文利用国家重大科技基础设施项目“中国大陆构造环境监测网络”积累的10年多的GPS连续观测资料以及中国大陆周边包含日本境内的IGS基准站数据，分析此次地震对中国大陆GPS基准站产生的同震影响，对地震前日本境内及周边IGS基准站基线时序结果的趋势性异常变化情况进行分析，并通过对中国大陆内部GPS基准站基线时间序列近期变化情况的总结分析，初步探讨日本 1GPS数据处理
本文使用的数据是1999年初至2011年9月25日的中国地壳运动观测网络基准站及中国大陆周边包含日本境内的IGS基准站记录的观测资料，GPS连续站观测值的数据处理采用GAMIT/GLOBK软件完成（Herring et al，2010）。数据处理的基本流程：首先利用GAMIT获
级地震对中国大陆构造形变场的影响。
龙源期刊网 http://www.qikan.com.cn
得网络工程33个连续站及中国大陆周围IGS测站的区域单日松弛解（王敏等，2005a，2005b，2007）。GAMIT计算中加入周边16个IGS站（IRKT、KIT3、USUD、SUWN、ULAB、TSKB、TWTF、DAEJ、KHAJ、IISC、MIZU、PIMO、POL2、SELE、YAKT和YSSK）相应时间段内记录的数据一并进行处理。再利用GLOBK将SOPAC给出的全球单日松弛解和计算所得的区域单日松弛解进行综合平差，在此基础上通过IGS核心站求解相对于全球参考框架ITRF2005的相似变换7参数，从而获得ITRF2005下的单日解，即GLOBK的NEU坐标值，不同于通常的站心坐标，其具体含义可参阅黄立人等（2006）的研究。要计算两点之间的大地线（武艳强等，2007，2009），须先获取两个端点精确经纬度，具体计算程序可参阅GAMIT/GLOBK程序包，最后根据大地主题解算公式（孔祥元等，2001）就可以计算得到两点间的大地线长度。
级地震前日本境内及周边IGS基准站基线时序结果分析2011年3月11日日本
级地震是自20世纪有仪器记录以来的第4大地震。日本列岛地处太平洋西北，是太平洋板块、北美板块、菲律宾板块和欧亚板块的交汇处。此次地震发生在太平洋板块与北美板块交界的俯冲带上，在该区域太平洋板块以80 mm/a的速率相对欧亚板块向西运动（DeMets et al，1994），插入北美板块下部。 为了探讨日本
级大震前GPS基线的变化情况，本文选择了部分IGS基准站进行分
析。包括中俄边界的基准站KHAJ，中国东北的基准站CHAN，朝鲜半岛上的基准站SUWN和DAEJ，日本境内的4个IGS基准站（YSSK、MIZU、USUD和TSKB），其中离此次地震最近的IGS基准站为MIZU，其次为TSKB和USUD（图1）。 由于所选择的部分IGS测站距离日本
级地震震中相对较近，因此GPS基线时间序列
的同震变化较大，尤其是与测站MIZU相连的基线最大同震变化达3 m。同震变化的巨大量级掩盖了震前的图1日本境内及周边IGS基准站和基线分布
eference stations and
baselines in Japan and its surrounding area形变变化，为了更好地分析震前GPS基线的变化情况，本文只给出了截止到地震前一天即2011年3月10日的基线时间序列结果（图2）。图2中MIZU-TSKB和MIZU-USUD两条基线时间序列在2008年6月13日的阶跃变化是由2008年6月13日日本本州岛北部岩手—宫城里氏 地震研究35卷第2期张风霜等：2011年日本
级大震引发的同震形变。
0级地震前后GPS基线时间序列分析本
文的分析侧重于GPS基线时间序列的趋势性转折变化，即时间序列的最低频部分，反映的是构造形变成分。而基线时间序列的非构造形变周期成分（如年周期、半年周期等）、随机漫步噪声、闪烁噪声、白噪声等中高频成分可以在时间序列的后期处理中通过各种滤波方法进行分解和剔除（黄立人，2006；黄立人，符养，2007；杨国华等，2007）。由图2可看出，由于长时间趋势变化幅度远大于高频信息（包括高频噪声和各种误差），在不剔除高频信息的情况
龙源期刊网 http://www.qikan.com.cn
下，同样能显示出构造活动引起的低频趋势性变化。因此本文只讨论GPS基线时间序列的最低频部分的构造形变信息，以观察其在大震发生前可能出现的异常。
本文实际分析的时序基线分布如图1所示，限于篇幅，图2只给出了其中部分基线时序曲线。图2中有4条基线（USUD-TSKB、YSSK-TSKB、SUWN-YSSK 和DAEJ-YSSK）震前的趋势性变化出现了转折。其中YSSK-TSKB从2006年下半年开始转折转平，而SUWN-YSSK和DAEJ-YSSK两条基线也都在2006年下半年发生转折，2006年下半年之前基线的变化为线性缩短，缩短速率约10 mm/a，之后基线缩短速率明显减小，约为3 mm/a。YSSK-TSKB、SUWN-YSSK 和DAEJ-YSSK这3条基线在2006年下半年的同步转折变化不能排除YSSK单点变化的可能性。但作为IGS基准站，选点及点位保护是比较完善的，甚至很可能是基岩点，假如在同一块基岩上再建一个点，这两个点的运动会比较一致，这里所谓的单点变化，很可能已不是通常意义上的因为点位不稳定引起的单点变化，也可能是代表局部一定范围内的地壳运动。
值得关注的是日本境内的基线USUD-TSKB的变化，USUD和TSKB这两个测站距离日本
级大震震中比YSSK测站近的多（YSSK是日本境内距震中最远的一个IGS基准站），基线USUD-TSKB呈东西向展布。在1999～2008年近10年的时间内这条基线的运动都是相对平稳的线性缩短，缩短速率约4 mm/a，2008年开始，该基线的运动出现转折。2008年至2011年3月11日地震前3年多的时间内基线只有小幅的波动，运动趋势既没有持续伸长也没有缩短，处于相对闭锁的状态。震后USUD-TSKB基线出现巨幅拉伸，同震伸长了约400 mm，由于同震形变巨大，3月11日后的时间序列图本文没有给出。距离地震震中较近的MIZU-TSKB和MIZU-USUD基线震前的趋势性运动并没有出现异常，与东西向展布的USUD-TSKB基线不同，MIZU-TSKB基线近南北向展布，而MIZU-USUD基线为北北东向展布。图2日本地震前部分GPS基线时间序列的变化
（a）USUD-TSKB基线；（b）MIZU-TSKB基线；（c）MIZU-USUD基线； （d）YSSK-TSKB基线；（e）SUWN-YSSK基线；（f）DAEJ-YSSK基线
表面上看以
级
上6条基线的变化不是很一致，但从统一的机理上还是能够给出较为圆满的解释。首先，6条基线都呈缩短的变化，基线的展布方向为北东向至东西向，这与该区域主要受来自北东向北美板块的挤压与来自东向太平洋板块的挤压有关。其次，相对长边缩短的量偏大，说明在统一地壳应力场的变化中基线的缩短有弹性变形的成分。当应变一定时，总弹性变形量与基线的长度成正比，但由于地壳的不均匀性，包括塑性、蠕变等非线性变形在空间域和时间域中均广泛存在，其正比关系并不严密。其中东西向展布的基线USUD-TSKB近10年的长趋势运动特征为线性缩短，尽管缩短量不是最大的，但考虑其边长短，其应变量是最大的，显示了地壳应力应变在空间上存在着不均匀性。且该基线离震中较近，是震前受太平洋板块东西向挤压较强烈的地方，因此东西向应变最大，同震巨幅拉伸达400 mm。地震破裂发生时，强大的东西向压应力突然释放，长期积累的压应变的弹性部分必然同时释放，导致USUD-TSKB基线巨幅拉伸。