Source time functions of the 1999, Jiji (Chi-Chi) earthquake from GDSN long period waveform data using aftershocks as empirical Green’s functions

A large earthquake (M _W=7.6) occurred in Jiji (Chi-Chi), Taiwan, China on September 20, 1999, and was followed by many moderate-size shocks in the following days. Two of the largest aftershocks with the magnitudes of M _W=6.1 and M _W=6.2, respectively, were used as empirical Green’s functions (EGFs) to obtain the source time functions (STFs) of the main shock from long-period waveform data of the Global Digital Seismograph Network (GDSN) including IRIS, GEOSCOPE and CDSN. For the M _W=6.1 aftershock of September 22, there were 97 pairs of phases clear enough from 78 recordings of 26 stations; for the M _W=6.2 aftershock of September 25, there were 81 pairs of phases clear enough from 72 recordings of 24 stations. For each station, 2 types of STFs were retrieved, which are called P-STF and S-STF due to being from P and S phases, respectively. Totally, 178 STF individuals were obtained for source-process analysis of the main shock. It was noticed that, in general, STFs from most of the stations had similarities except that those in special azimuths looked different or odd due to the mechanism difference between the main shock and the aftershocks; and in detail, the shapes of the STFs varied with azimuth. Both of them reflected the stability and reliability of the retrieved STFs. The comprehensive analysis of those STFs suggested that this event consisted of two sub-events, the total duration time was about 26 s, and on the average, the second event was about 7 s later than the first one, and the moment-rate amplitude of the first event was about 15% larger than that of the second one. Foundation item: State Natural Science Foundation of China (49904004) and IPGP of France. Contribution No. 02FE2007, Institute of Geophysics, China Seismological Bureau.     

Source time functions of the Gonghe,China earthquake retrieved from long-period digital waveform data using empirical Green’s function technique

An earthquake ofM _S=6.9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region,M _S=5.5 on May 7, 1990,M _S=6.0 on Jan. 3, 1994 andM _S=5.7 on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CD-SN) are deconvolved for the source time functions by the correspondent recordings of the three aftershocks as empirical Green’s functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) obtained are nearly identical. The RSTFs suggest theM _S=6.9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about 8 s. Comparing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from P-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-period wavform data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of theM _S=6.0 event on Jan. 3, 1994 and theM _S=5.7 event on Feb. 16, 1994 are quite simple, both RSTFs are single impulses. The RSTFs of theM _S=6.9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment of the three aftershocks deduced from the relative scalar seismic moments of theM _S=6.9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of theM _S=6.9 main shock calculated using the three aftershocks as EGF are 22 (theM _S=6.0 aftershock being EGF), 26 (theM _S=5.7 aftershock being EGF) and 66 (theM _S=5.5 aftershock being EGF), respectively. Deducing from those results, the relative scalar sesimic moments of theM _S=6.0 to theM _S=5.7 events, theM _S=6.0 to theM _S=5.5 events and theM _S=5.7 to theM _S=5.5 events are 1.18, 3.00 and 2.54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1.15, 3.43, and 3.05.     

Strong motion observations and recordings from the great Wenchuan Earthquake

The National Strong Motion Observation Network System （NSMONS） of China is briefly introduced in this paper. The NSMONS consists of permanent free-field stations, special observation arrays, mobile observatories and a network management system. During the Wenchuan Earthquake, over 1,400 components of acceleration records were obtained from 460 permanent free-field stations and three arrays for topographical effect and structural response observation in the network system from the main shock, and over 20,000 components of acceleration records from strong aftershocks occurred before August 1, 2008 were also obtained by permanent free-field stations of the NSMONS and 59 mobile instruments quickly deployed after the main shock. The strong motion recordings from the main shock and strong aftershocks are summarized in this paper. In the ground motion recordings, there are over 560 components with peak ground acceleration （PGA） over 10 Gal, the largest being 957.7 Gal. The largest PGA recorded during the aftershock exceeds 300 Gal.     

