Relationships between moment magnitude and fault parameters:theoretical and semi-empirical relationships

Fault parameters are important in earthquake hazard analysis.In this paper,theoretical relationships between moment magnitude and fault parameters including subsurface rupture length,downdip rupture width,rupture area,and average slip over the fault surface are deduced based on seismological theory.These theoretical relationships are further simplified by applying similarity conditions and an unique form is established.Then,combining the simplified theoretical relationships between moment magnitude and fault parameters with seismic source data selected in this study,a practical semi-empirical relationship is established.The seismic source data selected is also to used to derive empirical relationships between moment magnitude and fault parameters by the ordinary least square regression method.Comparisons between semi-empirical relationships and empirical relationships show that the former depict distribution trends of data better than the latter.It is also observed that downdip rupture widths of strike slip faults are saturated when moment magnitude is more than 7.0,but downdip rupture widths of dip slip faults are not saturated in the molnent magnitude rangcs of this study.

RESEARCH ON THE DEFORMATIONS IN QINGHAI_TIBET PLATEAU AND ITS MARGINS BY INVERTING SEISMIC MOMENT TENSORS AND GPS VELOCITIES

We have determined approximate average rates of deformation in the Qinghai_Tibet plateau and its margins from the GPS data for last 10 years and the moment tensors from earthquakes between 1900 and 1999.We also determined the strain rate (seismic strain rate) associated with the seismic deformation using 254 M w ≥5.0 earthquakes,and estimated the shortening and extension rates for every block in the area as well.We also estimated the strain rate (geodetic strain rate)by 80 GPS sites’ velocity vectors and analyzed characteristic of kinematics by two kinds of strain rates and discussed earthquake potential in the area.As a result,the deformation rates from seismic moment tensors and from GPS velocities are basically agreed with each other.It is feasible to analyze seismic risk by comparing geodetic strain rate with seismic strain rate based on the opinion that strain energy will be released through earthquake.It is concluded that there is no strong earthquake potential (>M7) in the Qinghai_Tibet plateau and its margins,but there is earthquake potential (>M5) in middle Tibet in a few years.

Scaling relations of moment magnitude, local magnitude,and duration magnitude for earthquakes originated in northeast India

In this study, we aim to improve the scaling between the moment magnitude （Mw）, local magnitude （ME）, and the duration magnitude （MD） for 162 earthquakes in Shillong-Mikir plateau and its adjoining region of northeast India by extending the Mw estimates to lower magnitude earthquakes using spectral analysis of P-waves from vertical component seismograms. The Mw-ML and Mw-MD relationships are determined by linear regression analysis. It is found that, Mw values can be considered consistent with ME and My, within 0.1 and 0.2 magnitude units respectively, in 90 % of the cases. The scaling relationships investigated comply well with similar relationships in other regions in the world and in other seismogenic areas in the northeast India region.

Study on the Moment Magnitude of Small and Moderate Earthquakes Located in the Inner Mongolia Region

Based on digital seismic waveform data from Inner Mongolia Digital Seismic Network,the source spectrum parameters of 182 small and moderate earthquakes from January,2009to September,2016 are derived,and the seismic momentmoment magnitude MW of the earthquakes are calculated.Theand the relationship between stress drop and magnitude are obtained using the linear regression method.It is clear that incorporating the moment magnitude into the seismic quick report catalog and the official earthquake catalog can enrich earthquake observation report content,thus providing better service for earthquake emergency and earthquake scientific research.

Control Parameters of Magnitude—Seismic Moment Correlation for the Crustal Earthquakes

