An SPH Framework for Earthquake-Induced Liquefaction Hazard Assessment of Geotechnical Structures

Earthquake-induced soil liquefaction poses significant risks to the stability of geotechnical structures worldwide.An understanding of the liquefaction triggering,and the post-failure large deformation behaviour is essential for designing resilient infrastructure.The present study develops a Smoothed Particle Hydrodynamics(SPH)framework for earthquake-induced liquefaction hazard assessment of geotechnical structures.The coupled flowdeformation behaviour of soils subjected to cyclic loading is described using the PM4Sand model implemented in a three-phase,single-layer SPH framework.A staggered discretisation scheme based on the stress particle SPH approach is adopted to minimise numerical inaccuracies caused by zero-energy modes and tensile instability.Further,non-reflecting boundary conditions for seismic analysis of semi-infinite soil domains using the SPH method are proposed.The numerical framework is employed for the analysis of cyclic direct simple shear test,seismic analysis of a level ground site,and liquefaction-induced failure of the Lower San Fernando Dam.Satisfactory agreement for liquefaction triggering and post-failure behaviour demonstrates that the SPH framework can be utilised to assess the effect of seismic loading on field-scale geotechnical structures.The present study also serves as the basis for future advancements of the SPH method for applications related to earthquake geotechnical engineering.

Large-scale shaking table test on unlined tunnel in fault zone under threedimensional earthquake

A fault is a geological structure characterized by significant displacement of rock masses along a fault plane within the Earth's crust.The Yunnan Tabaiyi Tunnel intersects multiple fault zones,making tunnel construction in fault-prone areas particularly vulnerable to the effects of fault activity due to the complexities of the surrounding geological environment.To investigate the dynamic response characteristics of tunnel structures under varying surrounding rock conditions,a three-dimensional large-scale shaking table physical model test was conducted.This study also aimed to explore the damage mechanisms associated with the Tabaiyi Tunnel under seismic loading.The results demonstrate that poor quality surrounding rock enhances the seismic response of the tunnel.This effect is primarily attributed to the distribution characteristics of acceleration,dynamic strain,and dynamic soil pressure.A comparison between unidirectional and multi-directional(including vertical)seismic motions reveals that vertical seismic motion has a more significant impact on specific tunnel locations.Specifically,the maximum tensile stress is observed at the arch shoulder,with values ranging from 60 to 100 k Pa.Moreover,NPR(Non-Prestressed Reinforced)anchor cables exhibit a substantial constant resistance effect under low-amplitude seismic waves.However,when the input earthquake amplitude reaches 0.8g,local sliding occurs at the arch shoulder region of the NPR anchor cable.These findings underscore the importance of focusing on seismic mitigation measures in fault zones and reinforcing critical areas,such as the arch shoulders,in practical engineering applications.

Geological risk assessment of traffic engineering construction among 7.0-8.5 magnitude earthquake areas:Practice from the Sichuan-Tibet transport corridor in the eastern Tibetan Plateau

At least 13 active fault zones have developed in the Ya'an-Linzhi section of the Sichuan-Tibet transport corridor,and there have been undergone 17 MS≥7.0 earthquakes,the largest earthquake is 1950 Chayu MS 8.5 earthquake,which has very strong seismic activity.Therefore,carrying out engineering construction in the Sichuan-Tibet transport corridor is a huge challenge for geological technological personnel.To determining the spatial geometric distribution,activity of active faults and geological safety risk in the Sichuan-Tibet transport corridor.Based on remote sensing images,ground surveys,and chronological tests,as well as the deep geophysical and current GPS data,we investigated the geometry,segmentation,and paleoearthquake history of five major active fault zones in the Ya'an-Linzhi section of the Sichuan-Tibet transport corridor,namely the Xianshuihe,Litang,Batang,Jiali-Chayu and Lulang-Yigong.The five major fault zones are all Holocene active faults,which contain strike-slip components as well as thrust or normal fault components,and contain multiple branch faults.The Selaha-Kangding segment of the Xianshuihe fault zone,the Maoyaba and Litang segment of the Litang fault zone,the middle segment(Yigong-Tongmai-Bomi)of Jiali-Chayu fault zone and Lulang-Yigong fault zone have the risk of experiencing strong earthquakes in the future,with a high possibility of the occurrence of MS≥7.0 earthquakes.The Jinsha River and the Palong-Zangbu River,which is a high-risk area for geological hazard chain risk in the Ya'an-Linzhi section of the Sichuan-Tibet transport corridor.Construction and safe operation Ya'an-Linzhi section of the Sichuan-Tibet transport corridor,need strengthen analysis the current crustal deformation,stress distribution and fault activity patterns,clarify active faults relationship with large earthquakes,and determine the potential maximum magnitude,epicenters,and risk range.This study provides basic data for understanding the activity,seismicity,and tectonic deformation patterns of the regional faults in the Sichuan-Tibet transport corridor.