A new technique for moment tensor inversion with applications to the 2008 Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>8.0 earthquake sequence

A M_S8.0 earthquake occurred in Wenchuan County, Sichuan Province, China, on May 12, 2008, and subsequently, numerous aftershocks followed. We obtained the moment tensor solutions and source time functions (STFs) for the Wenchuan earthquake and its seven larger aftershocks (M_S5.0～6.0) by a new technique of moment tensor inversion using the broadband and long-period seismic waveform data from the Global Seismic Network (GSN). Firstly, the theoretical background and technical flow of the new technique was briefly introduced, and an aftershock of the Wenchuan earthquake sequence was employed to illustrate the real procedure for inverting the moment tensor; secondly, the moment tensor solutions and STFs of the eight events, including the main shock, were presented, and finally, the interpretation of the results was made. The agreement of our results with the GCMT results indicates the new approach is efficient and feasible. By using this approach, not only the moment tensor solution can be obtained but also the STF can be retrieved; the inverted STFs indicate that the source rupture process may be complicated even for the moderate earthquakes. The inverted focal mechanisms of the Wenchuan earthquake sequence show that the most of the aftershocks occurred in the main faults of the Longmenshan fault zone with predominantly thrustingwith minor right-lateral strike-slip component, but some of them may have occurred in the subfaults with strike-slip faulting in the vicinity of the main faults.     

New Data on Aftershocks of the December 7, 1988, Spitak Earthquake

The paper addresses the collection and analysis of new data on aftershocks that occurred within 20 days of the main shock of the December 7, 1988, Spitak earthquake, Mw = 6.8. The data were used to improve the location of aftershock hypocenters and magnitudes. Available data concerning this 20-day period were the least reliable in terms of completeness, representativeness, and the accuracy of hypocenter location and, in particular, estimation of energy classes and magnitudes. New data were retrieved from the records and bulletins of the seismic stations of the regional and global networks. Hypocenter parameters were determined by means of the minimization of wave travel-time residuals and subsequent double-difference hypocenter relocation. Digital records of the Obninsk and Arti seismic stations (Δ = 15°–18°) and five more distant stations (Δ = 34°–53°) were used to more accurately estimate the surface-wave magnitude of the main shock and strongest aftershock. The aftershock catalog of the Spitak earthquake was substantially revised. First, the previous hypocenter locations (Aref’ev et al., 1991) were improved using the double-difference method; second, new data were retrieved from the bulletins of Caucasian seismic stations. The minimum magnitude of completeness (M _c = 1.9) of the new catalog for the first 20 days after the main shock (when there were no epicentral observations) is the same as that for the period from December 7, 1988, to December 31, 1989. The new catalog contains information on 2090 aftershocks with magnitude M = 1.9 and more for the period from December 7, 1988, to December 31, 1989. The double-difference method allowed the location of the epicenters of clustered earthquakes to be reliably estimated with a longitude error of no more than 4.6 km, a latitude error of 4 km, and a depth error of 5 km. The new spatial distribution of the aftershock hypocenters is better correlated with the tectonic setting than the old data. The new catalog can be used to assess seismic hazard after strong earthquakes in the region.     

张北地区震前、震时和震后的震源机制

对于张北地震之前出现在小震条带上的地震，我们收集了区域地震台网的P波初动和垂直向记录的直达P、S的最大振幅，由理论地震图反演出15次地震的震源机制；同样的方法反演得到145次余震的震源机制；连同刘瑞丰给出的6．2级主震、哈佛大学提供的5．6级强余震的矩张量解一起分析。发现，震前小震震源机制优势取向的平均解、震后余震震源机制优势取向的平均解和主震以及强余震的矩张量解，基本一致。震前、震时和震后的震源机制的优势取向，没有发生可以觉察的显著变化。表明震区始终受到构造应力场的控制，呈现出一定的稳定性。     

A new technique for moment tensor inversion with applications to the 2008 Wenchuan M S8.0 earthquake sequence

