芦山MS7.0地震前水平位移和同震水平位移研究

本文利用“中国地壳运动观测网络（二期）”多个GPS连续观测站观测数据处理结果,将2013年4月20日四川芦山MS7.0地震区域参考框架同震水平位移与全球参考框架同震水平位移进行比较,结果表明两组框架解一致,说明两种参考框架均可当作位错参考框架,也即全球参考框架同震水平位移也可视为区域参考框架同震水平位移．区域参考框架下GPS连续观测站地震前的水平位移和同震水平位移结果表明,震前数年,SCTQ站西侧的GPS站构造运动十分显著,而该站水平位移却很小,即出现反常的闭锁．但该站的同震水平位移使其弹性回跳至正常构造水平位移水平,因此SCTQ站震前的位移闭锁是水平位移空间分布中的异常,是芦山MS7.0地震的前兆．水平位移时空变化表明,该站震前和震时位移完全符合里德的弹性回跳理论．区域参考框架中位移时间系列和同震水平位移的综合研究有助于对芦山地震地壳运动前兆的认识和解释．尽管本文未能直接采用其它GPS连续观测站的资料,但结合本文和其它研究结果可以证实,震中附近其它站地震前后的变化与SCTQ站类似．基于芦山地震前水平位移和同震水平位移及其与前兆关系的研究,本文进一步讨论了GPS监测网的布设、 数据处理和分析等问题．      

2014年新疆于田MS7.3地震同震位移及位错反演研究

利用于田震中300 km范围内的1个GPS连续站和12个GPS流动站数据,解算得到了2014年新疆于田M_S7.3地震地表同震位移,并反演了发震断层滑动分布,探讨此次地震对周边断裂的影响.地表同震位移结果显示,GPS观测到的同震位移范围在平行发震断裂带的北东-南西向约210 km,垂直发震断裂带的北西-南东方向约为120 km,同震位移量大于10 mm的测站位于震中距约120 km以内;同震位移特征整体表现为北东-南西方向的左旋走滑和北西-南东方向的拉张特征,其中在北东-南西方向,I069测站位移最大,约为32.1 mm,在北西-南东方向,XJYT测站位移最大,约为28.1 mm;位错反演结果表明,最大滑动位于北纬36.05°,东经82.60°,位于深部约16.6 km,最大错动量为2.75 m,反演震级为M_W7.0,同震错动呈椭圆形分布,以左旋走滑为主并具有正倾滑分量,两者最大比值约为2.5：1,同震错动延伸至地表,并向北东方向延伸,总破裂长度约50 km,地表最大错动约1.0 m;同震水平位移场模拟结果显示贡嘎错断裂、康西瓦断裂和普鲁断裂等不同位置主应变特征具有差异性,这种差异特征是否影响断裂带以及周围区域的应力构造特征,值得关注.     

GPS测定的四川九寨沟MW6.6地震同震形变

2017年8月8日四川阿坝州九寨沟发生M_W6.6地震,震源机制解显示该地震为左旋走滑型地震。对震中周围的GPS连续站观测资料进行处理,获得高频GPS动态形变和静态同震水平位移。震中100km范围内四川松潘和甘肃武都站观测到1 Hz动态形变。距离震中约69km的松潘站观测的同震水平位移为7.4mm。根据少量的GPS静态同震位移反演的同震破裂模型显示本次地震的最大滑动量为376mm,地震矩为7.25×1018 N·m,等效矩震级为M_W6.6。正演计算的同震三维形变场显示本次地震的最大水平位移可达4~5cm,垂直位移呈四象限分布,最大可达1.5cm,区域内10个流动GPS站可观测到同震形变。     

近场位移数据约束的2013年芦山地震破裂模型及其构造意义

