The moment magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> and the energy magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>e</Subscript>: common roots and differences

Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M _w and energy magnitude M _e is outlined and critically discussed. The formulas for M _w and M _e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M ₀ and radiated seismic energy E _S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M _w and M _e are linked via the parameter Θ = log(E _S/M ₀), and the formula for M _e can be written as M _e = M _w + (Θ + 4.7)/1.5. This relationship directly links M _e with M _w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M _w and M _e can significantly differ. In fact, Θ is assumed to be constant when calculating M _w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V _R and seismic wave radiation efficiency η _R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M _e may sometimes differ by more than one magnitude unit from M _w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M _w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M _e is more suitable than M _w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M _w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E _S and M _e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E _S in characterizing the seismic source, the use of M _e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E _S estimations in order to allow M _e to be extensively used as an important complement to M _w in common seismological practice and its applications.     

Geotechnical engineering aspects of the catastrophic Tohoku (Japan) earthquake of March 11, 2011: the review of first results

Geotechnical engineering aspects of the catastrophic earthquake, which occurred in Japan on March, 11, 2011 and called Tohoku earthquake are discussed. A review is presented of the first results obtained mainly by Japanese scientists based on records of seismic networks of Japan K-NET, Kik-net and on GPS data. The basic concepts of seismic zoning in Japan and the location of the Tohoku-oki earthquake on the seismic zoning maps are described, as well as models of the source process obtained by various authors based on teleseismic data, strong motion data, GPS data, and tsunami observations. The recorded peak accelerations and velocities and their correspondence to the current empirical attenuation curves are discussed. The records of the Tohoku earthquake made by Japanese seismic networks K-NET, Kik-net and some others represent unique seismological material and the most complete seismic database (including vertical array records) in the near-source zone of a strongest earthquake with magnitude M _w = 9. These data will be studied by seismologists all over the world for many years and, probably, they will answer many questions of geotechnical seismic engineering.     

Critical evaluation of strong motion in Kocaeli and Düzce (Turkey) earthquakes

Turkey was struck by two major events on August 17th and November 12th, 1999. Named Kocaeli (M_w=7.4) and Düzce (M_w=7.2) earthquakes, respectively, the two earthquakes provided the most extensive strong ground motion data set ever recorded in Turkey. The strong motion stations operated by the General Directorate of Disaster Affairs, the Kandilli Observatory and Earthquake Research Institute of Bogazici University and Istanbul Technical University have produced at least 27 strong motion records for the Kocaeli earthquake within 200 km of the fault. Kocaeli earthquake has generated six motions within 20 km of the fault adding significantly to the near-field database of ground motions for M_w>=7.0 strike–slip earthquakes. The paper discusses available strong motion data, studies their attenuation characteristics, analyses time domain, as well as spectral properties such as spectral accelerations with special emphasis on fault normal and fault parallel components and the elastic attenuation parameter, kappa. A simulation of the Kocaeli earthquake using code FINSIM is also presented.     

2017年九寨沟7.0级地震前不同级别地震空区的演化过程分析

结合传统地震学方法与数字地震学方法,回顾总结了日常分析预报中针对九寨沟7.0级地震开展的相关工作,梳理总结了地震前提出的甘东南地震异常信息,进一步讨论了九寨沟地震前中期、短期及临震异常特征。通过分析九寨沟7.0级地震之前区域范围内不同震级活动图像,发现甘青川交界地区5级空区在震前向震中一侧收缩、4级地震超长平静458天、3级空区在震前4天被打破的现象,同时发现,震中附近区域震源机制一致性较高,反映了震源区的高应力水平。     

近断层脉冲型地震动作用下隔震结构地震反应分析

隔震结构在远震场地减震效果良好，但是近断层地震动的明显的长周期速度和位移脉冲运动可能对隔震建筑等长周期结构的抗震性能和设计带来不利影响，需要深入探讨。本文首先讨论近断层地震动的长周期脉冲运动特征，然后以台湾集集地震8条典型近震记录和其它4条常用近震记录以及4条远震记录作为地震动输入，对两幢安装铅芯橡胶隔震支座的钢筋混凝土框架隔震结构进行非线性地震反应时程分析，通过比较探讨了算例计算结果，定量说明隔震结构的近震脉冲效应显著，是隔震设计不容忽视的问题。     

Pulse-like versus non-pulse-like ground motion records: Spectral shape comparisons and record selection strategies