A MS8.0 earthquake occurred in Wenchuan County, Sichuan Province, China, on May 12, 2008, and subsequently, numerous aftershocks followed. We obtained the moment tensor solutions and source time functions (STFs) for the Wenchuan earthquake and its seven larger aftershocks (MS5.0~6.0) by a new technique of moment tensor inversion using the broadband and long-period seismic waveform data from the Global Seismic Network (GSN). Firstly, the theoretical background and technical flow of the new technique was briefly introduced, and an aftershock of the Wenchuan earthquake sequence was employed to illustrate the real procedure for inverting the moment tensor; secondly, the moment tensor solutions and STFs of the eight events, including the main shock, were presented, and finally, the interpretation of the results was made. The agreement of our results with the GCMT results indicates the new approach is efficient and feasible. By using this approach, not only the moment tensor solution can be obtained but also the STF can be retrieved; the inverted STFs indicate that the source rupture process may be com-plicated even for the moderate earthquakes. The inverted focal mechanisms of the Wenchuan earthquake sequence show that the most of the aftershocks occurred in the main faults of the Longmenshan fault zone with predomi-nantly thrustingwith minor right-lateral strike-slip component, but some of them may have occurred in the sub-faults with strike-slip faulting in the vicinity of the main faults.     

攀枝花M_S6.1地震余震强震动流动观测

对2008年8月30日攀枝花发生的6.1级地震强余震流动观测记录进行了分析处理,绘制了余震等值线图,讨论r相关问题.通过对P波分析看出:P波显示出明显的方向性,强余震的地震动加速度时程呈现出多普勒效应.为今后强震动观测资料的使用及强震动流动台站的布设,提供了有益的经验.     

宁洱地震序列的震源机制解分析

利用云南数字地震台网资料得到宁洱地震序列的主震、5.1级强余震和52个余震震源机制解.分析表明,该地震序列的发震断裂呈NW走向,倾角陡立.在接近水平的近南北向压应力作用下,断裂具有右旋走滑的错动性质.主震、强余震和众多的余震都发生在同一发震断裂上.大量的余震震源机制解结果与主震一致,是地震序列的主要破裂类型,但还存在与主要类型不一致的倾向滑动类型,这可能与余震破裂起始点的微构造控制作用有关,但是它们呈水平向的应力轴与主震的主应力轴一致.NW向断层作右旋走向错动,滑动断层推挤的象限都是逆冲类型的余震,而拉张的象限都是正断层类型的地震.宁洱地震序列的震源机制和周围4次5级以上地震的震源机制相同,表明震源区应力场和区域应力场完全一致,宁洱地震的孕育和发生受区域应力场的控制.     

2017年九寨沟MS7.0地震震源区横波分裂变化特征

对2017年九寨沟M_S7.0地震序列的横波分裂的时空变化特征进行了分析.通过横波窗内S波质点运动图的分析,从九寨沟地震震源区各个地震台站的近震横波记录中提取了横波分裂的快波偏振方向和慢波延迟时间.观测结果显示,震源区各台站的上地壳各向异性在空间上存在分区特征,时间上有随时间的趋势性变化特征.空间上,位于震源区北部余震区内的3个台站中,发震断层东面的台站L5112和L5111只有一个突出的快波偏振优势方向（NNE向）,而西面的L6202台有两个快波偏振优势方向（除了NNE向,还有一个近EW向）,体现了余震区剧烈调整的地壳应力和构造复杂断裂的综合作用;余震区外的3个台站中,震源区东部靠近塔藏断裂（东）附近的JZG台的快波偏振优势方向为NW向,与塔藏断裂（东）的走向一致,南部的台站L5110和L5113的快波偏振优势方向为近EW向,与区域主压应力方向一致;余震区内各台站的平均慢波延迟时间大于余震区外各台站,反映了九寨沟地震孕育过程中余震区的应力积累强于其周边区域.时间上,快波偏振方向在主震后前期离散度较大,随着时间的推移,离散度在后期有逐渐变小的趋势;慢波延迟时间在主震后较大,但随着时间的推移,也表现出逐渐减小,趋于稳定.横波分裂随时间逐渐减小和趋于稳定的变化特征反映了九寨沟地震在孕震中积累的应力,随着主震和余震的发生而导致的应力释放和调整,应力大小和调整幅度逐渐减小,后期趋于稳定.     

