用连续GPS与远震体波联合反演2015年尼泊尔中部M_S8.1地震破裂过程

由于印度-欧亚板块碰撞,位于板块边界带的喜马拉雅地区大震频繁,但对其活动性的认识仍十分有限.2015年4月25日尼泊尔中东部地区时隔80年再次发生8级地震,为研究板缘地震提供了一次难得机遇.本文用西藏和尼泊尔的GPS连续观测数据和全球分布的远震地震波记录联合反演此次特大地震的破裂过程,结果显示此次地震发生在印度板块与青藏高原接触边界面——喜马拉雅主滑脱断层上.北倾11°、近东西(295°)走向的断层面破裂约100km长(博卡拉到加德满都),130km宽(从加德满都深入我国西藏吉隆县),破裂以逆冲滑动为主,平均幅度达到2.4m,释放的地震矩高达9.4×1020 N·m.反演结果还显示,震源体主要破裂分布深度范围为5~25km,应无地表破裂,属于一次盲地震.基于GPS资料推测的地壳现今运动速率及1833年地震的震源位置,我们推测地震在此次地震破裂区域复发的周期可能为150~200a,而极震区以南的深部滑脱断层仍保持闭锁,未来仍有导致灾害性大震的可能性.     

远震体波有限破裂过程反演敏感性测试及其在汶川地震上的应用

本文利用远震P和SH波反演得到2008年5月12日汶川大地震(M_W=7.9)的一系列有限破裂模型。使用一种基于小波变换的模拟退火非线性反演方法, 我们将主断层划分成若干个子断层, 在反演时同时确定每个子断层上的滑移量、 滑动角、 上升时间(rise time)以及平均破裂速度。我们首先根据一个假定的破裂模型生成理论地震图, 将该理论地震数据作为输入进行反演, 对该有限破裂反演方法进行了一系列测试, 以验证反演对断层倾角、 平均破裂速度、 最大破裂深度等参数的敏感性。然后我们采用4个不同倾角的断层面来对汶川地震远震体波记录进行反演。结果表明, 若对只在一个断层面上模拟该地震, 30°倾角是个较为合适的值。反演的结果还表明, 此次地震有两个主要的能量释放区域, 并且主断层面存在倾角变化的可能性。在将来的研究中, 可以结合GPS, InSAR测地学以及强震等数据, 来对强震的破裂过程做更细致的研究。     

GEODETIC AND TELESEISMIC CONSTRAINTS ON SLIP DISTRIBUTION OF 2015 MW6.4 PISHAN EARTHQUAKE

On July 3^rd, 2015, a M_W6.4 earthquake occurred on Pishan County, Xinjiang, located in the front of western Kunlun thrust belt, which is the largest earthquake(M_W6.0~7.0)in the past 40 years in this region. In this study, we collected both the near-filed geodetic coseismic deformation observations including 4 GPS sites and one high-resolution ALOS-2 InSAR imagery, and far-field teleseismic P waveforms from 25 stations provided by IRIS/USGS, to invert the fault parameters(strike and dip)and coseismic rupture model of 2015 M_W6.4 Pishan earthquake. Using the finite fault theory, a non-linear simulated annealing algorithm was employed to resolve our joint inversion problem. The strike (120°~130°) and dip angle(35°~40°)of optimal models are different from that of some previous studies, and the dip change is strongly constrained by combined data than that of strike. In fixing the geometric parameters of optimal fault model, we also considered data weight(5)(geodetic data/teleseismic P waveforms)and constrained weight from moment and smooth factor(2.5). Clearly, our results indicate that the slip distribution mainly concentrates in the depth range from 9 to 16km and a length range of 20km along the strike direction, which is similar to the spatial distribution of the relocated aftershocks. The maximum slip is~95cm. The seismic moment release is 5.45×10¹⁸N·m, corresponding to M_W6.42. Compared with the single data set, geodetic data or teleseismic waveform, our joint inversion model could simultaneously constrain the seismic moment and slip distribution well, thus avoiding effectively a lower-resolution rupture distribution determined by teleseismic-only inversion and a bias released moment estimated by the geodetic-only inversion. Importantly, we should consider both the near-field geodetic data and far-field teleseismic data in retrieving the rupture model for accurately describing the seismogenic structure of active fault in western Kunlun region.     