Pulse-like records are well recognized for their potential to impose higher demands on structures when compared with ordinary records. The increased severity of the structural response usually caused by pulse-like records is commonly attributed to the spectral increment around the pulse period. By comparing the building response to sets of spectrally equivalent pulse-like and ordinary records, we show that there are characteristics of pulse-like records beyond the shape of the acceleration response spectrum that affect the results of nonlinear dynamic analysis. Nevertheless, spectral shape together with the ratio of pulse period to the first-mode structural period, T_p/T₁, are confirmed as “sufficient” predictors for deformation and acceleration response metrics in a building, conditioned on the seismic intensity. Furthermore, the average spectral acceleration over a period range, AvgSA, is shown to incorporate to a good proxy for spectral shape, and together with T_p/T₁, form an efficient and sufficient intensity measure for response prediction to pulse-like ground motions. Following this latter route, we propose a record selection scheme that maintains the consistency of T_p with the hazard of the site but uses AvgSA to account for the response sensitivity to spectral shape.     

Deterministic Seismic Scenarios for Imphal City

In this article, the spatial variation of ground motion in Imphal City has been estimated by the finite-fault seismological model coupled with site response analysis. The important seismic sources around Imphal City have been identified from the fault map and past seismicity data. The rock level acceleration time histories at Imphal City for the 1869 Cachar (M_w 7.5) earthquake and a hypothetical M_w 8.1 event in the Indo-Burma subduction zone have been estimated by a stochastic finite-fault model. Soil investigation data of 122 boreholes have been collected from several construction projects in Imphal City. Site response analysis has been carried out and the surface level ground motion has been determined for Imphal City for these two earthquake events. The results are presented in the form of peak ground acceleration (PGA) contour map. From the present study it has been ascertained that the maximum amplification for PGA over Imphal City is as high as 2.5. The obtained contour maps can serve as guidelines for identifying vulnerable areas and disaster mitigation in Imphal City.     

Response spectrum-oriented pulse identification and magnitude scaling of forward directivity pulses in near-fault ground motions

This paper presents a novel approach to identify the pulse-like motions in earthquake recordings that dominate the maximum structural responses over a wide period range. The identification method is based on the congruence relationship between the response spectrum and the dimensionless П-response spectrum established in this study through straightforward dimensional arguments of linear and bilinear SDOF oscillators subject to pulse-like ground motions. By evaluating the geometric match and dislocations of the П-response spectrum of a given waveform with the dimensional response spectrum in bi-logarithm plotting, one can identify the simple pulses and their parameters that match simultaneously the kinematic characteristics and the response spectrum of earthquake recordings that exhibit pulse-like features. The developed pulse identification method has been implemented in a computer program and applied successfully to detect the pulse-like motions in the PEER NGA strong motion database. Both velocity and acceleration pulses potentially due to forward directivity effects in near fault regions are identified. The identified velocity pulses show strong correlation with the seismological parameters. They are subsequently used in regression analysis to derive the empirical scaling laws that relate the directivity pulse parameters to the earthquake magnitude and rupture distance. The study confirms some magnitude scaling laws in literature and demonstrates the accuracy and efficiency of the proposed pulse identification method.     

Unified Local Magnitude Scale for Earthquakes of South Balkan Area

—A homogeneous earthquake catalog spanning 1964–1995 for the southern Balkan area is presented, by expressing the size of the earthquakes in a unified local magnitude scale. The strategy followed to produce this catalog is also presented. Local magnitudes calculated by six Balkan seismological centers (ATH, THE, ISK, TIR, TTG, SKO) have been used in order to obtain relations between the estimated ML values from the Greek seismological centers (ATH and THE) and the remaining four Balkan stations. Since it was found that local magnitudes estimated by ATH and THE are almost identical, they have been used as one data set so as to correlate with the data of each one of the remaining four seismological stations. Based on the proposed relations, a unified local magnitude, MLGR, is given for each earthquake of the regional catalog. A published scaling relation between the ML values from ATH and THE networks and the corresponding seismic moment magnitude, Mw, was used in combination with the above relations, in order to enable the conversion of any ML value from any station into Mw. The catalog completeness has been checked and the b-value has been calculated for the complete data sample.     

Pulse-like ground motion observed during the 6 February 2023 MW7.8 Pazarcık Earthquake (Kahramanmaraş, SE Türkiye)