Key technologies of earthquake early warning system for China’s high-speed railway

Purpose–The purpose of this study is to introduce the top-level design ideas and the overall architecture of earthquake early-warning system for high speed railways in China,which is based on P-wave earthquake early-warning and multiple ways of rapid treatment.Design/methodology/approach–The paper describes the key technologies that are involved in the development of the system,such as P-wave identification and earthquake early-warning,multi-source seismic information fusion and earthquake emergency treatment technologies.The paper also presents the test results of the system,which show that it has complete functions and its major performance indicators meet the design requirements.Findings–The study demonstrates that the high speed railways earthquake early-warning system serves as an important technical tool for high speed railways to cope with the threat of earthquake to the operation safety.The key technical indicators of the system have excellent performance:The first report time of the P-wave is less than three seconds.From the first arrival of P-wave to the beginning of train braking,the total delay of onboard emergency treatment is 3.63 seconds under 95%probability.The average total delay for power failures triggered by substations is 3.3 seconds.Originality/value–The paper provides a valuable reference for the research and development of earthquake early-warning system for high speed railways in other countries and regions.It also contributes to the earthquake prevention and disaster reduction efforts.

Three-dimensional fault model and features of chained hazards of the Luding M_(S)6.8 earthquake,Sichuan Province,China

The M_(S)6.8 Luding earthquake in 2022 is located on the NNW-trending Moxi segment of the Xianshuihe fault with left-lateral strike-slip behavior.This area is where the Xianshuihe,Anninghe,Daliangshan and Longmenshan faults intersect.China Earthquake Administration has identified that intersection area,among the Moxi segment of the Xianshuihe fault,the Anninghe fault,the Daliangshan fault and the southern part of the Longmenshan fault,as a high-magnitude earthquake hazard area.According to existing data on the Luding earthquake,including the focal parameters,the spatial distribution of re-located aftershocks,dominated azimuth of the earthquake intensities and earthquake-induced ground fissures,we built a 3D earthquake fault model.We found that two discontinuous NNW-trending vertical strike-slip faults with left stepping were the seismogenic faults of the Luding earthquake.Its coseismic left-lateral dislocation triggered transtensional slips and aftershocks on the NW-trending secondary faults at its northernmost tensile area.Meanwhile,local crustal coseismic shortening on the side of Mt.Gongga triggered the aftershocks on the NE-and NW-trending secondary conjugated strike-slip faults,which were confirmed by GNSS observations and In SAR deformation field around the epicenter.This earthquake rupturing pattern also controlled the spatial distribution of the earthquake intensity IX area and earthquake chain hazards.The Coulomb stress calculation shows that the Luding earthquake increases the risk of high-magnitude earthquake occurrence on the southernmost part of the Xianshuihe fault and the Anninghe fault.Finally,we suggested doing good monitoring of the Anninghe fault and the southernmost part of the Xianshuihe fault and avoiding active faults with seismogenic capacity and areas prone to earthquake-chained hazards during the site selection and planning of reconstruction.

The MW5.5 earthquake on August 6,2023,in Pingyuan,Shandong,China:A rupture on a buried fault

