由Hi—Net台阵成像得出的2004年苏门答腊-安达曼地震的范围、持续时间和速度

2004年损失惨重的苏门答腊-安达曼地震是已经记录到的最大地震之一。这种地震的破坏程度取决于断层滑动的范围和大小。在苏门答腊-安达曼地震后几小时得到的最初可靠的矩震级估计（Dziewonski and Woodhouse，1983）值为9．0级，但是，最近的长周期简正模式分析表明，     

The aftershock dynamics of the Sumatra-Andaman earthquake

The aftershocks of the catastrophic Sumatra-Andaman earthquake of December 26, 2004 (M = 9.0) are analyzed in the general context of the theory of critical phenomena. The analysis relies on the idea that, according to this theory, critical transitions have two key properties. The first is that the intensity of the fluctuations in a dynamical system monotonically increases with the approach of the bifurcation point, so that at a certain time instant, a sufficiently strong internal pulse initiates the catastrophe. This transition can be treated as spontaneous. The second property is that the reactance of the dynamical system drastically increases on the approach of the bifurcation. Even a weak external perturbation in the near-threshold interval can result in a catastrophe. In this case, it is reasonable to refer to the critical transition as an induced transition. The aftershocks of the Sumatra-Andaman earthquake are likely to demonstrate the typical features of induced seismicity. First, the strongest aftershock (M = 7.2) occurred 3 h 20 min after the main shock. It could have probably been induced by the round-trip seismic echo. Second, it was found that the spectral density of the aftershock sequence significantly increases at about ～0.3 mHz, which is close to the frequency of the spheroidal mode ₀S₂. This suggests that the spheroidal oscillations of the Earth, which are excited by the main seismic shock, modulate the aftershock activity. Both hypotheses are supported by the analysis of the aftershocks of the Tohoku earthquake of March 11, 2011 (M = 9.0).     

用10分钟P波窗估算2004年12月26日苏门答腊—安达曼地震的辐射能量

苏门答腊—安达曼MW9.1级地震的破裂过程持续了大约500秒,这一时间几乎是一般情况下用于计算远震辐射能量的P到时和PP到时之间时间窗长度的两倍。为了测量整个地震所辐射的P波,我们将时间窗扩展为从P波到时到S波到时,并用扩展窗对震中距大于60°的台站的地震记录进行分析。这些持续时间8～10分钟的窗内包含了PP,PPP,ScP震相和其他一些多次反射震相。为了测量包含这些附加震相的影响,我们计算了由扩展窗(P波到时和S波到时之间)得到的震源谱和由标准窗(P波到时和PP波到时之间)得到的震源谱的比值。对扩展窗的分析是在假设它只包含P-pP-sP波群的情况下进行的。我们分析了发生在MW9.1级主震附近具有相似深度和震源机制的4个相对较小的地震事件。这些地震的震级范围从2005年1月9日的MW6.0级余震到震中位于苏门答腊-安达曼地震南部的2005年3月28日MW8.6级的尼亚斯地震。将得到的这4个地震事件的震源谱比值取平均,就得到扩展窗的频变算子。然后对扩展窗得到的2004年12月26日主震的震源谱进行校正,就得到苏门答腊-安达曼地震整个破裂过程(～600秒)的完整的或校正的震源谱。我们计算的地震辐射总能量为1.4×1017J。经过对整个地震的校正震源谱与破裂过程的前～250秒的震源谱(由标准远震窗得到)比较,我们发现主震破裂过程的前半部分辐射的地震波能量多于破裂过程后半部分辐射能量,尤其对于周期从3秒到40秒的地震波,这种现象更加明显。     

从2004年苏门答腊—安达曼巨型逆冲断层破裂得到的经验教训

发生于2004年的苏门答腊—安达曼地震,因其规模巨大,并且造成了毁灭性的后果,而得到广泛的研究。研究者通过分析大量高质量的地震、大地测量和地质数据,提出很多构造模型,以进一步了解这次地震的破裂长度、持续时间、断层几何形状、破裂速度和滑动过程等信息。一些最新的模型尽管在细节上有所不同,但在大尺度特征上都是一致的,例如,沿着整个长1300～1500km的破裂区域都有明显的同震滑动。块体滑动的速度非常快,足以辐射出地震波。由于这次地震的规模很大,使用传统方法进行数据处理会面临很多困难和挑战,这也促使人们创造和发展新的分析和反演方法,如多源反演、高频体波成像,以及通过卫星观测海啸浪高和重力变化等。苏门答腊巨型逆冲断层地震是40年来发生的最大地震,也是到目前为止资料记载保存最完好的地震,但它的一些特性与其他大的俯冲带地震看起来没有根本上的不同。     

