Higher degree moment tensor inversion of Mani earthquake using far-field broadband recording

Breakthrough point source model, extended earthquake source model is used to calculate more seismic source parameters in this paper. We express seismic source using higher degree moment tensors, to reduce a large number terms originally presenting in higher degree moment tensor representation, Haskell rupture model is used. We inverted the source parameters of Mani earthquake in Tibet using broad-band body wave of 32 stations of Global Seismograph Network (GSN), the results show that it is a strike-slip fault, rupture direction is 75° , rupture duration is 19 s, the fault plan is f =77° , d =88° , l =0° , the auxiliary plane is f =347° , d =90° , l =178° , and the fault dimension is 47 km′ 28 km. These results will give new quantitative data for earth dynamics and have practical meaning for seismic source tomography research.     

Source parameters of bohai earthquake, July 18, 1969 and yongshan earthquake, May 11, 1974 determined by synthetic seismograms of teleseismic P waves

The source parameters of the Bohai Sea earthquake, July 18, 1969 and Yongshan, Yunnan earthquake, May 11, 1974 were determined by full — wave theory synthetic seismograms of teleseismic P waves. P+pP+sP wereform were calculated with WKBJ approximation and real integral paths. One — dimensional unilateral, finite propagation source was also considered. By trail — and — error in comparing the theoretical seismograms with the observational ones of WWSSN stations, the source parameters were obtained as follow: for Bohai earthquake, φ=195°, δ=85°, λ=65°,M _o=0.9×10¹⁹Nm,L=59.9km.V _R=3.5km/s, ∧ _R =160°; for Yongshan earthquake, φ=240°, δ=80°, ∧=150°,M _o=1.3×10¹⁸Nm,L=48.8km,V _R=3km/s, ∧ _R =−10°, where φ is strike, δ dip angle, λ slip angle,M _o seismic moment,L rupture length,V _R rupture propagation speed. As III type fractures the faulting propagated along the fault planes, and ∧ _R is the angle from the strike to the propagation direction. Yongshan earthquake showed complexity in its focal process, having four sub—ruptures during the first 60 seconds. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 1–8, 1991.     

基于恢复地震数据获取震级、震源机制及破裂过程的评价——以2013年四川芦山MW6.6地震为例

区域地震波形对于震源研究非常重要,但限幅问题限制了区域地震台网数据的运用,并影响到震源参数测定的准确度.本文利用恢复后的芦山地震区域地震波形,研究了芦山地震的震级、点源机制解以及破裂过程.基于震中距99~300 km恢复前与恢复后地震数据获取的面波震级分别为7.01与7.06级.分别利用7个震中距150~250 km宽频带台站的恢复前和恢复后的数据反演点源机制解,与参考机制解相比,滑动角偏差自13°减小到了4°.基于7个震中距81~134 km的区域地震波形联合远场数据获得的震源破裂过程结果,其主要参数（如滑动分布、破裂速度等）与强地面运动波形联合远场数据得到的结果具有很好的一致性.研究结果表明,本文所采用的数据恢复方法具有较高的可靠性,有效提高了震源参数测定的准确度.     

Rupture process of the M L=4.1 earthquake in Huailai Basin on July 20, 1995

On July 20, 1995, an earthquake of M _L=4.1 occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, seismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocentres were precisely located. About 2 hours after the occurrence of the main shock, a smaller event of M _L=2.0 took place at 40.323°N, 115.447°E with a focal depth of 5.0 km, which is very close to the main shock. Using the M _L=2.0 earthquake as an empirical Green’s function, a regularization method was applied to retrieve the far-field source-time function (STF) of the main shock. Considering the records of HDSN are the type of velocity, to depress high frequency noise, we removed instrument response from the records of the two events, then integrated them to get displacement seismogram before applying the regularization method. From the 5 field stations, P phases in vertical direction which mostly are about 0.5 s in length were used. The STFs obtained from each seismic phases are in good agreement, showing that the M _L=4.1 earthquake consisted of two events. STFs from each station demonstrate an obvious “seismic Doppler effect”. Assuming the nodal plane striking 37° and dipping 40°, determined by using P wave first motion data and aftershock distribution, is the fault plane, through a trial and error method, the following results were drawn: Both of the events lasted about 0.1 s, the rupture length of the first one is 0.5 km, longer than the second one which is 0.3 km, and the rupture velocity of the first event is 5.0 km/s, larger than that of the second one which is about 3.0 km/s; the second event took place 0.06 s later than the first one; on the fault plane, the first event ruptured in the direction γ=140° measured clockwise from the strike of the fault, while the second event ruptured at γ=80°, the initial point of the second one locates at γ=−100° and 0.52 km from the beginning point of the first one. Using far-field ground displacement spectrum measurement method, the following source parameters about the M _L=4.1 earthquake were also reached: the scalar earthquake moment is 3.3×10¹³ N·m, stress drop 4.6 MPa, rupture radius 0.16 km. Contribution No. 99FE2022, Institute of Geophysics, China Seismological Bureau. This study is supported by the Chinese Joint Seismological Science Foundation (95-07-411).     