In this study, we analyzed 100 three-component strong ground motion records observed within 200 km of the causative fault of the 6 February 2023 M_W7.8 Pazarcık (Kahramanmaraş) Earthquake in SE Türkiye. The wavelet method was utilized to identify and analyze the characteristics of pulse-like ground motions in the near-fault region, while considering the uncertainty of the pulse orientation during the analysis. Our investigation focused on the effects of the focal mechanism and rupture process on the spatial distribution, pulse orientation, and maximum pulse direction of the observed pulse-like ground motion. We also analyzed the amplitude and period of the observed ground pulses and the effect of long-period amplification on the ground motion response spectra. Our results indicated the following: (1) A total of 21 typical ground velocity pulses were observed during this earthquake, exhibiting complex characteristics due to the influence of the strike-slip mechanism and rupture directivity. Most ground pulses (17 out of 21) were recorded within 20 km of the fault, in a wide range of orientations, including normal and parallel to the fault direction. The waveforms exhibited unidirectional features, indicating the effects of left-lateral fault slip. Distinct pulses observed more than 20 km from the fault were mainly oriented normal to the fault. The waveforms were bidirectional with double- or multi-round trips as a result of rupture directivity. (2) The amplitudes of the observed pulses ranged from 30.5 to 220.0 cm/s, with the largest peak velocity of 220.0 cm/s observed at Station 3138. The pulse periods ranged from 2.3 to 14.5 s, with the longest pulse period of 14.5 s observed at Station 3116. The amplitude and period of the pulses observed during this earthquake were comparable to those of similar-magnitude global earthquakes. The amplitude of the pulses decreased significantly with increasing fault distance, whereas the pulse period was not significantly affected by the fault distance. (3) Compared with non-pulse records, the velocity pulse records had a pronounced amplification effect on the acceleration response spectra near the pulse period, with factors ranging from 2.1 to 5.8. The larger velocity pulses also significantly amplified the velocity response spectra, particularly over the long periods. This significant amplification effect of the pulses on the response spectra leads to empirical models underestimating the long-period earthquake ground motion.     

近断层速度脉冲地震动的三维有限差分模拟

根据台湾西部地质地貌特征和1999年集集M_W7.6地震的研究成果,建立三维速度结构模型和震源模型,并采用三维有限差分法对双冬断层可能产生的近断层脉冲型地震动进行数值模拟。结果表明,方向性效应引起的双向速度脉冲集中在垂直于断层滑动分量的方向上,而滑冲效应引起的单向速度脉冲则集中在平行于断层滑动分量的方向上。受方向性效应和上盘效应的共同调制,近断层脉冲型地震动反映出不对称带状分布的特征,速度脉冲主要分布在距离断层面约10 km的范围内。凹凸体的特性影响着地震动的时空分布,由地震波场显示南投和台中处于强地震动危险区。近场脉冲型地震动的研究对分析速度脉冲形成机理以及地震危险性有一定的参考意义。     

Empirical Global Relations Converting M S and m b to Moment Magnitude

The existence of several magnitude scales used by seismological centers all over the world and the compilation of earthquake catalogs by many authors have rendered globally valid relations connecting magnitude scales a necessity. This would allow the creation of a homogeneous global earthquake catalog, a useful tool for earthquake research. Of special interest is the definition of global relations converting different magnitude scales to the most reliable and useful scale of magnitude, the moment magnitude, M _W. In order to accomplish this, a very large sample of data from international seismological sources (ISC, NEIC, HRVD, etc.) has been collected and processed. The magnitude scales tested against M _W are the surface wave magnitude, M _S, the body wave magnitude, m _b, and the local magnitude, M _L. The moment magnitudes adopted have been taken from the CMT solutions of HRVD and USGS. The data set used in this study contains 20,407 earthquakes, which occurred all over the world during the time period 1.1.1976–31.5.2003, for which moment magnitudes are available. It is shown that well-defined relations hold between M _W and m _b and M _S and that these relations can be reliably used for compiling homogeneous, with respect to magnitude, earthquake catalogs.     

Validation of ground‐motion simulations for historical events using MDoF systems

