2009年3月19日Mw7.6级汤加大地震的“震前扰动”现象

汶川大地震爆发后,作者在武汉大地测量国家野外科学观测站的重力仪资料中发现了震前约48小时的“重力扰动”现象、并在随后对全国数十个台站宽带地震仪资料的分析中得到了证实.2009年3月19日,太平洋岛国汤加海域爆发了M_w7.6级（里氏7.9级）大地震,宽带地震仪资料中再次出现了显著的、信号周期为4~8秒的“震前扰动”现象：从大地震前约16天的3月4日开始一直持续到地震发生.     

IRIS台网地震仪资料出现2010年1月12日海地大地震"震前扰动现象"

2010年1月12日,加勒比岛国海地爆发了Mw7.0级大地震.震中附近的宽带地震仪记录在1月10—1月12日观测到了显著的“震前扰动”现象.扰动信号的周期为3~10秒.北美大陆的宽带地震仪也记录到了类似的扰动信号.海地大地震“震前扰动”发生期间,加勒比海域及周围海域没有发生任何海洋风暴.海地大地震的“震前扰动”是由某种非台风因素引发的.     

宽带地震仪资料证实汶川大地震“震前重力扰动”

在武汉大地测量国家野外科学观测站LacosteET重力仪资料中发现汶川大地震"震前扰动"现象之后,又在全国数十个台站宽带地震仪资料中发现了汶川大地震的"震前扰动".和LacosteET重力仪观测到的"震前扰动"一样,宽带地震仪资料的"震前扰动"也是在震前两天开始并逐渐增大,一直持续到汶川大地震发生,信号周期也是4～8 s.     

强台风对汶川大地震和昆仑山大地震“震前扰动”影响的分析

2008年8.0级汶川大地震发生前约40 h起,中国地震局测震台网中的宽带地震仪、重力仪和倾斜仪都记录到了明显的低频扰动现象.然而,在此期间西太平洋上空产生了一次强台风Rammasun（威马逊）;与此类似,2001年8.1级昆仑山大地震发生前约60 h起,国家地震局测震台网中的宽带地震仪也记录到了明显的低频扰动现象,在此期间中国南海上空产生了一次强台风LingLing（玲玲）.本文针对这两大地震的“震前扰动”现象进行了研究,旨在探讨其“震前扰动”现象是否与台风引起的巨浪产生的地脉动效应有关.我们的研究结果表明：在中国内陆观测到的异常扰动与台风路径以及台风的强度密切相关.近海台风可在中国大陆内地引起明显的异常扰动;而远离大陆海岸位于深水海域上空的台风,很难在中国大陆内地引起明显的异常地脉动变化.台风LingLing的路径距离中国海岸最近处仅有约700 km,而台风Rammasun距离中国大陆海岸最近处却达1300 km.本文的对比分析研究表明：路径和强度类似于LingLing的台风在中国内陆可引起明显的异常扰动.因此,中国内陆和沿海许多地震台记录到的昆仑山大地震的“震前扰动”源于LingLing地脉动效应的可能性较大;但位于新疆的地震台记录到的扰动信号频率却明显偏低,这一现象值得注意.路径和强度类似于Rammasun的台风很难在中国内陆引起明显的异常扰动,因此汶川大地震的“震前扰动”完全源于台风Rammasun地脉动效应的可能性较小.汶川大地震的“震前扰动”产生的原因和机理值得我们去进行深入的探讨和研究.     

“震前扰动”现象:从模糊走向清晰——与傅容珊教授商榷

本文介绍了作者近期研究发现的汶川大地震“震前扰动”信号的“时间偏移”现象和昆仑山大地震“震前扰动”信号的“频率偏移”现象,并对“震前扰动”与“台风影响”的关系进行了讨论.不同意《地震前兆还是其它因素?——与“汶川大地震宽带地震仪短临异常及成因初探”作者商榷》将汶川大地震“震前扰动”的原因简单地归咎于台风并完全排除其与地...     

汶川大地震前非台风扰动现象的研究