Estimating the size of an earthquake using short-period seismograms of the first three seconds: A simulated experiment using the 1999 Chi-Chi earthquake sequence

Most of the present earthquake early warning systems are based on broadband or strong motion recordings. How-ever, the short-period instruments are still deployed. It is well-known that short-period recordings have saturation problems for large earthquakes when estimating the size of an earthquake. Thus, it is necessary to make clear the magnitude at which saturation starts to occur for the commonly used τc and Pd measurements, respectively. To investigate the possibility of using short-period seismic recordings for earthquake early warning, we conducted a simulated experiment using the strong motion data of the 1999 Chi-Chi earthquake sequence including its main shock and 31 aftershocks, with magnitude span from 4 to 7.6. The strong motion acceleration recordings were convolved with the instrument response of short-period seismographs in northern China to simulate short-period seismograms. Parameters τc and Pd from the first-three-second seismograms were calculated for the simulated short-period recordings and compared with that obtained by the original strong ground motion recordings. The result showed that to some extent, short-period recordings can be used for threshold earthquake early warning, while the magnitude saturation of Pd estimation can be up to 6.5, better than τc estimation.     

GRACE detection of the medium- to far-field coseismic gravity changes caused by the 2004 MW9.3 Sumatra-Andaman earthquake

Large earthquakes cause observable changes in the Earth's gravity field, which have been detected by the Gravity Recovery and Climate Experiment (GRACE). Since most previous studies focus on the detection of near-field gravity effects, this study provides the results from the medium- to far-field gravity changes caused by the 2004 Sumatra-Andaman earthquake that are recorded within GRACE monthly solutions. Utilizing a spherical-earth dislocation model we documented that large-scale signals predominate in the global field of the coseismic gravity changes caused by the earthquake. After removing the near-field effects, the coseismic gravity changes show a negative anomaly feature with an average magnitude of -0.18×10^-8 m·s^-2 in the region ranging ~40° from the epicenter, which is considered as the "medium field" in this study. From the GRACE data released by Center for Space Research from August 2002 to December 2008, we retrieved the large-scale gravity changes smoothed with 3 000 km Gaussian filter. The results show that the coseismic gravity changes detected by GRACE in the medium field have an average of (-0.20€±0.06)×10^-8 m·s^-2, which agrees with the model prediction. The detection confirms that GRACE is sensitive to large-scale medium-field coseismic gravitational effects of mega earthquakes, and also validates the spherical-earth dislocation model in the medium field from the perspective of satellite gravimetry.     

Wave separation for the great Sumatra-Andaman earthquake with regional seismic array

We analyzed the seismic waveforms from the December 26,2004 Sumatra-Andaman earthquake recorded at broadband seismic stations in western Europe.Previous studies involving of the beam-forming technique and high frequency analysis suggest that the earthquake ruptured with a duration of around 500 s.This very long duration makes P wave overlap with later arrivals such as PP wave,which follows P in about 200 s.Since P waves are crucial for modeling earthquake processes,we propose an iterative method to separate...     

The response of water level in the YuZ-5 well, Kamchatka to the magnitude 9.3, Sumatra-Andaman earthquake of December 26, 2004

Forced and free oscillations of water level were recorded in the YuZ-5 well, Kamchatka due to the passage of seismic waves from the Sumatra-Andaman earthquake of December 26, 2004, M _w = 9.3, hypocentral distance 8250 km. The greatest amplitude of water level oscillations, at least 5 cm, was observed during the onset of seismic surface waves with a typical period of 20–50 s. The total duration of the forced and free water level oscillations was about ten hours. The available theoretical models that describe oscillations of water level in a well due to seismic waves and rapid injection of water were used to estimate the transmissivity of the aquifer. The values obtained exceed by at least two orders of magnitude the transmissivity derived from pumping test measurements. A hypothesis was proposed to explain the temporary increase in aquifer transmissivity during the passage of seismic waves by invoking disturbances in the structure of the crack-pore space and a sharp increase in aquifer rock permeability.     

Spheroidal oscillations of the Earth stimulated by the Sumatra-Andaman earthquake with CDSN data