以往的研究显示了2013年芦山MS7.0级地震发震断层的隐伏逆冲断层基本特征,但是破裂深部细节差异较大.本文以近场密集的同震形变数据约束芦山地震破裂面几何形状及滑动分布,结果显示芦山地震破裂面具有铲状结构,上部16km为43°~50°高角度断层,深部16~25km为小于27°的低角度断层,破裂深度与重定位的余震分布深度一致.破裂分布模型清楚显示上下两个断层上各有一个滑动幅度大于0.5m的峰值破裂区,最大滑动量1.5m位于13km深处.重定位的余震分布基本都落在最大滑动量等值线外部库仑应力增加的区域.芦山地震破裂面几何形状和滑动分布特征与2008年汶川MS8.0级地震映秀—北川破裂相似,支持龙门山冲断带发育大规模的近水平滑脱层,是青藏高原东缘地壳缩短增厚、龙门山挤压隆升的重要证据.     

Coseismic displacement estimate of the 2013 M_S7.0 Lushan,China earthquake based on the simulation of near-fault displacement field

Usually, GPS observation provides direct evidence to estimate coseismic displacement. However, GPS stations are scattered, sparse and cannot provide a detailed distribution of coseismic displacement. Strong ground motion records share the same disadvantages as GPS in estimating coseismic displacement. Estimations from In SAR data can provide displacement distributions; however, the resolution of such methods is limited by the analysis techniques. The paper focuses on estimating the coseismic displacement of the M_S7.0 Lushan earthquake on April 20, 2013 using a simulation of the wave field based on the elastic wave equation instead of a quasi-static equation. First, the media and source models were constructed by comparing the simulated velocity and the record velocity of the ground motion. Then simulated static displacements were compared with GPS records. Their agreement validates our results. Careful analysis of the distribution of simulated coseismic displacements near the fault reveals more details of the ground motion. For example, an uplift appears on the hanging wall of the fault,rotation is associated with the horizontal displacement, the fault strike and earthquake epicenter provide the main control on motion near the faults, and the motion on the hanging wall is stronger than that on the footwall. These results reveal additional characteristics of the ground motion of the Lushan earthquake.     

Coseismic displacement estimate of the 2013 <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>7.0 Lushan,China earthquake based on the simulation of near-fault displacement field

Usually, GPS observation provides direct evidence to estimate coseismic displacement. However, GPS stations are scattered, sparse and cannot provide a detailed distribution of coseismic displacement. Strong ground motion records share the same disadvantages as GPS in estimating coseismic displacement. Estimations from InSAR data can provide displacement distributions; however, the resolution of such methods is limited by the analysis techniques. The paper focuses on estimating the coseismic displacement of the M _S7.0 Lushan earthquake on April 20, 2013 using a simulation of the wave field based on the elastic wave equation instead of a quasi-static equation. First, the media and source models were constructed by comparing the simulated velocity and the record velocity of the ground motion. Then simulated static displacements were compared with GPS records. Their agreement validates our results. Careful analysis of the distribution of simulated coseismic displacements near the fault reveals more details of the ground motion. For example, an uplift appears on the hanging wall of the fault, rotation is associated with the horizontal displacement, the fault strike and earthquake epicenter provide the main control on motion near the faults, and the motion on the hanging wall is stronger than that on the footwall. These results reveal additional characteristics of the ground motion of the Lushan earthquake.     

2014年于田MS7.3级地震近场地表运动和同震位移的初步分析结果

基于2014年于田Ms7.3级地震周边地区的GPS连续观测站数据,得到了此次地震的近场地表运动和同震位移,分析了此次地震的地壳形变特征.结果表明,此次地震引起的地表运动显著,GPS记录到的最大位移振幅远超过最终的同震位移,其中距离震中约60km的XJYT站记录到了西向74.7mm和北向39.9mm的最大位移振幅,同震位移为西向20.0mm和北向9.8mm.2008年于田Ms7.3级地震在一定程度上触发了此次地震.     

芦山“4·20”7.0级地震同震位移场考证及地表变形模式初探