Interrelation among several important links in seismic risk analysis

Introduction In the probability analysis method of seismic risk considering time-space inhomogeneity of seismic activity and adopted commonly in China (State Seismological Bureau, 1996) (called in-homogeneous distribution model for short), the division of seismic statistical regions, delimitation of potential seismic sources and estimation of seismicity parameters are the main links that affect significantly the estimation of ground motion parameters of a site. HUANG and WU (2005) studied …     

Historical 1942 Ecuador and 1942 Peru subduction earthquakes and earthquake cycles along Colombia-Ecuador and Peru subduction segments

Two large shallow earthquakes occurred in 1942 along the South American subduction zone inclose proximity to subducting oceanic ridges: The 14 May event occurred near the subducting Carnegie ridge off the coast of Ecuador, and the 24 August event occurred off the coast of southwestern Peru near the southern flank of the subducting Nazca ridge. Source parameters for these for these two historic events have been determined using long-periodP waveforms,P-wave first motions, intensities and local tsunami data.We have analyzed theP waves for these two earthquakes to constrain the focal mechanism, depth, source complexity and seismic moment. Modeling of theP waveform for both events yields a range of acceptable focal mechanisms and depths, all of which are consistent with underthrusting of the Nazca plate beneath the South American plate. The source time function for the 1942 Ecuador event has one simple pulse of moment release with a duration of 22 suconds, suggesting that most of the moment release occurred near the epicenter. The seismic moment determined from theP waves is 6–8×10²⁰N·m, corresponding ot a moment magnitude of 7.8–7.9. The reported location of the maximum intensities (IX) for this event is south of the main shock epicenter. The relocated aftershcks are in an area that is approximately 200 km by 90 km (elongated parallel to the trench) with the majority of aftershocks north of the epicenter. In contrast, the 1942 Peru event has a much longer duration and higher degree of complexity than the Ecuador earthquake, suggesting a heterogeneous rupture. Seismic moment is released in three distinct pulses over approximately 74 seconds; the largest moment release occurs 32 seconds after rupture initiation. the seismic moment as determined from theP waves for the 1942 Peru event is 10–25×10²⁰N·m, corresponding to a moment magnitude of 7.9–8.2. Aftershock locations reported by the ISS occur over a broad area surrounding the main shock. The reported locations of the maximum intensities (IX) are concentrated south of the epicenter, suggesting that at least part of the rupture was to the south.We have also examined great historic earthquakes along the Colombia-Ecuador and Peru segments of the South American subduction zone. We find that the size and rupture length of the underthrusting earthquakes vary between successive earthquake cycles. This suggests that the segmentation of the plate boundary as defined by earthquakes this century is not constant.     

Mainshocks and aftershocks of the 2002 molise seismic sequence, southern Italy

In October and November 2002, the Molise region (southern Italy) was struck by two moderate magnitude earthquakes within 24 hours followed by an one month long aftershocks sequence. Soon after the first mainshock (October 31st, 10.32 UTC, M_w 5.7), we deployed a temporary network of 35 three-component seismic stations. At the time of occurrence of the second main event (November 1st, 15.08 UTC, M_w 5.7) the eight local stations already installed allowed us to well constrain the hypocentral parameters. We present the location of the two mainshocks and 1929 aftershocks with 2 < M_L < 4.2. Earthquake distribution reveals a E-trending 15 km long fault system composed by two main segments ruptured by the two mainshocks. Aftershocks define two sub-vertical dextral strike-slip fault segments in agreement with the mainshock fault plane solutions. P- and T-axes retrieved from 170 aftershocks focal mechanisms show a coherent kinematics: with a sub-horizontal NW and NE-trending P and T-axes, respectively. For a small percentage of focal mechanisms (∼ 10%) a rotation of T axes is observed, resulting in thrust solutions. The Apenninic active normal fault belt is located about 80 km westward of the 2002 epicentral area and significant seismicity occurs only 20-50 km to the east, in the Gargano promontory. Seismic hazard was thought to be small for this region because neither historical earthquake are reported in the Italian seismic catalogue or active faults were previously identified. In this context, the 2002 seismic sequence highlights the existence of trans-pressional active tectonics in between the extensional Apenninic belt and the Apulian foreland.     