In connection with conversion from energy class KR (KR = log10E R, where ER — seismic energy, J) to the universal magnitude estimation of the Tien Shan crustal earthquakes the development of the self-coordinated correlation of the magnitudes (mb , ML, Ms ) and KR with the seismic moment M0 as the base scale became necessary. To this purpose, the first attempt to develop functional correlations in the magnitude—seismic moment system subject to the previous studies has been done. It is assumed that in the expression M (mb , ML , Ms) = Ki + zi log10M0 , the coefficients ki? and zi? are controlled by the parameters of ratio?(where;f0 —corner frequency, Brune, 1970, 1971;M0, N×m). According to the new theoretical predictions common functional correlation of the advanced magnitudes Mm (mbm = mb , MLm = ML , MSm = MS ) from log10M0,? log10t0? and the elastic properties (Ci) can be presented as , where , and , for the averaged elastic properties of the Earth’s crust for thembmthe coefficients Ci= –11.30 and di = 1.0, for MLm: Ci = –14.12, di = 7/6;for MSm : Ci = –16.95 and di = 4/3. For theTien Shan earthquakes (1960-2012 years) it was obtained that , and on the basis of the above expressions we received that MSm = 1.59mbm – 3.06. According to the instrumental data the correlation Ms = 1.57mb – 3.05 was determined. Some other examples of comparison of the calculated and observed magnitude - seismic moment ratios for earthquakes of California, the Kuril Islands, Japan, Sumatra and South America are presented.

A Catalogue of Earthquakes With Moment Magnitude of the Center-Southern Asia Test Area of the Global Seismic Hazard Assessment Program

We compile a regional catalogue of earthquakes with moment magnitude of the Center-Southern Asia test area(20°-35°N,85°-105°E)for the Global Seismic Hazard Assessment Program(GSHAP).There are significant inhomogeneous and uncompleted data,and n uniformity of earthquake magnitudes in this test area because this region is situated on the boundaries of many countries,such as China,India,Nepal,Vietnam,etc.We establish a relationship between Gutenberg surface-wave magnitude and IASPEI surface-wave magnitude,which can be used for conversion of different magnitude scales into moment magnitude for this catalogue.A catalogue of events with Mw≥6.0 of this test area is given at the end of this paper.

Moment magnitudes of two large Turkish earthquakes on February 6,2023 from long-period coda

Two large earthquakes(an earthquake doublet)occurred in south-central Turkey on February 6,2023,causing massive damages and casualties.The magnitudes and the relative sizes of the two mainshocks are essential information for scientific research and public awareness.There are obvious discrepancies among the results that have been reported so far,which may be revised and updated later.Here we applied a novel and reliable long-period coda moment magnitude method to the two large earthquakes.The moment magnitudes(with one standard error)are 7.95±0.013 and 7.86±0.012,respectively,which are larger than all the previous reports.The first mainshock,which matches the largest recorded earthquakes in the Turkish history,is slightly larger than the second one by 0.11±0.035 in magnitude or by 0.04 to 0.18 at 95%confidence level.

Seismic safety assessment with non-Gaussian random processes for train-bridge coupled systems

Extensive high-speed railway(HSR)network resembled the intricate vascular system of the human body,crisscrossing mainlands.Seismic events,known for their unpredictability,pose a significant threat to both trains and bridges,given the HSR’s extended operational duration.Therefore,ensuring the running safety of train-bridge coupled(TBC)system,primarily composed of simply supported beam bridges,is paramount.Traditional methods like the Monte Carlo method fall short in analyzing this intricate system efficiently.Instead,efficient algorithm like the new point estimate method combined with moment expansion approximation(NPEM-MEA)is applied to study random responses of numerical simulation TBC systems.Validation of the NPEM-MEA’s feasibility is conducted using the Monte Carlo method.Comparative analysis confirms the accuracy and efficiency of the method,with a recommended truncation order of four to six for the NPEM-MEA.Additionally,the influences of seismic magnitude and epicentral distance are discussed based on the random dynamic responses in the TBC system.This methodology not only facilitates seismic safety assessments for TBC systems but also contributes to standard-setting for these systems under earthquake conditions.

Towards fast focal mechanism inversion of shallow crustal earthquakes in the Chinese mainland