The study presented in this paper addresses the issue of engineering validation of Graves and Pitarka's (2010) hybrid broadband ground motion simulation methodology with respect to some well‐recorded historical events and considering the response of multiple degrees of freedom (MDoF) systems. Herein, validation encompasses detailed assessment of how similar is, for a given event, the seismic response due to comparable hybrid broadband simulated records and real records. In the first part of this study, in order to investigate the dynamic response of a wide range of buildings, MDoF structures are modeled as elastic continuum systems consisting of a combination of a flexural cantilever beam coupled with a shear cantilever beam. A number of such continuum systems are selected including the following: (1) 16 oscillation periods between 0.1 and 6 s; (2) three shear to flexural deformation ratios to represent respectively shear‐wall structures, dual systems, and moment‐resisting frames; and (3) two stiffness distributions along the height of the systems, that is, uniform and linear. Demand spectra in terms of generalized maximum interstory drift ratio (IDR) and peak floor acceleration (PFA) are derived using simulations and actual recordings for four historical earthquakes, namely, the 1979 M_w 6.5 Imperial Valley earthquake, 1989 M_w 6.8 Loma Prieta earthquake, 1992 M_w 7.2 Landers earthquake, and 1994 M_w 6.7 Northridge earthquake. In the second part, for two nonlinear case study structures, the IDR and PFA distributions over the height and their statistics, are obtained and compared for both recorded and simulated time histories. These structures are steel moment frames designed for high seismic hazard, 20‐story high‐rise and 6‐story low‐rise buildings. The results from this study highlight the similarities and differences between simulated and real records in terms of median and intra‐event standard deviation of logs of seismic demands for MDoF building systems. This general agreement, in a broad range of moderate and long periods, may provide confidence in the use of the simulation methodology for engineering applications, whereas the discrepancies, statistically significant only at short periods, may help in addressing improvements in generation of synthetic records. Copyright © 2013 John Wiley & Sons, Ltd.     

Evidence for a Tang-Song Dynasty great earthquake along the Longmen Shan Thrust Belt prior to the 2008 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.9 Wenchuan earthquake,China

The magnitude (M _w) 7.9 Wenchuan earthquake occurred on 12 May 2008 in the Longmen Shan region of China, the transition zone between the Tibetan Plateau and the Sichuan Basin, resulting in widespread damage throughout central and western China. The steep, high-relief eastern margin of the Tibetan Plateau has undergone rapid Cenozoic uplift and denudation accompanied by folding and thrusting, yet no large thrust earthquakes are known prior to the 2008 M _w 7.9 Wenchuan earthquake. Field and excavation investigations reveal that a great historical earthquake occurred in the Sichuan region that ruptured a >200-km-long thrust fault within the Longmen Shan Thrust Belt, China, which also triggered the 2008 M _w 7.9 Wenchuan earthquake. The average co-seismic slip amount produced by this historical earthquake is estimated to be 2–3 m, comparable with that caused by the 2008 Wenchuan earthquake. Paleoseismic and archaeological evidence and radiocarbon dating results show that the penultimate great earthquake occurred in the Sichuan region during the late Tang-Song Dynasty, between AD 800 and 1000, suggesting a recurrence interval of ~1,000–1,200 years for Wenchuan-magnitude (M = ~8) earthquakes in the late Holocene within the Longmen Shan Thrust Belt. This finding is in contrast with previous estimates of 2,000–10,000 years for the recurrence interval of large earthquakes within the Longmen Shan Thrust Belt, as obtained from long-term slip rates based on the Global Positioning System and geological data, thereby necessitating substantial modifications to existing seismic-hazard models for the densely populated region at the eastern marginal zone of the Tibetan Plateau.     

Probabilistic Assessment of Earthquake Recurrence in Northeast India and Adjoining Regions

Northeast India and adjoining regions (20°–32° N and 87°–100° E) are highly vulnerable to earthquake hazard in the Indian sub-continent, which fall under seismic zones V, IV and III in the seismic zoning map of India with magnitudes M exceeding 8, 7 and 6, respectively. It has experienced two devastating earthquakes, namely, the Shillong Plateau earthquake of June 12, 1897 (M _w 8.1) and the Assam earthquake of August 15, 1950 (M _w 8.5) that caused huge loss of lives and property in the Indian sub-continent. In the present study, the probabilities of the occurrences of earthquakes with magnitude M ≥ 7.0 during a specified interval of time has been estimated on the basis of three probabilistic models, namely, Weibull, Gamma and Lognormal, with the help of the earthquake catalogue spanning the period 1846 to 1995. The method of maximum likelihood has been used to estimate the earthquake hazard parameters. The logarithmic probability of likelihood function (ln L) is estimated and used to compare the suitability of models and it was found that the Gamma model fits best with the actual data. The sample mean interval of occurrence of such earthquakes is estimated as 7.82 years in the northeast India region and the expected mean values for Weibull, Gamma and Lognormal distributions are estimated as 7.837, 7.820 and 8.269 years, respectively. The estimated cumulative probability for an earthquake M ≥ 7.0 reaches 0.8 after about 15–16 (2010–2011) years and 0.9 after about 18–20 (2013–2015) years from the occurrence of the last earthquake (1995) in the region. The estimated conditional probability also reaches 0.8 to 0.9 after about 13–17 (2008–2012) years in the considered region for an earthquake M ≥ 7.0 when the elapsed time is zero years. However, the conditional probability reaches 0.8 to 0.9 after about 9–13 (2018–2022) years for earthquake M ≥ 7.0 when the elapsed time is 14 years (i.e. 2009).     