The 1986 Dharamsala earthquake of Himachal Himalaya – estimates of source parameters,average intrinsic attenuation and site amplification functions

The 26th April 1986 Dharamsala earthquake (m_b 5.5) occurred in the Kangra region of Himachal Himalaya, which lies in the rupture zone of great Kangra earthquake of 1905. This was the first moderate sized earthquake to be recorded at a few sites of the strong ground motion array in the NW Himalaya. The accelerograms of this earthquake have been used to estimate its source parameters, site amplification functions and to estimate the effective shear wave attenuation factor Q_β in the frontal region of Himachal Himalaya. A double couple fault plane solution for the earthquake has been obtained based on the spectra of the transverse component of the accelerograms. The estimated values of the source parameters are seismic moment: 2.1×10²⁴ dyne-cm, static stress drop (Δσ): 36 bars, source radius (r): 2.8 km and moment magnitude (M_w): 5.4. The estimated average values of effective shear wave attenuation factor Q_β for various sites are in the range of 125 to 300 with an overall spatial average of 239. The influence of local site effects on the observed PGA values have been examined on the basis of site amplification functions.     

Analysis and interpretation of the October 21, 1766 earthquake in the Southeastern Caribbean

The October 21, 1766 earthquake is the most widely felt event in the seismic history of Trinidad and Venezuela. Previous works diverged on the interpretation of the historical data available for this event. They associated the earthquake either with the Lesser Antilles subduction zone, with strike-slip motion along El Pilar fault, or with intraplate deformation at the edge of Guyana shield. Isoseismal areas are proposed after a new search and analysis of primary and secondary sources of historical information. Two of the largest earthquakes of the twentieth century which occurred in the region, the 1968 (M _S 6.4, h = 103 km), and the 1997 (M _W 6.9, h = 25 km) events, for which both intensity data and instrumentally determined source parameters are available, are used to calibrate the isoseismal areas and to interpret them in terms of source depth and magnitude. It is concluded that the large extent of intensity values higher than V is diagnostic of the depth (85 ± 20 km) of the 1766 source, and of local amplifications of ground motion due to soft soil conditions and to strong contrasts of impedance at the edge of Guyana shield. It is proposed that the event occurred either in slab, or close to the bottom lithospheric interface between the Caribbean and South American plates (∼11°N; ∼62.5°W). The value of the magnitude is estimated at 6.5 < M _S < 7.5 depending on the source depth and on the decay of ground motion as a function of distance. Deep and intermediate depth earthquakes can induce important casualties in Trinidad, Venezuela, and Guyana, possibly more damaging than those induced by shallower earthquakes along the strike of El Pilar Fault.     

CPO-induced seismic anisotropy in UHP eclogites