We have developed an automatic regional focal mechanism inversion system based on the Earthquake Rapid Report(ERR) system and the real-time three-component seismic waveform stream of 1 000 broadband seismic stations provided by the China Earthquake Networks Center(CENC). The system can rapidly provide a double couple solution and centroid depth within 5–15 min after receiving earthquake information from the ERR system.The data processing is triggered by earthquake information obtained from the ERR system. The system is capable of determining the focal mechanism of all shallow-depth earthquakes in the Chinese mainland with a magnitude of 5.5–6.5. It utilizes waveform data recorded by seismic stations located within 500 km from the epicenter,enabling the reporting of a focal mechanism solution within 5–15 min of an earthquake occurrence. Additionally,the system can assign a corresponding grade(A B C) to the focal mechanism solution. We processed a total of 301earthquakes that occurred from 2021 to June 2022, and after the quality control, 166 of them were selected.These selected solutions were manually checked, and 160 of them were compiled in a focal mechanism catalog.This catalog can be conveniently downloaded online via the Internet. The automatic focal mechanism solution of earthquakes in eastern China exhibits a good agreement with that provided by the Global Centroid Moment Tensor(GCMT), when available. The average Kagan angle between this catalog and GCMT is 22°, and the average difference in MWis 0.17. Furthermore, compared with GCMT, the minimum magnitude of our catalog has been reduced from approximately 5.0 to 4.0. The correlation between the centroid depth and crustal thickness in the Chinese mainland confirms the distribution of the centroid depth.

GPU-acceleration 3D rotated-staggered-grid solutions to microseismic anisotropic wave equation with moment tensor implementation

To improve the accuracy of microseismic inversion,seismic anisotropy and moment tensor source should be carefully considered in the forward modelling stage.In this study,3D microseismic anisotropy wave forward modelling with a moment tensor source was proposed.The modelling was carried out based on a rotated-staggered-grid(RSG)scheme.In contrast to staggered-grids,the RSG scheme defines the velocity components and densities at the same grid,as do the stress components and elastic parameters.Therefore,the elastic moduli do not need to be interpolated.In addition,the detailed formulation and implementation of moment-tensor source loaded on the RSG was presented by equating the source to the stress increments.Meanwhile,the RSG-based 3D wave equation forward modelling was performed in parallel using compute unified device architecture(CUDA)programming on a graphics processing unit(GPU)to improve its efficiency.Numerical simulations including homogeneous and anisotropic models were carried out using the method proposed in this paper,and compared with other methods to prove the reliability of this method.Furthermore,the high efficiency of the proposed approach was evaluated.The results show that the computational efficiency of proposed method can be improved by about two orders of magnitude compared with traditional central processing unit(CPU)computing methods.It could not only help the analysis of microseismic full wavefield records,but also provide support for passive source inversion,including location and focal mechanism inversion,and velocities inversion.

Moment tensor and stress inversion solutions of acoustic emissions during compression and tensile fracturing in crystalline rocks

We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in the specimens introduce a complex stress field,resulting in a spatial and temporal variation of focal mechanisms.Specifically,we consider two experimental setups:(1)where the rock is loaded in compression to generate primarily shear-type fractures and(2)where the material is loaded in indirect tension to generate predominantly tensile-type fractures.In each test,we first decompose AE moment tensors into double-couple(DC)and non-DC terms and then derive unambiguous normal and slip vectors using k-means clustering and an unstructured damped stress inversion algorithm.We explore temporal and spatial distributions of DC and non-DC events at different loading levels.The majority of the DC and the tensile non-DC events cluster around the pre-cut flaws,where macro-cracks later develop.Results of stress inversion are verified against the stress field from finite element(FE)modeling.A good agreement is found between the experimentally derived and numerically simulated stress orientations.To the best of the authors’knowledge,this work presents the first case where stress inversion methodologies are validated by numerical simulations at laboratory scale and under highly heterogeneous stress distributions.