为了区分汶川大地震的震前扰动现象中的台风因素和非台风因素,本文研究了中国大陆宽带地震仪在汶川大地震前记录到的异常扰动信号的时频特征.研究结果表明汶川大地震的震前扰动主要由两种扰动构成,二者动态特征完全不同.其中优势频率为0.2～0.25 Hz的扰动主要与台风Rammasun有关.这种台风扰动在沿海地区较强,在内陆地区较弱,其震动源在靠近台风运动路径的海底.另一种优势频率为0.1～0.18 Hz的扰动与台风无关,这种非台风扰动在地震发生前约10 h突然急剧增强,其最大值出现在地震爆发时刻.非台风扰动在靠近震中的地区较强,在沿海和西部地区较弱.震源扰动扫描算法计算初步定位的结果显示其震动源不在海底,而是分散在震中附近的内陆地区.汶川大地震前的非台风扰动是否与汶川大地震有关,值得进行更深入的研究.     

汶川大地震“震前扰动”存在“第三类脉动”吗？

采用国家地震台网数十台地震仪的观测资料研究了汶川大地震和昆仑山大地震的“震前扰动”现象,结果表明,“震前扰动”现象与强台风的影响密切相关,适用于英国海洋学家Higgins的海浪非线性干涉理论中的“第二类脉动”.然而,通过对信号时频特征的细致分析后发现：“震前扰动”信号存在着明显的“时频偏移”现象,不完全符合“第二类脉动”的特征.因此本文认为,汶川大地震的“震前扰动”除了包含“第二类脉动”信号外,还可能存在着“第三类脉动”,这对强震短临前兆的研究有重要意义.     

2001年11月14日昆仑山大地震震前扰动信号的再研究

研究了2001年11月14日昆仑山M_S8.1大地震前新疆地区宽频带地震仪记录到的异常扰动事件.新疆地震台网在发生强震的前三天记录到了强烈的扰动信号.由于新疆地震台距离昆仑山大地震的震中很近,因此该震前扰动的成因备受关注,一直被认为可能是与大地震相关的慢地震事件,可能是大地震的某种前兆现象.本研究扩大了观测异常扰动事件的范围,分析了欧亚大陆187个宽频带地震仪在2001年11月的观测记录,发现这些地震仪在11月10—12日期间都记录到了时频特征类似的强扰动信号,其中最强的扰动信号出现在挪威台站KONO.调查研究发现挪威海域在11月10—12日期间出现了一次超强北极风暴,我们对扰动源位置初步定位后发现大地震前欧亚大陆各台站记录的最大扰动信号主要来源于挪威西南海岸.本文的研究结果确认昆仑山大地震前新疆台网记录到的强烈震前扰动信号不是与大地震相关的慢地震事件,而是挪威海域北极风暴激发的微地震（地脉动）.本文进一步研究表明北极风暴激发的地脉动可显著影响中国大陆的地震和重力观测.在分析中国大陆地震台网记录的强扰动信号时,北极风暴的影响不能忽视.     

超导重力仪检测2011年日本东北Mw9.0级地震前的重力扰动信号

对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 为检验日本M_w9.0级地震是否存在震前扰动现象,本研究利用全球超导重力仪记录到的地震前后7天内20组秒采样数据进行分析. 经潮汐、大气改正等处理去除仪器的漂移及残余潮汐效应,得到非潮汐重力变化曲线.结果表明大部分振幅大于30×10^-8 m·s^-2的曲线反映了全球M_w≥6级地震引起的高频波动信号,其中11组数据在3月9日M_w7.3级前震之前出现了扰动现象.震前扰动可分解为三个频段,其中,低于0.1 Hz和高于0.18 Hz的分量分别反映了地震波动信号及非构造信息,中间频段(0.118~0.18 Hz)信号能够较大程度地压制地震波动信号、并同时保留异常扰动信息.它的振幅在3月7日10时之前基本保持约1×10^-8 m·s^-2,之后开始逐渐增大,到3月9日7.3级前震前后达到最大,此后振荡衰减,振幅保持约(5~10)×10^-8 m·s^-2,直至主震发生.中间频段信号的变化特征与主震前的应力迁移过程以及实验记录到的地震成核过程有许多相似之处;不过,震前重力异常是否与主震前的应力加速积累有关,仍待进一步研究.     

