2009年3月19日汤加地震破裂过程快速反演

2009年3月19日汤加MS7.9地震发生后,我们迅速利用地震破裂过程快速反演技术反演了全球地震台网(GSN)的宽频带波形资料,获得了这次地震的破裂过程．随后,收集了更多的资料重新进行了反演,反演结果的主要特征与快速反演结果基本一致．这次地震的震级为MW7.8,破裂持续时间约为70s,断层面上主要有3个滑动量较大的区域．      

2009年9月29日萨摩亚群岛地区地震破裂过程快速反演

运用地震破裂过程快速反演方法,在2009年9月29日萨摩亚群岛地区地震发生后,采用全球地震台网(GSN)的宽频带地震资料,快速反演了这次地震的破裂过程,并于震后3.5小时内得出了这次地震破裂过程的反演结果.结果表明,这次萨摩亚群岛地区地震的破裂过程具有如下基本特征:①矩震级约为MW8.0;②地震主要破裂持续时间约为138s;③滑动量在断层面上的分布比较复杂,整个地震破裂包含至少2个滑动量较大的区域;④这次地震基本上是一次单侧破裂事件,破裂主要朝向西北方向.     

覆盖土层场地地震断裂反应分析方法

基于有限元和拟静力的根本原理，提出了一种基于断层位错的覆盖土层场地地震断裂反应弹塑性有限元分析方法，计算中考虑了土层厚度、断层种类和倾角、软弱夹层以及土的非线性等因素，通过编制的程序计算了地表断裂位移，模拟了土层的破裂过程。     

Site effects in Mexico City eight years after the September 1985 Michoacan earthquakes

The purpose of this paper is to take a comprehensive look at site effects in Mexico City for the 1985 Michoacan earthquake. We examine, successively, 1D and 2D models. For the latter, we consider in detail both large scale and small scale heterogeneities, using extensively the Aki-Larner wave propagation method, in the version given by Bard and Gariel. In particular, we make a critical review of the different explanations proposed for the large duration of strong ground motion in the lake zone. Our purpose is two-sided. We first outline the difference between what is well established and what remains still unexplained regarding the seismic response of Mexico City basin. On the other hand, we wish to make explicit the conditions that the proposed models require to explain strong motion duration. Our results allow us to qualify the models proposed to date and to point out what could be the experiments and the new data required to find a truly satisfactory explanation of strong ground motion at Mexico City.     

Spatio-temporal rupture process of the 2008 great Wenchuan earthquake

Focal mechanism and dynamic rupture process of the Wenchaun M _s8.0 earthquake in Sichuan province on 12 May 2008 were obtained by inverting long period seismic data from the Global Seismic Network (GSN), and characteristics of the co-seismic displacement field near the fault were quantitatively analyzed based on the inverted results to investigate the mechanism causing disaster. A finite fault model with given focal mechanism and vertical components of the long period P-waves from 21 stations with evenly azimuthal coverage were adopted in the inversion. From the inverted results as well as aftershock distribution, the causative fault of the great Wenchuan earthquake was confirmed to be a fault of strike 225°/dip 39°/rake 120°, indicating that the earthquake was mainly a thrust event with right-lateral strike-slip component. The released scalar seismic moment was estimated to be about 9.4×10²⁰-2.0×10²¹ Nm, yielding moment magnitude of M _w7.9–8.1. The great Wenchuan earthquake occurred on a fault more than 300 km long, and had a complicated rupture process of about 90 s duration time. The slip distribution was highly inhomogeneous with the average slip of about 2.4 m. Four slip-patches broke the ground surface. Two of them were underneath the regions of Wenchuan-Yingxiu and Beichuan, respectively, with the first being around the hypocenter (rupture initiation point), where the largest slip was about 7.3 m, and the second being underneath Beichuan and extending to Pingwu, where the largest slip was about 5.6 m. The other two slip-patches had smaller sizes, one having the maximum slip of 1.8 m and lying underneath the north of Kangding, and the other having the maximum slip of 0.7 m and lying underneath the northeast of Qingchuan. Average and maximum stress drops over the whole fault plane were estimated to be 18 MPa and 53 MPa, respectively. In addition, the co-seismic displacement field near the fault was analyzed. The results indicate that the features of the co-seismic displacement field were coincident with those of the intensity distribution in the meizoseismal area, implying that the large-scale, large-amplitude and surface-broken thrust dislocation should be responsible for the serious disaster in the near fault area. Supported by the National Basic Research Program of China (Grant No. 2004CB418404-4) and the National Natural Science Foundation of China (Grant Nos. 40574025 and 40874026)     

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.     

The source rupture feature of the southern Taiwan Straits earthquake of September 16，1994 （Ms7．3） ＆ the analysis of earthquake circumstance in southeastern coast of China

ＴｈｅｓｏｕｒｃｅｒｕｐｔｕｒｅｆｅａｔｕｒｅｏｆｔｈｅｓｏｕｔｈｅｒｎＴａｉｗａｎＳｔｒａｉｔｓｅａｒｔｈｑｕａｋｅｏｆＳｅｐｔｅｍｂｅｒ１６，１９９４（Ｍｓ７．３）＆ｔｈｅａｎａｌｙｓｉｓｏｆｅａｒｔｈｑｕａｋｅｃｉｒｃｕｍｓｔａｎｃｅｉｎｓｏｕ...     

2010年4月14日青海玉树地震破裂过程快速反演

基于远震资料的破裂过程反演方法,通过反演2010年4月14日玉树地震的全球宽频带地震垂直向P波波形记录,先后4次反演得到了玉树地震的破裂过程,并比较和讨论了这4次结果.结果表明,玉树地震的破裂过程具有如下基本特征:①地震主要由两次子事件组成,分别对应于震中附近以及震中东南方向上的两块滑动量集中的破裂区域,其中与第2次子事件对应的震中东南方向上滑动集中的区域破裂贯穿至地表;②最大滑动量和最大滑动速率分别为2.1m和1.1m/s,断层滑动速率较大;③玉树地震总体上是一次单侧破裂事件,破裂从初始破裂点(即仪器测定的震源位置)开始,主要向震中东南方向扩展,由地震多普勒效应导致在震中东南方向上产生强烈能量聚焦,是玉树城区遭受严重破坏在震源方面的主要原因.     

The source rupture feature of the southern Taiwan Straits earthquake of September 16, 1994 ( M S7.3) & the analysis of earthquake circumstance in southeastern coast of China

According to the source mechanism of the main shock and the distribution feature of the aftershocks occurring in the southern Taiwan Straits on Sept. 16, 1994, in this paper the authors analysed the source rupture feature of the major earthquake, demonstrated that this seismic sequence possessed the characteristics of a large intraplate earthquake. And according to the seismotectonic background and the historical seismicity in the area, the authors clarified the active characteristics of the seismically active belts along northwestern direction and analysed preliminarily the earthquake circumstance in the southeastern coast of China.     

Rock magnetic properties and ore microscopy of the iron ore deposit of Las Truchas, Michoacan, Mexico

Iron ore and host rocks have been sampled (90 oriented samples from 19 sites) from the Las Truchas mine, western Mexico. A broad range of magnetic parameters have been studied to characterize the samples: saturation magnetization, Curie temperature, density, susceptibility, remanence intensity, Koenigsberger ratio, and hysteresis parameters. Magnetic properties are controlled by variations in titanomagnetite content, deuteric oxidation, and hydrothermal alteration. Las Truchas deposit formed by contact metasomatism in a Mesozoic volcano-sedimentary sequence intruded by a batholith, and titanomagnetites underwent intermediate degrees of deuteric oxidation. Post-mineralization hydrothermal alteration, evidenced by pyrite, epidote, sericite, and kaolin, seems to be the major event that affected the minerals and magnetic properties. Magnetite grain sizes in iron ores range from 5 to >200 μm, which suggest dominance of multidomain (MD) states. Curie temperatures are 580±5°C, characteristic of magnetite. Hysteresis parameters indicate that most samples have MD magnetite, some samples pseudo-single domain (PSD), and just a few single domain (SD) particles. AF demagnetization and IRM acquisition indicate that NRM and laboratory remanences are carried by MD magnetite in iron ores and PSD–SD magnetite in host rocks. The Koenigsberger ratio falls in a narrow range between 0.1 and 10, indicating the significance of MD and PSD magnetites.     

Surface rupture zone of the 1303 Hongtong M=8 earthquake, Shanxi Province