1989—1999大同地震序列的隐伏断层研究:库仑应力分析和余震JHD重定位

1989年到1999年,大同—阳高地区发生了一系列MS≥5的中强地震.本文基于前人对1989年三次MS≥5地震的震源机制反演的结果,通过建立不同断层模型,利用库仑应力方法,计算前震对于主震,以及前震和主震对于余震的库仑应力触发关系,提出了一种可能的破裂模型,即1989年前震沿北西西方向发生左旋破裂,之后主震和余震沿北北东方向发生右旋破裂.根据这种破裂模式计算得出,前震发生后,主震震源处的库仑应力增加了约2×105 Pa,余震震源处的库仑应力出现下降;主震发生后,余震处的库仑应力出现回升,最后余震处的库仑应力几乎没有变化.基于大同地震台网的近场观测数据,用JHD(Joint Hypocentral determination)定位方法,对1999年11月1日MS=5.6地震后一个月的余震进行重定位,得到一条走向118°,倾角85°的左旋走滑断层,余震的深度分布在5km至20km范围内,显示该断层是隐伏断层.另外提出对主震震中位置约10km的修正.本文对1989年三次MS≥5地震序列和1999年MS=5.6地震余震空间分布的研究揭示该地区存在两条活跃的共轭隐伏走滑断层(1989年主震的北北东方向和1999年地震的北西西方向),并且推断已知的大王村断裂和团堡断裂是地下这两条共轭的隐伏走滑断层构造/地震活动在地表的响应.     

1970年通海7.7级大地震强余震触发

1970年1月5日云南通海发生了M_S7.7地震, 震后发生了多次M_S＞5.0的强余震。 文中计算了1970年通海7.7级大地震后, 主震分别在5次强余震破裂面上诱发的库仑破裂应力变化(ΔCFS)。 结果表明, 有4次强余震发生在库仑破裂应力增加(ΔCFS＞0)的地区, 增加的范围为10^-2～10^-1 MPa; 有1次强余震按2种震源机制解结果给出的破裂面计算, 得到2种结果, 分别发生在库仑破裂应力变化为正和在库仑破裂应力变化为负的地区。 研究结果表明, 主震位错产生的库仑破裂应力变化可能是1970年通海7.7级大地震强余震活动的重要原因。     

唐山、海城、盐源-宁蒗地震强余震特征

对大地震发生后强余震的震源机制,基于主要地震断层上小余震的时空分布进行了详细研究.近期在中国境内发生的强烈板内地震,如1976年唐山、1975年海城、1976年盐源——宁蒗等地震曾被中国地震台网很好记录.主震刚发生后的小余震以及强余震前后的小余震的震源用 S——P 时间的收敛法进行测定.强余震的断层是依靠这些小余震的震中分布来描绘的.可以看出有三种情况.第一,在主震断层区内,一次强余震发生在较小余震的比较不密集的部位;第二,一次强余震沿断层发生并与主震断层是共轭的;第三,一次强余震在沿平行并离开主震断层的断层上发生.可以认为,地质条件控制着强余震的发生.作为强余震的可能前兆,则有唐山地震的一次强余震,震级为7.1.其发生前,出现一些前震以及在接近主震断层方向上的某些地震台所记录的初动方向有变化.可以设想,这里的应力场在强余震发生前是有所变化的.      

利用余震震中分析芦山MS7.0 地震发震构造