针对2013年四川芦山7.0级地震震源机制与GPS同震位移场不完全一致的问题,利用高精度水平仪对LS05、LS06等6个GPS连续观测临时站观测墩进行了倾斜测量,对SCTQ墩面水准点高差进行了复测,并对该区域地表变形进行了调查。倾斜测量结果表明,LS06等5个站观测墩相对建站时均无明显倾斜。LS05站观测墩标的确相对建站时发生了倾斜,估计其倾斜量在0.4°～1.2°,倾斜方向以NW为主。无论该点倾斜发生于哪个时段,均不影响同震位移以S向运动为主的特征。对SCTQ基准站观测墩面水准测量结果表明,该站点观测墩非常稳定,未发生倾斜。野外地表调查表明,地壳挤压缩短变形特征较为明显,LS05测点周边确有SW向加速运动和剪切变形现象。综合分析认为,LS05、SCTQ和LS06等GPS测点表现出的左旋应变释放特征真实可靠。     

Strong motion displacement waveforms using 10‐Hz precise point positioning GPS: an assessment based on free oscillation experiments

The recently developed precise point positioning (PPP) technique permits to compute instantaneous coordinates of a GPS station relatively to distant reference stations and waveforms of ground displacements during strong motions at the 1 Hz level. This is another application of GPS, different from the computation of static coseismic movements or of accurate monitoring of dynamic displacements of structures using a static receiver and a nearby moving receiver (DGPS). Recently, earthquake ground displacement waveforms using 10‐Hz GPS data have also been calculated, but no independent evidence to assess their quality exists. To overcome this problem, we evaluated the output of 10‐Hz PPP results on the basis of supervised learning experiments. Semistatic and dynamic displacements (damped harmonic oscillations) of known characteristics of the order of a few centimeter were produced and were recorded by GPS, an accelerometer, and a robotic total station. Time series of instantaneous displacements were analyzed using different PPP techniques and were compared with reference (true) values derived from DGPS and the other sensors. Our analysis revealed that the PPP‐derived coordinates are contaminated by long‐period noise but they can display the details of semistatic displacements, while their short‐period component describes well the pattern of waveforms and spectra (at least up to 4 Hz) of dynamic displacements, with up to 20 mm accuracy for isolated points. These results indicate that 10‐Hz PPP‐GPS is useful for earthquake engineering and can safely be used to reconstruct waveforms of deflections of the ground and of various points on structures during strong motions. Copyright © 2014 John Wiley & Sons, Ltd.     

GNSS观测的九寨沟7.0级地震同震位移初步结果

2017年8月8日四川省九寨沟县发生了7.0级地震,中国大陆构造环境网络与北斗地基增强系统的GNSS连续观测共同监测到了此次地震的同震位移（坐标：东向为正,北向为正）,结果显示：3个站点记录到了明显的同震位移,距离震中43 km的九寨沟台站（SCJZ）在东西向的位移为-9.8±1.5 mm,在南北向的位移为3.3±0.7 mm;距离震中65 km的松潘站（SCSP）在东西向的位移为-1.8±0.7 mm,在南北向的位移为-7.7±0.6 mm;距离震中77 km的舟曲站（GSZQ）在东西向的位移为0.4±1.2 mm,在南北向的位移为3.6±0.8 mm.通过同震位移分布特征,可以推测此次地震为一次左旋走滑型事件,引起水平向同震位移大致不超过150 km范围,地震对东南侧的龙门山断裂带影响非常小,对北侧的塔藏断裂和西侧的岷江断裂处引起的同震位移为厘米级.同震位移的反演结果显示：断层面上滑动量主要集中在7 km深度,最大量值约为0.4 m,平均滑动角为-15°,利用滑动分布计算的相应矩震级为M_W6.4,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.     

2016年日本本州东岸近海MS7.2地震前后的地壳运动