Introduction The 1303 Shanxi Hongtong M=8 earthquake is the earliest M=8 event determined in histori-cal records in China and the largest recorded in Shanxi fault-depression system in history. Some researchers have discussed the tectonic environment of this earthquake (DENG, et al, 1973; DENG, 1984; DENG, XU, 1994, 1995; Seismo-geological Brigade, State Seismological Bureau, Depart-ment of Geology and Geography, Peking University, 1979; LIU, XIAO, 1982; ZHANG, JIA, 1986; SU, …     

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

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

Fast inversion of rupture process of the 14 April 2010 Yushu,Qinghai, earthquake

Four results of the rupture process of 14 April 2010 Yushu, Qinghai, earthquake, obtained by inverting the broadband seismic data of Global Seismographic Network (GSN) based on the inversion method of earthquake rupture process, were compared and discussed. It is found that the Yushu earthquake has several basic characteristics as follows: ① There exist two principal sub-events which correspond to two slip-concentrated patches being located near the hypocenter and to the southeast of the epicenter. The rupture of the slip-concentrated patch to the southeast of the epicenter broke though the ground surface; ② The peak slip and peak slip-rate are about 2.1 m and 1.1 m/s, respectively, indicating that the Yushu earthquake is an event with large slip-rate on the fault plane; ③ Overall the Yushu earthquake is a unilateral rupture event with the rupture mainly propagating southeastward. The strong focusing of the seismic energy in the southeast of the epicenter due to the "seismic Doppler effect" reasonably accounts for the tremendous damage in the Yushu city.     

反投影全球子台网P波记录研究2010年4月14日玉树地震破裂过程

我们运用反投影远震P波研究了2010年4月14日玉树地震的破裂过程.通过台站响应函数对从IRIS下载的经过波形筛选的206个震中距在30°~95°间的台站进行选取.29个能够得到最好台站响应函数且波形相似度高的台站的记录被用于最终的反投影研究.在研究中我们运用4阶方根叠加法和2 s窗长、1 s滑动的滑动窗对总时长为100 s,始于P波前20 s的波形记录进行反投影.反投影结果显示玉树地震震源破裂时间为20 s,破裂尺度为60 km.在玉树地震的整个破裂过程中能量主要在两个时间和空间上释放.一个能量释放点是震后6 s,位于震中附近;另一个是震后12 s,位于震中位置东南,靠近结古镇的区域,且这次能量释放是玉树地震破裂过程中能量最大的.最大能量释放地点和实际震害最重的区域相吻合.反投影远震P波得到的震源破裂模型是软件自动生成的震动图的有效补充,能为震后应急救援服务.     

Dynamic rupture process of the 1999 Chi-Chi,Taiwan,earthquake

In this study, we preliminarily investigated the dynamic rupture process of the 1999 Chi-Chi, Taiwan, earthquake by using an extended boundary integral equation method, in which the effect of ground surface can be exactly included. Parameters for numerical modeling were carefully assigned based on previous studies. Numerical results indicated that, although many simplifications are assumed, such as the fault plane is planar and all heterogeneities are neglected, distribution of slip is still consistent roug...     

Dynamic rupture process of the 1999 Chi-Chi,Taiwan, earthquake

In this study, we preliminarily investigated the dynamic rupture process of the 1999 Chi-Chi, Taiwan, earthquake by using an extended boundary integral equation method, in which the effect of ground surface can be exactly included. Parameters for numerical modeling were carefully assigned based on previous studies. Numerical results indicated that, although many simplifications are assumed, such as the fault plane is planar and all heterogeneities are neglected, distribution of slip is still consistent roughly with the results of kinematic inversion, implying that for earthquakes in which ruptures run up directly to the ground surface, the dynamic processes are controlled by geometry of the fault to a great extent. By taking the common feature inferred by various kinematic inversion studies as a restriction, we found that the critical slip-weakening distance <i<D</i<<sub<c</sub< should locate in a narrow region 60 cm, 70 cm, and supershear rupture might occur during this earthquake, if the initial shear stress before the mainshock is close to the local shear strength.     

Application of a parameter-shifted grey wolf optimizer for earthquake dynamic rupture inversion

Global optimization is an essential approach to any inversion problem. Recently, the grey wolf optimizer (GWO) has been proposed to optimize the global minimum, which has been quickly used in a variety of inversion problems. In this study, we proposed a parameter-shifted grey wolf optimizer (psGWO) based on the conventional GWO algorithm to obtain the global minimum. Compared with GWO, the novel psGWO can effectively search targets toward objects without being trapped within the local minimum of the zero value. We confirmed the effectiveness of the new method in searching for uniform and random objectives by using mathematical functions released by the Congress on Evolutionary Computation. The psGWO algorithm was validated using up to 10,000 parameters to demonstrate its robustness in a large-scale optimization problem. We successfully applied psGWO in two-dimensional (2D) synthetic earthquake dynamic rupture inversion to obtain the frictional coefficients of the fault and critical slip-weakening distance using a homogeneous model. Furthermore, this algorithm was applied in inversions with heterogeneous distributions of dynamic rupture parameters. This implementation can be efficiently applied in 3D cases and even in actual earthquake inversion and would deepen the understanding of the physics of natural earthquakes in the future.     

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.     

汶川M_W7.9地震起始破裂断层几何结构

本文收集使用紫坪铺水库台网记录到的汶川地震主震P波波形资料,利用P波反投影叠加法获取了2008年5月12日汶川M_W7.9地震起始破裂的时空演化过程.通过分析本次大地震起始破裂阶段(0～1s)破裂点在三维空间内的分布特征,确定了本次大地震起始破裂位置及起始破裂断层几何结构模型.得到以下结果:汶川地震起始破裂点位于31.013±0.002°N、103.392±0.002°E,震源深度为8.2±0.4km,发震时刻为2008年5月12日14∶27∶58.80±0.4.汶川地震起始破裂的最佳断层面走向为NE48°,倾向NW35°,起始阶段破裂的深度范围为地下7.5～9km.