Source process of the 14 November 2001 western Kunlun Mountain M_S=8.1 earth-quake

Introduction On November 14, 2001, a great earthquake occurred in the western Kunlun Mountain area(Figure 1). The original time is 09h26min10.0s (UTC); the hypocentral location is 35.95°N,90.54°E; focal depth is 10 km from USGS National Earthquake Information Center (NEIC);MS=8.1 from China Seismic Network and Mw=7.8 from Harvard and Earthquake Research Insti-tute (ERI), University of Tokyo. This earthquake, known as the western Kunlun Mountain earth-quake, is an extraord…     

A detailed source study of the Orta (Çankırı) earthquake of June 6, 2000 (MS = 6.1): An intraplate earthquake in central Anatolia

The devastating zmit and Düzce earthquakes were followed by theOrta intra-plate earthquake (M_S = 6.1) occurred in the central Anatolianblock on June 6, 2000. The focal mechanism, aftershock distribution andthe field studies (Emre et al., 2000) suggest a movement on a 21-km longDodurga fault striking nearly N-S where the sense of motion is left-lateralstrike-slip with considerable amount of normal component. We applied theconstrained linear finite-fault inversion method of Hartzell and Heaton(1983) to the teleseismic P and SH waveforms to derive a coseismic slipdistribution model for the earthquake. Time windows approach is appliedallowing variable rise times and rupture velocities. The source-rise timefunction is discretized into consecutive time intervals that stand for slipcontribution of individual subfaults. Although no clear surface ruptureswere associated with the earthquake, the resulting slip model suggestscoseismic slip in the order of several tens of centimetres. Our coseismicslip distribution model identifies two slip patches with the followingmaximum slip values: (1) the larger one (42 cm) is located to the southof the hypocenter at depth range of 4–8 km and (2) the smaller one(31 cm) is located just above and north of the hypocenter. Theslip-model yield a seismic moment of 1.0 × 10¹⁸ Nm, most of whichis released from the rupture over the depth of 8 km.     

基于InSAR和远场地震波联合反演2008年MW6.3大柴旦地震震源破裂过程

利用InSAR同震形变升、降轨数据和远场地震波数据,基于均方根最小与标量地震矩最小双重约束下的模拟退火方法,联合反演2008年11月10日大柴旦M_W6.3地震震源破裂过程.结果表明,2008年大柴旦地震震源破裂过程整体表现为沿倾向方向从深部破裂起始点处向上往地表传播,且破裂未到达地表;在前7 s内,滑动沿西北和东南两个方向传播,7 s后主要沿东南方向传播,破裂过程时间持续约为11 s,同震滑动分布主要集中在地下10~20 km范围内,最大滑动量达-0.71 m;反演结果揭示本次地震为西南倾高角度兼具微量走滑分量的逆冲破裂事件,反演矩张量为3.96×10¹⁸N·m,矩震级约M_W6.37.通过大柴旦地震发震断层和破裂机制综合分析,初步判断发震断层为西南倾向的大柴旦—宗务隆山断裂.     

2021年日本福岛县冲地震的震源破裂过程分析基于采用经验格林函数方法的波形反演

利用中小震作为经验格林函数,选取0.2~2.0 Hz频段的强震数据进行波形反演,获得了2021年福岛县冲地震的破裂过程。结果表明:该地震的破裂主要集中在断层面距离震源约25 km的区域内,沿震源向东北延伸约5 km,向西南延伸约20 km;在该区域内,识别出两个滑动量集中的区域,均分布在距离震源约15 km西南侧,主要滑动量集中区域最大滑动量约3.2 m,深度基本与震源一致;次要滑动量集中区域略比主要滑动量集中区域浅约18 km。该震源模型对应矩震级M_w7.3,破裂速度为2.4 km/s。通过选择不同的中小震组合进行波形反演,结果对该震源模型特性无显著影响,揭示了该震源模型的稳健性;基于该震源模型合成反演分析中未使用的强震观测台站强震动,获得的合成波形与观测波形有很好的相关性,充分证明了该震源模型时空特征的可靠性。     