Ultrahigh-pressure (UHP) eclogites often show strong plastic deformation and anisotropy of seismic properties. We report in this paper the seismic velocity and anisotropy of eclogite calculated from the crystallographic preferred orientations (CPOs) of constituent minerals (garnet, omphacite, quartz and rutile) and single crystal elastic properties. We also compared the calculated results with the measured results in similar eclogites. Our results suggest that (1) Except that garnet is a seismically quasi-isotropic mineral, omphacite, quartz, coesite and rutile all have strong seismic anisotropies (AVp = 23.0%―40.9%, Max. AVs = 18.5%―47.1%). They are the major sources for anisotropy in eclogite. The average seismic velocities are fast in garnet and rutile, moderate in omphacite and coesite, and slow in quartz. (2) The deformed eclogites have the maximum Vp (8.33―8.75 km/s) approximately parallel to foliation and lineation, the minimum Vp (8.25―8.62 km/s) approximately normal to foliation and lineation and the Vp anisotropies of 1.0―1.7%. Their Vs are 4.93―4.97 km/s. The corresponding maximum anisotropies (0.73%―1.78%) of Vs are at 45° to both foliation and lineation and the minimum anisotropies at positions normal to lineation on the foliation plane. The Vs1 polarization planes are approximately parallel to foliation. The mean Vp and Vs of eclogite under UHP peak metamorphism conditions (P = 3―5 GPa, T = 900―1100℃) are estimated to be 3.4%―7.2% and 6.3%―12.1% higher than those at ambient pressure and temperature conditions, respectively. (3) Omphacite component dominates the anisotropy of eclogite while garnet component reduces the anisotropy and increases the seismic velocities. Quartz component has a small effect on the anisotropy but reduces the seismic velocities of eclogite. The effect of rutile component is negligible on seismic properties of eclogite due to its trivial volume fraction. (4) The increase of volume fraction of omphacite in eclogite will reduce the seismic velocities and increase the anisotropy. Omphacitite has seismic velocities reduced by 6%―8% and anisotropies increased to 3%―4% compared to those of garnetite. Our results suggest that the seismic properties calculated with single crystal elastic properties and CPOs are equivalent to those measured in laboratory. Moreover, it provides insights into the mineral physical interpretations of eclogite seismic properties.     

Source mechanism and source parameters of May 28, 1998 earthquake,Egypt

On May 28, 1998, a moderate size earthquake of mb 5.5 occurred offshore the northwestern part of Egypt (latitude 31.45°N and longitude 27.64°E). It was widely felt in the northern part of Egypt. Being the largest well-recorded event in the area for which seismic data from the global digital network are available, it provides an excellent opportunity to study the tectonic process and present day stress field occurring along the offshore Egyptian coast. The source parameters of this event are determined using three different techniques: modeling of surface wave spectral amplitudes, regional waveform inversion, and teleseismic body waveform inversion. The results show a high-angle reverse fault mechanism generally trending NNW–SSE. The P-axis trends ENE–WSW consistently with the prevailed compression stress along the southeastern Hellenic arc and southwestern part of the Cyprean arc. This unexpected mechanism is most probably related to a positive inversion of the NW trending offshore normal faults and confirms an extension of the back thrusting effects towards the African margin. The estimated focal depth ranges from 22 to 25 km, indicating a lower crustal origin earthquake owing to deep-seated tectonics. The source time function indicates a single source with rise time and total rupture duration of 2 and 5 s, respectively. The seismic moment (M _o) and the moment magnitude (M _w) determined by the three techniques are 1.03 × 10¹⁷ Nm, 5.28; 1.24 × 10¹⁷ Nm, 5.33; and 1.68 × 10¹⁷ Nm, 5.42; respectively. The calculated fault radius, stress drop, and the average dislocation assuming a circular fault model are 7.2 km, 0.63 Mpa, and 0.11 m, respectively.     

2015年11月23日青海祁连MS5.2地震震源机制解及发震构造初步探讨

采用CAP（Cut and Paste）方法反演了2015年11月23日青海祁连M_S5.2主震的震源机制解,其最佳双力偶解：节面Ⅰ走向109°、倾角58°、滑动角21°,节面Ⅱ走向8°、倾角72°、滑动角146°,矩震级M_W5.16,矩心震源深度约为9 km。结合震区的活动构造,判定发震断层面为节面Ⅰ,推测托勒山北缘活动断裂中段为此次地震的发震断裂。     

The 1992 Tafilalt seismic crisis (Anti-Atlas, Morocco)

