Relationship between earthquake and volcanic eruption inferred from historical records

Introduction From the records of Wudalianchi volcanic group eruption in 1720~1721 obtained from the Man ethnic group files of Heilongjiang General Yamen in Qing Dynasty (WU, 1998; CHEN, WU, 2003), we have discovered the eruption time, state, material and scale of Laoheishan and Hu-oshaoshan volcanoes, as well as numerous seismic records. These historical materials are discov-ered for the first time although the study on Wudalianchi volcanic group has a long history. These earthquakes co…     

Attenuation of peak ground accelerations from the great Wenchuan earthquake

Over 800 accelerograms recorded by 272 ground-level stations during the Wenchuan earthquake are used to analyze the influence of rupture distance, local site conditions and azimuth on peak ground accelerations (PGAs). To achieve a better understanding of the characteristics of ground motions, the spatial distributions of the EW, NS and UD components of PGAs arc obtained. Comparisons between the EW and NS components, the fault-normal and fault-parallel components,and the vertical and horizontal components of PGAs are performed, and the regression formula of the vertical-to-horizontal ratio of PGAs is developed. The attenuation relationship of peak horizontal accelerations (PHAs) is compared with several contemporary attenuation relationships. In addition, an analysis of residuals is conducted to identify the potential effects of rupture distance, azimuth and site conditions on the observed values of PHAs. The analysis focuses on medium-hard soil site conditions, as they provided most of the data used in this study.     

Attenuation in the uppermost inner core from PKP recordings at African seismological stations

The attenuation factor Q_P at the top of the inner core is evaluated by using the amplitude spectral ratio of PKPdf and PKPbc phases observed at African stations (BGCA mostly), from strong deep earthquakes in the Pacific Ocean area. The maximum depth of penetration of the PKPdf phase into the inner core (IC) is roughly 377 km, and the sampled region of IC is centered beneath the Southern Indian Ocean. The derived mean value of Q_P is 249 ± 31 (95% confidence level) in the frequency range 0.2–2 Hz, where no frequency dependence of attenuation has been reliably observed. By using Student’s t-test, we show that the value is statistically significantly different (with a probability greater than 95%) from other mean values of Q derived by using the same method, for both the western (180 °W to 40 °E) and eastern (40 °E to 180 °E) hemispheres of the IC. The decrease of Q with the radius of the turning point (denoted by r^TP), according to Q_P = 840 − 0.62 r^TP, has a moderate statistical support (the R-squared value is 38%). A slightly increase of Q as a function of the angle of the PKPdf path within the inner core with respect to the Earth’s spin axis is observed, in agreement with various investigations performed in the time domain. However, the value of the anisotropy, if any, is suggested to be around 3%.     

由精确震源机制推断的日本海岸地震带的局部应力集中:板内地震断层应力积累过程的意义

在日本西南的板内地区发现了最大主应力轴（σ1）方位角的空间变化。根据大量精确震源机制进行应力反演推测,沿日本海岸地震带的口。轴方位角取向为N110°E～N130°E,而在周边地区几乎为东西方向,即N90°E～N100°E。沿南海海沟内陆板块的地壳浅部广泛观测到东西方向的最大水平压应力。然而,在该地震带上仅观测到WNW-ESE方向的最大水平应力。用该地震带之下下地壳中的无震断层或延性断层带的形变解释了该地震带内及其周围应力场的这种空间变化。     

Fractal analysis of earthquake swarms of Vogtland/NW-Bohemia intraplate seismicity

The special type of intraplate microseismicity with swarm-like occurrence of earthquakes within the Vogtland/NW-Bohemian Region is analysed to reveal the nature and the origin of the seismogenic regime. The long-term data set of continuous seismic monitoring since 1962, including more than 26000 events within a range of about 5 units of local magnitude, provides an unique database for statistical investigations. Most earthquakes occur in narrow hypocentral volumes (clusters) within the lower part of the upper crust, but also single event occurrence outside of spatial clusters is observed. Temporal distribution of events is concentrated in clusters (swarms), which last some days until few month in dependence of intensity. Since 1962 three strong swarms occurred (1962, 1985/86, 2000), including two seismic cycles. Spatial clusters are distributed along a fault system of regional extension (Leipzig-Regensburger Störung), which is supposed to act as the joint tectonic fracture zone for the whole seismogenic region. Seismicity is analysed by fractal analysis, suggesting a unifractal behaviour of seismicity and uniform character of seismotectonic regime for the whole region. A tendency of decreasing fractal dimension values is observed for temporal distribution of earthquakes, indicating an increasing degree of temporal clustering from swarm to swarm. Following the idea of earthquake triggering by magma intrusions and related fluid and gas release into the tectonically pre-stressed parts of the crust, a steady increased intensity of intrusion and/or fluid and gas release might account for that observation. Additionally, seismic parameters for Vogtland/NW-Bohemia intraplate seismicity are compared with an adequate data set of mining-induced seismicity in a nearby mine of Lubin/Poland and with synthetic data sets to evaluate parameter estimation. Due to different seismogenic regime of tectonic and induced seismicity, significant differences between b-values and temporal dimension values are observed. Most significant for intraplate seismicity are relatively low fractal dimension values for temporal distribution. That observation reflects the strong degree of temporal earthquake clustering, which might explain the episodic character of earthquake swarms and support the idea of push-like triggering of earthquake avalanches by intruding magma.     