利用近场高频GPS、 强地面运动和远场地震波形数据联合反演2008年汶川MS8.0地震的震源时空破裂过程

联合近场GPS测站1-Hz运动学位移、 强震仪加速度波形和全球台站P震相波形作为约束, 以时空滑动分布约束条件和ABIC模型参数选择方法, 结合先验的滑动方向变化范围, 反演2008年汶川M_S8.0地震的震源时空破裂过程, 给出了能够综合反映震源破裂过程的统一模型。 结果表明, 汶川地震总体上存在4个主要的破裂区, 最主要的一个破裂区位于震源东北40~120 km, 断层面上的最大位错量约为10 m, 主体滑动分布在2~20 km深度范围, 破裂达到地表; 第二个主体破裂区位于断层破裂带南段, 最大滑动量达到6 m; 另外2个主体滑动区位于断层破裂带北段, 但滑动破裂量小于断层南段破裂区的滑动量, 滑动破裂值最大值为4 m, 超过1 m的区域在走向上超过70 km。 反演得到的断层滑动模型的地震矩为9.5×10²¹ Nm, 相应的矩震级为M_W7.95。 汶川地震破裂表现为单侧破裂, 起始破裂在汶川下方16 km深度, 向东北方向一致性地传播, 过程持续~120 s。 在地震发生后0~10 s内, 破裂集中在震源起始破裂区, 滑动破裂值为~1.0 m, 之后破裂向东北方向扩展, 震后20~40 s是主要的破裂时段。 在40~60 s, 破裂跨越断层南段和北段。 在80~90 s破裂最大值开始下降, 在100~110 s时, 下降为~0.5 m, 在110~120 s时, 下降为~0.1 m。 加入近场GPS测站1-Hz 波形数据与近场强震仪波形和远场长周期体波联合反演, 提高了震源破裂模型的空间分辨率, 特别是浅部滑动破裂区的分辨率, 反演的最大滑动破裂值比不用1-Hz 波形数据反演的结果增大, 表明近场1-Hz GPS波形数据对于揭示汶川地震的时空破裂过程具有重要的作用。     

Moment tensor inversion and source rupture process of the September 27, 2003 M_S=7.9 earthquake occurred in the border area of China, Russia and Mongolia

IntroductionOn September 27, 2003, an earthquake of MS=7.9 struck the border area of China, Russia and Mongolia. According to the field investigation from the Earthquake Administration of XinjiangAutonomous Region, the whole northern Tianshan region felt the hit. Buildings and structures within six counties and one city in Altay region, which is total about 0.11×106 km2 area, were damaged to different extent and caused certain economic losses. The epicenter determined by China National …     

Source process of the 14 November 2001 western Kunlun Mountain M S=8.1 earthquake