On August 6,2023,a magnitude MW5.5 earthquake struck Pingyuan County,Dezhou City,Shandong Province,China.This event was significant as no large earthquakes had been recorded in the region for over a century,and no active fault had been previously identified.This study collects 1309 P-wave arrival times and 866 S-wave arrival times from 74 seismic stations less than 200 km to the epicenter to constrain the spatial distribution of the mainshock and its 125 early aftershocks by the double difference earthquake relocation method,and selects 864 P-waveforms from 288 stations located within 800 km of the epicenter to constrain the focal mechanism solution of the mainshock through centroid moment tensor inversion.The relocation and the inversion indicate,the Pingyuan MW5.5 earthquake was caused by a rupture on a buried fault,likely an extensive segment of the Gaotang fault.This buried fault exhibited a dip of approximately 75°to the northwest,with a strike of 222°,similar to the Gaotang fault.The rupture initiated at the depth of 18.6 km and propagated upward and northeastward.However,the ground surface was not broken.The total duration of the rupture was~6.0 s,releasing the scalar moment of 2.5895×1017 N·m,equivalent to MW5.54.The moment rate reached the maximum only 1.4 seconds after the rupture initiation,and the 90%scalar moment was released in the first 4.6 s.In the first 1.4 seconds of the rupture process,the rupture velocity was estimated to be 2.6 km/s,slower than the local S-wave velocity.As the rupture neared its end,the rupture velocity decreased significantly.This study provides valuable insights into the seismic characteristics of the Pingyuan MW5.5 earthquake,shedding light on the previously unidentified buried fault responsible for the seismic activity in the region.Understanding the behavior of such faults is crucial for assessing seismic hazards and enhancing earthquake preparedness in the future.

Rapid report of source parameters of 2023 M6.2 Jishishan,Gansu earthquake sequence

The M6.2 earthquake in Jishishan,Gansu Province,on December 18,2023,caused extraordinary earthquake disasters.It was located in the northern part of the north−south seismic zone,which is a key area for earthquake monitoring in China.The newly built dense strong motion stations in this area provide unprecedented conditions for high-precision earthquake relocation,especially the earthquake focal depth.This paper uses the newly built strong motion and traditional broadband seismic networks to relocate the source locations of the M3.0 and above aftershocks and to invert their focal mechanisms.The horizontal error of earthquake location is estimated to be 0.5−1 km,and the vertical error is 1−2 km.The focal depth range of aftershocks is 9.6−14.6 km,distributed in a 12-km-long strip with SSE direction.Aftershocks in the south are more concentrated horizontally and vertically,while aftershocks in the north are more scattered.The focal mechanisms of the main shock and aftershocks are relatively consistent,and the P-axis orientation is consistent with the regional strain direction.There is a seismic blank area of M3.0 and above,about 3−5 km between the main shock and aftershocks.It is suggested that the energy released by the main shock rupture is concentrated in this area.Based on the earthquake location and focal mechanism of the main shock,it is inferred that the Northern Lajishan fault zone is the seismogenic structure of the main shock,and the main shock did not occur on the main fault,but on a secondary fault.The initial rupture depth and centroid depth of the main shock were 12.8 and 14.0 km,respectively.The source rupture depth may not be the main reason for the severe earthquake disaster.

Characterization and spatial analysis of coseismic landslides triggered by the Luding Ms 6.8 earthquake in the Xianshuihe fault zone, Southwest China

On September 5, 2022, a magnitude Ms 6.8 earthquake occurred along the Moxi fault in the southern part of the Xianshuihe fault zone located in the southeastern margin of the Tibetan Plateau,resulting in severe damage and substantial economic loss. In this study, we established a coseismic landslide database triggered by Luding Ms 6.8 earthquake, which includes 4794 landslides with a total area of 46.79 km^(2). The coseismic landslides primarily consisted of medium and small-sized landslides, characterized by shallow surface sliding. Some exhibited characteristics of high-position initiation resulted in the obstruction or partial obstruction of rivers, leading to the formation of dammed lakes. Our research found that the coseismic landslides were predominantly observed on slopes ranging from 30° to 50°, occurring at between 1000 m and 2500 m, with slope aspects varying from 90° to 180°. Landslides were also highly developed in granitic bodies that had experienced structural fracturing and strong-tomoderate weathering. Coseismic landslides concentrated within a 6 km range on both sides of the Xianshuihe and Daduhe fault zones. The area and number of coseismic landslides exhibited a negative correlation with the distance to fault lines, road networks, and river systems, as they were influenced by fault activity, road excavation, and river erosion. The coseismic landslides were mainly distributed in the southeastern region of the epicenter, exhibiting relatively concentrated patterns within the IX-degree zones such as Moxi Town, Wandong River basin, Detuo Town to Wanggangping Township. Our research findings provide important data on the coseismic landslides triggered by the Luding Ms 6.8 earthquake and reveal the spatial distribution patterns of these landslides. These findings can serve as important references for risk mitigation, reconstruction planning, and regional earthquake disaster research in the earthquake-affected area.