对2013年4月20日芦山M_S7.0地震主震震中周边29个地震台记录到的震后一年多的微、 小余震,利用Hypo71绝对定位方法进行定位,获得了约1960次地震的震源位置. 结果显示,芦山地震余震在平面上主要沿双石—大川分支断裂及其周边分布,在垂向上主要集中在大约5—20 km深度之间的两条余震交叉带上. 其中一条余震带倾向NW,倾角在12 km左右深度发生变化,浅部倾角较陡,该余震带延伸至地表与双石—大川分支断裂和新开店断裂之间的推测隐伏断裂位置相重合; 另一条余震带倾向SE,其延伸至地表的位置与双石—大川断裂非常接近,但与该断裂倾向相反. 主震震源位置与两条余震带相交的位置接近,且芦山地震主震的两个节面产状与这两条余震带的深部几何形态正好对应, 表明芦山地震主震可能是两条余震带所对应的两条断裂同时活动的结果.      

Variations of shear wave splitting in the 2008 Wenchuan earthquake region

Through the analysis of S-wave particle motion of local events in the shear wave window, the polariza-tion directions of the faster shear wave and the delay times between the faster and the slower shear waves were derived from seismic recordings at the stations near the fault zones. The shear wave split-ting results of seven stations in the area of Longmenshan fault zone reveal spatial variation of the po-larization directions of the fast shear wave. The directions at stations in the southeastern side of the Longmenshan fault zone (in the Sichuan Basin area) are in the NE direction, whereas the direction at station PWU (in the Plateau), which is in the northwestern side of the faults, is in the EW direction. Systematic changes of the time delays between two split shear waves were also observed. At station L5501 in the southern end of the aftershock zone, the delay times of the slower shear wave decrease systematically after the main shock. After the main shock, the delay times at station PWU were longer than those before the earthquake. Seismic shear wave splitting is caused mostly by stress-aligned microcracks in the rock below the stations. The results demonstrate changes of local stress field dur-ing the main-shock and the aftershocks. The stress in the southern part of Wenchuan seismogenic zone was released by the main-shock and the aftershocks. The crustal stresses were transferred to the northeastern part of the zone, resulting in stress increase at station PWU after the main-shock.     

The Al Hoceima Mw 6.4 earthquake of 24 February 2004 and its aftershocks sequence

The Al Hoceima Mw 6.4 earthquake of 24 February 2004 that occurred in the eastern Rif region of Morocco already hit by a large event in May 1994 (Mw 5.9) has been followed by numerous aftershocks in the months following the event. The aftershock sequence has been monitored by a temporary network of 17 autonomous seismic stations during 15 days (28 March–10 April) in addition to 5 permanent stations of the Moroccan seismic network (CNRST, SPG, Rabat). This network allowed locating accurately about 650 aftershocks that are aligned in two directions, about N10-20E and N110-120E, in rough agreement with the two nodal planes of the focal mechanism (Harvard). The aftershock alignments are long enough, about 20 km or more, to correspond both to the main rupture plane. To further constrain the source of the earthquake main shock and aftershocks (mb > 3.5) have been relocated thanks to regional seismic data from Morocco and Spain. While the main shock is located at the intersection of the aftershock clouds, most of the aftershocks are aligned along the N10-20E direction. This direction together with normal sinistral slip implied by the focal mechanism is similar with the direction and mechanisms of active faults in the region, particularly the N10E Trougout oblique normal fault. Indeed, the Al Hoceima region is dominated by an approximate ENE-SSW direction of extension, with oblique normal faults. Three major 10–30 km-long faults, oriented NNE-SSW to NW-SE are particularly clear in the morphology, the Ajdir and Trougout faults, west and east of the Al Hoceima basin, respectively, and the NS Rouadi fault 20 km to the west. These faults show clear evidence of recent vertical displacements during the late Quaternary such as uplifted alluvial terraces along Oued Rihs, offset fan surfaces by the Rouadi fault and also uplifted and tilted abandoned marine terraces on both sides of the Al Hoceima bay.However, the N20E direction is in contrast with seismic sources identified from geodetic inversions, which favour but not exclusively the N110-120E rupture directions, suggesting that the 1994 and 2004 events occurred on conjugate faults. In any event, the recent seismicity is thus concentrated on sinistral N10-20E or N110-120E dextral strike-slip faults, which surface expressions remain hidden below the 3–5 km-thick Rif nappes, as shown by the tomographic images build from the aftershock sequence and the concentration of the seismicity below 3 km. These observations may suggest that strain decoupling between the thrusted cover and the underlying bedrock and highlights the difficulty to determine the source properties of moderate events with blind faults even in the case of good quality recorded data.     