2016年11月22日日本本州东岸近海发生东日本M_W9.0大地震的M_S7.2强余震。本文利用美国内华达大学内华达大地测量实验室网站获得此次大地震周围共30多个GPS连续观测站及其它台站的IGS08全球参考框架坐标时间序列,采用北京附近的GPS连续观测站BJSH作为区域位移参考框架的核心站,获取了此次M_S7.2强余震的同震水平位移和区域参考框架位移时间序列,得到了此次强余震前后的位移时空变化图像。结果显示:尽管本州东岸近海M_S7.2地震作为东日本M_W9.0地震的强余震,受其震后形变的影响强烈,但其地壳水平形变的前兆规律与已观测到的大地震一致;不同的是东日本M_W9.0地震前的垂直位移无积累,而本州东岸近海M_S7.2地震前后的垂直位移保持M_W9.0地震后均匀而缓慢的衰减变化;临震前震中附近的GPS连续观测站的东西分量明显减速,有的甚至减至零,是明显的短期前兆异常。此外,本文还进一步讨论了两种同震位移及其成因,并推断地壳水平运动挤压是此次地震的成因。      

芦山7级地震的同震位移估计和震源滑动模型反演尝试

利用自动经验基线校正方法,分析2013年4月20日芦山M_S7.0级地震13个近场强震动台的观测资料,以估算同震位移场分布,并据此反演了震源滑动模型.经与GPS结果比较,两种不同方法给出的芦山7级地震的水平近场同震位移场幅度都不超过cm级,均显示为典型的逆冲型地震（兼有少量左旋走滑错动）.强震最大水平和垂直永久位移分别为4.9 cm和4.4 cm,分别出现在51YAM台和51QLY 台.两种资料反演的震源滑动模型虽显示多事件特征,但主要滑动均集中在第一次事件,即初始滑动点两侧的走向长约30 km、倾向长约25 km的相对集中的较小范围内,强震和GPS模型的最大滑动量分别为1.14 m和1.09 m,较为一致.其余子事件滑动量小且分布零散,不能排除其数值效应的因素.反演矩震级均在M_w6.7左右,地表破裂应该不明显.文章还讨论了目前在我国利用近场强震动记录估计M_w6~7级地震同震位移场存在的困难和问题,为今后类似工作提供参考.     

2011年日本东北地区太平洋近海地震对亚洲板块及韩国大地控制网的影响分析

本研究分析了2011年3月11日发生的M_w9.0日本东北地区太平洋近海地震对亚洲地区和韩国国内GPS卫星常年跟踪站的位移影响.为此,利用了日本东北地区太平洋近海地震发生前后两周(2011年3月4日到3月18日)的GPS站点数据,包括震中附近地区(韩国,中国,中国台湾地区,日本和俄罗斯)55个GPS卫星常年跟踪站和284个IGS 全球跟踪站,并采用GAMIT/GLOBK软件进行处理和平差,估算出所有GPS站点的同震形变.结果显示,日本东北地区太平洋近海地震引起的同震形变影响在亚洲地区比较明显,包括日本和附近国家,距离震中2702 km的中国武汉(WUHN)站也观测到同震形变.为精确分析日本东北地区太平洋近海地震对韩国国家大地控制网的影响,通过GAMIT/GLOBK软件计算出韩国GPS卫星常年跟踪站之间的基线长度变形,并分析出弹性变形量.结果表明:大部分GPS站点均向震中方向膨胀,且向震中的垂直方向收缩.由日本东北地区太平洋近海地震导致的最大剪应变达到韩国国家大地控制网年均变形率的约7倍,对韩国的地壳产生14.5~57.7 mm的水平位移,并导致韩国国家大地控制网产生弹性变形.因此,在不及时更新维护韩国国家大地控制网的情况下,GPS测量成果将会发生最大20 mm的位置误差.     

Modification of fault slip models of the Mw 9.0 Tohoku Earthquake by far field GPS observations

