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

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

The October 11, 1999 and November 08, 2006 Beni Suef Earthquakes

Two felt moderate-sized earthquakes with local magnitudes of 4.9 on October 11, 1999 and 4.3 on November 08, 2006 occurred southeast of Beni Suef and Cairo cities. Being well recorded by the digital Egyptian National Seismic Network (ENSN) and some regional broadband stations, they provided us with a unique opportunity to study the tectonic process and present-day stress field acting on the northern part of the Eastern Desert of Egypt. In this study, we analyze the main shocks of these earthquakes and present 15 well recorded aftershocks (0.9 ≤ ML ≤ 3.3) which have small errors on both horizontal and vertical axes. The relocation analysis using the double difference algorithm clearly reveals a NW trending fault for the 1999 earthquake. The spatial distribution of its aftershocks indicates a propagation of rupture from the SW towards the NW along a fault length ~5 km dipping nearly ~40°SW. We also determined the focal mechanisms of the two main shocks by two methods (polarities and amplitudes ratios of P, S_V and S_H and regional waveform inversion). Our results indicate a normal faulting mechanism with a slight shear component for the first event, while pure normal faulting for the second one. The spatial distribution of the 1999 aftershocks sequence along with the retrieved focal mechanism confirmed the NW plane as the true fault plane. While for the 2006 event, the few aftershocks do not reveal any fault geometry; its focal mechanism indicated a pure normal fault nearly trending WNW-ESE that corresponds more likely to the extension of the 1999 earthquake fault. The seismicity distribution between the two earthquake sequences reveals a noticeable gap that may be a site of a future event. The NNE-SSW extensional stress indicated by the mechanisms of these events is in agreement with the regional stress field and the rifting of the northern Red Sea in its northern branches (Gulf of Suez and Gulf of Aqaba). The source parameters (seismic moment, moment magnitude, fault radius, stress drop and displacement across the fault) were also estimated and compared based on both the regional waveform inversion and the displacement spectra and interpreted in the context of the tectonic setting. The obtained results imply a reactivation of the pre-exiting NW-SE faults as a result of extensional deformation from the northern Red Sea-Gulf of Suez rifts.     

Present-day seismicity,stress field and crustal deformation of Egypt

In this study we investigate present-day seismicity and crustal deformation of Egypt based on a comprehensive earthquake catalog from 1900 to 2004 by focal mechanism stress inversion and by recent GPS observations. Spatial distribution of earthquake epicenters indicates that Egypt has been suffered from both interplate and intraplate earthquakes. Most earthquake activity (more than 70%) has been concentrated in northern Egypt along the geologically documented borders of Sinai subplate (northern Red Sea and its two branches Suez rift and Aqaba–Dead Sea transform). The majority of inland earthquake focal mechanisms in Egypt are normal with strike-slip component or strike-slip faulting events. Only a small minority, namely four events, exhibits reverse faulting. The inversion method of Gephart and Forsyth (1984) was applied to calculate the orientation of the principle stress axes and the shape of the stress tensor. The best fitting tensor in Egypt shows homogeneity stress field. The tension stress regime is dominant in northern Egypt. The stress directions are well resolved by the 95% confidence limits, the relative stress magnitude has a value of about 0.3. However, along southern Egypt the strike-slip regime is dominant. The shape factor (R-value) is 0.5, which means that the deviatoric components of σ₁ and σ₃ are of the same magnitude, but of opposite signs. The average horizontal velocity of GPS stations in Egypt is 5.15± 1.1 mm/year in mostly NNW direction. The results of deformation analysis indicate that the northern Egypt is deformed more than the southern part. Only the Egyptian-Mediterranean coastal–Nile Delta zone dominates as a compression deformation area. However, an extensional deformation has been observed throughout the rest of country. This means that the relative motion of African plate with respect to both Eurasian and Arabian has highly controlled the deformation processes in Egypt.     

Coupling Between Seismic Activity and Hydrogeochemistry at the Shillong Plateau, Northeastern India