The 26 December 2004 Sumatra Tsunami: Analysis of Tide Gauge Data from the World Ocean Part 1. Indian Ocean and South Africa

The M_w = 9.3 megathrust earthquake of December 26, 2004 off the northwest coast of Sumatra in the Indian Ocean generated a catastrophic tsunami that was recorded by a large number of tide gauges throughout the World Ocean. Part 1 of our study of this event examines tide gauge measurements from the Indian Ocean region, at sites located from a few hundred to several thousand kilometers from the source area. Statistical characteristics of the tsunami waves, including wave height, duration, and arrival time, are determined, along with spectral properties of the tsunami records.     

Estimation of strong ground motions from hypothetical earthquakes on the Cascadia subduction zone,Pacific Northwest

Strong ground motions are estimated for the Pacific Northwest assuming that large shallow earthquakes, similar to those experienced in southern Chile, southwestern Japan, and Colombia, may also occur on the Cascadia subduction zone. Fifty-six strong motion recordings for twenty-five subduction earthquakes ofM _s7.0 are used to estimate the response spectra that may result from earthquakesM _w<81/4. Large variations in observed ground motion levels are noted for a given site distance and earthquake magnitude. When compared with motions that have been observed in the western United States, large subduction zone earthquakes produce relatively large ground motions at surprisingly large distances. An earthquake similar to the 22 May 1960 Chilean earthquake (M _w 9.5) is the largest event that is considered to be plausible for the Cascadia subduction zone. This event has a moment which is two orders of magnitude larger than the largest earthquake for which we have strong motion records. The empirical Green's function technique is used to synthesize strong ground motions for such giant earthquakes. Observed teleseismicP-waveforms from giant earthquakes are also modeled using the empirical Green's function technique in order to constrain model parameters. The teleseismic modeling in the period range of 1.0 to 50 sec strongly suggests that fewer Green's functions should be randomly summed than is required to match the long-period moments of giant earthquakes. It appears that a large portion of the moment associated with giant earthquakes occurs at very long periods that are outside the frequency band of interest for strong ground motions. Nevertheless, the occurrence of a giant earthquake in the Pacific Northwest may produce quite strong shaking over a very large region.     

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Source Parameters of the 8 October, 2005 M<Subscript>w</Subscript>7.6 Kashmir Earthquake

During the last six years, the National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5 broadband seismographs and 10 accelerographs in the Kachchh seismic zone, Gujarat, with the prime objective to monitor the continued aftershock activity of the 2001 M_w7.7 Bhuj mainshock. The reliable and accurate broadband data for the M_w 7.6 (8 Oct., 2005) Kashmir earthquake and its aftershocks from this network, as well as from the Hyderabad Geoscope station, enabled us to estimate the group velocity dispersion characteristics and the one-dimensional regional shear-velocity structure of peninsular India. Firstly, we measure Rayleigh- and Love-wave group velocity dispersion curves in the range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8-km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Observation and Modeling of the January 2009 West Papua, Indonesia Tsunami

We modeled a tsunami from the West Papua, Indonesia earthquakes on January 3, 2009 (M _w?=?7.7). After the first earthquake, tsunami alerts were issued in Indonesia and Japan. The tsunami was recorded at many stations located in and around the Pacific Ocean. In particular, at Kushimoto on Kii Peninsula, the maximum amplitude was 43?cm, larger than that at Manokwari on New Guinea Island, near the epicenter. The tsunami was recorded on near-shore wave gauges, offshore GPS sensors and deep-sea bottom pressure sensors. We have collected more than 150 records and used 72 stations?? data with clear tsunami signals for the tsunami source modeling. We assumed two fault models (single fault and five subfaults) which are located to cover the aftershock area. The estimated average slip on the single fault model (80?×?40?km) is 0.64?m, which yields a seismic moment of 1.02?×?10²⁰?Nm (M _w?=?7.3). The observed tsunami waveforms at most stations are well explained by this model.