Microearthquake reveals the lithospheric structure at midocean ridges and oceanic transform faults

Mid-ocean ridge and oceanic transforms are among the most prominent features on the seafloor surface and are crucial for understanding seafloor spreading and plate tectonic dynamics,but the deep structure of the oceanic lithosphere remains poorly understood.The large number of microearthquakes occurring along ridges and transforms provide valuable information for gaining an indepth view of the underlying detailed seismic structures,contributing to understanding geodynamic processes within the oceanic lithosphere.Previous studies have indicated that the maximum depth of microseismicity is controlled by the 600-℃isotherm.However,this perspective is being challenged due to increasing observations of deep earthquakes that far exceed this suggested isotherm along mid-ocean ridges and oceanic transform faults.Several mechanisms have been proposed to explain these deep events,and we suggest that local geodynamic processes(e.g.,magma supply,mylonite shear zone,longlived faults,hydrothermal vents,etc.)likely play a more important role than previously thought.

Myths about Earthquakes:Quo vadis?

We are living in a world of numbers and calculations with enormous amount of pretty fast user-friendly software ready for an automatic output that may lead to a discovery or,alternatively,mislead to a deceptive conclusion,erroneous claims and predictions.As a matter of fact,nowadays,Science can disclose Natural Hazards,assess Risks,and deliver the state-of-the-art Knowledge of looming disaster in advance catastrophes along with useful Recommendations on the level of risks for decision making regarding engineering design,insurance,and emergency management.

Research on the Extraction of OLR Anomaly Prior to Ms 7.5 Sand Point,Alaska Earthquake based on IPI Method

Based on the understanding that the seismic fault system is a nonlinear complex system,Rundle(1995)introduced the nonlinear threshold system used in meteorology to analyze the ocean-atmosphere interface and the El Ni?o Southern Oscillation into the study of seismic activity changes,and then proposed the PI method(Rundle et al.,2000a,b).Wu et al.(2011)modified the Pattern Informatics Method named MPI to extract the ionospheric anomaly by using data from DEMETER satellites which is suitable for 1–3 months short-term prediction.

The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

On February 6,2023,two earthquakes with magnitudes of M_(W) 7.8 and M_(W) 7.5 struck southeastern Turkey,causing significant casualties and economic losses.These seismic events occurred along the East Anatolian Fault Zone,a convergent boundary between the Arabian Plate and the Anatolian Subplate.In this study,we analyze the M_(W) 7.8 and M_(W) 7.5 earthquakes by comparing their aftershock relocations,tomographic images,and stress field inversions.The earthquakes were localized in the upper crust and exhibited steep dip angles.Furthermore,the aftershocks occurred either close to the boundaries of low and high P-wave velocity anomaly zones or within the low P-wave velocity anomaly zones.The East Anatolia Fault,associated with the M_(W) 7.8 earthquake,and the SürgüFault,related to the M_(W) 7.5 earthquake,predominantly experienced shear stress.However,their western sections experienced a combination of strike-slip and tensile stresses in addition to shear stress.The ruptures of the M_(W) 7.8 and M_(W) 7.5 earthquakes appear to have bridged a seismic gap that had seen sparse seismicity over the past 200 years prior to the 2023 Turkey earthquake sequence.

Strong Earthquake Model of North China Craton: A Case of 2023 Mw 5.5 Earthquake in Pingyuan County, Shandong Province, China

North China is one of the high-risk areas for destructive strong earthquakes in China's Mainland, with a history of numerous significant seismic events. On August 6, 2023, an Mw5.5 earthquake struck Pingyuan County, Dezhou City, in Shandong Province, China. This earthquake was the largest in the eastern North China Craton(NCC) since the Tangshan earthquake of 1976. Due to the absence of surface ruptures, the fault responsible for the Pingyuan Mw5.5 earthquake remains unclear. To reveal the subsurface geological structure near the earthquake epicenter, this study utilized highresolution two-dimensional seismic reflection profiles to interpret pre-existing faults.

Identifying drivers of urban landuse changes in the Wenchuan earthquake- affected area by using night-time light data