基于SAR影像偏移量获取汶川地震二维形变场

本文以ALOS卫星PALSAR影像为数据源,采用强度图像偏移量方法获得的整个汶川地震地表二维形变场显示,整个映秀-北川地表破裂带全长约240 km,从西南端的虹口往北东方向一直延伸到青川县附近,在虹口及北川县城所在地为两个形变量最大区域,偏移量可达4~6 m,局部更是达到了6~8 m.在高川乡附近出现一斜列拉分阶区,宽约8~10 km.在映秀-北川断层的地表破裂迹线南侧约12 km处还有一条汉旺-白鹿次级破裂带,从漩口镇一直延伸到秀水镇,长度大约100 km,在白鹿附近形变量较大,可达3~4 m.另外在小鱼洞附近可见一个NW 走向、长宽约10×5 km、形变幅度达3~4 m的连接以上两条破裂带的地表破裂带,性质为逆冲兼具左旋走滑.研究表明,利用SAR影像偏移量法能够获取近场几米量级的大形变量及客观揭示断层破裂迹线的真实形态和分段特征,可以成为野外观测、InSAR等手段的有益补充,综合以上几种观测手段,优势互补,我们可以构建更为真实的断层模型,进而对汶川地震的复杂破裂过程有更深入的了解.     

Source characteristics of the 1992 Nicaragua tsunami earthquake inferred from teleseismic body waves

We analyzed the broadband body waves of the 1992 Nicaragua earthquake to determine the nature of rupture. The rupture propagation was represented by the distribution of point sources with moment-rate functions at 9 grid points with uniform spacing of 20 km along the fault strike. The moment-rate functions were then parameterized, and the parameters were determined with the least squares method with some constraints. The centroid times of the individual moment-rate functions indicate slow and smooth rupture propagation at a velocity of 1.5 km/s toward NW and 1.0 km/s toward SE. Including a small initial break which precedes the main rupture by about 10 s, we obtained a total source duration of 110 s. The total seismic moment isM _o =3.4×10²⁰ Nm, which is consistent with the value determined from long-period surface waves,M _o =3.7×10²⁰ Nm. The average rise time of dislocation is determined to be 10 s. The major moment release occurred along a fault length of 160 km. With the assumption of a fault widthW=50 km, we obtained the dislocationD=1.3 m. From andD the dislocation velocity isD=D/0.1 m/s, significantly smaller than the typical value for ordinary earthquakes. The stress drop =1.1 MPa is also less than the typical value for subduction zone earthquakes by a factor of 2–3. On the other hand, the apparent stress defined by 2E _s /M _o , where andE _s are respectively the rigidity and the seismic wave energy, is 0.037 MPa, more than an order of magnitude smaller than . The Nicaragua tsunami earthquake is characterized by the following three properties: 1) slow rupture propagation; 2) smooth rupture; 3) slow dislocation motion.     

使用小震震源机制解研究山东地区背景应力场

基于山东地区1970—2012年5月387次中小地震震源机制解数据,使用Gephart和Forsyth的 FMSI方法计算了该区的平均构造应力方向; 使用区域应力张量阻尼反演方法,结合共轭梯度法求解了应力张量的阻尼最小二乘最优化问题,并根据Lund和Townend推导给出的由构造应力张量计算最大水平主应力方向的公式,反演计算了山东地区构造应力场的空间分布特征,并讨论了该区构造应力场的变化及其与郯庐断裂带的关系．     

Attenuation models of the earth