Transient hydrogeochemical anomalies were detected in a granite-hosted aquifer, which is located at a depth of 110 m, north of the Shillong Plateau, Assam, India, where groundwater chemistry is mainly buffered by feldspar alteration to kaolinite. Their onsets preceded moderate earthquakes on December 9, 2004 (M_W = 5.3) and February 15, 2005 (M_W = 5.0), respectively, 206 and 213 km from the aquifer. The ratios [Na+K]/Si, Na/K and [Na+K]/Ca, conductivity, alkalinity and chloride concentration began increasing 3–5 weeks before the M_W = 5.3 earthquake. By comparison with field, experimental and theoretical studies, we interpret a transient switchover between source aquifers, which induced an influx of groundwater from a second aquifer, where groundwater chemistry was dominantly buffered by the alteration of feldspar to smectite. This could have occurred in response to fracturing of a hydrological barrier. The ratio Ba/Sr began decreasing 3–6 days before the M_W = 5.0 earthquake. We interpret a transient switchover to anorthite dissolution caused by exposure of fresh plagioclase to groundwater interaction. This could have been induced by microfracturing, locally within the main aquifer. By comparison with experimental studies of feldspar dissolution, we interpret that hydrogeochemical recovery was facilitated by groundwater interaction and clay mineralization, which could have been coupled with fracture sealing. The coincidence in timing of these two hydrogeochemical events with the only two M_W ≥ 5 earthquakes in the study area argues in favor of cause-and-effect seismic-hydrogeochemical coupling. However, reasons for ambiguity include the lack of similar hydrogeochemical anomalies coupled with smaller seismic events near the monitoring station, the >200 km length scale of inferred seismic-hydrogeochemical coupling, and the potential for far-field effects related to the Great Sumatra–Andaman Islands Earthquake of December 26, 2004.     

Relocation of the MS5.7 Earthquake Sequence in Songyuan, Jilin Province, 2018 and the Analysis of Its Seismogenic Structure

The 2018,Songyuan,Jilin M_S5. 7 earthquake occurred at the intersection of the FuyuZhaodong fault and the Second Songhua River fault. The moment magnitude of this earthquake is M_W5. 3,the centroid depth by the waveform fitting is 12 km,and it is a strike-slip type event. In this paper,with the seismic phase data provided by the China Earthquake Network, the double-difference location method is used to relocate the earthquake sequence,finally the relocation results of 60 earthquakes are obtained. The results show that the aftershock zone is about 4. 3km long and 3. 1km wide,which is distributed in the NE direction. The depth distribution of the seismic sequence is 9km-10 km. 1-2 days after the main shock,the aftershocks were scattered throughout the aftershock zone,and the largest aftershock occurred in the northeastern part of the aftershock zone. After 3-8 days,the aftershocks mainly occurred in the southwestern part of the aftershock zone. The profile distribution of the earthquake sequence shows that the fault plane dips to the southeast with the dip angle of about 75°. Combined with the regional tectonic setting,focal mechanism solution and intensity distribution,we conclude that the concealed fault of the Fuyu-Zhaodong fault is the seismogenic fault of the Songyuan M_S5. 7 earthquake. This paper also relocates the earthquake sequence of the previous magnitude 5. 0 earthquake in 2017. Combined with the results of the focal mechanism solution,we believe that the two earthquakes have the same seismogenic structure,and the earthquake sequence generally develops to the southwest. The historical seismic activity since 2009 shows that after the magnitude 5. 0 earthquake in 2017,the frequency and intensity of earthquakes in the earthquake zone are obviously enhanced,and attention should be paid to the development of seismic activity in the southwest direction of the earthquake zone.     

Relative seismic moment tensor determination for Vrancea intermediate depth earthquakes

The method of relative seismic moment tensor determination proposed byStrelitz (1980) is extended a) from an interactive time domain analysis to an automated frequency domain procedure, and b) from an analysis of subevents of complex deep-focus earthquakes to the study of individual source mechanism of small events recorded at few stations.The method was applied to the recovery of seismic moment tensor components of 95 intermediate depth earthquakes withM _L=2.6–4.9 from the Vrancea region, Romania. The main feature of the obtained fault plane solutions is the horizontality ofP axes and the nonhorizontal orienaation ofT axes (inverse faulting). Those events with high fracture energy per unit area of the fault can be grouped unambiguously into three depth intervals: 102–106 km, 124–135 km and 141–152 km. Moreover, their fault plane solutions are similar to ones of all strong and most moderate events from this region and the last two damaging earthquakes (November 10, 1940 withM _W=7.8 and March 4 1977 withM _W=7.5) occurred within the third and first depth interval, respectively. This suggests a possible correlation at these depths between fresh fracture of rocks and the occurrence of strong earthquakes.     

