Parameters of Tsunami Source Versus Earthquake Magnitude

This study employs an empirical-analytical approach in combination with Monte-Carlo method to establish relationships between earthquake moment magnitude and upper limits of tsunami source parameters: double-amplitude of vertical bottom deformation, displaced water volume, and potential energy of initial elevation. The approach is based on the Okada solution for a finite rectangular fault and empirical scaling laws for earthquake sources. Results are compared to empirical and theoretical relationships published previously. Parameters of some recent tsunami sources, those for which USGS provides slip distribution data, are considered in light of the established relationships.     

远场地震电磁辐射强度与孕震震级关系

基于垂直电偶极子激发源模型,首先研究了远场地震电磁场辐射场及其分布特征,其次分析了远场地震电磁辐射强度与孕震震级之间的关系.结果表明垂直电偶极子激发源可以解释孕震区激发电磁辐射现象,是激发地震电磁辐射最有可能的激发源之一;远场地震电磁辐射强度的对数与地震震级呈线性关系,直线的斜率与孕震区裂隙特性分布及其电磁特性有关.研究结果对进一步认识地震电磁辐射现象和拓宽地震监测预报思路均有积极意义.     

利用烈度数据点估算中强历史地震震级与震中的方法研究

本文基于地震烈度数据点,用椭圆参数方程对烈度点空间分布进行最小二乘拟合,得到各烈度区椭圆烈度估计线,然后对烈度估计线进行统计回归,得到一组适用于青、甘、川、滇4省中强地震的椭圆烈度分布模型I=f（M,R）.基于所建烈度分布模型,联立考虑中心点和方向性的椭圆数学方程,代入全部烈度数据点估算地震震级和宏观震中.本文的试算震例表明了此方法的可行性,并进行了不确定性分析.利用该方法对明清时期4次烈度点较少的中强历史地震参数进行了估算,所得结果表明了此方法对历史地震有效.     

四川地震震级与烈度区面积和震源深度的关系

本文按照地震地质分区,同时考虑7级以上地震的空间分布特征,分四川为鲜水河、松潘—龙门山及马边地震带等几个区带,对四川1900年以来既有仪器测定参数、又有宏观烈度考查的33个地震资料,应用多元线性回归分析方法,统计拟合出了震级与不同烈度区面积、烈度值、震源深度之间的多组关系式,其一般形式为:■相关系数R大都在0.90左右,用部分地震检验的结果表明,震级误差一般小于0.3,震源深度误差一般小于5公里,基本消除了离散情况,笔者用求得的关系式对1900年以前的有三条以上等震线的部分历史地震的震级、震源深度进行了标定,结果也比较令人满意。     

基于历史地震数据的地震烈度与人员死亡率相关性研究

基于历史震例数据,通过对不同烈度的人员死亡率的分析发现,人员死亡率与地震烈度存在着正相关关系,即人员死亡率随烈度的增加而增大,不同烈度的人员死亡率有各自的分布区间范围,而通过对相邻烈度的人员死亡率的分析发现,相邻烈度人员死亡率存在着1个数量级的倍数关系,一般在10倍左右,在低烈度区域,倍数关系集中于偏向大于10倍的区间范围;在高烈度区域,倍数关系集中于小于10倍的区间范围。在此结果上,基于地震烈度、震级等因素构建了烈度人员死亡率模型,决定系数R²值为0.8667,拟合结果相对较好,能够为后续基于分烈度人员死亡率的地震人员死亡评估模型提供参考。     

Intensity Retrieved from Public Monitoring Video——Cases from the May 12,2008,Wenchuan Earthquake

Macroseismic investigation is presented with a new potential means of investigation in the ＂digital age＂. In this paper we studied several cases in the May 12, 2008, Wenchuan Ms8.0 earthquake with a view to exploring the feasibility of retrieving intensity, or even strong ground motion parameters, from public monitoring video which was originally deployed for security purposes. We used public monitoring video records from 44 sites across the meizoseismal region and its surroundings to estimate the intensity. By using the video from a bank in Dujiangyan, Sichuan Province and making the simple assumption that ground vibration is in the form of a propagating harmonic wave, we estimated that the local ground motion acceleration was 0.14g.     

Tsunami Source of the 2010 Mentawai, Indonesia Earthquake Inferred from Tsunami Field Survey and Waveform Modeling

