利用流动GPS测定2011年日本MW9.0地震远场同震位移
赵静旸 武艳强 梁洪宝 杜凯夫 时爽爽
摘要： 采用1999—2017年我国东北和华北地区的流动GPS观测资料，对比坐标时间序列直接解算、连续GPS同震位移插值结果约束解算、位错模型模拟同震位移约束解算和由上述2种外部数据同时约束解算等4种解算策略，获得了2011年日本MW9.0地震远场同震位移，并由该远场同震位移约束反演了此次地震的静态位错。结果表明：在利用流动GPS解算量级较小的同震位移时，由外部可靠数据作为约束可获得较为理想的结果；流动GPS测定的2011年日本MW9.0地震在我国东北和华北地区造成的远场同震位移最大可达33 mm，以东向位移为主；由远场同震水平位移和近场同震垂直位移约束得到的静态位错对应的矩震级MW为8.77，最大同震位错量为21.09 m。
关键词：流动GPS；约束解算；同震位移；日本MW9.0地震
中图分类号：P315.725 文献标识码：A 文章编号：1000-0666（2018）03-0451-06
0 引言
精确的震区形变场，有助于构建精细的地震破裂模型、认识和理解地震发生机理。GPS、InSAR等现代大地测量技术，具有观测连续性好、空间覆盖率高等特点，已被广泛应用于多次地震的同震变形特征分析工作中（万永革等，2004；郝明等，2012；单新建等，2014；杨国华等，2015；占伟等，2015；Wu et al，2016；季灵运等，2017）。2011年3月11日，日本海沟地区（38.3°N，142.4°E）发生MW9.0地震，并引起了强烈的海啸。得益于日本良好的地震监测基础，此次地震在近场引起的同震地壳形变，通过GPS、InSAR、海底形变观测等手段得到了很好的解读，相关结果已应用于此次地震同震位错、破裂过程反演等研究（刁法启等，2011，2012；Wang et al，2012；Bletery et al，2015）。
大地震除了会在震中区形成位移之外，还会导致大区域的变形，如Plafker（1965）发现1964年阿拉斯加地震的同震位移发生在俯冲带400～600 km范围内，王敏等（2006）的研究结果表明2004年印尼苏门答腊MW9.3地震的同震位移东西、南北的范围均达到6 000～7 000 km。相关研究表明，日本地震同样在中国大陆产生了永久的变形影响：王敏等（2011）的研究结果显示此次地震造成我国东北和华北地区产生毫米至厘米级的同震水平位移，最大值为35 mm；杨少敏等（2011）认为此次地震水平位移影响范围达2 000 km，造成中国东北东移1～3 cm，华北东移3～8 mm，朝鲜半岛东移2 cm；蔡华等（2012）解算分析了整个中国大陆地区的远场同震位移，认为东北地区的同震位移接近20 cm，华北地区同震位移在10 cm左右，鄂尔多斯及周围、阿尔金—祁连山—阿拉善、新疆等地区同震位移为5 mm左右。以上研究成果均是通过GPS连续观测数据获得，由于测站分布稀疏，上述分析仅为区域概略结果，难以了解此次地震在中国大陆远场同震形变细部的大小与空间的变化，而远场同震位移的获取对研究地壳构造运动、地震动力学特征及判定未来地震趋势均具有重要意义。
随着中国地壳运动观测网络和中国大陆构造环境监测网络的建立，中国大陆目前建立了相对密集的GPS流动监测网络，这些测站自1999年开始观测，目前已在2011年前后积累了相对丰富的观测数据，为利用流动GPS资料解算此次日本巨震在中国大陆的远场同震形变提供了条件。本文充分利用已有密集流动GPS观测资料，对比相关数据处理方法，提取受此次日本巨震影响显著的我国华北和东北地区的同震形变信息，并利用该远场形变作为约束反演此次地震的静态位错分布。
1 GPS数据处理
根据前人的研究成果，日本MW9.0地震在中国大陆产生的同震变形在甘肃、四川边界以及华南块体附近不显著，因此，本文选定的流动GPS数据分布在（110°～135°E，30°～55°N）范围内。该区域内流动GPS一般逐年或隔年观测1期，每期至少观测4 d，为了保证解算结果的可靠性，选取的测站自1999年开始观测，观测截止时间为2017年，在2011年3月11日前后至少有3期观测数据。经统计，有318个测站符合选定标准。基础数据处理由GAMIT/GLOBK软件完成，先以单天观测为一解算时段，采用双差模式处理得到单日松弛解，继而选定全球80个IGS站作为框架点，处理得到各测站在ITRF2014框架下的坐标时间序列。
式中：t0为任意选定的初始时刻；y（t0）为t0时刻的位置；v为构造运动速度；三角函数表示非构造运动引起的周期信号；aj，fj，φj分别表示周期信号的振幅、频率和初相位；ri为未能模型化的观测噪声；teq表示地震发生时刻；c为地震的同震位移；对数函数表示震后变形；p为震后位移系数；τ为震后弛豫时间。由于采用的流动GPS观测间隔较长，且远场变形量级较小，本文忽略震后变形和周期性的非构造变形，假设地震发生前后的构造運动速度不同，将y（t0），v1，v2和c作为待定参数通过最小二乘进行求解。由于流动观测点获得的垂向位移缺乏可靠性，本文的求解仅针对水平坐标分量。
2 流动GPS同震位移解算策略对比分析 为了克服流动GPS数据观测间隔时间长、观测误差不易剔除等缺陷，获取可靠的同震形变信息，对比分析4种解算策略，分别为：利用坐标时间序列直接解算；利用连续GPS观测计算的同震位移在流动点位上的插值进行约束解算；利用位错模型模拟的流动点位上的同震位移约束解算；同时采用上述2种外部同震位移约束解算。
利用坐标时间序列直接解算时，由各坐标值的误差确定其权重，同震位移解算结果如图1a所示。图1a显示2011年日本MW9.0地震引起的研究区内的同震位移以东向运动为主，位移整体上稍显杂乱，反映了不同区域对同一构造活动的敏感度不同。图2a为直接解算获得的东向位移分量与经度分布的关系图，显示随着震中距的加大，位移量不断减小，个别测站出现表示西向位移的负值。Wang等（2014）的研究表明，弹性半空间中逆冲破裂造成的同震水平位移符合随着震中距的加大而指数下降的变化特征，因此我们利用如式（3）所示的指数函数对东向位移和经度进行拟合，以拟合值和实测值的差异水平作为评价所获同震位移精度的标准：
式中：D表示东向位移分量；lon表示经度；a和b表示指数函数系数。上述直接解算的同震位移拟合得到拟合值和实测值差异的均方根值为5.78 mm。
地震发生后，短则数小时之内，长则数天之内，即可利用连续GPS观测获取较高精度的同震形变信息，但受限于连续测站的分布，此类信息往往空间分辨率不足。将Wang等（2014）解算的分布于中国大陆的陆态网络基准站的同震位移，通过最小二乘配置法插值到流动观测点位上，对于每一个单个点位，将插值结果作为一个观测方程参与观测方程组的解算，对流动点位的同震位移进行约束解算，解算时除流动观测本身的误差外，还考虑连续站同震位移本身的精度。获取的同震位移如图1b所示，显示在此约束条件下获得的同震位移表现出有序的特性，整体上仍然以东向位移为主，位移量值从3.3 mm到33.1 mm不等。同样用指数函数拟合东向位移分量与经度的关系，如图2b所示，拟合值和实测值差异的均方根值为2.90 mm。
除实测的同震位移外，地震发生之后可以利用地震波等手段快速获取初步的断层运动信息，继而可以通过位错模型模拟获得地表点位的同震位移，该模拟信息也可作为外部数据对流动GPS观测进行约束解算。本文的模拟同震位移通过球状分层地球位错理论（Sun et al，2010）计算获得，其中断层滑动模型采用Wei等（2012）的研究结果，地球模型采用Wang等（2014）的PREM2模型。由于模拟同震位移没有误差信息，解算时对于相应测站结果赋予2倍连续GPS插值结果的误差。获取的同震位移如图1c所示，显示由模拟数据约束解算的结果有序性优于直接解算的结果，劣于连续GPS观测约束解算的结果，个别测站仍有如东西向位移与其它测站相反的现象。指数函数拟合东向位移分量与经度的关系（图2c），拟合值和实测值差异的均方根值为3.99 mm。
除上述3种解算策略外，利用连续GPS结果和位错模型模拟结果同时约束进行流动GPS同震位移的解算，外部数据的精度影响因素设置同上，解算结果如图1d所示。该结果与仅使用连续GPS信息进行约束的解算结果相近，同震位移显示出较好的有序性，位移方向均指向震中，位移量值从3.4 mm到33.3 mm不等，最大位移出现在辽宁和吉林东部地区。利用指数函数拟合东向位移分量与经度的关系（图2d），拟合值和实测值差异的均方根值为2.81 mm。
3 静态位错反演
为了验证远场同震位移在地震静态位错反演中的作用，以由2种外部数据同时约束获得的远场水平同震位移场和由日本岛陆地连续GPS获得的垂直同震位移场为约束，反演日本MW9.0地震断层上的滑动分布，其中日本岛垂直位移场由美国ARIA研究组提供（ftp：//sideshow.jpl.nasa.gov/pub/usrs/ARIA/AIRA_coseismic_offset.v0.3.table）。反演采用SDM程序（Wang et al，2011，2013）进行，断层模型参考震中位置，北端边界设定在40°N，上边缘沿日本以东的海沟展布，深度至60 km，断层倾角在上边界设定为7°并随深度线性变化。反演时，断层被离散为10 km×10 km的若干子断层，滑动角约束在65°～115°变化，分层介质模型参数与刁法启等（2011）一致，取自CRUST2.0。具体反演过程附加断层面上的应力降平滑作为约束条件，平滑因子根据模型粗糙度和数据拟合残差之间的折中曲线确定为0.2。反演结果如图3所示，结果显示最大同震位错为21.09 m，位错以逆冲性滑动为主，平均滑动角为98.71°，相应的矩震级MW为8.77，这与刁法启等（2011，2012）利用近场陆地同震位移约束反演获得的结果相近，但滑动量和矩震级略有偏小，这是由于近场形变观测的约束能力比远场形变观测的约束能力强造成的。
4 讨论与结论
本文利用我国东北和华北地区1999—2017年的流动GPS观测资料，对比了利用坐标时序直接解算、附加连续GPS同震位移插值结果约束解算、附加位错模型模拟同震位移约束解算以及同时附加上述2种外部同震位移约束解算4种同震位移解算策略，获取2011年日本MW9.0地震在中国大陆引起的远场同震水平位移，并利用获得的远场水平位移结合近场垂直位移对震中区的静态滑动分布进行了反演。研究结果表明：
（1）直接通过流动GPS坐标时间序列解算获得的远场同震位移能够显示局部的变形特征，但由于观测间隔时间较长，流动观测易受观测环境变化等偶然因素的影响，且偶然影响引起的观测质量下降不易被有效识别，直接解算的结果中难免含有一定程度的误差，特别是对于华北等本身构造运动量级不大、人类活动丰富的地区。
（2）附加可靠的外部約束条件，可以有效改善流动GPS资料解算同震位移的精度，如利用连续GPS同震位移插值结果约束和利用位错模型模拟同震位移结果约束等解算方案，其中由于连续GPS解算的同震位移的精度较高，4种解算策略中，有该资料参与约束的两种解算方案得到的同震位移结果相近、精度最高。解算结果显示2011年日本MW9.0地震在我国东北和华北地区引起的同震变形以东向位移为主，位移方向指向震中，位移量值从3 mm到33 mm不等，最大位移出现在辽宁和吉林东部地区。
（3）远场同震形变可以有效地辅助震源区同震位错特征信息的获取，但以远场形变约束为主的反演结果较以近场形变约束的反演结果的滑动量级和震级均较小，本文反演的静态位错对应的矩阵级MW为8.77，最大同震位错量为21.09 m。
参考文献：
蔡华，孙汉荣，赵齐乐，等.2012.GPS测定的2011年日本9.0级地震的中国大陆地区同震位移场[J].武汉大学学报（信息科学版），37（8）：953-955.
单新建，屈春燕，郭利民，等.2014.基于InSAR与GPS观测的汶川同震垂直形变场的获取[J].地震地质，36（3）：718-730.
刁法啟，熊熊，倪四道，等.2011.利用GPS位移反演日本MW9.0仙台地震及MW7.9强余震静态位错模型[J].科学通报，56（24）：1999-2005.
刁法启，熊熊，郑勇.2012.MW9.0日本Tohoku大地震静态位错模型：陆地GPS资料和海底GPS/Acoustic资料联合反演的结果[J].科学通报，57（18）：1676-1683.
国家重大科学工程“中国地壳运动观测网络”项目组.2008.GPS测定的2008年汶川MS8.0级地震的同震位移场[J].中国科学：地球科学，38（10）：1195-1206.
郝明，王庆良，崔笃信，等.2012.根据GPS和水准资料反演2008年汶川地震同震破裂模型[J].地震研究，35（3）：323-329.
季灵运，刘传金，徐晶，等.2017.九寨沟MS7.0地震的InSAR观测及发震构造分析[J].地球物理学报，60（10）：4069-4082.
万永革，王敏，沈正康，等.2004.利用GPS和水准测量资料反演2001年昆仑山口西8.1级地震的同震滑动分布[J].地震地质，26（3）：393-404.
王敏，李强，王凡，等.2011.全球定位系统测定的2011年日本宫城MW9.0级地震远场同震位移[J].科学通报，56（20）：1593-1596.
王敏，张培震，沈正康，等.2006.全球定位系统（GPS）测定的印尼苏门达腊巨震的远场同震地表位移[J].科学通报，51（3）：365-368.
杨国华，朱爽，梁洪宝，等.2015.芦山MS7.0级地震震前及同震地表形变[J].武汉大学学报（信息科学版），40（1）：121-127.
杨少敏，聂兆生，贾志革，等.2011.GPS解算的日本MW9.0级地震的远场同震地表位移[J].武汉大学学报（信息科学版），36（11）：1336-1339.
占伟，武艳强，梁洪宝，等.2015.GPS观测结果反映的尼泊尔MW7.8地震孕震特征[J].地球物理学报，58（5）：1818-1826.
Bletery Q，Sladen A，Delouis B，et al.2015.A detailed source model for the MW9.0 Tohoku‐Oki earthquake reconciling geodesy，seismology，and tsunami records[J].Journal of Geophysical Research Solid Earth，119（10）：7636-7653.
Plafker G.1965.Tectonic Deformation Associated with the 1964 Alaska Earthquake：The earthquake of 27 March 1964 resulted in observable crustal deformation of unprecedented areal extent[J].Science，148（3678）：1675-1687.
Sun W，Okubo S，Fu G，et al.2010.General formulations of global co-seismic deformations caused by an arbitrary dislocation in a spherically symmetric earth model—applicable to deformed earth surface and space‐fixed point[J].Geophysical Journal of the Royal Astronomical Society，177（3）：817-833.
Wang C，Ding X，Shan X，et al.2012.Slip distribution of the 2011 Tohoku earthquake derived from joint inversion of GPS，InSAR and seafloor GPS/acoustic measurements[J].Journal of Asian Earth Sciences，57（6）：128-136.
Wang R，Schurr B，Milkereit C，et al.2011.An Improved Automatic Scheme for Empirical Baseline Correction of Digital Strong-Motion Records[J].Bulletin of the Seismological Society of America，101（5）：2029-2044.
Wang R，Parolai S，Ge M，et al.2013.The 2011 MW9.0 Tohoku Earthquake：Comparison of GPS and Strong-Motion Data[J].Bulletin of the Seismological Society of America，103（2B）：1336-1347.
Wang W，Sun W，Wu Y，et al.2014.Modification of fault slip models of the MW9.0 Tohoku Earthquake by far field GPS observations[J].Journal of Geodynamics，75（4）：22-33.
Wei S，Graves R，Helmberger D，et al.2012.Sources of shaking and flooding during the Tohoku-Oki earthquake：A mixture of rupture styles[J].Earth & Planetary Science Letters，333-334（6）：91-100.
Wu Y，Jiang Z，Liang H，et al.2016.Coseismic deformations of the 2015 MW7.8 Gorkha earthquake and interseismic strain accumulation in the Himalayan tectonic belt and Tibetan plateau[J].Tectonophysics，670：144-154.
Abstract The campaign GPS observations in the northeastern China from 1999 to 2017 are used to derive the far-field coseismic displacements of the 2011 Tohoku-oki MW9.0 earthquake，Japan.The strategies include direct calculation with the coordinate time series，constrained calculation with the interpolated coseismic displacements from continuous GPS observations，constrained calculation with the simulated coseismic displacements from dislocation model and constrained calculation with both the above two external data are compared.Then the best far-field coseismic horizontal displacements as well as the near field coseismic vertical displacements are used to invert the static dislocation model of the earthquake.The results show that when the campaign GPS data are used to inverse the coseismic displacements of small magnitude，using external reliable data as constraint is conducive to getting more satisfactory results.And the maximum far-field coseismic displacements in northeastern China can reach 33 mm and the dominated displacements is eastward.The resultant moment magnitude corresponding to the coseismic dislocation inversion is MW8.77 and the largest dislocation of this earthquake is 21.09 m.
Keywords：campaign GPS observation；constrained calculation；coseismic displacement；the Tohoku-oki MW9.0 earthquake
摘要： 采用1999—2017年我国东北和华北地区的流动GPS观测资料，对比坐标时间序列直接解算、连续GPS同震位移插值结果约束解算、位错模型模拟同震位移约束解算和由上述2种外部数据同时约束解算等4种解算策略，获得了2011年日本MW9.0地震远场同震位移，并由该远场同震位移约束反演了此次地震的静态位错。结果表明：在利用流动GPS解算量级较小的同震位移时，由外部可靠数据作为约束可获得较为理想的结果；流动GPS测定的2011年日本MW9.0地震在我国东北和华北地区造成的远场同震位移最大可达33 mm，以东向位移为主；由远场同震水平位移和近场同震垂直位移约束得到的静态位错对应的矩震级MW为8.77，最大同震位错量为21.09 m。
关键词：流动GPS；约束解算；同震位移；日本MW9.0地震
中图分类号：P315.725 文献标识码：A 文章编号：1000-0666（2018）03-0451-06
0 引言
精确的震区形变场，有助于构建精细的地震破裂模型、认识和理解地震发生机理。GPS、InSAR等现代大地测量技术，具有观测连续性好、空间覆盖率高等特点，已被广泛应用于多次地震的同震变形特征分析工作中（万永革等，2004；郝明等，2012；单新建等，2014；杨国华等，2015；占伟等，2015；Wu et al，2016；季灵运等，2017）。2011年3月11日，日本海沟地区（38.3°N，142.4°E）发生MW9.0地震，并引起了强烈的海啸。得益于日本良好的地震监测基础，此次地震在近场引起的同震地壳形变，通过GPS、InSAR、海底形变观测等手段得到了很好的解读，相关结果已应用于此次地震同震位错、破裂过程反演等研究（刁法启等，2011，2012；Wang et al，2012；Bletery et al，2015）。
大地震除了会在震中区形成位移之外，还会导致大区域的变形，如Plafker（1965）发现1964年阿拉斯加地震的同震位移发生在俯冲带400～600 km范围内，王敏等（2006）的研究结果表明2004年印尼苏门答腊MW9.3地震的同震位移东西、南北的范围均达到6 000～7 000 km。相关研究表明，日本地震同样在中国大陆产生了永久的变形影响：王敏等（2011）的研究结果显示此次地震造成我国东北和华北地区产生毫米至厘米级的同震水平位移，最大值为35 mm；杨少敏等（2011）认为此次地震水平位移影响范围达2 000 km，造成中国东北东移1～3 cm，华北东移3～8 mm，朝鲜半岛东移2 cm；蔡华等（2012）解算分析了整个中国大陆地区的远场同震位移，认为东北地区的同震位移接近20 cm，华北地区同震位移在10 cm左右，鄂尔多斯及周围、阿尔金—祁连山—阿拉善、新疆等地区同震位移为5 mm左右。以上研究成果均是通过GPS连续观测数据获得，由于测站分布稀疏，上述分析仅为区域概略结果，难以了解此次地震在中国大陆远场同震形变细部的大小与空间的变化，而远场同震位移的获取对研究地壳构造运动、地震动力学特征及判定未来地震趋势均具有重要意义。
随着中国地壳运动观测网络和中国大陆构造环境监测网络的建立，中国大陆目前建立了相对密集的GPS流动监测网络，这些测站自1999年开始观测，目前已在2011年前后积累了相对丰富的观测数据，为利用流动GPS资料解算此次日本巨震在中国大陆的远场同震形变提供了条件。本文充分利用已有密集流动GPS观测资料，对比相关数据处理方法，提取受此次日本巨震影响显著的我国华北和东北地区的同震形变信息，并利用该远场形变作为约束反演此次地震的静态位错分布。
1 GPS数据处理
根据前人的研究成果，日本MW9.0地震在中国大陆产生的同震变形在甘肃、四川边界以及华南块体附近不显著，因此，本文选定的流动GPS数据分布在（110°～135°E，30°～55°N）范围内。该区域内流动GPS一般逐年或隔年观测1期，每期至少观测4 d，为了保证解算结果的可靠性，选取的测站自1999年开始观测，观测截止时间为2017年，在2011年3月11日前后至少有3期观测数据。经统计，有318个测站符合选定标准。基础数据处理由GAMIT/GLOBK软件完成，先以单天观测为一解算时段，采用双差模式处理得到单日松弛解，继而选定全球80个IGS站作为框架点，处理得到各测站在ITRF2014框架下的坐标时间序列。
式中：t0为任意选定的初始时刻；y（t0）为t0时刻的位置；v为构造运动速度；三角函数表示非构造运动引起的周期信号；aj，fj，φj分别表示周期信号的振幅、频率和初相位；ri为未能模型化的观测噪声；teq表示地震发生时刻；c为地震的同震位移；对数函数表示震后变形；p为震后位移系数；τ为震后弛豫时间。由于采用的流动GPS观测间隔较长，且远场变形量级较小，本文忽略震后变形和周期性的非构造变形，假设地震发生前后的构造運动速度不同，将y（t0），v1，v2和c作为待定参数通过最小二乘进行求解。由于流动观测点获得的垂向位移缺乏可靠性，本文的求解仅针对水平坐标分量。
2 流动GPS同震位移解算策略对比分析 为了克服流动GPS数据观测间隔时间长、观测误差不易剔除等缺陷，获取可靠的同震形变信息，对比分析4种解算策略，分别为：利用坐标时间序列直接解算；利用连续GPS观测计算的同震位移在流动点位上的插值进行约束解算；利用位错模型模拟的流动点位上的同震位移约束解算；同时采用上述2种外部同震位移约束解算。
利用坐标时间序列直接解算时，由各坐标值的误差确定其权重，同震位移解算结果如图1a所示。图1a显示2011年日本MW9.0地震引起的研究区内的同震位移以东向运动为主，位移整体上稍显杂乱，反映了不同区域对同一构造活动的敏感度不同。图2a为直接解算获得的东向位移分量与经度分布的关系图，显示随着震中距的加大，位移量不断减小，个别测站出现表示西向位移的负值。Wang等（2014）的研究表明，弹性半空间中逆冲破裂造成的同震水平位移符合随着震中距的加大而指数下降的变化特征，因此我们利用如式（3）所示的指数函数对东向位移和经度进行拟合，以拟合值和实测值的差异水平作为评价所获同震位移精度的标准：
式中：D表示东向位移分量；lon表示经度；a和b表示指数函数系数。上述直接解算的同震位移拟合得到拟合值和实测值差异的均方根值为5.78 mm。
地震发生后，短则数小时之内，长则数天之内，即可利用连续GPS观测获取较高精度的同震形变信息，但受限于连续测站的分布，此类信息往往空间分辨率不足。将Wang等（2014）解算的分布于中国大陆的陆态网络基准站的同震位移，通过最小二乘配置法插值到流动观测点位上，对于每一个单个点位，将插值结果作为一个观测方程参与观测方程组的解算，对流动点位的同震位移进行约束解算，解算时除流动观测本身的误差外，还考虑连续站同震位移本身的精度。获取的同震位移如图1b所示，显示在此约束条件下获得的同震位移表现出有序的特性，整体上仍然以东向位移为主，位移量值从3.3 mm到33.1 mm不等。同样用指数函数拟合东向位移分量与经度的关系，如图2b所示，拟合值和实测值差异的均方根值为2.90 mm。
除实测的同震位移外，地震发生之后可以利用地震波等手段快速获取初步的断层运动信息，继而可以通过位错模型模拟获得地表点位的同震位移，该模拟信息也可作为外部数据对流动GPS观测进行约束解算。本文的模拟同震位移通过球状分层地球位错理论（Sun et al，2010）计算获得，其中断层滑动模型采用Wei等（2012）的研究结果，地球模型采用Wang等（2014）的PREM2模型。由于模拟同震位移没有误差信息，解算时对于相应测站结果赋予2倍连续GPS插值结果的误差。获取的同震位移如图1c所示，显示由模拟数据约束解算的结果有序性优于直接解算的结果，劣于连续GPS观测约束解算的结果，个别测站仍有如东西向位移与其它测站相反的现象。指数函数拟合东向位移分量与经度的关系（图2c），拟合值和实测值差异的均方根值为3.99 mm。
除上述3种解算策略外，利用连续GPS结果和位错模型模拟结果同时约束进行流动GPS同震位移的解算，外部数据的精度影响因素设置同上，解算结果如图1d所示。该结果与仅使用连续GPS信息进行约束的解算结果相近，同震位移显示出较好的有序性，位移方向均指向震中，位移量值从3.4 mm到33.3 mm不等，最大位移出现在辽宁和吉林东部地区。利用指数函数拟合东向位移分量与经度的关系（图2d），拟合值和实测值差异的均方根值为2.81 mm。
3 静态位错反演
为了验证远场同震位移在地震静态位错反演中的作用，以由2种外部数据同时约束获得的远场水平同震位移场和由日本岛陆地连续GPS获得的垂直同震位移场为约束，反演日本MW9.0地震断层上的滑动分布，其中日本岛垂直位移场由美国ARIA研究组提供（ftp：//sideshow.jpl.nasa.gov/pub/usrs/ARIA/AIRA_coseismic_offset.v0.3.table）。反演采用SDM程序（Wang et al，2011，2013）进行，断层模型参考震中位置，北端边界设定在40°N，上边缘沿日本以东的海沟展布，深度至60 km，断层倾角在上边界设定为7°并随深度线性变化。反演时，断层被离散为10 km×10 km的若干子断层，滑动角约束在65°～115°变化，分层介质模型参数与刁法启等（2011）一致，取自CRUST2.0。具体反演过程附加断层面上的应力降平滑作为约束条件，平滑因子根据模型粗糙度和数据拟合残差之间的折中曲线确定为0.2。反演结果如图3所示，结果显示最大同震位错为21.09 m，位错以逆冲性滑动为主，平均滑动角为98.71°，相应的矩震级MW为8.77，这与刁法启等（2011，2012）利用近场陆地同震位移约束反演获得的结果相近，但滑动量和矩震级略有偏小，这是由于近场形变观测的约束能力比远场形变观测的约束能力强造成的。
4 讨论与结论
本文利用我国东北和华北地区1999—2017年的流动GPS观测资料，对比了利用坐标时序直接解算、附加连续GPS同震位移插值结果约束解算、附加位错模型模拟同震位移约束解算以及同时附加上述2种外部同震位移约束解算4种同震位移解算策略，获取2011年日本MW9.0地震在中国大陆引起的远场同震水平位移，并利用获得的远场水平位移结合近场垂直位移对震中区的静态滑动分布进行了反演。研究结果表明：
（1）直接通过流动GPS坐标时间序列解算获得的远场同震位移能够显示局部的变形特征，但由于观测间隔时间较长，流动观测易受观测环境变化等偶然因素的影响，且偶然影响引起的观测质量下降不易被有效识别，直接解算的结果中难免含有一定程度的误差，特别是对于华北等本身构造运动量级不大、人类活动丰富的地区。
（2）附加可靠的外部約束条件，可以有效改善流动GPS资料解算同震位移的精度，如利用连续GPS同震位移插值结果约束和利用位错模型模拟同震位移结果约束等解算方案，其中由于连续GPS解算的同震位移的精度较高，4种解算策略中，有该资料参与约束的两种解算方案得到的同震位移结果相近、精度最高。解算结果显示2011年日本MW9.0地震在我国东北和华北地区引起的同震变形以东向位移为主，位移方向指向震中，位移量值从3 mm到33 mm不等，最大位移出现在辽宁和吉林东部地区。
（3）远场同震形变可以有效地辅助震源区同震位错特征信息的获取，但以远场形变约束为主的反演结果较以近场形变约束的反演结果的滑动量级和震级均较小，本文反演的静态位错对应的矩阵级MW为8.77，最大同震位错量为21.09 m。
参考文献：
蔡华，孙汉荣，赵齐乐，等.2012.GPS测定的2011年日本9.0级地震的中国大陆地区同震位移场[J].武汉大学学报（信息科学版），37（8）：953-955.
单新建，屈春燕，郭利民，等.2014.基于InSAR与GPS观测的汶川同震垂直形变场的获取[J].地震地质，36（3）：718-730.
刁法啟，熊熊，倪四道，等.2011.利用GPS位移反演日本MW9.0仙台地震及MW7.9强余震静态位错模型[J].科学通报，56（24）：1999-2005.
刁法启，熊熊，郑勇.2012.MW9.0日本Tohoku大地震静态位错模型：陆地GPS资料和海底GPS/Acoustic资料联合反演的结果[J].科学通报，57（18）：1676-1683.
国家重大科学工程“中国地壳运动观测网络”项目组.2008.GPS测定的2008年汶川MS8.0级地震的同震位移场[J].中国科学：地球科学，38（10）：1195-1206.
郝明，王庆良，崔笃信，等.2012.根据GPS和水准资料反演2008年汶川地震同震破裂模型[J].地震研究，35（3）：323-329.
季灵运，刘传金，徐晶，等.2017.九寨沟MS7.0地震的InSAR观测及发震构造分析[J].地球物理学报，60（10）：4069-4082.
万永革，王敏，沈正康，等.2004.利用GPS和水准测量资料反演2001年昆仑山口西8.1级地震的同震滑动分布[J].地震地质，26（3）：393-404.
王敏，李强，王凡，等.2011.全球定位系统测定的2011年日本宫城MW9.0级地震远场同震位移[J].科学通报，56（20）：1593-1596.
王敏，张培震，沈正康，等.2006.全球定位系统（GPS）测定的印尼苏门达腊巨震的远场同震地表位移[J].科学通报，51（3）：365-368.
杨国华，朱爽，梁洪宝，等.2015.芦山MS7.0级地震震前及同震地表形变[J].武汉大学学报（信息科学版），40（1）：121-127.
杨少敏，聂兆生，贾志革，等.2011.GPS解算的日本MW9.0级地震的远场同震地表位移[J].武汉大学学报（信息科学版），36（11）：1336-1339.
占伟，武艳强，梁洪宝，等.2015.GPS观测结果反映的尼泊尔MW7.8地震孕震特征[J].地球物理学报，58（5）：1818-1826.
Bletery Q，Sladen A，Delouis B，et al.2015.A detailed source model for the MW9.0 Tohoku‐Oki earthquake reconciling geodesy，seismology，and tsunami records[J].Journal of Geophysical Research Solid Earth，119（10）：7636-7653.
Plafker G.1965.Tectonic Deformation Associated with the 1964 Alaska Earthquake：The earthquake of 27 March 1964 resulted in observable crustal deformation of unprecedented areal extent[J].Science，148（3678）：1675-1687.
Sun W，Okubo S，Fu G，et al.2010.General formulations of global co-seismic deformations caused by an arbitrary dislocation in a spherically symmetric earth model—applicable to deformed earth surface and space‐fixed point[J].Geophysical Journal of the Royal Astronomical Society，177（3）：817-833.
Wang C，Ding X，Shan X，et al.2012.Slip distribution of the 2011 Tohoku earthquake derived from joint inversion of GPS，InSAR and seafloor GPS/acoustic measurements[J].Journal of Asian Earth Sciences，57（6）：128-136.
Wang R，Schurr B，Milkereit C，et al.2011.An Improved Automatic Scheme for Empirical Baseline Correction of Digital Strong-Motion Records[J].Bulletin of the Seismological Society of America，101（5）：2029-2044.
Wang R，Parolai S，Ge M，et al.2013.The 2011 MW9.0 Tohoku Earthquake：Comparison of GPS and Strong-Motion Data[J].Bulletin of the Seismological Society of America，103（2B）：1336-1347.
Wang W，Sun W，Wu Y，et al.2014.Modification of fault slip models of the MW9.0 Tohoku Earthquake by far field GPS observations[J].Journal of Geodynamics，75（4）：22-33.
Wei S，Graves R，Helmberger D，et al.2012.Sources of shaking and flooding during the Tohoku-Oki earthquake：A mixture of rupture styles[J].Earth & Planetary Science Letters，333-334（6）：91-100.
Wu Y，Jiang Z，Liang H，et al.2016.Coseismic deformations of the 2015 MW7.8 Gorkha earthquake and interseismic strain accumulation in the Himalayan tectonic belt and Tibetan plateau[J].Tectonophysics，670：144-154.
Abstract The campaign GPS observations in the northeastern China from 1999 to 2017 are used to derive the far-field coseismic displacements of the 2011 Tohoku-oki MW9.0 earthquake，Japan.The strategies include direct calculation with the coordinate time series，constrained calculation with the interpolated coseismic displacements from continuous GPS observations，constrained calculation with the simulated coseismic displacements from dislocation model and constrained calculation with both the above two external data are compared.Then the best far-field coseismic horizontal displacements as well as the near field coseismic vertical displacements are used to invert the static dislocation model of the earthquake.The results show that when the campaign GPS data are used to inverse the coseismic displacements of small magnitude，using external reliable data as constraint is conducive to getting more satisfactory results.And the maximum far-field coseismic displacements in northeastern China can reach 33 mm and the dominated displacements is eastward.The resultant moment magnitude corresponding to the coseismic dislocation inversion is MW8.77 and the largest dislocation of this earthquake is 21.09 m.
Keywords：campaign GPS observation；constrained calculation；coseismic displacement；the Tohoku-oki MW9.0 earthquake