Rapid report of June 1, 2022 MW 5.9 Lushan earthquake,China with geodetic and teleseismic data

Timely response to earthquake characterization can facilitate earthquake emergency rescue and further scientific investigations. On June 1, 2022, M_W 5.9 earthquake occurred in the southern area of the Longmenshan fault zone. This event also happened at the south end of the Dayi seismic gap and is the largest earthquake that has occurred in this seismic gap since the 1970 M 6.2 event. The slip-distribution model constrained by the seismic waveforms suggests a thrust-dominated faulting mechanism. The main slip occurs at a depth of ～14 ?km, and the cumulative energy is released in the first 6 ?s. The variations of Coulomb stress caused by the mainshock show a positive change in the southwest area of the Dayi seismic gap, indicating possible activation of future earthquakes. In addition, we emphasize the importance of rapid estimation of deformation for near-field hazard delineation, especially when interferometric radar fails to image coseismic deformation in a high relief terrain.     

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

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

Occurrence probability of moderate to large earthquakes in Italy based on new geophysical methods

We develop new approaches to calculating 30-year probabilities for occurrence of moderate-to-large earthquakes in Italy. Geodetic techniques and finite-element modelling, aimed to reproduce a large amount of neotectonic data using thin-shell finite element, are used to separately calculate the expected seismicity rates inside seismogenic areas (polygons containing mapped faults and/or suspected or modelled faults). Thirty-year earthquake probabilities obtained from the two approaches show similarities in most of Italy: the largest probabilities are found in the southern Apennines, where they reach values between 10% and 20% for earthquakes of M _W ≥ 6.0, and lower than 10% for events with an M _W ≥ 6.5.     

Analysis on the master event method and precise location of 1997 Jiashi strong earthquake swarm in western China

IntroductionEarthquakelocationisoneoftheoldestinverseproblemsinseismology.Preciseearthquakelocationisthebasisoftheseismicresearch.Theprecisionoftheearthquakelocationdependsonmanyfactors,suchasthequalityofthestationnetwork,theprecisionoftheseismicwavearrivaltimemeasured,andtheknowledgeaboutthecrustalvelocitystructureandsoon.AccordingtotheresearchmadebyZhao(1983),thehorizontallocationerroroftheeanhquakelocatedbythetraditionalabsolutemethodinthecentralareaofNorthChinaisabouttZkin,intheedgeregio…     

Multifractal Analysis of Earthquakes in the Southeastern Iran-Bam Region

Earthquakes in Iran and neighbouring regions are closely connected to their position within the geologically active Alpine-Himalayan belt. Modern tectonic activity is forced by the convergent movements between two plates: The Arabian plate, including Saudi Arabia, the Persian Gulf and the Zagros Ranges of Iran, and the Eurasian plate. The intensive seismic activity in this region is recorded with shallow focal depth and magnitude rising as high as M_w = 7.8. The study region can be attributed to a highly complex geodynamic process and therefore is well suited for multifractal seismicity analysis. Multifractal analysis of earthquakes (m_b ≥ 3) occurring during 1973 – 2006 led to the detection of a clustering pattern in the narrow time span prior to all the large earthquakes: M_w = 7.8 on 16.9.1978; M_w = 6.8 on 26.12.2003; M_w = 7.7 on 10.5.97. Based on the spatio-temporal clustering pattern of events, the potential for future large events can be assessed. Spatio-temporal clustering of events apparently indicates a highly stressed region, an asperity or weak zone from which the rupture propagation eventually nucleates, causing large earthquakes. This clustering pattern analysis done on a well-constrained catalogue for most of the fault systems of known seismicity may eventually aid in the preparedness and earthquake disaster mitigation.     

Fracture characteristics of the 1997 Jiashi, Xinjiang, China, earthquake swarm inferred from source spectra