Free oscillation and body wave data are used to construct average Q models for the earth. The data set includes fundamental and overtone observations of the radial, spheroidal and toroidal modes, ScS observations and amplitudes of body waves as a function of distance. The preferred model includes a low-Q zone at both the top and the bottom of the mantle. In these regions the seismic velocities are likely to be frequency dependent in the “seismic” band. Absorption in the mantle is predominantly due to losses in shear. Compressional absorption may be important in the inner core.A grain-boundary relaxation model is proposed that explains the dominance of shear over compressional dissipation, the roughly frequency independent average values for Q and the variation of Q with depth. In the high-Q regions, the lithosphere and the midmantle (200–2000 km), Q is predicted to be frequency dependent. However, the low-Q regions of the earth, where Q is roughly frequency independent, dominate the observations of attenuation.     

Attenuation of high-frequency earthquake waves in South America

Earlier attenuation studies for the South American continent indicate that for Sn and Lg waves there is low attenuation in the shield region east of the Andes, whereas in the west of South America, with some exceptions, there is high attenuation. Those studies, however, were nonquantitative. In this study theQ ₀ (1-Hz values) for Lg waves for South America are presented, based on a scattering model ofAki (1969) as extended byHerrmann (1980) for the coda waves of shallow local and near-regional earthquakes. The results of the codaQ method are compared with those obtained byNuttli's method (1973). TheseQ ₀ values are in good agreement with the apparentQ ₀ of Lg waves obtained by the latter method. The data were obtained from over 100 local and regional earthquakes recorded by 12 WWSSN stations throughout continental South America. They provided a range of frequencies from 0.4 to 1.4 Hz, for which the frequency dependence ofQ was investigated by assumingQ=Q ₀(f/f ₀). The observed data indicate that the tectonic region of western South America is characterized by lowQ ₀ and a large frequency-dependent factor , the values ranging from 150 to 350 and 0.4 to 0.7, respectively. TheQ ₀ values increase in the shield region east of the Andes, but frequency dependence decreases. The average crustalQ ₀ for north and central Argentina ranges from 420 to 580, and ranges from 0.2 to 0.3. TheQ ₀ is larger in the Brazil region, ranging from 580 to 980, and varies from 0.0 to 0.2. In the lower-attenuation region of eastern South America higher values of attenuation correlate with greater thickness of the sedimentary layers.     

近震源地震波波形资料反演震源破裂过程的可靠性分析

用近震源波形资料拟合反演地震的震源破裂过程,所包含的一些不确定因素将对反演结果的精度及可靠性产生影响,文中的数值实验分析了所假定的反演断层模型参数的某些不确定性对反演结果的影响程度,并对观测波形的截取长度对反演精度的影响进行了讨论．结果表明：(1)近震源地震波形资料能较好地分辨断层浅部的破裂过程.然而对断层深部的位错分布的约束和反演能力较差．联合使用近、远场地震波资料进行反演,能反演出一个更为完全的整个断层破裂过程的图像.(2)用近震源地震波资料反演时,反演结果对所假定的反演断层的走向和倾角非常敏感.断层走向偏离真实值2°或倾角偏离真实值5°都会导致一个虚假的反演结果.(3)反演中所使用的介质速度结构模型的不确定性,也会对反演结果产生影响.     

利用InSAR资料反演缅甸Mw6.8地震断层滑动分布

2011年3月24日缅甸东北部发生M_w6.8级地震.本文利用覆盖该地区的升降轨ALOS PALSAR数据,获取了该次地震的同震形变场,并采用灰度配准技术获取了其地表破裂位移.针对影像中因轨道不精确造成的非线性长波长误差,本文采用二次多项式曲面法予以去除,获取了更为精确的同震形变场.最后,基于弹性半空间形变模型反演了该地震的断层滑动分布.结果表明,该地震断层滑动以左旋走滑为主,兼具少许的倾滑运动分量,断层滑动主要集中分布在断层面0~10 km深度范围,最大滑动量达5 m,位于地表以下5 km深处.反演获得的地震标量矩为2.49×10¹⁹N·m,震级约为M_w6.8级.     

利用GPS和InSAR数据反演2011年日本东北MW9.0地震断层的同震滑动分布

本文首先对Envisat/ASAR数据进行干涉处理, 获取2011年日本东北M_W9.0地震的地表InSAR同震形变场; 然后通过对InSAR同震形变数据重采样方法的深入分析, 选择条纹率法结合干涉图的空间相干性对InSAR同震形变数据进行重采样; 最后基于弹性半空间位错模型, 联合InSAR与GPS形变数据, 采用最小二乘法反演发震断层的滑动分布. 研究结果表明: 考虑相干性的条纹率重采样方法, 更适用于形变场中存在除断层外的有限边界、 且形变场范围较大的InSAR数据重采样处理; 断层滑动主要发生在地表以下50 km范围内, 最大滑动量为49.9 m, 矩张量为4.89×1022 N·m, 所对应的矩震级为M_W9.1, 与地震学反演的结果比较吻合.      