广西苍梧M_S5.4地震前华南地区连续重力异常分析

正实践表明,大震前的异常扰动现象已被高精度重力仪和宽频带地震仪监测到。如Lacaste ET重力仪和宽频带地震仪监测和证实了汶川地震前1—2天内存在周期为4—8 s的重力扰动,Gphone重力仪监测到日本M_S9.0大震前的连续观测重力异常和高频扰动,以及宽频带地震仪监测到的汤加M_W7.6地震和海地M_W7.0地震前的"震前扰动"现象。到目前为止,有关重力仪记录到的震前"重力扰动"及宽频带地震仪等仪器记录到的"震前扰动"的报道大都来自7.0     

汶川地震前后波速比变化特征的再研究

收集整理了四川省地震台网1990—2012年产出的地震直达波震相数据, 利用单台多震和达法对2008年汶川M_S8.0地震前后四川地区的波速比变化特征进行了分析研究(共筛选出13个数据量较为丰富、 连续性较好的台站). 结果显示: 其中位于龙门山断裂中北段西侧的4个台站, 震前出现长达7年左右的中长期波速比低值异常; 其它9个台站的波速比震前变化基本稳定. 震前波速比出现异常的4个台站的分布与汶川M_S8.0地震孕震区范围大体一致, 从而为研究汶川M_S8.0地震前是否存在地壳介质特性的时间变化过程, 提供了有力的判定依据.      

Low CodaQcin the Epicentral Region of the 2001 Bhuj Earthquake ofMw7.7

On 26 January, 2001 (03:46:55,UT) a devastating intraplate earthquake of M_w 7.7 occurred in a region about 5 km NW of Bhachau, Gujarat (23.42°N, 70.23°E). The epicentral distribution of aftershocks defines a marked concentration along an E-W trending and southerly dipping (45°) zone covering an area of (60 × 40) km². The presence of high seismicity including two earthquakes of magnitudes exceeding 7.7 in the 200 years is presumed to have caused a higher level of shallow crustal heterogeneity in the Kutch area; a site lying in the seismic zone V (zone of the highest seismicity for potentially M8 earthquakes) on the seismic zoning map of India. Attenuation property of the medium around the epicentral area of the Bhuj earthquake covering a circular area of 61,500 km² with a radius of 140 km is studied by estimating the coda-Q_c from 200 local earthquakes of magnitudes varying from 3.0–4.6. The estimated Q₀ values at locations in the aftershock zone (high seismicity) are found to be low in comparison to areas at a distance from it. This can be attributed to the fact that seismic waves are highly scattered for paths through the seismically active and fractured zone but they are well behaved outside the aftershock zone. Distribution of Q₀ values suggests that the local variation in Q₀ values is probably controlled by local geology. The estimated Q₀ values at different stations suggest a low value of Q=(102 ± 0.80)^*f^{(0.98 ± 0.02)} indicating an attenuative crust beneath the entire region. The frequency-dependent relation indicates a relatively low Q_c at lower frequencies (1–3 Hz) that can be attributed to the loss of energy due to scattering attenuation associated with heterogeneities and/or intrinsic attenuation due to fluid movement in the fault zone and fluid-filled cracks. The large Q_c at higher frequencies may be related to the propagation of backscattered body waves through deeper parts of the lithosphere where less heterogeneity is expected. Based on the attenuation curve estimated for Q₀=102, the ground acceleration at 240 km distance is 13% of 1 g i.e., 0.13 g agreeing well with the ground acceleration recorded by an accelerograph at Ahmedabad (0.11 g). Hence, it is inferred that the Q₀ value obtained from this study seems to be apt for prediction of ground motion for the region.     

九寨沟MS 7.0地震前单台波速比变化特征