Effectiveness of two conventional methods for seismic retrofit of steel and RC moment resisting frames based on damage control criteria

This study investigates the efficiency of two types of rehabilitation methods based on economic justification that can lead to logical decision making between the retrofitting schemes. Among various rehabilitation methods, concentric chevron bracing（CCB） and cylindrical friction damper（CFD） were selected. The performance assessment procedure of the frames is divided into two distinct phases. First, the limit state probabilities of the structures before and after rehabilitation are investigated. In the second phase, the seismic risk of structures in terms of life safety and financial losses（decision variables） using the recently published FEMA P58 methodology is evaluated. The results show that the proposed retrofitting methods improve the serviceability and life safety performance levels of steel and RC structures at different rates when subjected to earthquake loads. Moreover, these procedures reveal that financial losses are greatly decreased, and were more tangible by the application of CFD rather than using CCB. Although using both retrofitting methods reduced damage state probabilities, incorporation of a site-specific seismic hazard curve to evaluate mean annual occurrence frequency at the collapse prevention limit state caused unexpected results to be obtained. Contrary to CFD, the collapse probability of the structures retrofitted with CCB increased when compared with the primary structures.

Moment Magnitude and Its Calculation

In this paper,we give a brief introduction to the proposal and development history of the earthquake magnitude concept. Moment magnitude MWis the best physical quantity for measuring earthquakes. Compared with other magnitude scales used traditionally,moment magnitude is not saturated for all earthquakes,regardless of big and small earthquakes,deep and shallow earthquakes,far field and near field seismic data,geodetic and geological data,moment magnitude can be measured,and can be connected with wellknown magnitude scales such as surface wave magnitude MS. Moment magnitude is a uniform magnitude scale,which is suitable for statistics with wide magnitude range.Moment magnitude is the preferred magnitude selected by the International Seismological community,and it is preferred by the departments responsible for publishing seismic information to the public. Moment magnitude is a uniform magnitude scale,which is suitable for statistics with wide magnitude range. Moment magnitude is a preferred magnitude for international seismology,it is preferred by the agency responsible for providing information about earthquakes to the public. We provide all formulas used in the calculation of moment magnitude,and the calculation steps in detail. We also analyzed some problems and rules to solve these problems by using different formulas and numerical value calculation steps.

Evaluation of the seismic moment of the April 20,2013 Lushan earthquake

The April 20, 2013 Lushan earthquake which occurred in Sichuan, China had only moderate thrust. However, the computed seismic moments （M0） for the Lushan earthquake calculated by several institutions differ significantly from 0.4 × 1019 to 1.69 ×019 Nm, up to four times difference. We evaluate ten computed Mos by using normal mode observations from superconducting gravimeters in China's Mainland. We compute synthetic normal modes on the basis of moment tensor solutions and fit them to the observed normal modes. Comparison of our results indicates that Mo is the main cause for some large differences between observations and synthetics. We sug- gest that a moment magnitude of Mw6.6, corresponding to a Mo of 0.97-1.08 × 1019 Nm, characterizes the size and strength of the seismic source of the Lushan earthquake.

Estimating seismic moments and Lg Q using Lg spectra