汶川Ms8．0级地震成因三“层次”分析——基于深部电性结构

2008年5月12日汶川Ms8级地震的发生不是局部地区孤立的构造事件,研究汶川地震的孕震机制,应该把局部分析和区域分析相结合,关注地壳上地幔直至地幔过渡带的深部结构.基于近年来在东北、华北和汶川地震附近地区进行的深部结构电磁探测结果,结合地震学等其他资料,从太平洋板块的俯冲、印度板块的碰撞和松潘甘孜地块的推挤三个"层次"探讨分析汶川特大地震的成因.太平洋板块向亚洲大陆的俯冲作用,导致中国大陆东部地幔过渡带深度较普遍地存在着停滞的板片,它对汶川地震的影响不可忽视.印度板块与青藏高原的碰撞,使组成高原的各地块发生向北和向东的运动,各地块向东的运动作用于南北地震带中南段,影响到该区域的地震活动.松潘甘孜地块向四川地块的推挤.使松潘甘孜地块运动方向和龙门山断裂带形成"丁"字形结构,龙门山断裂带显示为较陡直的电性边界,加剧了汶川地震前的应力积累,可能是汶川地震发生的最直接的诱因.     

Renewal models for earthquake predictability

The new Database of Italy’s Seismogenic Sources (Basili et al. 2008) identifies areas with a degree of homogeneity in earthquake generation mechanism judged sufficiently high. Nevertheless, their seismic sequences show rather long and regular interoccurrence times mixed with irregularly distributed short interoccurrence times. Accordingly, the following question could naturally arise: do sequences consist of nearly periodic events perturbed by a kind of noise; are they Poissonian; or short interoccurrence times predominate like in a cluster model? The relative reliability of these hypotheses is at present a matter of discussion (Faenza et al., Geophys J Int 155:521–531, 2003; Corral, Proc Geoph 12:89–100, 2005, Tectonophysics 424:177–193, 2006). In our regions, a statistical validation is not feasible because of the paucity of data. Moreover, the classical tests do not clearly suggest which one among different proposed models must be favoured. In this paper, we adopt a model of interoccurrence times able to interpret the three different hypotheses, ranging from exponential to Weibull distributions, in a scenario of increasing degree of predictability. In order to judge which one of these hypotheses is favoured, we adopt, instead of the classical tests, a more selective indicator measuring the error in respect to the chosen panorama of possible truths. The earthquake prediction is here simply defined and calculated through the conditional probability of occurrence depending on the elapsed time t₀ since the last earthquake. Short-term and medium-term predictions are performed for all the Italian seismic zones on the basis of datasets built in the context of the National Projects INGV-DPC 2004–2006, in the frame of which this research was developed. The mathematical model of interoccurrence times (mixture of exponential and Weibull distributions) is justified in its analytical structure. A dimensionless procedure is used in order to reduce the number of parameters and to make comparisons easier. Three different procedures are taken into consideration for the estimation of the parameter values; in most of the cases, they give comparable results. The degree of credibility of the proposed methods is evaluated. Their robustness as well as their sensitivity are discussed. The comparison of the probability of occurrence of a Maw >5.3 event in the next 5 and 30 years from January 1, 2003, conditional to the time elapsed since the last event, shows that the relative ranking of impending rupture in 5 years is roughly maintained in a 30-year perspective with higher probabilities and large fluctuations between sources belonging to the same macro region.     

2020年10月22日四川北川MS4.7地震强震动特征分析

2020年10月22日11时03分37秒四川省绵阳市北川县发生MS4.7地震,四川强震动台网与预警烈度速报台网在震区建成较密集的台站,获取了532组三分量加速度记录,有助于开展区域地震动衰减和地震动特征研究.本文对强震记录进行常规处理后计算出强震动记录的相关参数,利用克里金插值方法得到地震动峰值加速度PGA和峰值速度PGV的空间分布图,长轴呈北西—南东方向.分析强震动记录PGA、PGV随距离的衰减规律,与常用衰减关系预测值进行对比,此次地震PGA的衰减特性与俞言祥和汪素云(2006)提出的中国西部地区水平向基岩加速度衰减关系有较好的一致性.北川MS4.7地震获得的密集强震动记录为建立区域衰减关系,以及开展基于经验格林函数方法(EGFM)再现大震强地震动场展布等研究提供了重要的数据支撑.     