GPS data from Crustal Movement Observation Network of China (CMONOC) are used to derive far-field co-seismic displacements induced by the M_w 9.0 Tohoku Earthquake. Significant horizontal displacements about 30 mm, 10 mm, and 20 mm were caused by this large event in northeast China, north China, and on the Korean peninsula respectively. Vectors of relatively large horizontal displacements with dominant east components pointed to the epicenter of this earthquake. The east components show an exponential decay with the longitude, which is characteristic of the decay of the co-seismic horizontal displacements associated with earthquakes of thrust rupture. The exponential fit of the east components shows that the influence of the co-seismic displacements can be detected by GPS at a distance of about 3200 km from the epicenter of the earthquake. By considering the capability of the far field displacements for constraining the inversion of the fault slip model of the earthquake, we use spherically stratified Earth models to simulate the co-seismic displacements induced by this event. Using computations and comparisons, we discuss the effects of parameters of layered Earth models on the results of dislocation modeling. Comparisons of the modeled and observed displacements show that far field GPS observations are effective for constraining the fault slip model. The far field horizontal displacements observed by GPS are used to modify the slips and seismic moments of fault slip models. The result of this work is applicable as a reference for other researchers to study seismic source rupture and crustal deformation.     

A Possible Detection of the 26 December 2004 Great Sumatra-Andaman Islands Earthquake with Solution Products of the International GNSS Service

The main goal of this work is to critically review the IGS solution products and Precise Point Positioning (PPP) in order to demonstrate their potential to contribute to studies of large earthquakes such as the one that devastated Southeast Asia on December 26th, 2004. In view of a possible detection of the Mw 9.0 Sumatra-Andaman Islands Earthquake of December 26, 2004, position solutions, ranging from intervals of years to one second, of four International GNSS Service (IGS) stations within 3000 km of the epicenter were examined. The IGS combined, cumulative solution product (IGS04P51), consisting of epoch and station velocity solutions and based on data spans of several years prior to the earthquake, was used as a reference. Four IGS combined weekly position solutions (igs04P1301-4), two weeks before and after the earthquake, were utilized for the weekly solution resolution. PPP static and kinematic solutions with IGS Final combined orbits and clocks were used for the mean daily and instantaneous 5-min and 1-sec epoch solutions, respectively. The most significant changes, detected by both weekly and daily solutions occurred in longitude. The nearest IGS station ntus, about 1000 km east of the epicenter, moved westward about 15 mm, while the more distant Indian station iisc (∼ 2300 km NW from the epicenter), shifted about 15 mm eastward. In spite of position errors caused by interpolation of the 5-min IGS clocks, the 1-sec solutions, based on separate data sets, available only for two stations (iisc, dgar), still showed seismic surface waves, in particular at the Indian station iisc. Precise daily IGS combined polar motion and length-of-day products, after correcting for the atmospheric effects, also likely detected, statistically significant, anomalistic excitations on December 26, 2004 that could be caused by this great earthquake.     

四川地震预警网外玛多7.4级地震预警分析研究

2021年5月22日2时4分青海省玛多县发生M7.4地震。四川地震预警台网成功地处理并产出了本次地震预警结果。本文针对四川地震预警系统产出的8次地震预警处理结果及预警减灾效能进行了分析。此次地震震中距四川地震预警网55km,首台触发后12.19s发布首次预警处理结果,前5次预警处理结果用时距离首台触发时间均小于35s。与正式地震目录对比,第1次预警处理结果震中位置偏差为17.1km,震级偏差为-1.8,盲区半径为121km。随着参与计算台站数量增多,震中位置偏差和震级偏差不断减小,盲区半径不断增大。由于地震破裂时间长和台站限幅等原因,8次预警震级均偏小。本次地震预警的有效获益区为震中距121～324km的环形区域,此区域内减灾效能显著。本次地震预警结果表明,四川地震预警系统在本次地震中预警产出效果较好,对于网外大地震预警有重要参考意义。     

地震台站空隙角对地震定位精度的影响

利用四川测震台网记录的芦山余震资料,选取20个M≥2.5地震,以Loc3D(川滇)定位方法定位结果为标准。采用MSDP中常用定位方法(HypoSAT、HYP2000、LocSAT、单纯型),在人为增减台站情况下,使得台站空隙角在25°—320°范围内任意变化,采用相同震相和相同台站对选定的地震事件进行定位,并对定位结果与标准结果对比,结果表明:震中位置在台站空隙角25°—130°任意变化时,4种定位方法效果较好,差别不大;在空隙角达到130°时,震中位移具有随着空隙角增大而逐步变大趋势;在空隙角达到190°时,单纯型法震中位移明显加大,出现随机突跳现象。总体而言,LocSAT法定位效果较好,平均震中位移2.6 km,HYP2000法平均震中位移2.8 km,HypoSAT平均震中位移4.1 km,单纯型定位效果较差,平均震中位移10.3 km。     