A multi-event and multi-station inverse method is presented in the paper to simultaneously estimate the seismic moments (M0) and source corner frequencies (fc) of several Jiashi (Xinjiang, China) earthquakes, as well as the apparent Lg Q models for the paths from Jiashi to eight seismic stations (WMQ, AAK, TLG, MAKZ, KUR, VOS, ZRN and CHK) in Central Asia. The resultant seismic moments correlate well with the M0 values obtained by Harvard University using the centroid moment tensor (CMT) inversion and the surface-wave magnitudes as well. After the correction by a typical value of average radiation coefficient for regional SV waves, the M0 values from Lg spectral inversion are still close to the corresponding values obtained from CMT inversion. The obtained ap- parent Q0Lg values (Lg Q at 1 Hz) are consistent with the tectonic features of corresponding propagation paths. The Q0Lg values are 351±87, 349±86 and 300±27 for the paths from Jiashi to AAK, TLG and MAKZ, respectively. They are smaller than Q0Lg values for the paths to KUR, VOS, ZRN and CHK, which are 553±72, 569±58, 550±57 and 603±65, respectively. These results agree with the condition that the paths to AAK, TLG and MAKZ mainly propagate through the mountainous Tianshan area where relatively strong seismic activities and large variations of topography are exhibited, while the paths to KUR, VOS, ZRN and CHK mainly propagate through the stable area of Kazak platform. The Q0Lg value for the path to WMQ is 462±56. This is also in agreement with the condition that the path to WMQ is basically along the border area between Tianshan Mountain and Tarim Basin, and along this path the variations of topography and crustal thickness are moderate in comparison with that along the path to MAKZ.

Seismic moment tensor representations and radiation patterns in unbounded media with elipsoidal cavities driven by low frequency pressure

Seismic waves generated in elastic media by ellipsoidal cavities driven by low frequency pressure can be solved by Eshelby′s method. In this paper we prove that this method, originally presented for elastostatic ellipsoidal inclusion problems, can also be applied to dynamic problems on the assumption that wavelengths are much longer than the dimensions of cavities. Comparison between the approximate solution for spherical cavities and its relevant precise solution shows that this method may be used for radii of cavities smaller than 100 m and for frequencies below 10 Hz. By this approach we show that, for low frequency seismic waves in the far field an explosive load on an ellipsoidal cavity can be equivalent to three dipoles oriented along the principal axes of that cavity. Seismic radiation patterns are given analytically and results thus obtained show that an explosion detonated in an ellipsoidal cavity may radiate significant S waves.

On the Consistency of Large Earthquake Moment and Strain Rate Inferred from GPS Data in North China

The new GPS data can map crustal strain rates over large areas with a useful degree of precision. Stable strain measurement results open the door for improved estimates of earthquake occurrence. The Kostrov’s formula (1974) translates the smoothed strain rates in North China into geodetic moment rates. In North China, the ratio of seismic moment released to moment accumulated from GPS measurement is 60.6% in NS direction, 68.9% in EW direction, and 104.1% in NE shear direction. The near unit ratio points to the reliability of GPS measurements there. The combination of historical seismicity and GPS measurement offers a powerful attack on earthquake hazard.

Using Long-Period Body Wave to Inverse Moment Tensors of the M_S6.8 Jiashi,Xinjiang Earthquake in 2003 and the Moderate and Small Earthquakes before and after It

An M_S6.8 strong earthquake took place in Jiashi,Xinjiang on February 24 of 2003.The digital wave form data recorded in Kashi and Wushi stations are selected to inverse the moment tensor solutions for the strong earthquake and the moderate and small earthquakes before and after it(108 earthquakes in 2001～2004).67 focal mechanism solutions have been calculated,and the results agree with those from Harvard University and USGS.The analysis reveals that before the strong earthquake,the moderate and small earthquake distribution was dispersed,and after the event the distribution was mainly concentrated around the strong earthquake.Before the strong earthquake,the seismic faults of the mid and small events had the character of strike-slip and normal faulting,and after the event,they exhibit strike-slip and thrust faulting.The region is dominated by near-NS horizontal compression from the southern block after the strong earthquake.

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

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

Applications of seismic moment tensor inversion in fast response to earthquakes

Using the technique of seismic moment tensor inversion, the source mechanisms of 10 earthquakes with Ms5.2that occurred in China from November 1996 to January 1998 were determined rapidly. The determined resultswere sent as 'Bulletins of Source Mechanism Parameters of Earthquakes' to the Seismic Regime Guards' Office,China Seismological Bureau, and the relevant provincial seismological bureaus. These bulletins have played rolein the fast response to large earthquakes.