To learn the process of urban land evolution before and after an earthquake is vital to formulate the urban reconstruction control policies and recovery measures in the earthquake-stricken areas.However,spatiotemporal evolution and its driving factors of urban land in earthquake-prone areas remains limited due to the scarcity of ground observation data.This research,leveraging night-time light remote sensing imagery and land cover data,conducted a comprehensive analysis of the long-term evolution characteristics of urban land in earthquake-prone areas.It introduced methodologies for assessing the socio-economic impact and the primary natural environmental factors driving urban land evolution in these regions.To validate the proposed methods,the 2008 Wenchuan earthquake-affected area in China was selected as a representative study area.The results indicated that the average Digital Number(DN)values in socio-economically impacted areas showed a trend of rising,falling,and then rising again after the earthquake.DN values in three types of damaged areas including Type Ⅱ,Type Ⅲ,and Type Ⅳ exceeded pre-earthquake levels.The analysis of determinative factors influencing urban land evolution revealed that slope and elevation were key elements in controlling urban land expansion before the earthquake,whereas factors such as slope,elevation,lithology,and faults had a stronger influence on urban land expansion after the earthquake.It can be seen that,in view of the differences in the natural conditions of regions for post-disaster reconstruction,the local government need to actively adjust and adapt to urban spatial planning,so as to leverage the scale effect of large-scale inputs of funds,facilities,human resources and other factors after the disaster,thus enhancing resilience and recovery efficiency in response to disaster impacts.

Stress triggering effect on the 2022 Honghe M_(S)5.0 earthquake with historical strong earthquakes

The 2022 Honghe M_(S)5.0 seismic event is intriguing due to its occurrence in the south of the Red River Fault,an area historically lacking seismic activities greater than M_(S)5.0.To elucidate the seismogenic mechanism and scrutinize stress-triggered interactions,we calculated co-seismic and post-seismic Coulomb stress alterations induced by nine historical seismic events(M≥6.0).The analysis reveals that these substantial seismic events provoked co-seismic stress augmentations of 1.409 bar and postseismic stress increments of 0.159 bar.Noteworthy seismic events,such as the 1833 Songming,1877Shiping,1913 Eshan,and 1970 Tonghai earthquakes,catalyzed the occurrence of the Honghe earthquake.Areas of heightened future seismic risk include the southern region of the Red River Fault and the eastern segments of the Shiping-Jianshui and Qujiang faults.Additionally,we assessed the correlation between the spatial distribution of aftershocks and the Coulomb stress shift triggered by the mainshock,taking into account the influence of calculation parameter settings.

Seismicity patterns before the 2021 Fin (Iran) doublet earthquakes using the region-time-length and time-to-failure methods

Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied. Furthermore, the estimation of the region-time-length (RTL) parameter has been proposed to detect seismic quiescence before the occurrence of a large earthquake. In addition, the time-to-failure method has been used to estimate the time occurrence of large earthquakes. Hence, in this study, to gain deeper insight into seismic activity in the southern Zagros region, we utilized the RTL algorithm to identify the quiescence and activation phases leading to the Fin doublet earthquakes. Temporal variations in the RTL parameter showed two significant anomalies. One corresponded to the occurrence time of the first earthquake (2017-12-12);the other anomaly was associated with the occurrence time of the second event (2021-11-14). Based on a negative value of the RTL parameter observed in the vicinity of the Fin epicenters (2021), seismic quiescence (a decrease in seismicity compared to the preceding background rate) was identified. The spatial distribution of the RTL prognostic parameters confirms the appearance of seismic quiescence surrounding the epicenter of the Fin doublet earthquakes (2021). The time-to-failure method was designed using precursory events that describe the acceleration of the seismic energy release before the mainshock. Using the time-to-failure method for the earthquake catalog, it was possible to estimate both the magnitude and time of failure of the Fin doublet. Hence, the time-tofailure technique can be a useful supplementary method to the RTL algorithm for determining the characteristics of impending earthquakes.

The 2024 M_(j) 7.6 Noto Peninsula, Japan earthquake caused by the fluid flow in the crust