印尼8.7级大震前后GPS观测站的地壳水平与垂直位移时间序列结果

利用IGS及中国地壳运动观测网络GPS连续观测站6年多时间的观测结果, 以中国大陆东部的稳定点组作为水平位移与垂直位移解的参考基准, 得到了2004年12月26日印尼苏门答腊8.7级与2005年3月28日8.5级大震前后GPS观测站的地壳水平及垂直位移时间序列的结果。 尽管所采用的GPS站点少, 且分布范围大, 但仍可看到在如此大的地震震前、 同震与震后地壳运动的特征, 为今后观测研究提供了十分有意义的结果。 文中所采用稳定点组基准有效消除水平位移场的平移与旋转, 而局部椭球面的不平性对计算区域地壳运动结果, 特别是对垂直位移的影响很小, 说明了位移解所描述的大范围区域位移场是合理的。 这两次地震的同震水平位移及震后的垂直位移影响量级达数厘米, 范围达远离8.7级地震震中4500 km之外, 甚至更远。     

新疆精河M_W6.3地震产生的静态应力变化研究

利用精河M_W6.3地震有限断层破裂模型,计算了精河地震产生的位移场、应力场、周围主要断层上的静态库仑应力变化以及主震对余震的触发作用。结果表明:(1)精河地震产生的地表隆升最大值约为6.6cm,沉降最大值约为1.8cm;水平位移方向呈现震中南北侧向震中汇聚、震中东西侧向外"流出"的特点。(2)精河地震产生的水平面应力场展布南北侧物质主要受到指向震中的拉张力作用,东西两侧物质主要受到因震中过剩物质东西向排出而导致的东西向挤压力作用。(3)震中西侧距震中约20km的库松木契克山前断裂中段和震中东北部距震中约50km的四棵树-古尔图南断裂西段的库仑应力加载均大于0.01MPa,即2处为地震危险区。(4)在震源深度为8~12km的余震事件中,约有85.5%处于库仑应力加载区,即受到主震的的触发作用;在深度为4~8km的余震事件中,约有87%受到主震的应力触发作用。     

A Rupture Model for the 1999 Chi-Chi Earthquake from Inversion of Teleseismic Data Using the Hybrid Homomorphic Deconvolution Method

This study investigates the kinematics of the rupture process of the M _L 7.3 Chi–Chi, Taiwan, earthquake on September 21, 1999. By applying the proposed hybrid homomorphic deconvolution method to deconvolve teleseismic broadband P-wave displacement recordings of the earthquake, this study derives the apparent source time functions (ASTFs) at ten stations located around the epicenter. To further characterize the fault, the kinematic history of the rupture was inverted from ASTFs using a genetic algorithm, coupled with nonlinear iterative technique. The calculated ASFTs reveal that the total rupture event lasted for approximately 27 s. Static slip distribution images indicate that most slip occurred at shallower portions of the fault plane, especially 20–55 km north of the epicenter. The maximum slip reached 20 m at 45 km north of the epicenter, and the average slip throughout the observed rupture area was approximately 2 m. Large asperities on the fault appeared at 25–35 km and 40–50 km north of the hypocenter, and coincided with relatively high rupture velocity. This suggests that the earthquake’s energy may have been released quickly. The rupture velocity decreased upon encountering an asperity, and increased again after passing the asperity. This implies that the rupture required more time to overcome the resistances of the asperities. The maximum rupture velocity was 3.8 km/s, while the average rupture velocity was approximately 2.2 km/s. The rise time distribution suggests that larger slip amplitudes generally correspond to shorter rise times on the subfaults.