The 2010 Mentawai earthquake (magnitude 7.7) generated a destructive tsunami that caused more than 500 casualties in the Mentawai Islands, west of Sumatra, Indonesia. Seismological analyses indicate that this earthquake was an unusual “tsunami earthquake,” which produces much larger tsunamis than expected from the seismic magnitude. We carried out a field survey to measure tsunami heights and inundation distances, an inversion of tsunami waveforms to estimate the slip distribution on the fault, and inundation modeling to compare the measured and simulated tsunami heights. The measured tsunami heights at eight locations on the west coasts of North and South Pagai Island ranged from 2.5 to 9.3 m, but were mostly in the 4–7 m range. At three villages, the tsunami inundation extended more than 300 m. Interviews of local residents indicated that the earthquake ground shaking was less intense than during previous large earthquakes and did not cause any damage. Inversion of tsunami waveforms recorded at nine coastal tide gauges, a nearby GPS buoy, and a DART station indicated a large slip (maximum 6.1 m) on a shallower part of the fault near the trench axis, a distribution similar to other tsunami earthquakes. The total seismic moment estimated from tsunami waveform inversion was 1.0 × 10²¹ Nm, which corresponded to M_w 7.9. Computed coastal tsunami heights from this tsunami source model using linear equations are similar to the measured tsunami heights. The inundation heights computed by using detailed bathymetry and topography data and nonlinear equations including inundation were smaller than the measured ones. This may have been partly due to the limited resolution and accuracy of publically available bathymetry and topography data. One-dimensional run-up computations using our surveyed topography profiles showed that the computed heights were roughly similar to the measured ones.     

南海北缘海啸波高对潜在海啸源震级偏差的敏感性

本文假设马尼拉海沟北段为潜在海啸源,基于中国地震台网对马尼拉海沟地区震级测定偏差,采用COMCOT(comell Multi-grid Coupled Tsunami Model)海啸数值模型,模拟南海海啸波传播.选取南海北缘3个特定地点,其中两个位于华南近海区域,另一个位于台湾岛南端近海区域,此外还在临近马尼拉海沟北段的深海地区选取了1个特定地点.分析这些特定地点最大海啸波以及最大海啸波到时对于震级测定偏差的敏感性.结果表明:马尼拉海沟北段地震如触发海啸,华南近海区域以及台湾岛南部近海区域最大海啸波振幅对震级偏差敏感,但最大海啸波振幅到时对于震级测定偏差不敏感;振幅最大的海啸波,二十几分钟即可波及台湾岛南端近岸区域,大约1小时后波及大陆华南近海北部区域.     

关于1856年黔江-咸丰地震震中烈度与震级的探讨

1856年黔江-咸丰地震在不同版本的目录中,震级和震中烈度均存在不小的差异。本文通过对1856年黔江-咸丰地震史料的收集,查阅更大范围的史志、地方志,详细分析了史料中此次地震造成的震害情况。对黔江-咸丰地震进行实地考查,并通过对西南山区近年来几次地震的对比分析,从震害角度、地震烈度图不同烈度面积的角度,进一步讨论震中烈度的可能性。采用多种方法对地震震级进行估算。但黔江-咸丰地震的震中烈度和震级确定依然是个难题,还需进行进一步研究。     

印度洋地震海啸引发的思考

从2004年12月印度洋发生地震海啸造成的灾难中得到启示，从完善我国的应急救援、防灾救灾体系、法律法规体系，预测系统、快速反应预警系统以及启动应急预案、组建应急指挥体系、组织应急救援队伍、实施应急救援等方面，阐述了建立应急救援指挥机构和各种应急队伍的必要性，绘制了灾害处理应急流程图，分析了应急队伍所应承担的任务，提出在我国应尽快建立地震海啸预警机制的建议。     

中国西部地区地震烈度衰减关系

本文收集了1991年之后我国西部地区71个地震的烈度等震线资料,并以此对1918—1989年间176次地震的烈度资料进行了补充,采用长、短轴椭圆模型重新拟合了西部地区的分区地震烈度衰减关系。结果表明,新疆区和川藏区2个统计单元内地震烈度衰减关系有显著不同,应作为不同的分区对待。同时与其他研究者给出的该地区地震烈度衰减关系的对比结果显示,本文结果合理可靠,体现了近、远场地震烈度分布的特点,适合于该地区工程地震研究和应用。     

5·12地震震级-累积频度关系的研究

古登堡-李希特提出的震级-累积频度关系,是地震分析统计预报和工程场地地震安全性评价的重要基础理论之一,其合理性被大量实际地震和岩石力学试验所证实。在汶川大地震一周年时,震级-累积频度关系仍是值得深入研究的根本性地震活动问题,本文对其进行了较为深入的研讨。     

历史地震资料与地震的中长期预测

故障诊断模型的基本内容是根据动态系统的外部特征来判断系统内部是否发生故障及确定故障发生的部位、时间和大小。由于故障诊断技术在监测及诊断故障的思路上和地震预报有很多相似之处，因此，把故障诊断技术应用于地震预报是可行的。由于该模型与其他数学模型一样，需要有较多的学习过程，所以历史地震资料在该模型中有重要作用。而其中的鲁棒性故障诊断模型在抑制各子模型的个性，凸现其共性方面有其特有的性质，所以把它作为一种综合模型，能抑制各子模型的个性，突出在地震预测方面的共性，从而提高预测的精确性。本文根据一个实例，说明了这种综合性模型的可行性。在文章的最后，由信息量的分析，说明了模型的鲁棒性特征。     