Broadband P and S waves source spectra of 12 MS5.0 earthquakes of the 1997 Jiashi, Xinjiang, China, earthquake swarm recorded at 13 GDSN stations have been analyzed. Rupture size and static stress drop of these earthquakes have been estimated through measuring the corner frequency of the source spectra. Direction of rupture propagation of the earthquake faulting has also been inferred from the azimuthal variation of the corner frequency. The main results are as follows: ①The rupture size of MS6.0 strong earthquakes is in the range of 10～20 km, while that of MS=5.0～5.5 earthquakes is 6～10 km.② The static stress drop of the swarm earthquakes is rather low, being of the order of 0.1 MPa. This implies that the deformation release rate in the source region may be low. ③ Stress drop of the earthquakes appears to be proportional to their seismic moment, and also to be dependent on their focal mechanism. The stress drop of normal faulting earthquakes is usually lower than that of strike-slip type earthquakes. ④ For each MS6.0 earthquake there exists an apparent azimuthal variation of the corner frequencies. Azimuthally variation pattern of corner frequencies of different earthquakes shows that the source rupture pattern of the Jiashi earthquake swarm is complex and no uniform rupture expanding direction exists.     

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

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

Possible strain events reflected in water levels in wells along San Jacinto fault zone,southern California

Water levels have been monitored in wells along the San Jacinto fault zone since 1977. The three largest earthquakes to occur within 30 km of the segment of the San Jacinto fault zone being monitored with continuous recorders showed magnitudesM of 4.5, 4.8, and 5.5. Two wells in Borrego Valley, 31 to 32 km southeast of theM=5.5 earthquake on 25 February 1980, showed anomalous spikes recording a probable strain event 88 hours before the earthquake. Two other wells 12 km northwest of the epicenter showed no water-level anomalies. No water-level anomalies preceded theM=4.8 earthquake near Anza on 15 June 1982. Anomalous water-level fluctuations occurred in a well near Ocotillo Wells, 13 km northeast of theM=4.5 earthquake on 22 March 1982, 19 to 23 days prior to the earthquake. Similar fluctuations in other wells have not been followed by sizable earthquakes. A simultaneous drop in water level occurred in four wells on 8 September 1982; this possible strain event was not associated with a sizable earthquake. The presumed strain events occur only in wells that show earth tides and may have been the result of creep on strands of the San Jacinto fault zone. Although water-level anomalies have occurred in only one or two wells prior to two out of three moderate (M=4.5–5.5) earthquakes, the simultaneous drop in water level on 8 September 1982 and the spikes in two wells before theM=5.5 earthquake on 25 February 1980 suggest that wells responsive to earth tides may detect strain events.     

Source characteristics of the deep Philippine earthquake cluster, 23 and 24 July 2010

We determine the rupture velocity, rupture area, stress drop and duration of four strong deep-focus earthquakes in the Philippines by back-projecting the teleseismic P waves. Four deep-focus earthquakes occurred in a totally consumed Molucca microplate; their focal depths were greater than 550 km and their moment magnitudes were between M _w 6.6 and M _w 7.6. By studying this deep-focus cluster, we are able to estimate the rupture velocity, rupture area and stress drop which would assist in constraining the physical mechanism for earthquakes deeper than 500 km. Since the Molucca microplate is totally consumed, little evidence is left on the surface for us to do research. This deep-focus cluster provides us the opportunity to reveal the properties of this totally consumed microplate by using seismic method for the first time. Four earthquakes in this deep-focus cluster all have multiple rupture subevents. The M _w 7.3 event ruptures in two subevents, the M _w 7.6 and M _w 7.4 events both have three subevents. The M _w 6.6 event has single peak on the amplitude as a function of time; however, its energy releases at two spatially separated areas. Our results show that this deep-focus cluster has a slow rupture velocity which is about 0.27 to 0.43 of the shear wave velocity, long-scaled duration, concentrated energy release area, and high stress drop. These source properties are similar to those of other deep earthquakes occurring in warm slabs and indicate that the totally consumed Molucca microplate possibly is a warm plate.     

Rift tectonics development in the light of geophysical data: Red sea region

Summary The paper discusses geophysical phenomena observed over tholeiitic dykes detected in the Eastern Desert of Egypt (Nubian Block). The dykes show specific physical properties and are either normally or reversely polarized along their whole length over 120 km. Dykes of similar properties may also be interpreted from airborne geophysical maps in Saudi Arabia (Arabian Block). Magnetic fields of the tholeiitic dykes were interpreted quantitatively: they fit the magnetic pattern over the Red Sea. The study emphasizes those geophysical phenomena which may serve as criteria for recognition of regions with rift tectonics development.     

Performance of the Italian strong motion network during the 2009, L��Aquila seismic sequence (central Italy)