The Tafilalt region, located at the eastern end of the Anti-Atlas chain in Morocco, was shaken on 23 and 30 October 1992 by two moderate earthquakes of magnitude mb ∼ 5 and intensity ∼ VI MSK64, which caused two deaths and great damage in the area between Erfoud and Rissani. The review of data available on the seismic crisis allowed us to improve the knowledge on the macroseismic, instrumental and source parameters of the earthquakes. The main results of the present study are: (1) location of the epicentres with the help of data from a close portable network allowed us to propose new epicentral coordinates at 31.361° N, 4.182° W (23 October) and 31.286° N, 4.347° W (30 October); both events have focal depths of 2 km; (2) the shock of 30 October was followed by a series of 305 aftershocks, most of which were located west of Rissani; the 61 best-constrained events had focal depths of 5 to 19 km and magnitudes 0.7 to 3; (3) the largest damage was located in an area between the two epicentres within the Tafilalt valley and was probably amplified by site effects due to the proximity of the water table within the Quaternary sediments; (4) focal mechanisms of the main events correspond to strike-slip faulting with fault planes oriented N–S (left lateral) and E–W (right lateral); the only mechanism available for the aftershocks also corresponds to strike-slip faulting; (5) spectral analysis shows that the scalar seismic moment (Mo) of the first event is slightly larger than the second; the corresponding values of Mw are 5.1 and 5.0, respectively; (6) the dimensions of the faults for a circular fault model are 7.7 ± 1.4 and 7.4 ± 1.2 km, respectively; the average displacement is 4 cm for the first event and 3.7 cm for the second; the stress drop is 0.4 and 0.3 MPa, respectively, in agreement with standard values; (7) the Coulomb Stress test performed for both earthquakes suggests a relationship between both events only when the used location is at the limit of the horizontal uncertainty; (8) finally, the occurrence of these shocks suggests that the Anti-Atlas is undergoing tectonic deformation in addition to thermal uplift as suggested by recent publications.     

Precise determination of focal parameters for July 20, 1995 M_L=4.1 earthquake sequence in the Huailai basin

IntroductionSince the late 1970s, the quickly developed global digital seismograph network has been providing high quality recordings of large earthquakes in global scale, based on which digital seismology has made great progress. Compared with large earthquakes, moderate and small sized shocks have more frequent occurrence, and comprise clues to geological tectonics and tectonic stress field in a region. Preceding and following a large earthquake, usually occur numbers of small events that im…     

A low magnitude seismic sequence near Isernia (Molise,Central Italy) in January 1986

Following the increase in seismic activity which occurred near Isernia (Molise, Central Italy) in January 1986, a digital seismic network of four stations with three-component, short-period seismometers, was installed in the area by the Osservatorio Vesuviano. The temporary network had an average station spacing of about 8–10 km and, in combination with permanent local seismic stations, allowed the accurate determination of earthquake locations during an operating period of about one month. Among the 1517 detected earthquakes, 170 events were selected with standard errors on epicentre and depth not greater than respectively 0.5 and 1.5 km. The most frequent focal depths ranged between 4 and 8 km, while the epicentres distribution covered a small area NE of Isernia of about 10 km². A main rupture zone could not be clearly identified from the spatial distribution of the earthquakes, suggesting a rupture mechanism involving heterogeneous materials. The activity was characterized by low energy levels, the largest earthquake, on January 18, 1986, havingM _D =4.0. The time sequence of events and pattern of seismic energy release revealed a strong temporal clustering of events, similar to the behaviour commonly associated with seismic swarms.     

青藏工程走廊地震危险性对比分析及区划

潜在震源区地震活动性参数、地震动衰减关系对地震危险性分析结果至关重要。以中国第五代地震动区划图潜在震源区划分方案为基础,采用2类震级分档分别建立自编及五代图潜在震源区空间分布函数,收集4组青藏高原及周缘地区地震动衰减关系,采用不同组合对青藏工程走廊沿线的81个场点进行概率地震危险性分析计算,得到50年超越概率10%(地震重现期475年)的各场点基岩地震动峰值加速度(PGA),并转换为一般场地(Ⅱ类)PGA,对计算结果进行对比分析,并与第五代地震动区划图归档上下限值进行比较。结果显示:采用我国西部地区地震动衰减关系计算得到的PGA最大,采用云南地区地震动衰减关系得到的PGA最小,采用川藏地区及青藏高原东北缘地震动衰减关系时居中;在同一地震动衰减关系下,采用自编空间分布函数计算得到的PGA普遍略大于采用第五代图空间分布函数时;在震级上限为8.5的潜在震源区及附近地区,潜在震源区空间分布函数震级分档对计算结果有显著影响。综合分析表明,采用自编Ⅱ型震级分档空间分布函数方案与川藏地区地震动衰减关系组合方案的计算结果最为理想。最后,采用该组合方案对青藏工程走廊50年超越概率10%的基岩场地PGA及一般场地PGA进行了区划。     

山东长岛地区地震波非弹性衰减Q值、场地响应及震源参数研究

