The Hualien M 7.3 earthquake on April 3,2024,was a significant and strong earthquake in Taiwan,China in the past two decades.The rupture process of the main shock and strong aftershocks is of great significance to the...The Hualien M 7.3 earthquake on April 3,2024,was a significant and strong earthquake in Taiwan,China in the past two decades.The rupture process of the main shock and strong aftershocks is of great significance to the subsequent seismic activity and seismogenic tectonic research.Based on local strong-motion data,we used the IDS(Iterative Deconvolution and Stacking)method to obtain the rupture process of the mainshock and two strong aftershocks on the 23rd.The rupture of the mainshock was mainly unilateral,lasting 31 s,with a maximum slip of 2m,and the depth of the large slip zone is about 41–49 km.There is a clear difference between the rupture depth of the main shock and the two strong aftershocks.The depths of the large slip zones of the latter two are 3–9 km and 8–10 km,respectively.There is also a significant difference in the seismogenic fault between the mainshock and the aftershocks,and we believe that there are two seismogenic fault zones in the study area,the deep and the shallow fault zone.The slip of the deep faults activates the shallow faults.展开更多
A great earthquake struck central Myanmar on March 28,2025,causing extensive and severe damage in Myanmar and neighboring countries such as Thailand.Observed waveforms indicated the possibility of supershear rupture a...A great earthquake struck central Myanmar on March 28,2025,causing extensive and severe damage in Myanmar and neighboring countries such as Thailand.Observed waveforms indicated the possibility of supershear rupture at the source of this event.To investigate this possibility,a joint source inversion using near-field and teleseismic waveforms was performed.The fault model used,with a total length of 580 km,comprises five segments with varying strike directions and dip angles.Within the first 10 s after initiation at the hypocenter,the rupture propagated southward at a high speed of 5–6 km/s.while rupture in the northern direction proceeded at a speed lower than the local S-wave velocity of 3.5 km/s.Subsequently,the rupture propagated bilaterally in both northern and southern directions at high speeds of 5–6 km/s.The largest asperity,with a maximum slip of 5–6 m,occurred in the shallow part of the fault model,approximately 240 km south of the hypocenter.Large slips of 2–3 m were also found near the southern end of the fault model.Slips on the fault plane were dominated by strike-slip components,though slips south of the largest asperity included significant dip-slip components.展开更多
On December 18,2023,an M_(s)6.2 earthquake jolted Jishishan County in the Linxia Hui Autonomous Prefecture in Northwest China's Gansu Province,causing substantial casualties and building collapses.The earthquake o...On December 18,2023,an M_(s)6.2 earthquake jolted Jishishan County in the Linxia Hui Autonomous Prefecture in Northwest China's Gansu Province,causing substantial casualties and building collapses.The earthquake occurred in the Qilian Block on the northeastern border of the Qinghai-Tibet Plateau,where faults are highly active and the geological structure is complex.In this study,we utilized methods such as relocation,focal mechanism solutions,and earthquake rupture processes to describe seismogenic faults.The results indicated that the majority of aftershocks occurred at a depth of 12 km.The centroid depth of the main shock and the depth of the maximum rupture point during the rupture process were also 12 km.Various geophysical methods exhibited a high degree of consistency in depth exploration.Aftershocks were distributed mainly to the west and north of the main shock and extended in the NNW direction,primarily through unilateral rupture.The main shock was a reverse thrust event with a small dextral strike-slip component.In this study,more regional data,such as previous GPS observations,field geological observations,and the distributions of the primary stress states in the region,were also incorporated.We inferred that the main shock was triggered by the main fault at the northern margin of the Lajishan Fault and that the movement of the main fault also activated some secondary faults.The compressive forces on both sides of the Lajishan Fault Zone led to the uplift of mountain areas,accompanied by some landslides,leading to this catastrophic earthquake event.In this article,the activity relationships among the 2022 M_(s)6.9 Menyuan earthquake,the 2019 M_(s)5.7 Xiahe earthquake,and the Jishishan earthquake under the action of regional stress are also discussed.This study provides additional evidence and new ideas for exploring the seismogenic process of the Lajishan Fault Zone and has implications for future in-depth research on underground activity in this region.展开更多
Focal mechanism and dynamic rupture process of the Wenchaun Ms8.0 earthquake in Sichuan province on 12 May 2008 were obtained by inverting long period seismic data from the Global Seismic Network (GSN), and characteri...Focal mechanism and dynamic rupture process of the Wenchaun Ms8.0 earthquake in Sichuan province on 12 May 2008 were obtained by inverting long period seismic data from the Global Seismic Network (GSN), and characteristics of the co-seismic displacement field near the fault were quantitatively ana-lyzed based on the inverted results to investigate the mechanism causing disaster. A finite fault model with given focal mechanism and vertical components of the long period P-waves from 21 stations with evenly azimuthal coverage were adopted in the inversion. From the inverted results as well as after-shock distribution, the causative fault of the great Wenchuan earthquake was confirmed to be a fault of strike 225°/dip 39°/rake 120°, indicating that the earthquake was mainly a thrust event with right-lateral strike-slip component. The released scalar seismic moment was estimated to be about 9.4×1020―2.0×1021 Nm, yielding moment magnitude of Mw7.9―8.1. The great Wenchuan earthquake occurred on a fault more than 300 km long, and had a complicated rupture process of about 90 s duration time. The slip distribution was highly inhomogeneous with the average slip of about 2.4 m. Four slip-patches broke the ground surface. Two of them were underneath the regions of Wenchuan-Yingxiu and Beichuan, respectively, with the first being around the hypocenter (rupture initiation point), where the largest slip was about 7.3 m, and the second being underneath Beichuan and extending to Pingwu, where the largest slip was about 5.6 m. The other two slip-patches had smaller sizes, one having the maximum slip of 1.8 m and lying underneath the north of Kangding, and the other having the maximum slip of 0.7 m and lying underneath the northeast of Qingchuan. Average and maximum stress drops over the whole fault plane were estimated to be 18 MPa and 53 MPa, respectively. In addition, the co-seismic displacement field near the fault was analyzed. The results indicate that the features of the co-seismic displacement field were coincident with those of the intensity distribution in the meizo-seismal area, implying that the large-scale, large-amplitude and surface-broken thrust dislocation should be responsible for the serious disaster in the near fault area.展开更多
This paper introduces techniques in Gaussian process regression model for spatiotemporal data collected from complex systems.This study focuses on extracting local structures and then constructing surrogate models bas...This paper introduces techniques in Gaussian process regression model for spatiotemporal data collected from complex systems.This study focuses on extracting local structures and then constructing surrogate models based on Gaussian process assumptions.The proposed Dynamic Gaussian Process Regression(DGPR)consists of a sequence of local surrogate models related to each other.In DGPR,the time-based spatial clustering is carried out to divide the systems into sub-spatio-temporal parts whose interior has similar variation patterns,where the temporal information is used as the prior information for training the spatial-surrogate model.The DGPR is robust and especially suitable for the loosely coupled model structure,also allowing for parallel computation.The numerical results of the test function show the effectiveness of DGPR.Furthermore,the shock tube problem is successfully approximated under different phenomenon complexity.展开更多
In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back proj...In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.展开更多
An earthquake of Ms=7.4 occurred in Mani, Xizang (Tibet), China on November 8, 1997. The moment tensor ofthis earthquake was inverted using the long period body wave form data from China Digital Seismograph Network(CD...An earthquake of Ms=7.4 occurred in Mani, Xizang (Tibet), China on November 8, 1997. The moment tensor ofthis earthquake was inverted using the long period body wave form data from China Digital Seismograph Network(CDSN). The apparent source time functions (AS TFs) were retrieved from P and S waves, respectively, using thedeconvolution technique in frequency domain, and the tempo-spatial rupture process on the fault plane was imagedby inverting the azimuth dependent AS TFs from different stations. The result of the moment tensor inversionindicates that the P and T axes of earthquake-generating stress field were nearly horizontal, with the P axis in theNNE direction (29), the T axis in the SEE direction (122) and that the NEE-SWW striking nodal plane andNNW-SSE striking nodal plane are mainly left-lateral and right-lateral strike-slip, respectively; that this earthquakehad a scalar seismic moment of 3.4xl02o N. .m, and a moment magnitude of Mw=7.6. Taking the aftershock distribution into account, we proposed that the earthquake rupture occurred in the fault plane with the strike of 250,the dip of 88 and the rake of 19. On the basis of the result of the moment tensor inversion, the theoretical seismograms were synthesized, and then the AS T Fs were retrieved by deconvoving the synthetic seismograms fromthe observed seismograms. The A S T Fs retrieved from the P and S waves of different stations identically suggestedthat this earthquake was of a simple time history, whose ASTF can be approximated with a sine function with thehalf period of about 10 s. Inverting the azimuth dependent A S T Fs from P and S waveforms led to the imageshowing the tempo-spatial distribution of the rupture on the fault plane. From the 'remembering' snap-shots, therupture initiated at the western end of the fault, and then propagated eastward and downward, indicating an overallunilateral rupture. However, the slip distribution is non-uniform, being made up of three sub-areas, one in thewestern end, about 10 km deep ('western area'), another about 55 kin away from the western end and about 35 Iondeep ('eastern area'), the third about 30 km away from the western end and around 40 km deep ('central area').The total rupture area was around 70 km long and 60 km wide. From the 'forgetting' snap-shots, the rupturingappeared quite complex, with the slip occurring in different position at different time, and the earthquake being ofthe characteristics of 'healing pulse'. Another point we have to stress is that the locations in which the ruptureinitiated and terminated were not where the main rupture took place. Eventually, the static slip distribution wascalculated, and the largest slip values of the three sub-areas were 956 cm, 743 cm and 1 060 cm, for the western.eastern and central areas, respectively. From the slip distribution, the rupture mainly distributed in the fault about70 km eastern to the epicenter; from the aftershock distribution. however, the aftershocks were very sparse in thewest to the epicenter while densely clustered in the east to the epicenter It indicated that the Maul Ms=7.9 earthquake was resulted from the nearly eastward extension of the NEE-SWW to nearly E-W striking fault in thenorthwestern Tibetan plateau.展开更多
Fast inversion of source rupture process for significant earthquake is introduced and applied to earthquake emergency response. We retrieve and use the apparent source time functions （ASTFs） to stabilize the rupture...Fast inversion of source rupture process for significant earthquake is introduced and applied to earthquake emergency response. We retrieve and use the apparent source time functions （ASTFs） to stabilize the rupture process inversion, aiming to get the rupture process results correctly and quickly. The applications to 22 significant earthquakes occurred between January 2009 and July 2011 show an average inversion time consumed of about four hours, suggesting a high efficiency of this method, which much improves the inversion speed of the rupture process and is helpful in earthquake disaster mitigation and relief.展开更多
Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s...Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s along the fault strike direction. The source rupture process consists of three sub-events, the first oc- curred near the hypocenter and the rest two ruptured along the up-dip direction and broke the sea bed, causing a maximum slip of about 30 m. The large-scale sea bed breakage may account for the tremendous tsunami disaster which resulted in most of the death and missing in this mega earthquake.展开更多
According to the China Earthquake Networks Center,a strong earthquake of M6.8 occurred in Luding County,Ganzi Tibetan Autonomous Prefecture,Sichuan Province,China(102.08°E,29.59°N),on September 5,2022,with a...According to the China Earthquake Networks Center,a strong earthquake of M6.8 occurred in Luding County,Ganzi Tibetan Autonomous Prefecture,Sichuan Province,China(102.08°E,29.59°N),on September 5,2022,with a focal depth of 16 km.Rapid determination of the source parameters of the earthquake sequence is vital for post-earthquake rescue,disaster assessment,and scientific research.Near-field seismic observations play a key role in the fast and reliable determination of earthquake source parameters.The numerous broadband seismic stations and strong-motion stations recently deployed by the National Earthquake Intensity Rapid Report and Early Warning project have provided valuable real-time near-field observation data.Using these near-field observations and conventional mid-and far-field seismic waveform records,we obtained the focal mechanism solutions of the mainshock and M≥3.0 aftershocks through the waveform fitting method.We were further able to rapidly invert the rupture process of the mainshock.Based on the evaluation of the focal mechanism solution of the mainshock and the regional tectonic setting,we speculate that the Xianshuihe fault formed the seismogenic structure of the M6.8 strong earthquake.The aftershocks formed three spatially separated clusters with distinctly different focal mechanisms,reflecting the segmented nature of the Xianshuihe fault.As more high-frequency information has been applied in this study,the absolute location of the fault rupture is better constrained by the near-field strong-motion data.The rupture process of the mainshock correlates well with the spatial distribution of aftershocks,i.e.,aftershock activities were relatively weak in the maximum slip area,and strong aftershock activities were distributed in the peripheral regions.展开更多
The spatial and temporal slip distribution of the Lushan earthquake was estimated using teleseismic body wave data. To perform a stable inversion, we applied smoothing constraints and determined their optimal relative...The spatial and temporal slip distribution of the Lushan earthquake was estimated using teleseismic body wave data. To perform a stable inversion, we applied smoothing constraints and determined their optimal relative weights on the observed data using an optimized Akaike＇ s Bayesian Information Criterion （ABIC）. The inversion generated the source parameters. Strike, dip and slip were 218°, 39° and 100. 8° ,respectively. A seismic moment （M0） was 2. 1 × 10^20 Nm with a moment magnitude （Mw） of 6. 8, and a source duration was approximately 30 second. The rupture propagated along the dip direction, and the maximum slip occurred at the hypocenter. The maximum slip was approximately 2. 1 m, although this earthquake did not cause an apparent surface rupture. The energy was mainly released within 10 second. In addition, the Lushan earthquake was apparently related to the 2008 Wenchuan earthquake. However, the question of whether it was an aftershock of the Wenchuan earthquake requires further study.展开更多
In this study we performed a classical spectrum analysis of seismic waveforms recorded at far field stations of the great MW7.9 Wenchuan earthquake to observe the shifts of the corner frequency with azimuth due to the...In this study we performed a classical spectrum analysis of seismic waveforms recorded at far field stations of the great MW7.9 Wenchuan earthquake to observe the shifts of the corner frequency with azimuth due to the Doppler effect.Our results show that this damaging great earthquake had a dominating rupture propagation direction of 64.0°.The equivalent radius of the fault rupture surface was estimated to be 33 km,yielding the rupture area of about 3 500 km2.Thus the length of the rupture fault surface is about 230 km if the depth（or width） extent is 15 km.The computer program developed in this study can quickly provide the information about the source of a future large（damaging） earthquake,which could be very useful for predicting aftershocks and planning the rescue operations.展开更多
Four results of the rupture process of 14 April 2010 Yushu, Qinghai, earthquake, obtained by inverting the broadband seismic data of Global Seismographic Network （GSN） based on the inversion method of earthquake rup...Four results of the rupture process of 14 April 2010 Yushu, Qinghai, earthquake, obtained by inverting the broadband seismic data of Global Seismographic Network （GSN） based on the inversion method of earthquake rupture process, were compared and discussed. It is found that the Yushu earthquake has several basic characteristics as follows： 1 There exist two principal sub-events which correspond to two slip-concentrated patches being located near the hypocenter and to the southeast of the epicenter. The rupture of the slip-concentrated patch to the southeast of the epicenter broke though the ground surface; 2 The peak slip and peak slip-rate are about 2.1 m and 1.1 m/s, respectively, indicating that the Yushu earthquake is an event with large slip-rate on the fault plane; 3 Overall the Yushu earthquake is a unilateral rupture event with the rupture mainly propagating southeastward. The strong focusing of the seismic energy in the southeast of the epicenter due to the ＂seismic Doppler effect＂ reasonably accounts for the tremendous damage in the Yushu city.展开更多
Joint inversion of teleseismic body-wave data and strong ground motion waveforms was applied to determine the rupture process of the 2010 Mentawai earthquake. To obtain stable solutions, smoothing and non-negative con...Joint inversion of teleseismic body-wave data and strong ground motion waveforms was applied to determine the rupture process of the 2010 Mentawai earthquake. To obtain stable solutions, smoothing and non-negative constraints were introduced. A total of 33 teleseismic stations and 5 strong ground motion stations supplied data. The teleseismic and strong ground motion data were separately windowed for 150 s and 250 s and bandpass filtered with frequencies of 0.001e1.0 Hz and 0.005e0.5 Hz, respectively. The finitefault model was established with length and width of 190 km and 70 km, and the initial seismic source parameters were set by referring to centroid moment tensor(CMT) solutions. Joint inversion results indicate that the focal mechanism of this earthquake is thrust fault type, and the strike, dip, and rake angles are generally in accordance with CMT results. The seismic moment was determined as 5.814 1020Nm(Mw7.8) and source duration was about 102 s, which is greater than those of other earthquakes of similar magnitude. The rupture nucleated near the hypocenter and then propagated along the strike direction to the northwest, with a maximum slip of 3.9 m. Large uncertainties regarding the amount of slip retrieved using different inversion methods still exist; however, the conclusion that the majority of slip occurred far from the islands at very shallow depths was found to be robust. The 2010 Mentawai earthquake was categorized as a tsunami earthquake because of the long rupture duration and the generation of a tsunami much larger than was expected for an earthquake of its magnitude.展开更多
On April 20, 2013, an M_s 7.0 earthquake struck Lushan County in Sichuan Province, China, and caused serious damage to the source region. We investigated the rupture process of the M_s7.0 Lushan earthquake by jointly ...On April 20, 2013, an M_s 7.0 earthquake struck Lushan County in Sichuan Province, China, and caused serious damage to the source region. We investigated the rupture process of the M_s7.0 Lushan earthquake by jointly inverting waveforms of teleseismic P waveforms and local strong motion records as well as static GPS observations. The inverted results indicate that the rupture of this earthquake was dominated by the failure of an asperity with a triangular shape and that the main shock was dominated by thrust slip. The earthquake released a total seismic moment of 1.01× 10^（19）Nm, with 92% of it being released during the first 11 s. The rupture had an average slip of 0.9 m and produced an average stress drop of 1.8 MPa. Compared with our previous work that was based mainly on a unique dataset, this joint inversion result is more consistent with field observations and the distribution of aftershock zones.展开更多
On 25 April, 2015, an Mw7.9 earthquake occurred in Nepal, which caused great economic loss and casualties. However, almost no surface ruptures were observed. Therefore, in order to interpret the phenomenon, we study t...On 25 April, 2015, an Mw7.9 earthquake occurred in Nepal, which caused great economic loss and casualties. However, almost no surface ruptures were observed. Therefore, in order to interpret the phenomenon, we study the rupture process of the earthquake to seek answers. Inversion of teleseismic body-wave data is applied to estimate the rupture process of the 2015 Nepal earthquake. To obtain stable solutions, smoothing and non-negative constraints are introduced. 48 teleseismic stations with good coverage are chosen. Finite fault model is established with length and width of 195 km and 150 km, and we set the initial seismic source parameters referring to CMT solutions. Inversion results indicate that the focal mechanism of this earthquake is a thrust fault type, and the strike, dip and rake angle are in accordance with CMT results. The seismic moment is 0.9195 ×10^（21）Nm（Mw7.9）, and source duration is about 70s. The rupture nucleated near the hypocenter and then propagated along the dip direction to the southeast, and the maximum slip amounts to 5.2 m. Uncertainties on the amount of slip retrieved by different inversion methods still exist, the overall characteristics are inconsistent. The lack of shallow slip during the 2015 Gorkha earthquake implies future seismic hazard and this region should be paid more attention to.展开更多
The source parameters, such as moment tensor, focal mechanism, source time function (STF) and temporal-spatial rupture process, were obtained for the January 26, 2001, India, MS=7.8 earthquake by inverting waveform da...The source parameters, such as moment tensor, focal mechanism, source time function (STF) and temporal-spatial rupture process, were obtained for the January 26, 2001, India, MS=7.8 earthquake by inverting waveform data of 27 GDSN stations with epicentral distances less than 90? Firstly, combining the moment tensor inversion, the spatial distribution of intensity, disaster and aftershocks and the orientation of the fault where the earthquake lies, the strike, dip and rake of the seismogenic fault were determined to be 92? 58?and 62? respectively. That is, this earthquake was a mainly thrust faulting with the strike of near west-east and the dipping direction to south. The seismic moment released was 3.51020 Nm, accordingly, the moment magnitude MW was calculated to be 7.6. And then, 27 P-STFs, 22 S-STFs and the averaged STFs of them were determined respectively using the technique of spectra division in frequency domain and the synthetic seismogram as Greens functions. The analysis of the STFs suggested that the earthquake was a continuous event with the duration time of 19 s, starting rapidly and ending slowly. Finally, the temporal-spatial distribution of the slip on the fault plane was imaged from the obtained P-STFs and S-STFs using an time domain inversion technique. The maximum slip amplitude on the fault plane was about 7 m. The maximum stress drop was 30 MPa, and the average one over the whole rupture area was 7 MPa. The rupture area was about 85 km long in the strike direction and about 60 km wide in the down-dip direction, which, equally, was 51 km deep in the depth direction. The rupture propagated 50 km eastwards and 35 km westwards. The main portion of the rupture area, which has the slip amplitude greater than 0.5 m, was of the shape of an ellipse, its major axis oriented in the slip direction of the fault, which indicated that the rupture propagation direction was in accordance with the fault slip direction. This phenomenon is popular for strike-slip faulting, but rather rare for thrust faulting. The eastern portion of the rupture area above the initiation point was larger than the western portion below the initiation point, which was indicative of the asymmetrical rupture. In other words, the rupturing was kind of unilateral from west to east and from down to up. From the snapshots of the slip-rate variation with time and space, the slip rate reached the largest at the 4th second, that was 0.2 m/s, and the rupture in this period occurred only around the initiation point. At the 6th second, the rupture around the initiation point nearly stopped, and started moving outwards. The velocity of the westward rupture was smaller than that of the eastward rupture. Such rupture behavior like a circle mostly stopped near the 15th second. After the 16th second, only some patches of rupture distributed in the outer region. From the snapshots of the slip variation with time and space, the rupture started at the initiation point and propagated outwards. The main rupture on the area with the slip amplitude greater than 5 m extended unilaterally from west to east and from down to up between the 6th and the 10th seconds, and the western segment extended a bit westwards and downwards between the 11th and the 13th seconds. The whole process lasted about 19 s. The rupture velocity over the whole rupture process was estimated to be 3.3 km/s.展开更多
On July 20, 1995, an earthquake of ML=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 ...On July 20, 1995, an earthquake of ML=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 WL=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 ML=2.0 earthquake as an empirical Green's function, a regularization method was applied toretrieve the far-field source-time fonction (STF) of the main shock. Considering the records of HDSN are the typeof 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 fieldstations, P phases in vertical direction which mostly are about 0.5 s in length were used. The STFs obtained fromeach seismic phases are in good agreement, showing that the Mt=4.1 eedhquake consisted of two events. STFsfrom each station demonstrate an obvious 'seismic Doppler effect'. Assuming the nodal plane striking 37° anddipping 40°, determined by using P wave first motion data and aftershock distribution, is the fault plane, through atrial and error method, the following results were drawn：Both of the events lasted about 0. 1s, the ruptUre length ofthe 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.0km/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 thefirst 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 ML=4.1 earthquake were also reached： the scalar earthquake moment is 3.3×1013 N·m, stress drop 4.6 MPa, rupture radius 0. 16 km.展开更多
The 10 January 2018 Mw7.5 Swan island,Honduras earthquake occurred on the Swan island fault,which is a transform plate boundary between the North American and Caribbean plates.Here we back-project the rupture process ...The 10 January 2018 Mw7.5 Swan island,Honduras earthquake occurred on the Swan island fault,which is a transform plate boundary between the North American and Caribbean plates.Here we back-project the rupture process of the earthquake using dense seismic stations in Alaska,and find that the earthquake ruptured at least three faults(three stages)for a duration of~40 s.The rupture speed for the longest fault(stage 3)is as fast as 5 km/s,which is much faster than the local shear wave velocity of~4 km/s.Supershear rupture was incidentally observed on long and straight strike-slip faults.This study shows a supershear rupture that occured on a strike-slip fault with moderate length,implying that supershear rupture might commonly occur on large strike-slip earthquakes.The common occurrence of supershear rupture on strike-slip earthquakes will challenge present understanding of crack physics,as well as strong ground motion evaluation in earthquake engineering.展开更多
基金sponsored by the Earthquake Spark Technology Project(XH23051B)。
文摘The Hualien M 7.3 earthquake on April 3,2024,was a significant and strong earthquake in Taiwan,China in the past two decades.The rupture process of the main shock and strong aftershocks is of great significance to the subsequent seismic activity and seismogenic tectonic research.Based on local strong-motion data,we used the IDS(Iterative Deconvolution and Stacking)method to obtain the rupture process of the mainshock and two strong aftershocks on the 23rd.The rupture of the mainshock was mainly unilateral,lasting 31 s,with a maximum slip of 2m,and the depth of the large slip zone is about 41–49 km.There is a clear difference between the rupture depth of the main shock and the two strong aftershocks.The depths of the large slip zones of the latter two are 3–9 km and 8–10 km,respectively.There is also a significant difference in the seismogenic fault between the mainshock and the aftershocks,and we believe that there are two seismogenic fault zones in the study area,the deep and the shallow fault zone.The slip of the deep faults activates the shallow faults.
基金Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration under Grant No.2024B15National Natural Science Fundation of China under Grant Nos.U2239252 and 52478567Natural Science Foundation of Heilongjiang Province under Grant No.JQ2023E002。
文摘A great earthquake struck central Myanmar on March 28,2025,causing extensive and severe damage in Myanmar and neighboring countries such as Thailand.Observed waveforms indicated the possibility of supershear rupture at the source of this event.To investigate this possibility,a joint source inversion using near-field and teleseismic waveforms was performed.The fault model used,with a total length of 580 km,comprises five segments with varying strike directions and dip angles.Within the first 10 s after initiation at the hypocenter,the rupture propagated southward at a high speed of 5–6 km/s.while rupture in the northern direction proceeded at a speed lower than the local S-wave velocity of 3.5 km/s.Subsequently,the rupture propagated bilaterally in both northern and southern directions at high speeds of 5–6 km/s.The largest asperity,with a maximum slip of 5–6 m,occurred in the shallow part of the fault model,approximately 240 km south of the hypocenter.Large slips of 2–3 m were also found near the southern end of the fault model.Slips on the fault plane were dominated by strike-slip components,though slips south of the largest asperity included significant dip-slip components.
基金funded by the National Natural Science Foundation of China(Grant No.42304072)。
文摘On December 18,2023,an M_(s)6.2 earthquake jolted Jishishan County in the Linxia Hui Autonomous Prefecture in Northwest China's Gansu Province,causing substantial casualties and building collapses.The earthquake occurred in the Qilian Block on the northeastern border of the Qinghai-Tibet Plateau,where faults are highly active and the geological structure is complex.In this study,we utilized methods such as relocation,focal mechanism solutions,and earthquake rupture processes to describe seismogenic faults.The results indicated that the majority of aftershocks occurred at a depth of 12 km.The centroid depth of the main shock and the depth of the maximum rupture point during the rupture process were also 12 km.Various geophysical methods exhibited a high degree of consistency in depth exploration.Aftershocks were distributed mainly to the west and north of the main shock and extended in the NNW direction,primarily through unilateral rupture.The main shock was a reverse thrust event with a small dextral strike-slip component.In this study,more regional data,such as previous GPS observations,field geological observations,and the distributions of the primary stress states in the region,were also incorporated.We inferred that the main shock was triggered by the main fault at the northern margin of the Lajishan Fault and that the movement of the main fault also activated some secondary faults.The compressive forces on both sides of the Lajishan Fault Zone led to the uplift of mountain areas,accompanied by some landslides,leading to this catastrophic earthquake event.In this article,the activity relationships among the 2022 M_(s)6.9 Menyuan earthquake,the 2019 M_(s)5.7 Xiahe earthquake,and the Jishishan earthquake under the action of regional stress are also discussed.This study provides additional evidence and new ideas for exploring the seismogenic process of the Lajishan Fault Zone and has implications for future in-depth research on underground activity in this region.
基金Supported by the National Basic Research Program of China (Grant No. 2004CB418404-4)the National Natural Science Foundation of China (Grant Nos. 40574025 and 40874026)
文摘Focal mechanism and dynamic rupture process of the Wenchaun Ms8.0 earthquake in Sichuan province on 12 May 2008 were obtained by inverting long period seismic data from the Global Seismic Network (GSN), and characteristics of the co-seismic displacement field near the fault were quantitatively ana-lyzed based on the inverted results to investigate the mechanism causing disaster. A finite fault model with given focal mechanism and vertical components of the long period P-waves from 21 stations with evenly azimuthal coverage were adopted in the inversion. From the inverted results as well as after-shock distribution, the causative fault of the great Wenchuan earthquake was confirmed to be a fault of strike 225°/dip 39°/rake 120°, indicating that the earthquake was mainly a thrust event with right-lateral strike-slip component. The released scalar seismic moment was estimated to be about 9.4×1020―2.0×1021 Nm, yielding moment magnitude of Mw7.9―8.1. The great Wenchuan earthquake occurred on a fault more than 300 km long, and had a complicated rupture process of about 90 s duration time. The slip distribution was highly inhomogeneous with the average slip of about 2.4 m. Four slip-patches broke the ground surface. Two of them were underneath the regions of Wenchuan-Yingxiu and Beichuan, respectively, with the first being around the hypocenter (rupture initiation point), where the largest slip was about 7.3 m, and the second being underneath Beichuan and extending to Pingwu, where the largest slip was about 5.6 m. The other two slip-patches had smaller sizes, one having the maximum slip of 1.8 m and lying underneath the north of Kangding, and the other having the maximum slip of 0.7 m and lying underneath the northeast of Qingchuan. Average and maximum stress drops over the whole fault plane were estimated to be 18 MPa and 53 MPa, respectively. In addition, the co-seismic displacement field near the fault was analyzed. The results indicate that the features of the co-seismic displacement field were coincident with those of the intensity distribution in the meizo-seismal area, implying that the large-scale, large-amplitude and surface-broken thrust dislocation should be responsible for the serious disaster in the near fault area.
基金co-supported by the National Natural Science Foundation of China(No.12101608)the NSAF(No.U2230208)the Hunan Provincial Innovation Foundation for Postgraduate,China(No.CX20220034).
文摘This paper introduces techniques in Gaussian process regression model for spatiotemporal data collected from complex systems.This study focuses on extracting local structures and then constructing surrogate models based on Gaussian process assumptions.The proposed Dynamic Gaussian Process Regression(DGPR)consists of a sequence of local surrogate models related to each other.In DGPR,the time-based spatial clustering is carried out to divide the systems into sub-spatio-temporal parts whose interior has similar variation patterns,where the temporal information is used as the prior information for training the spatial-surrogate model.The DGPR is robust and especially suitable for the loosely coupled model structure,also allowing for parallel computation.The numerical results of the test function show the effectiveness of DGPR.Furthermore,the shock tube problem is successfully approximated under different phenomenon complexity.
基金supported by the National Key R&D Program of China(No.2022YFF0800601)National Scientific Foundation of China(Nos.41930103 and 41774047).
文摘In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.
文摘An earthquake of Ms=7.4 occurred in Mani, Xizang (Tibet), China on November 8, 1997. The moment tensor ofthis earthquake was inverted using the long period body wave form data from China Digital Seismograph Network(CDSN). The apparent source time functions (AS TFs) were retrieved from P and S waves, respectively, using thedeconvolution technique in frequency domain, and the tempo-spatial rupture process on the fault plane was imagedby inverting the azimuth dependent AS TFs from different stations. The result of the moment tensor inversionindicates that the P and T axes of earthquake-generating stress field were nearly horizontal, with the P axis in theNNE direction (29), the T axis in the SEE direction (122) and that the NEE-SWW striking nodal plane andNNW-SSE striking nodal plane are mainly left-lateral and right-lateral strike-slip, respectively; that this earthquakehad a scalar seismic moment of 3.4xl02o N. .m, and a moment magnitude of Mw=7.6. Taking the aftershock distribution into account, we proposed that the earthquake rupture occurred in the fault plane with the strike of 250,the dip of 88 and the rake of 19. On the basis of the result of the moment tensor inversion, the theoretical seismograms were synthesized, and then the AS T Fs were retrieved by deconvoving the synthetic seismograms fromthe observed seismograms. The A S T Fs retrieved from the P and S waves of different stations identically suggestedthat this earthquake was of a simple time history, whose ASTF can be approximated with a sine function with thehalf period of about 10 s. Inverting the azimuth dependent A S T Fs from P and S waveforms led to the imageshowing the tempo-spatial distribution of the rupture on the fault plane. From the 'remembering' snap-shots, therupture initiated at the western end of the fault, and then propagated eastward and downward, indicating an overallunilateral rupture. However, the slip distribution is non-uniform, being made up of three sub-areas, one in thewestern end, about 10 km deep ('western area'), another about 55 kin away from the western end and about 35 Iondeep ('eastern area'), the third about 30 km away from the western end and around 40 km deep ('central area').The total rupture area was around 70 km long and 60 km wide. From the 'forgetting' snap-shots, the rupturingappeared quite complex, with the slip occurring in different position at different time, and the earthquake being ofthe characteristics of 'healing pulse'. Another point we have to stress is that the locations in which the ruptureinitiated and terminated were not where the main rupture took place. Eventually, the static slip distribution wascalculated, and the largest slip values of the three sub-areas were 956 cm, 743 cm and 1 060 cm, for the western.eastern and central areas, respectively. From the slip distribution, the rupture mainly distributed in the fault about70 km eastern to the epicenter; from the aftershock distribution. however, the aftershocks were very sparse in thewest to the epicenter while densely clustered in the east to the epicenter It indicated that the Maul Ms=7.9 earthquake was resulted from the nearly eastward extension of the NEE-SWW to nearly E-W striking fault in thenorthwestern Tibetan plateau.
基金supported by the National Natural Science Foundation of China (No. 41090291)the Research Project in Earthquake Science, CEA (Nos. 201108002 and 200808068)
文摘Fast inversion of source rupture process for significant earthquake is introduced and applied to earthquake emergency response. We retrieve and use the apparent source time functions （ASTFs） to stabilize the rupture process inversion, aiming to get the rupture process results correctly and quickly. The applications to 22 significant earthquakes occurred between January 2009 and July 2011 show an average inversion time consumed of about four hours, suggesting a high efficiency of this method, which much improves the inversion speed of the rupture process and is helpful in earthquake disaster mitigation and relief.
基金financially supported by the National Natural Science Foundation of China (Nos. 90915012 and 41090291)the Research Project in Earthquake Science, CEA (No.201108002)
文摘Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s along the fault strike direction. The source rupture process consists of three sub-events, the first oc- curred near the hypocenter and the rest two ruptured along the up-dip direction and broke the sea bed, causing a maximum slip of about 30 m. The large-scale sea bed breakage may account for the tremendous tsunami disaster which resulted in most of the death and missing in this mega earthquake.
基金supported by the China Spark Program of Earthquake Science and Technology(No.XH23051B)National Key R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(No.2017YFC1500304)。
文摘According to the China Earthquake Networks Center,a strong earthquake of M6.8 occurred in Luding County,Ganzi Tibetan Autonomous Prefecture,Sichuan Province,China(102.08°E,29.59°N),on September 5,2022,with a focal depth of 16 km.Rapid determination of the source parameters of the earthquake sequence is vital for post-earthquake rescue,disaster assessment,and scientific research.Near-field seismic observations play a key role in the fast and reliable determination of earthquake source parameters.The numerous broadband seismic stations and strong-motion stations recently deployed by the National Earthquake Intensity Rapid Report and Early Warning project have provided valuable real-time near-field observation data.Using these near-field observations and conventional mid-and far-field seismic waveform records,we obtained the focal mechanism solutions of the mainshock and M≥3.0 aftershocks through the waveform fitting method.We were further able to rapidly invert the rupture process of the mainshock.Based on the evaluation of the focal mechanism solution of the mainshock and the regional tectonic setting,we speculate that the Xianshuihe fault formed the seismogenic structure of the M6.8 strong earthquake.The aftershocks formed three spatially separated clusters with distinctly different focal mechanisms,reflecting the segmented nature of the Xianshuihe fault.As more high-frequency information has been applied in this study,the absolute location of the fault rupture is better constrained by the near-field strong-motion data.The rupture process of the mainshock correlates well with the spatial distribution of aftershocks,i.e.,aftershock activities were relatively weak in the maximum slip area,and strong aftershock activities were distributed in the peripheral regions.
基金jointly supported by the Director of the Foundation of the Institute of Seismology,China Earthquake Administration(IS201102643)the National Natural Science Foundation of China(41004020)
文摘The spatial and temporal slip distribution of the Lushan earthquake was estimated using teleseismic body wave data. To perform a stable inversion, we applied smoothing constraints and determined their optimal relative weights on the observed data using an optimized Akaike＇ s Bayesian Information Criterion （ABIC）. The inversion generated the source parameters. Strike, dip and slip were 218°, 39° and 100. 8° ,respectively. A seismic moment （M0） was 2. 1 × 10^20 Nm with a moment magnitude （Mw） of 6. 8, and a source duration was approximately 30 second. The rupture propagated along the dip direction, and the maximum slip occurred at the hypocenter. The maximum slip was approximately 2. 1 m, although this earthquake did not cause an apparent surface rupture. The energy was mainly released within 10 second. In addition, the Lushan earthquake was apparently related to the 2008 Wenchuan earthquake. However, the question of whether it was an aftershock of the Wenchuan earthquake requires further study.
文摘In this study we performed a classical spectrum analysis of seismic waveforms recorded at far field stations of the great MW7.9 Wenchuan earthquake to observe the shifts of the corner frequency with azimuth due to the Doppler effect.Our results show that this damaging great earthquake had a dominating rupture propagation direction of 64.0°.The equivalent radius of the fault rupture surface was estimated to be 33 km,yielding the rupture area of about 3 500 km2.Thus the length of the rupture fault surface is about 230 km if the depth（or width） extent is 15 km.The computer program developed in this study can quickly provide the information about the source of a future large（damaging） earthquake,which could be very useful for predicting aftershocks and planning the rescue operations.
基金supported by China Postdoctoral Science Foundation funded project (20080440435)the project (DQJB09B06) from Institute of Geophysics (IGP),China Earthquake Administration (CEA). Contribution No. is 10FE3002, IGP-CEA
文摘Four results of the rupture process of 14 April 2010 Yushu, Qinghai, earthquake, obtained by inverting the broadband seismic data of Global Seismographic Network （GSN） based on the inversion method of earthquake rupture process, were compared and discussed. It is found that the Yushu earthquake has several basic characteristics as follows： 1 There exist two principal sub-events which correspond to two slip-concentrated patches being located near the hypocenter and to the southeast of the epicenter. The rupture of the slip-concentrated patch to the southeast of the epicenter broke though the ground surface; 2 The peak slip and peak slip-rate are about 2.1 m and 1.1 m/s, respectively, indicating that the Yushu earthquake is an event with large slip-rate on the fault plane; 3 Overall the Yushu earthquake is a unilateral rupture event with the rupture mainly propagating southeastward. The strong focusing of the seismic energy in the southeast of the epicenter due to the ＂seismic Doppler effect＂ reasonably accounts for the tremendous damage in the Yushu city.
基金supported by National Natural Science Foundation of China (41304046)
文摘Joint inversion of teleseismic body-wave data and strong ground motion waveforms was applied to determine the rupture process of the 2010 Mentawai earthquake. To obtain stable solutions, smoothing and non-negative constraints were introduced. A total of 33 teleseismic stations and 5 strong ground motion stations supplied data. The teleseismic and strong ground motion data were separately windowed for 150 s and 250 s and bandpass filtered with frequencies of 0.001e1.0 Hz and 0.005e0.5 Hz, respectively. The finitefault model was established with length and width of 190 km and 70 km, and the initial seismic source parameters were set by referring to centroid moment tensor(CMT) solutions. Joint inversion results indicate that the focal mechanism of this earthquake is thrust fault type, and the strike, dip, and rake angles are generally in accordance with CMT results. The seismic moment was determined as 5.814 1020Nm(Mw7.8) and source duration was about 102 s, which is greater than those of other earthquakes of similar magnitude. The rupture nucleated near the hypocenter and then propagated along the strike direction to the northwest, with a maximum slip of 3.9 m. Large uncertainties regarding the amount of slip retrieved using different inversion methods still exist; however, the conclusion that the majority of slip occurred far from the islands at very shallow depths was found to be robust. The 2010 Mentawai earthquake was categorized as a tsunami earthquake because of the long rupture duration and the generation of a tsunami much larger than was expected for an earthquake of its magnitude.
基金supported by a grant from the Chinese Earthquake Administration (No.201308013)the National Natural Science Foundation of China (Nos.41604057, 40974034, and 41021003)a key project from the Institute of Geodesy and Geophysics
文摘On April 20, 2013, an M_s 7.0 earthquake struck Lushan County in Sichuan Province, China, and caused serious damage to the source region. We investigated the rupture process of the M_s7.0 Lushan earthquake by jointly inverting waveforms of teleseismic P waveforms and local strong motion records as well as static GPS observations. The inverted results indicate that the rupture of this earthquake was dominated by the failure of an asperity with a triangular shape and that the main shock was dominated by thrust slip. The earthquake released a total seismic moment of 1.01× 10^（19）Nm, with 92% of it being released during the first 11 s. The rupture had an average slip of 0.9 m and produced an average stress drop of 1.8 MPa. Compared with our previous work that was based mainly on a unique dataset, this joint inversion result is more consistent with field observations and the distribution of aftershock zones.
基金supported by National Natural Science Foundation of China (41304046)
文摘On 25 April, 2015, an Mw7.9 earthquake occurred in Nepal, which caused great economic loss and casualties. However, almost no surface ruptures were observed. Therefore, in order to interpret the phenomenon, we study the rupture process of the earthquake to seek answers. Inversion of teleseismic body-wave data is applied to estimate the rupture process of the 2015 Nepal earthquake. To obtain stable solutions, smoothing and non-negative constraints are introduced. 48 teleseismic stations with good coverage are chosen. Finite fault model is established with length and width of 195 km and 150 km, and we set the initial seismic source parameters referring to CMT solutions. Inversion results indicate that the focal mechanism of this earthquake is a thrust fault type, and the strike, dip and rake angle are in accordance with CMT results. The seismic moment is 0.9195 ×10^（21）Nm（Mw7.9）, and source duration is about 70s. The rupture nucleated near the hypocenter and then propagated along the dip direction to the southeast, and the maximum slip amounts to 5.2 m. Uncertainties on the amount of slip retrieved by different inversion methods still exist, the overall characteristics are inconsistent. The lack of shallow slip during the 2015 Gorkha earthquake implies future seismic hazard and this region should be paid more attention to.
基金973 Project (G1998040705) from Ministry of Science and Technology, P. R. China the National Science Foundation of China under grant No.49904004.
文摘The source parameters, such as moment tensor, focal mechanism, source time function (STF) and temporal-spatial rupture process, were obtained for the January 26, 2001, India, MS=7.8 earthquake by inverting waveform data of 27 GDSN stations with epicentral distances less than 90? Firstly, combining the moment tensor inversion, the spatial distribution of intensity, disaster and aftershocks and the orientation of the fault where the earthquake lies, the strike, dip and rake of the seismogenic fault were determined to be 92? 58?and 62? respectively. That is, this earthquake was a mainly thrust faulting with the strike of near west-east and the dipping direction to south. The seismic moment released was 3.51020 Nm, accordingly, the moment magnitude MW was calculated to be 7.6. And then, 27 P-STFs, 22 S-STFs and the averaged STFs of them were determined respectively using the technique of spectra division in frequency domain and the synthetic seismogram as Greens functions. The analysis of the STFs suggested that the earthquake was a continuous event with the duration time of 19 s, starting rapidly and ending slowly. Finally, the temporal-spatial distribution of the slip on the fault plane was imaged from the obtained P-STFs and S-STFs using an time domain inversion technique. The maximum slip amplitude on the fault plane was about 7 m. The maximum stress drop was 30 MPa, and the average one over the whole rupture area was 7 MPa. The rupture area was about 85 km long in the strike direction and about 60 km wide in the down-dip direction, which, equally, was 51 km deep in the depth direction. The rupture propagated 50 km eastwards and 35 km westwards. The main portion of the rupture area, which has the slip amplitude greater than 0.5 m, was of the shape of an ellipse, its major axis oriented in the slip direction of the fault, which indicated that the rupture propagation direction was in accordance with the fault slip direction. This phenomenon is popular for strike-slip faulting, but rather rare for thrust faulting. The eastern portion of the rupture area above the initiation point was larger than the western portion below the initiation point, which was indicative of the asymmetrical rupture. In other words, the rupturing was kind of unilateral from west to east and from down to up. From the snapshots of the slip-rate variation with time and space, the slip rate reached the largest at the 4th second, that was 0.2 m/s, and the rupture in this period occurred only around the initiation point. At the 6th second, the rupture around the initiation point nearly stopped, and started moving outwards. The velocity of the westward rupture was smaller than that of the eastward rupture. Such rupture behavior like a circle mostly stopped near the 15th second. After the 16th second, only some patches of rupture distributed in the outer region. From the snapshots of the slip variation with time and space, the rupture started at the initiation point and propagated outwards. The main rupture on the area with the slip amplitude greater than 5 m extended unilaterally from west to east and from down to up between the 6th and the 10th seconds, and the western segment extended a bit westwards and downwards between the 11th and the 13th seconds. The whole process lasted about 19 s. The rupture velocity over the whole rupture process was estimated to be 3.3 km/s.
文摘On July 20, 1995, an earthquake of ML=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 WL=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 ML=2.0 earthquake as an empirical Green's function, a regularization method was applied toretrieve the far-field source-time fonction (STF) of the main shock. Considering the records of HDSN are the typeof 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 fieldstations, P phases in vertical direction which mostly are about 0.5 s in length were used. The STFs obtained fromeach seismic phases are in good agreement, showing that the Mt=4.1 eedhquake consisted of two events. STFsfrom each station demonstrate an obvious 'seismic Doppler effect'. Assuming the nodal plane striking 37° anddipping 40°, determined by using P wave first motion data and aftershock distribution, is the fault plane, through atrial and error method, the following results were drawn：Both of the events lasted about 0. 1s, the ruptUre length ofthe 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.0km/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 thefirst 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 ML=4.1 earthquake were also reached： the scalar earthquake moment is 3.3×1013 N·m, stress drop 4.6 MPa, rupture radius 0. 16 km.
基金This work is supported by the National Key R&D Program of China(No.2018YFC0603500)Programme on Global Change and Air-Sea Interaction(GASI-GEOGE-02)+2 种基金NSFC(Nos.41474050,41425012,41874062 and 41922025)111 project(No.BP0719022)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)CUG170602(D.W.).
文摘The 10 January 2018 Mw7.5 Swan island,Honduras earthquake occurred on the Swan island fault,which is a transform plate boundary between the North American and Caribbean plates.Here we back-project the rupture process of the earthquake using dense seismic stations in Alaska,and find that the earthquake ruptured at least three faults(three stages)for a duration of~40 s.The rupture speed for the longest fault(stage 3)is as fast as 5 km/s,which is much faster than the local shear wave velocity of~4 km/s.Supershear rupture was incidentally observed on long and straight strike-slip faults.This study shows a supershear rupture that occured on a strike-slip fault with moderate length,implying that supershear rupture might commonly occur on large strike-slip earthquakes.The common occurrence of supershear rupture on strike-slip earthquakes will challenge present understanding of crack physics,as well as strong ground motion evaluation in earthquake engineering.