Great earthquake at 7.3 ka inferred from tsunami deposits in the Sukumo Bay area,Southwestern Japan

Tsunami deposits in Kyushu Island, Southwestern Japan, have been attributed to the 7.3 ka Kikai caldera eruption, but their origin has not been confirmed. We analyzed an 83-cm-thick Holocene event deposit in the SKM core, obtained from incised valley fill in the coastal lowlands near Sukumo Bay, Southwestern Shikoku Island. We confirmed that the event deposit contains K-Ah volcanic ash from the 7.3 ka eruption. The base of the event deposit erodes the underlying inner-bay mud, and the deposit contains material from outside the local terrestrial and marine environment, including angular quartz porphyry from a small inland exposure, oyster shell debris, and a coral fragment. Benthic foraminifers and ostracods in the deposit indicate various habitats, some of which are outside Sukumo Bay. The sand matrix contains low-silica volcanic glass from the late stage of the Kikai caldera eruption. We also documented the same glass in an event deposit in the MIK1 core, from the incised Oyodo River valley in the Miyazaki Plain on Southeastern Kyushu. These two 7.3 ka tsunami deposits join other documented examples that are widely distributed in Southwestern Japan including the Bungo Channel and Beppu Bay in Eastern Kyushu, Tachibana Bay in Western Kyushu, and Zasa Pond on the Kii Peninsula as well as around the caldera itself. The tsunami deposits near the caldera have been divided into older and younger 7.3 ka tsunami deposits, the younger ones matching the set of widespread deposits. We attribute the younger 7.3 ka tsunami deposits to a large tsunami generated by a great interplate earthquake in the Northern part of the Ryukyu Trench and (or) the Western Nankai Trough just after the late stage of the Kikai caldera eruption and the older 7.3 ka tsunami deposits to a small tsunami generated by an interplate earthquake or Kikai caldera eruption.     

Fracture characteristics of the 1997 Jiashi, Xinjiang, China, earthquake swarm inferred from source spectra

Broadband P and S waves source spectra of 12 MS5.0 earthquakes of the 1997 Jiashi, Xinjiang, China, earthquake swarm recorded at 13 GDSN stations have been analyzed. Rupture size and static stress drop of these earthquakes have been estimated through measuring the corner frequency of the source spectra. Direction of rupture propagation of the earthquake faulting has also been inferred from the azimuthal variation of the corner frequency. The main results are as follows: ①The rupture size of MS6.0 strong earthquakes is in the range of 10～20 km, while that of MS=5.0～5.5 earthquakes is 6～10 km.② The static stress drop of the swarm earthquakes is rather low, being of the order of 0.1 MPa. This implies that the deformation release rate in the source region may be low. ③ Stress drop of the earthquakes appears to be proportional to their seismic moment, and also to be dependent on their focal mechanism. The stress drop of normal faulting earthquakes is usually lower than that of strike-slip type earthquakes. ④ For each MS6.0 earthquake there exists an apparent azimuthal variation of the corner frequencies. Azimuthally variation pattern of corner frequencies of different earthquakes shows that the source rupture pattern of the Jiashi earthquake swarm is complex and no uniform rupture expanding direction exists.     

由震源谱推断1997年新疆伽师强震群破裂特性

分析了1997年伽师强震群中12次MS5.0地震的Ｐ波和Ｓ波位移谱,用遗传算法和目测法联合计算了这些地震Ｐ,Ｓ波震源谱的拐角频率,并由此推算出这些地震的震源尺度和静态应力降（以下简称应力降）．根据观测震源谱拐角频率的方位性变化特征,对相关地震的破裂方向进行了估算．主要结果如下:① 伽师强震群中６级强震的震源破裂尺度在10～16ｋｍ,５级地震的破裂尺度在6～10ｋｍ,个别５级地震的破裂尺度超过10ｋｍ;② 伽师强震群中,中、强地震的应力降均在0.1ＭＰａ左右,这种显著偏低的应力降特性表明,该强震群震源区积累的应力释放缓慢,这意味着伽师强震群活动还将持续相当长时间;③ 应力降大小与相应地震震级大小存在统计正相关关系．同时还可清楚地看出,正断层型地震的应力降整体低于走滑型地震的应力降;④伽师强震群中６ 级强震的震源谱拐角频率有一定的方向性,震源谱拐角频率方位分布图显示出这组强震的破裂机制较为复杂,无一致的破裂方向．