利用山东台网记录的长岛震群2017年2月14日—9月1日期间的波形与震相资料研究长岛地区非弹性衰减系数,得到该地区介质平均Q值与频率f的关系式为Q(f)=363.9f1.374 1。采用Moya等[1]提出的利用遗传算法联合反演得到长岛周边台站的场地响应,根据Brune模型震源参数计算公式求解长岛震群序列地震震源参数。结果显示,各个震源参数之间均存在一定的相关关系,地震矩随ML震级的增大而增加,地震矩与破裂半径R之间存在半对数关系,拐角频率fc随地震矩的增大而减少;长岛地震序列的应力降数值普遍偏小,最大不超过0.9MPa,这意味着长岛震源区整体构造应力较低,也可能指示长岛震群为低摩擦应力的断层作用;震源参数随时间的变化方面,整体而言,长岛震群地震应力降变化起伏很大,在M4.1地震发生前,拐角频率与应力降均发生快速下降后随即翻转上升的现象,证明在M4.1地震发生前震源区整体应力的挤压逐渐增强。     

Active segment of the 12 November 2003 Mw 5.6 earthquake at Salsipuedes oceanic basin,Gulf of California. Mexico

We analyzed records of eight seismic stations of the autonomous broadband seismograph network of a joint project between Utrecht University (the Netherlands), California Institute of Technology, and Centro de Investigación Científica y de Estudios Superiores de Ensenada (CICESE). These stations recorded the Mw 5.6 earthquake that occurred on 12 November 2003 at Salsipuedes basin in the middle of the Gulf of California 2 km west of the island Angel de la Guarda. This event was located at 29.16º N and 113.37º W, 30 km northeast of Bahia de los Angeles. A foreshock and hundreds of aftershocks were recorded in the 48 hours after its origin time. With the location of 29 earthquakes we identified the active segment, perpendicular to the main transform fault NW–SE of Canal de Ballenas, representing the transtensional boundary between the Pacific and North American plates. The direction of the active fault described is consistent with the normal fault mechanism reported by the National Earthquake Information Center (strike=39º, dip=34º, slip=–44º).From the duration magnitude of 456 aftershocks, we calculated a b-value of 1.14±0.28; furthermore, we calculated a seismic moment of (3.5 ±3.3) X10¹⁷Nm, a source radius of 3.7 ± 2.63 km, and a static stress drop of 3.94 ± 1.15 MPa (39.4 ± 11.5 bar.).     

The Yecheng, Xinjiang, China, earthquake of February 14, 1980 —A destructive intermediate earthquake

In this study, we have checked the location and focal depth of the Yecheng earthquake (m _b = 6.0, maximum intensity VII) of February 14, 1980. The result shows that this is an intermediate event with a focal depth of 90 km. The microepicenter is located at 36. 4°N, 76, 9°E, while the macroepicenter is at 37.3°N, 76.9°E, 90 km to the north of the microepicenter. This is the first destructive intermediate event in China which led to a damage as severe as of intensity VII. The focal mechanism of the event is determined to be of thrust type. Combined with the analysis of seismological and geological data in surrounding area, the possible relation between the event and plate movement has been discussed. The result in this paper indicates that in some particular place, the destructive effect of intermediate event should be considered in seismic hazard assessment. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 137–143, 1992.     

The Ain Temouchent (Algeria) Earthquake of December 22<Superscript>nd</Superscript>, 1999

—On December 22^nd, 1999 an earthquake of Magnitude M_w : 5.7 occurred at Ain Temouchent (northwest Algeria). This moderate seismic event was located in a region characterized by a low seismic activity where few historical events have been observed. The earthquake, with a maximum intensity of VII (MSK scale), caused serious damages to the Ain Temouchent city and its surroundings. In the epicentral area, 25 people died and about 25,000 people were made homeless. Some minor breaks have been observed in several areas in the field. They were mainly related to minor collapses in the landscape or in volcanic cavities. The focal mechanism has been studied by using broadband data at regional and teleseismic distances, and different methods. The fault-plane solution has been estimated from first motions of P wave. Depth and source time function have been estimated from the modeling of body waveforms. Scalar seismic moment and source dimension have been obtained from spectral analysis. Results show thrust motion, with a horizontal pressure axis oriented in a NW-SE direction, a depth of 4 km and a simple source time function with time duration of 5 s. Scalar seismic moment estimated from waveform modeling is 4.7 × 10¹⁷ Nm, and spectral analysis gives a value of 1.7 × 10¹⁷ Nm and a source radius of 7.5 km.