On April 6, 2009, the town of L’Aquila in the Abruzzo region (central Italy) was struck by a seismic event at 01:32 (UTC), of magnitude M_W = 6.3. The mainshock was followed by a long period of intense seismic activity and within seven days after the mainshock there were seven events of magnitude M_W ≥ 5 that occurred from April 6 to April 13. This long seismic sequence was characterized by a complex rupture mechanism that involved two major normal faults of the central Apennines: the Paganica and the Gorzano faults. The strong-motions of the mainshock were recorded by 64 stations of the Italian Strong-motion Network (RAN) operated by the National Civil Protection Department (DPC). Six stations of a local strong-motion array were working in NW L’Aquila suburb area. One of them, located at about 6 km from the Paganica fault surface tip-line, set up in trigger mode, recorded continuously for more than 20 min the mainshock and the aftershocks. Besides the mainshock, the RAN stations recorded in total 78 foreshocks and aftershocks of M_L ≥ 3.5, during the period from January to December 2009. The corresponding waveforms provide the most extensive digital strong ground motion data set ever recorded in Italy. Moreover, the 48 three-component observations of events of magnitude M_W ≥ 5, recorded at a distance less than 15 km from each of the major involved faults, provide a significant increasing of near-field records available for the Italian territory. Six days after the mainshock, the strong-motion dataset, referred to preliminary locations of the events with M_L ≥ 4.0, was made available on the DPC web site () and at the same time it was delivered to the ITACA database (). This dataset has been used by many authors in scientific papers and by engineers, geophysicists and geologists for professional technical works. In this paper, the present-day available strong-motion signals from the L’Aquila sequence and the performance of the Italian strong-motion network in terms of the number and quality of recorded data, the geometry and data transmission system are described. In addition the role of the temporary network that represents an extension of the permanent Italian strong-motion network, supporting the emergency response by civil protection authorities and improving the network coverage has been evaluated.     

Seismicity in the block mountains between Halle and Leipzig,Central Germany: centroid moment tensors,ground motion simulation,and felt intensities of two <Emphasis Type="Italic">M</Emphasis> ≈ 3 earthquakes in 2015 and 2017

On April 29, 2017 at 0:56 UTC (2:56 local time), an M_W =?2.8 earthquake struck the metropolitan area between Leipzig and Halle, Germany, near the small town of Markranstädt. The earthquake was felt within 50 km from the epicenter and reached a local intensity of I₀ = IV. Already in 2015 and only 15 km northwest of the epicenter, a M_W =?3.2 earthquake struck the area with a similar large felt radius and I₀ = IV. More than 1.1 million people live in the region, and the unusual occurrence of the two earthquakes led to public attention, because the tectonic activity is unclear and induced earthquakes have occurred in neighboring regions. Historical earthquakes south of Leipzig had estimated magnitudes up to M_W ≈?5 and coincide with NW-SE striking crustal basement faults. We use different seismological methods to analyze the two recent earthquakes and discuss them in the context of the known tectonic structures and historical seismicity. Novel stochastic full waveform simulation and inversion approaches are adapted for the application to weak, local earthquakes, to analyze mechanisms and ground motions and their relation to observed intensities. We find NW-SE striking normal faulting mechanisms for both earthquakes and centroid depths of 26 and 29 km. The earthquakes are located where faults with large vertical offsets of several hundred meters and Hercynian strike have developed since the Mesozoic. We use a stochastic full waveform simulation to explain the local peak ground velocities and calibrate the method to simulate intensities. Since the area is densely populated and has sensitive infrastructure, we simulate scenarios assuming that a 12-km long fault segment between the two recent earthquakes is ruptured and study the impact of rupture parameters on ground motions and expected damage.     

Macroseismic evaluation of the impact of the Olyutorskii earthquake and its larger aftershocks

This paper presents the results from a macroseismic survey of the impact and consequences of the M _W = 7.6 April 20(21), 2006 Olyutorskii earthquake in the area of the Koryak Autonomous Okrug and the adjacent areas in Kamchatka and Magadan regions. The earthquake was felt over an area of about 400000 sq. km with intensities of II to IX–X on the MSK-64 scale. Information was gathered from 37 population centers situated in this area and was used to present a summary of felt effects, to construct an isoseismals map, and to determine the macroseismic magnitude.     

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

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