Based on digital teleseismic P-wave seismograms recorded by 28 long-period seismograph stations of the global seismic network, source process of the November 14, 2001 western Kunlun Mountain M _S=8.1 (M _W=7.8) earthquake is estimated by a new inversion method. The result shows that the earthquake is a very complex rupture event. The source rupture initiated at the hypocenter (35.95°N, 90.54°E, focal depth 10 km, by USGS NEIC), and propagated to the west at first. Then, in several minutes to a hundred minutes and over a large spatial range, several rupture growth points emerged in succession at the eastern end and in the central part of the finite fault. And then the source rupture propagated from these rupture growth points successively and, finally, stopped in the area within 50 km to the east of the centroid position (35.80°N, 92.91°E, focal depth 15 km, by Harvard CMT). The entire rupture lasted for 142 s, and the source process could be roughly separated into three stages: The first stage started at the 0 s and ended at the 52 s, lasting for 52 s and releasing approximately 24.4% of the total moment; The second stage started at the 55 s and ended at the 113 s, lasting for 58 s and releasing approximately 56.5% of the total moment; The third stage started at the 122 s and ended at the 142 s, lasting for 20 s and releasing approximately 19.1% of the total moment. The length of the ruptured fault plane is about 490 km. The maximum width of the ruptured fault plane is about 45 km. The rupture mainly occurred within 30 km in depth under the surface of the Earth. The average static slip in the underground rocky crust is about 1.2 m with the maximum static slip 3.6 m. The average static stress drop is about 5 MPa with the maximum static stress drop 18 MPa. The maximum static slip and the maximum stress drop occurred in an area within 50 km to the east of the centroid position. Foundation item: Joint Seismological Science Foundation of China (103066) and Foundation of the Seismic Pattern and Digital Seismic Data Application Research Office of Institute of Earthquake Science of the China Earthquake Administration.     

2010年2月27日智利地震破裂过程快速反演

运用地震破裂过程快速反演方法,在2010年2月27日智利地震发生后,采用全球地震台网(GSN)的宽频带地震资料,反演了这次地震的破裂过程,在震后约3.5小时得到了这次地震破裂过程的反演结果.结果表明,这次智利地震的破裂过程具有如下基本特征:①矩震级为MW8.6;②地震破裂持续时间约为150s;③包括4个滑动量集中分布区域,最大滑动量约为8m,最大滑动速率约为0.8m/s;④这次地震总体上是一次不对称的双侧破裂事件,破裂从破裂起始点(震源)开始,同时向南北两个方向扩展,但以向北扩展的破裂为主.     

THE SOURCE PARAMETERS AND RUPTURE PROCESS OF THE MW7.0 EARTHQUAKE IN ALASKA,USA ON DECEMBER 1, 2018

A magnitude M_W7.0 earthquake struck north of Anchorage, Alaska, USA on 1 December 2018. This earthquake occurred in the Alaska-Aleutian subduction zone, on a fault within the subducting Pacific slab rather than on the shallower boundary between the Pacific and North American plates. In order to better understand the earthquake source characteristics and slip distribution of source rupture process as well as to explore the effect of tectonic environment on dynamic triggering of earthquake, the faulting geometry, slip distribution, seismic moment, source time function are estimated from broadband waveforms downloaded from IRIS Data Management Center. We use the regional broadband waveforms to infer the source parameters with ISOLA package and the teleseismic body wave recorded by stations of the Global Seismic Network is employed to conduct slip distribution inversion with iterative deconvolution method. The focal mechanism solution indicates that the Alaska earthquake occurred as the result of tensile-type normal faulting, the estimated centroid depth from waveform inversion shows that the earthquake occurred at the depth of 56.5km, and the centroid location is 10km far away in northeast direction relative to the location of initial epicenter. We use the aftershock distribution to constrain the fault-plane strike of a normal fault to set up the finite fault model, the finite fault inversion shows that the earthquake slip distribution is concentrated mainly on a rectangular area with 30km×20km, and the maximum slip is up to 3.6m. In addition, the slip distribution shows an asymmetrical distribution and the range of possible rupture direction, the direction of rupture extends to the northeast direction, which is same as that of aftershock distribution for a period of ten days after the mainshock. It is interesting to note that a seismic gap appears in the southwest of the seismogenic fault, we initially determined that the earthquake was a typical normal fault-type earthquake that occurred in the back-arc extensional environment of the subduction collision zone between the Pacific plate and the North American plate, this earthquake was not related to tectonic movement of faults near the Earth's surface. Due to the influence of high temperature and pressure during the subduction of the Pacific plate toward to the north, the subduction angle of the Pacific plate becomes steep, causing consequently the backward bending deformation, thus forming to a tensile environment at the trailing edge of the collision zone and generating the M_W7.0 earthquake in Alaska.     