The 15 August 2007 Peru Earthquake and Tsunami: Influence of the Source Characteristics on the Tsunami Heights

The tsunami caused by the 2007 Peru earthquake (M_w 8.0) provoked less damage than by the seismic shaking itself (numerous casualties due to the earthquake in the vicinity of Pisco). However, it propagated across the Pacific Ocean and small waves were observed on one tide gauge in Taiohae Bay (Nuku Hiva, Marquesas, French Polynesia). We invert seismological data to recover the rupture pattern in two steps. The first step uses surface waves to find a solution for the moment tensor, and the second step uses body waves to compute the slip distribution in the source area. We find the slip distribution to consist of two main slip patches in the source area. The inversion of surface waves yields a scalar moment of 8.9 10²⁰ Nm, and body-wave inversion gives 1.4 10²¹ Nm. The inversion of tsunami data recorded on a single deep ocean sensor also can be used to compute a fault slip pattern (yielding a scalar moment of 1.1 10²¹ Nm). We then use these different sources to model the tsunami propagation across the Pacific Ocean, especially towards Nuku Hiva. While the source model taken from the body-wave inversion yields computed tsunami waves systematically too low with respect to observations (on the central Pacific Ocean DART buoy as on the Polynesian tide gauge), the source model established from the surface-wave inversion is more efficient to fit the observations, confirming that the tsunami is sensitive to the low frequency component of the source. Finally we also discuss the modeling of the late tsunami arrivals in Taiohae Bay using several friction coefficients for the sea bottom.     

Re-examination of the Source Mechanism of the 1998 Papua New Guinea Earthquake and Tsunami

— Simulation of tsunami propagation and runup of the 1998 Papua New Guinea (PNG) earthquake tsunami using the detailed bathymetry measured by JAMSTEC and adding bathymetric data at depths less than 60 m is carried out, reproducing the tsunami energy focus into Warapu and Arop along the Sissano Lagoon. However, the computed runup heights in the lagoon are still lower than those measured. Even if the error in estimating the fault parameters is taken into consideration, computational results are similar. Analysis by the wave ray method using several scenarios of the source size of the tsunami and location by the wave ray method suggests that a source characterized by small size in water 1,000-m deep approximately 25 km offshore the lagoon, best fits the arrival determined from the interviews with eyewitnesses. A two-layer numerical model simulating the interaction of the tsunami with a landslide is employed to study the behavior of a landslide-generated tsunami with different size sand depths of the initial slide just outside the lagoon. A landslide model with a volume of 4–8 × 10⁹ m³ is selected as the best in order to reproduce the distribution of the measured tsunami runup in the lagoon. The simulation of a tsunami generated in two stages, fault and landslide, could show good agreement with the runup heights and distribution of the arrival time, but a time gap of around 10 minutes remains, suggesting that a tsunami generated by the mainshock at 6:49 PM local time is too small for people to notice, and the following tsunami triggered by landslide or mass movement near the lagoon about ten minutes after the mainshock attacked the coast and caused the huge damage.     

古登堡与里克特震级-频度关系的系数研究

众所周知,大小地震之间的著名统计关系 logn(Mi)=a-bMi logN(Mi)=A-BMi 式中i=0,1,2,…,k,其中,M.为最小统计震级档,Mk为最大震级档.古登堡在统计时,是将地震按震级分档的,即n是震级为Mi+△M的地震次数,为(1)式.里克特采用按震级由强至弱累计,即N是震级大于等于Mi的地震次数,为(2)式.     

Discussion on the Relationship between Different Earthquake Magnitude Scales and the Effect of Seismic Station Sites on Magnitude Estimation

Based on the earthquake catalog reported by the Chinese digital seismic network in recent years, we select the earthquakes with both surface wave magnitude and local magnitude and fit them into a relationship between the two magnitudes. The systematic difference is found from the formula which has been used for 30 years. Because of a large dynamic range and wide frequency range of the current digital observation system, in addition to a larger number of stations and earthquakes being used compared to before, the relation obtained in this paper seems more reliable. Our calculation shows that there is no significant difference before and after magnitude conversion so we suggest the abandonment of magnitude conversion. The site response of a station consists of amplification at different frequencies. The amplification is equal to about 1 and changes little with frequency at stations located on basement rock, and it is greater than 1 at low frequency ranges and less than 1 at high frequency ranges at stations located on sediment layers. The difference between magnitudes from single station located on sediment layer and the average magnitude from the whole network increases from negative to positive with period. It seems that there is no fixed station correction factor and the station correction method does not work to improve the accuracy and magnitude estimates.