2013年4月20日四川芦山M7.0级地震与余震精确定位及发震构造初探

使用汇集在四川台网中心的固定台站、震后架设的流动台站、周边水库台站等震中距150 km以内的震相数据,选用分层速度模型,对芦山7.0级地震及震后9天内的余震利用双差定位法进行了重新定位.给出了芦山7.0级地震的发震时刻为2013-04-20 08:02:46.8,震中位置30.278°N,102.989°E,震源深度16.67 km,给出了3324次余震的双差定位结果,并对发震构造进行了分析.结果表明:芦山地震主破裂长度约40 km,下倾宽度约20 km,破裂视面积约800 km²,主破裂沿南西走向,倾角约40°.余震震源优势深度为10~22 km.余震沿南西走向,主要集中于大邑-名山断裂上盘.     

Stochastic procedure for the simulation of synthetic main shock‐aftershock ground motion sequences

According to the current seismic codes, structures are designed to resist the first damaging earthquake during their service life. However, after a strong main shock, a structure may still face damaging aftershocks. The main shock‐aftershock sequence may result in major damage and eventually the collapse of a structure. Current studies on seismic hazard mainly focus on the modeling and simulation of main shocks. This paper proposes a 3‐step procedure to generate main shock‐aftershock sequences of pairs of horizontal components of a ground motion at a site of interest. The first step generates ground motions for the main shock using either a source‐based or site‐based model. The second step generates sequences of aftershocks' magnitudes, locations, and times of occurrence using either a fault‐based or seismicity‐based model. The third step simulates pairs of ground motion components using a new empirical model proposed in this paper. We develop prediction equations for the controlling parameters of a ground motion model, where the predictors are the site condition and the aftershock characteristics from the second step. The coefficients in the prediction equations and the correlation between the model parameters (of the 2 horizontal components of 1 record and of several records in 1 sequence) are estimated using a database of aftershock accelerograms. A backward stepwise deletion method is used to simplify the initial candidate prediction equations and avoid overfitting the data. The procedure, based on easily identifiable engineering parameters, is a useful tool to incorporate effects of aftershocks into seismic analysis and design.     

The May 13, 1995, Kozani-Grevena (NW Greece) earthquake: Source study and its tectonic implications

At 08:47 GMT, on May 13, 1995, a strong earthquake of M_s = 6.6 occurred in the NW part of Greece (Western Macedonia) and caused serious damage in the Kozani and Grevena prefectures, but fortunately no fatalities. The maximum observed macroseismic intensity was IX + of the Modified Mercalli scale. The main shock was preceded by several foreshocks and followed by intense aftershock activity lasting several months.The Institute of Geodynamics of the National Observatory of Athens, in order to monitor and study the aftershock activity, installed a seismic network of nine (9) stations operated for a period of 50 days. Thousands of aftershocks were recorded. Based on the analysis of recorded data, a NE-SW trending zone dipping NW is defined.In the field a surface rupture of normal slip was observed, following a NE-SW direction for a length of 8 km with a 4 cm down throw of the NW area. This break was located along a pre-existing minor normal fault, while a main fault system exists 10 km to the SE.The focal mechanism of the main shock shows normal faulting, which is in agreement with the field observations. Moreover focal mechanisms of several well defined aftershocks were computed, showing various types of faulting.