四川汶川8.0级地震震源过程

2008年5月12日在青藏高原与四川盆地交界的龙门山山脉发生了M_s8.0级强烈地震,引发山体滑坡等地质灾害,造成了巨大的人员伤亡和经济损失.本文利用远场体波波形记录结合近场同震位移数据,根据地质资料和地震形成的地表破裂轨迹,构造了一个双“铲状”有限地震断层模型,利用反演技术重建地震的破裂过程.结果显示汶川大地震主要是沿龙门山构造带的映秀-北川断裂和灌县-江油断裂发生的逆冲兼右旋走滑破裂事件.断层面上的滑动分布显示两个高滑动区先后发生在地震破坏最为严重的映秀和北川地区,最大滑动量高达1200～1250 cm,且破裂过程也显示一定的复杂性.地震破裂的平均走滑量略大于平均倾滑量,与多种观测资料获得的震前龙门山断裂带构造变形相一致,推断是由于长期区域应力场作用和龙门山地区特殊的物质组成和结构孕育了这次千年尺度的强烈地震.     

2009年10月7日瓦努阿图地震破裂过程快速反演

运用地震破裂过程快速反演方法,在2009年10月7日瓦努阿图地震发生后,采用全球地震台网(GSN)的宽频带地震资料,快速反演了这次地震的破裂过程,在震后3小时内得出了这次地震破裂过程的反演结果.结果表明,这次瓦努阿图地震的破裂过程具有如下基本特征:①矩震级约为MW7.8;②地震主要破裂持续时间约为80s;③滑动量在断层面上的分布比较简单,整个地震破裂只包含一个滑动量较大的区域;④这次地震的破裂模式接近圆盘破裂,其中心位于震源处.     

Rupture process of the 2011 Tohoku earthquake from the joint inversion of teleseismic and GPS data

Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s along the fault strike direction. The source rupture process consists of three sub-events, the first occurred near the hypocenter and the rest two ruptured along the up-dip direction and broke the sea bed, causing a maximum slip of about 30 m. The large-scale sea bed breakage may account for the tremendous tsunami disaster which resulted in most of the death and missing in this mega earthquake.     

芦山4.20地震破裂过程及其致灾特征初步分析

采用全球地震台网的远震地震波数据快速反演了2013年4月20日雅安芦山地震(芦山4.20地震)的破裂过程,在震后3小时得到并发布了相关测定结果.结果显示这次地震的震级约M_W6.8,包括两次破裂子事件,但都发生在断层面上震源(破裂起始点)附近,整个地震没有表现出明显的破裂方向性.雅安芦山地震的破裂没有大规模出露地表,主要地震灾害集中在断层上盘,位于芦山县和宝兴县一带.     

Inversion of near-field waveform data for earthquake source rupture process ( II ): Inversion of Chi Chi earthquake source, 20 September, 1999, Taiwan

The new inversion algorithm developed based on the recent progress in the nonlinear programming study by us is used to invert the earthquake source process of Chi Chi earthquake M _w7.6, 20 Semptember, 1999, Taiwan. A curve fault model is constructed in our inversion to make the fault model close to the real rupturing fault to reduce the influence from the discrepancy between the constructed fault model and the real rupturing fault. The results show that (1) the rupture process of the Chi Chi earthquake source lasted about 32 seconds and the main faulting occurred between 6th to 21st second after the start of the ruptures and the high slip area were mainly located at the northern segment of the fault. (2) The slip was dominated by thrust faulting. The average rake angle was 64.5°, which was very consistent with those inverted by USGS, Harvard and CWB (Central Weather Bureau of Taiwan). The amount of the moment inverted in this paper was 7.76×10²⁰ NM, which was a slightly bigger than those inverted by USGS and Harvard. (3) A clear nucleation step existed in the source faulting process and it lasted about 6 seconds. The moment release rate accelerated obviously at the end of the nucleation step. The faulting started from the southern segment and mainly occurred at the northern segment after 10 seconds. At the end of this paper, we analyzed the reliability of the inversion result via comparing with the GPS observations and discussed its scientific signification.     