收集、整理四川省地震台网2013年5月至2017年9月产出的震相观测报告数据,采用单台多震和达法,重点分析2017年8月8日九寨沟M_S 7.0强震前龙门山断裂带中北段地区波速比变化特征。结果表明：龙门山断裂带中段西侧的茂县地震台,在发震前2年左右,波速比出现一次持续时间约一年半的下降异常;发震前半年左右,波速比出现明显回返特征,波速比回返至均值上下时即发生九寨沟M_S 7.0地震;震源区周边其他6个地震台波速比在发震前变化较平稳。茂县地震台波速比低值异常对判定九寨沟M_S 7.0强震前龙门山断裂带中北段可能存在孕震体具有重要参考价值。     

Ultrasonic measurements of Vp and Qp: relaxation spectrum of complex modulus on basalt melts

Ultrasonic compressional wave velocity V_p and quality factor Q_p have been measured in alkali basalt, olivine basalt and basic andesite melts in the frequency range of 3.4–22 MHz and in the temperature range of 1100–1400°C. Velocity and attenuation of the melts depend on frequency and temperature, showing that there are relaxation mechanisms in the melts. Complex moduli are calculated from the ultrasonic data. The results fit well a complex modulus of Arrhenius temperature dependence with log-normal Gaussian distribution in relaxation times of attenuation. The analysis yields average relaxation time, its activation energy, relaxed modulus, unrelaxed modulus and width of Gaussian distribution in relaxation times. Relaxed modulus is smaller (17.5 GPa) for basic andesite melt of high silica and high alumina contents than for the other two basalt melts (18.1–18.4 GPa). The most probable relaxation times decrease from ～ 3 × 10^?10 s for basic andesite to ～ 10^?11 s for alkali basalt at 1400°C. Activation energies of attenuation, ranging from 270 to 340 kJ mol^?1 in the three melts, are highest in basic andesite. Longitudinal viscosity values and their temperature dependences are also calculated from V_p and Q_p data. The volume viscosity values are estimated from the data using the shear viscosity values. Longitudinal, volume and shear viscosities and their activation energies are highest in the basic andesite melt of the most polymerized structure.     

Vp and Vp/Vs structures in the crust and upper mantle of the Taiwan region, China

A tomographic study of the V _p and V _p/V _s structures in the crust and upper mantle beneath the Taiwan region of China is conducted by simultaneous inversion of P and S arrival times. Compared with the previous tomographic results, the spherical finite difference technique is suitable for the strong heterogeneous velocity structure, and may improve the accuracy in the travel time and three-dimensional ray tracing calculations. The V _p and V _p/V _s structures derived from joint inversion and the relocated earthquakes can provide better constraints for analyzing the lateral heterogeneity and deep tectonic characters in the crust and upper mantle. Our tomographic results reveal significant relations between the seismic wavespeed structure and the tectonic characters. In the shallow depth, sedimentary basins and orogen show distinct wavespeed anomalies, with low V _p, high V _p/V _s in basins and high V _p, low V _p/V _s in orogen. As the suture zone of Eurasian Plate and Philippine Sea Plate, Longitudinal Valley is characterized by a significant high V _p/V _s anomaly extending to the middle-lower crust and upper mantle, which reflects the impact of rock cracking, partial melting, and the presence of fluids. In the northeast Taiwan, the V _p, V _p/V _s anomalies and relocated earthquakes depict the subducting Philippine Sea Plate under the Eurasian Plate. The high V _p of oceanic plate and the low V _p, high V _p/V _s atop the subducted oceanic plate extend to 80 km depth. Along the east-west profiles, the thickness of crust reaches 60 km at the east of Central Range with eastward dipping trend, which reveals the eastward subduction of the thickened and deformed crust of the Eurasian continental plate. Supported by Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KZCX3-SW-234-2), National Basic Research Program of China (Grant No. 2007CB411701), National High Technology Research and Development Program of China (Grant No. 2006AA09A101-0201) and National Natural Science Foundation of China (Grant Nos. 40804016, 40704013)     

兴蒙造山带诺敏河火山群地壳厚度与波速比研究