The 0S2~0S54 spheroidal modes of Earth’s free oscillations, triggered by the great Sumatra-Andaman earthquake of 26 December 2004 are retrieved from VHZ data recorded by seven upgraded stations of China Digital Seismograph Network (CDSN). We compare these spheroidal modes with theoretical free oscillation spectra calculated from the Preliminary Reference Earth Model (PREM) and find a coincidence between their periods. Spectral splitting phenomenon is observed obviously in 0S2, 0S3, 0S4, 2S1 and 1S2 free oscillation modes. What is most noticeable is that the oscillation mode 2S1 is reported for the second time (the first time by Rosat et al) without any data stacking. We simulated the split singlet of 0S2 mode on seven CDSN stations based on general focal mechanism and seismic moment of the earthquake. The result shows that seismic moment of the earthquake can reach 1023 N.m. We also find that the recording of Earth’s free oscillations carries abundant information of source mechanism and earthquake location, which is applicable to the detailed study of source rupture parameters.     

用前3秒P波估计地震大小的早期预警方法是否适用于短周期记录？-1999年集集地震序列的模拟实验

现有的地震早期预警系统大多采用宽频带记录或强震记录, 然而很多地震台站仍然使用短周期地震仪. 用短周期地震记录估计大地震震级时会出现饱和的现象早已众所周知. 但需要研究的问题是tau;c和Pd方法估计地震震级时出现饱和的震级到底是多少.为研究使用短周期记录进行地震早期预警的可能性,本文使用1999年集集地震序列的强震数据（包括主震和31次余震）进行模拟实验,所选地震序列的震级在4级到7.6级范围内. 用强震记录与短周期地震仪响应函数的卷积模拟短周期记录,利用其前3秒波形计算参数tau;c和Pd,并与用原始的强地面运动记录计算得到的结果进行比较. 结果表明, 在一定程度上,短周期地震记录仍可用作地震的早期阈值预警,并且Pd方法直到6.5级才出现饱和,好于tau;c方法.      

全球地震台网记录苏门答腊-安达曼地震的地幔波

介绍印度尼西亚苏门答腊-安达曼Ms8．7地震在全球地震台网（GSN）及所含“新一代中国数字地震台网”（NCDSN）记录的地幔面波。主要讨论本次大地震地幔瑞利波Ri,j（环球瑞利面波）无频散的3种可能原因,同时提出地幔妹泽波Mi,j的问题。     

Source parameters inversion of the 2013 Lushan earthquake by combining teleseismic waveforms and local seismograms

On April 20, 2013, a magnitude M s 7.0 earthquake occurred in Lushan, Sichuan Province, China, and caused heavy casualties and economic losses. Based on the local broadband waveforms in Sichuan and adjacent provinces regional networks and teleseismic broadband records from IRIS stations, the focal mechanism and the focal depth are determined by the CAP (Cut And Paste) and its upgraded methods, CAPtele and CAPjoint, respectively. The results show that the focal mechanisms and depth from different methods are steady, and the best double couple solution derived from the joint inversion is 210°, 44°, and 91° for strike, dip, and rake angles respectively for one nodal plane and 29°, 46°, and 89° for another with 16 km focal depth and M w 6.66 moment magnitude. In order to verify the reliability of the results, a number of tests are performed based on local seismograms with different velocity models. They indicate that there is about 10 degree’s fluctuation in focal mechanisms and about 2 km variation in focal depth with a complex velocity structure. Furthermore, inverted by re-sampling the teleseismic waveforms on the basis of epicentral distance, the solutions are consistent with each other, which manifests that the teleseismic records are effective for constraining source parameters of the Lushan earthquake.     

Coherency analysis of accelerograms recorded by the UPSAR array during the 2004 Parkfield earthquake

Spatial variability of near‐fault strong motions recorded by the US Geological Survey Parkfield Seismograph Array (UPSAR) during the 2004 Parkfield (California) earthquake is investigated. Behavior of the lagged coherency for two horizontal and the vertical components is analyzed by separately examining the decay of coherency with frequency and distance. Assumptions, approximations, and challenges that are involved in estimation of the coherency from recorded data are presented in detail. Comparison of the UPSAR coherency estimates with coherency models that are commonly used in engineering practice sheds light on the advantages and limitations of different approaches to modeling the coherency, as well as on similarities and differences in the spatial variability exhibited by seismic ground motion arrays at different sites. Copyright © 2013 John Wiley & Sons, Ltd.     

Field Survey of Tsunami Effects in Sri Lanka due to the Sumatra-Andaman Earthquake of December 26, 2004

The December 26, 2004 Sumatra-Andaman earthquake that registered a moment magnitude (M_w) of 9.1 was one of the largest earthquakes in the world since 1900. The devastating tsunami that resulted from this earthquake caused more casualties than any previously reported tsunami. The number of fatalities and missing persons in the most seriously affected countries were Indonesia - 167,736, Sri Lanka - 35,322, India - 18,045 and Thailand - 8,212. This paper describes two field visits to assess tsunami effects in Sri Lanka by a combined team of Japanese and Sri Lankan researchers. The first field visit from December 30, 2004 – January 04, 2005 covered the western and southern coasts of Sri Lanka including the cities of Moratuwa, Beruwala, Bentota, Pereliya, Hikkaduwa, Galle, Talpe, Matara, Tangalla and Hambantota. The objectives of the first field visit were to investigate the damage caused by the tsunami and to obtain eyewitness information about wave arrival times. The second field visit from March 10–18, 2005 covered the eastern and southern coasts of Sri Lanka and included Trincomalee, Batticaloa, Arugam Bay, Yala National Park and Kirinda. The objectives of the second visit were mainly to obtain eyewitness information about wave arrival times and inundation data, and to take relevant measurements using GPS instruments.     

利用区域宽频带数据反演鲁甸MS6.5级地震震源破裂过程

基于有限断层模型反演方法,我们利用区域宽频带数据反演得到了2014年8月3日鲁甸M_S6.5级地震的震源破裂过程.反演结果显示：此次地震的发震断层走向为北北西向,破裂主要以左旋走滑为主,位移主要发生在震源左上方,最大滑动量为0.7 m,模型显示断层破裂可能接近地表,破裂长度约10 km.此次地震释放的标量地震矩为1.97×10¹⁸ N·m,相当于矩震级为M_w 6.1,地震能量主要在前15 s释放.鲁甸地震有四个显著的特点：（1）位移主要集中在浅部,从11 km起破点开始迅速向上传播,大部分位于10 km以上且最大位移位于深度3 km处,从模型来看,破裂可能接近地表,因此地表震动较为强烈;（2）应力降比较大,计算显示释放的同震静态应力降约为2.8 MPa;（3）破裂速度较快,在地表附近超过了2.5 km·s^-1;（4）主震可能发生在一个共轭断层系上.这四个特点可能是导致此次地震造成如此重大人员伤亡和财产损失的最重要的原因.     

Bi-normalized response spectral characteristics of the 1999 Chi-Chi earthquake

To develop uniform and seismic environment-dependent design spectrum, common acceleration response spectral characteristics need to be identified. In this paper, a bi-normalized response spectrum (BNRS) is proposed, which is defined as a spectrum of peak response acceleration normalized with respect to peak acceleration of the excitation plotted vs. the natural period of the system normalized with respect to the spectrum predominant period, Tp. Based on a statistical analysis of records from the 1999 Chi-Chi earthquake, the conventionally normalized response spectrum(NRS) and the BNRS are examined to account for the effects of soil conditions, epicentral distance, hanging wall and damping. It is found that compared to the NRS, the BNRS is much less dependent on these factors. Finally, some simple relationships between the BNRS for a specified damping ratio and that for a damping ratio of 5%, and between the spectra predominant period and epicentral distance for different soil types are provided.     

采用体波谱振幅相关系数方法研究地震序列的震源机制变化过程

发生在同一震源区的两个地震, 当地震之间的距离远小于地震到台站之间的距离时, 同一个台站的直达P波、 S波零频观测谱值之比与其辐射图型因子之比线性相关, 而辐射图型因子取决于震源机制解的几何特性。 利用同一震源区两次事件的相同台站的波形记录计算谱振幅相关系数, 该系数可认为是震源机制异同的归一化的数值描述; 据此对地震序列事件进行聚类分组以及计算序列事件的滑动平均相关系数时序曲线。 对2001年雅江和2003年大姚强震序列, 以及2002年漾濞和2003年洱源小震序列进行分析的研究结果表明: ① 强震序列与小震序列的滑动平均相关系数时序曲线形态存在明显差异; ② 某些强震序列的初期震源机制有从散乱到一致的过程。 分析认为, 孕震应力场的应力强度对孕震区内地震的破裂方向有明显的约束作用。     

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.