Moment tensor inversion and source rupture process of the September 27, 2003 M S=7.9 earthquake occurred in the border area of China, Russia and Mongolia

We conducted moment tensor inversion and studied source rupture process for M _S=7.9 earthquake occurred in the border area of China, Russia and Mongolia on September 27 2003, by using digital teleseismic P-wave seismograms recorded by long-period seismograph stations of the global seismic network. Considering the aftershock distribution and the tectonic settings around the epicentral area, we propose that the M _S=7.9 earthquake occurred on a fault plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of M _S=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tensor is M ₀=0.97×10²⁰ N·m. Rupture mainly occurred on the shallow area with 110 km long and 30 km wide, the location in which the rupture initiated is not where the main rupture took place, and the area with slip greater than 0.5 m basically lies within 35 km deep middle-crust under the earth surface. The maximum static slip is 3.6 m. There are two distinct areas with slip larger than 2.0 m. We noticed that when the rupture propagated towards northwest and closed to the area around the M _S=7.3 hypocenter, the slip decreased rapidly, which may indicate that the rupture process was stopped by barriers. The consistence of spatial distribution of slip on the fault plane with the distribution of aftershocks also supports that the rupture is a heterogeneous process owing to the presence of barriers.     

Inversion of near-field waveform data for earthquake source rupture process (Ⅱ): Inversion of Chi Chiearthquake source, 20 September, 1999, Taiwan

The new inversion algorithm developed based on the recent progress in the nonlinear programming study by us is used to invert the earthquake source process of Chi Chi earthquake Mw7.6, 20 Semptember,1999, Taiwan. A curve fault model is constructed in our inversion to make the fault model close to the real rupturing fault to reduce the influence from the discrepancy between the constructed fault model and the real rupturing fault. The results show that (1) the rupture process of the Chi Chi earthquake source lasted about 32 seconds and the main faulting occurred between 6th to 21st second after the start of the ruptures and the high slip area were mainly located at the northern segment of the fault. (2) The slip was dominated by thrust faulting. The average rake angle was 64.5°, which was very consistent with those inverted by USGS, Harvard and CWB (Central Weather Bureau of Taiwan). The amount of the moment inverted in this paper was 7.76×1020 NM, which was a slightly bigger than those inverted by USGS and Harvard. (3) A clear nucleation step existed in the source faulting process and it lasted about 6 seconds. The moment release rate accelerated obviously at the end of the nucleation step. The faulting started from the southern segment and mainly occurred at the northern segment after 10 seconds. At the end of this paper, we analyzed the reliability of the inversion result via comparing with the GPS observations and discussed its scientific signification.     

2014年云南鲁甸M_W6.1地震:一次共轭破裂地震

烈度与余震分布显示2014年云南鲁甸MW6.1(MS6.5)地震的发震构造较复杂.为深入了解鲁甸地震的发震断层与破裂特征,本文考虑了单一断层破裂和共轭断层破裂的情况,对震中距250km范围内的近震资料(宽频带资料和强震资料)和远震体波资料进行了反演,得到了鲁甸地震的破裂过程,探讨了滑动分布与余震分布之间的关系.根据反演得到的滑动分布、震源时间函数和波形拟合,认为鲁甸地震是一次在北西向主压应力与北东向主张应力的统一应力场下发生的两条共轭断层先后破裂的一次复杂地震事件.在破裂开始后0~2s,破裂主要发生在ENE—WSW向(近东西向)的断层上,随后NNW—SSE向(近南北向)断层开始破裂,释放了大部分的地震矩.由于近南北向断层南段(即震中以南)的破裂规模较大,且以左旋走滑为主,对近东西向断层的西段起到了一定程度的解锁作用,可能是震中以西无明显主震破裂但存在密集余震分布的主要原因.