利用布设于兴蒙造山带诺敏河火山群地区的宽频带流动地震台站资料,基于接收函数方法,获取了该地区的地壳厚度与波速比值.研究结果显示,该地区的地壳厚度介于32~38 km,莫霍面深度在空间上分布特征与五大连池为中心的火山带分布具有较好的一致性：沿着火山带延展方向地壳较薄.该地区的波速比介于1.74~1.84,波速比在空间上与地壳厚度变化具有一致性：高波速比主要集中于靠近五大连池火山带地区,向诺敏河火山和小古里河火山延展.研究认为：诺敏河火山与五大连池火山带可能具有相同的岩浆来源,可能与富钾岩石圈地幔拆沉作用造成的地幔热物质上涌有关.研究区地壳厚度与波速比呈现负相关关系,暗示该地区可能发生过岩浆底侵作用.     

Earthquake prediction on the basis of Vp/Vs variations — a case history

On November 7, 1976, an earthquake of the strike-slip fault type, and of magnitude 6.7, occurred at two mountainous localities, Yanyuan and Ninglang, in the border region of the Szechuan and Yunnan Provinces, China. One year before the earthquake, a prediction was made as to the magnitude and location of the impending earthquake by the present authors, on the basis of the results of a general survey of an area of about 20,000 km² for V_p/V_s ratio variations. The prediction of time of occurrence was made afterwards by the combined analysis of results of observations of some other precursory phenomena. The actual occurrence of the event was generally considered as being in agreement with the prediction.The present study may be taken as a new example for the detection of V_p/V_s ratio variations prior to an earthquake of the strike-slip fault type. By considering the difference between the shear rupture and stick-slip motion and the anisotropy induced by dilatancy, a preliminary discussion is made concerning the related results.     

永清MS4.3地震和廊坊MS3.0地震的发震构造研究

基于首都圈数字地震台网的宽频带资料,首先采用CAP方法确定了永清M_S4.3地震和廊坊M_S3.0地震的震源机制解:永清地震节面Ⅰ的走向、倾角和滑动角分别为52°,62°和?140°,节面Ⅱ的走向、倾角和滑动角分别为300°,55°和?35°;廊坊地震节面I的走向、倾角和滑动角分别为48°,57°和?147°,节面Ⅱ的走向、倾角和滑动角分别为299°,63°和?38°。两次地震的震源机制解较为一致,推测它们可能具有相同的发震断层。利用近震转换波获得两次地震的震源深度,分别为19 km和13 km。利用双差法对两次地震的主余震进行重新定位,结果显示:两个地震序列的震中均呈NE向分布,余震震源深度均浅于主震震源深度,震源深度分别集中在17—20 km和12—13 km范围内,两个序列的短轴剖面揭示了震源分布均呈现倾向SE,倾角陡立的特点。将地震序列的分布与震源机制解的结果进行对比,认为两个序列的水平展布方向与其对应的主震震源机制解中节面Ⅰ的走向比较接近,深度分布的高倾角特征也与节面Ⅰ比较相似,因此认为发震断层面均为节面Ⅰ。通过将震源机制解中节面Ⅰ的参数和地震序列的分布与区域活动断层的产状性质进行比较,取得了一些关于发震构造和地震成因的重要认识:① 永清M_S4.3地震和廊坊M_S3.0地震的发震构造不是上地壳的先存正断裂?河西务断裂,不排除与中下地壳的新生构造或深大断裂有关;② 永清、廊坊地震发生在13—19 km深度上,结合地壳结构、断裂构造以及区域流变结构等资料,推测该深度范围可能是廊固凹陷的壳内脆性?韧性转换区域,是地震孕育和发生的有利构造部位。      

87Rb/87Sr study of diogenites

The five diogenites, Johnstown, Roda, Ellemeet, Shalka and Tatahouine, give scattered data in the⁸⁷Rb/⁸⁶Sr,⁸⁷Sr/⁸⁶Sr diagram. This can result from a disturbance which occurred later than 4.45 Ga ago. However, it is shown that if samples of sufficient size were analyzed, there meteorites could plot on the eucrite isochron and are thereby in agreement with a genetic relation between eucrites, howardites and diogenites. The age of eucrite differentiation from diogenites has been computed using data from the two families yielding an age of 4.47±0.1Ga(2σ) (λ=1.42×10^?11a^?1), the initial⁸⁷Sr/⁸⁶Sr ratio being BABI.