On January 1, 2024 at 16:10:09 JST, an M_(j) 7.6 earthquake struck the Noto Peninsula in the southern part of the Sea of Japan. This location has been experiencing an earthquake swarm for more than three years. Here, we provide an overview of this earthquake, focusing on the slip distribution of the mainshock and its relationship with the preceding swarm. We also reexamined the source areas of other large earthquakes that occurred around the Sea of Japan in the past and compared them with the Matsushiro earthquake swarm in central Japan from1964 to 1968. The difference between the Matsushiro earthquake swarm and the Noto earthquake swarm is the surrounding stress field. The Matsushiro earthquake swarm was a strike-slip stress field, so the cracks in the crust were oriented vertically. This allowed fluids seeped from the depths to rise and flow out to the surface. On the other hand, the Noto area was a reverse fault stress field. Therefore, the cracks in the earth's crust were oriented horizontally. Fluids flowing underground in deep areas could not rise and spread over a wide area in the horizontal plane. This may have caused a large amount of fluid to accumulate underground, triggering a large earthquake. Although our proposed mechanism does not take into account other complex geological conditions into consideration, it may provide a simple way to explain why the Noto swarm is followed by a large earthquake while other swarms are not.

Rapid determination of source parameters of the M 6.2 Jishishan earthquake in Gansu Province and its application in emergency response

In this study, we swiftly determined the focal parameters(focal mechanism, seismic imaging process, magnitude)of the Jishishan earthquake, leveraging a solved fault model to assess the intensity field and casualties promptly.The investigation began by retrieving the source mechanism through the P-wave initial motion and W-phase method. This enabled us to chart the spatial and temporal distribution of energy release in the source area via the back-projection technique. Following this, we estimated the earthquake's intensity field by merging the source inversion findings with the ground motion prediction equation. This analysis facilitated the evaluation of earthquake casualties, utilizing the theoretical intensity field and a casualty assessment model. Our findings indicate that the fault type is a thrust fault, characterized by a unilateral rupture in the direction of NW, with a rupture length spanning approximately 10–15 km and a duration ranging between 8 and 10 s. The earthquake's magnitude varied from M 5.9 to M 6.2. The demarcated high-intensity areas, as per our intensity assessment, align closely with the actual survey results. Furthermore, the predicted total casualties and identified critical rescue zones closely match the real-world casualty figures. These insights offer crucial technical support for governmental emergency command and rescue operations.

Overview of China high-speed railway earthquake early-warning system

Amidst China's aggressive expansion of its high-speed rail network,the intersection of these lines with seismic fault zones has elevated the risk profile for high-speed rail travel.To counteract the potential dangers posed by seismic disturbances,China has introduced a comprehensive high-speed railway earthquake early-warning system.This article presents an in-depth examination of this system,encompassing aspects such as its developmental evolution,architectural design,and pivotal technologies.Furthermore,it ventures into the realm of future enhancements and developmental pathways for the system,fusing emergent findings from earthquake early warning research with advancements in artificial intelligence.

Rapid report of the December 18,2023 M_(S)6.2 Jishishan earthquake,Gansu,China

On December 18,2023,the Jishishan area in Gansu Province was jolted by a M_(S) 6.2 earthquake,which is the most powerful seismic event that occurred throughout the year in China.The earthquake occurred along the NWtrending Lajishan fault(LJSF),a large tectonic transformation zone.After this event,China Earthquake Networks Center(CENC)has timely published several reports about seismic sources for emergency responses.The earthquake early warning system issued the first alert 4.9 s after the earthquake occurrence,providing prompt notification that effectively mitigated panics,injuries,and deaths of residents.The near real-time focal mechanism solution indicates that this earthquake is associated with a thrust fault.The distribution of aftershocks,the rupture process,and the recorded amplitudes from seismic monitoring and GNSS stations,all suggest that the mainshock rupture predominately propagates to the northwest direction.The duration of the rupture process is~12 s,and the largest slip is located at approximately 6.3 km to the NNW from the epicenter,with a peak slip of 0.12 m at~8 km depth.Seismic station N0028 recorded the highest instrumental intensity,which is 9.4 on the Mercalli scale.The estimated intensity map shows a seismic intensity reaching up to IX near the rupture area,consistent with field survey results.The aftershocks(up to December 22,2023)are mostly distributed in the northwest direction within~20 km of the epicenter.This earthquake caused serious casualties and house collapses,which requires further investigations into the impact of this earthquake.