Understanding Micropolar Theory in the Earth Sciences II: The Seismic Moment Tensor

SUMMARY Seismic moment tensors are an important tool and input variable for many studies in the geosciences. The theory behind the determination of moment tensors is well established. They are routinely and (semi-) automatically calculated on a global scale. However, on regional and local scales, there are still several difficulties hampering the reliable retrieval of the full seismic moment tensor. In an earlier study, we showed that the waveform inversion for seismic moment tensors can benefit significantly when incorporating rotational ground motion in addition to the commonly used translational ground motion. In this study, we test, what is the best processing strategy with respect to the resolvability of the seismic moment tensor components: inverting three-component data with Green’s functions (GFs) based on a 3-D structural model, six-component data with GFs based on a 1-D model, or unleashing the full force of six-component data and GFs based on a 3-D model? As a reference case, we use the inversion based on three-component data and 1-D structure, which has been the most common practice in waveform inversion for moment tensors so far. Building on the same Bayesian approach as in our previous study, we invert synthetic waveforms for two test cases from the Korean Peninsula: one is the 2013 nuclear test of the Democratic People’s Republic of Korea and the other is an Mw 5.4 tectonic event of 2016 in the Republic of Korea using waveform data recorded on stations in Korea, China and Japan. For the Korean Peninsula, a very detailed 3-D velocity model is available. We show that for the tectonic event both, the 3-D structural model and the rotational ground motion, contribute strongly to the improved resolution of the seismic moment tensor. The higher the frequencies used for inversion, the higher is the influence of rotational ground motions. This is an important effect to consider when inverting waveforms from smaller magnitude events. The explosive source benefits more from the 3-D structural model than from the rotational ground motion. Nevertheless, the rotational ground motion can help to better constraint the isotropic part of the source in the higher frequency range.

Download Full-text

Seismic Moment Tensor Solutions from GeoNet Data to Provide a Moment Magnitude Scale for New Zealand

10.26686/wgtn.16947268.v1 ◽

2021 ◽

Author(s):

◽

Elizabeth de Joux Robertson

Keyword(s):

New Zealand ◽

Seismic Moment ◽

Moment Tensor ◽

Velocity Model ◽

Moment Tensor Inversion ◽

Moment Magnitude ◽

Seismic Moment Tensor ◽

Waveform Data ◽

Moment Tensors ◽

The Moment

The aim of this project is to enable accurate earthquake magnitudes (moment magnitude, MW) to be calculated routinely and in near real-time for New Zealand earthquakes. This would be done by inversion of waveform data to obtain seismic moment tensors. Seismic moment tensors also provide information on fault-type. I use a well-established seismic moment tensor inversion method, the Time-Domain [seismic] Moment Tensor Inversion algorithm (TDMT_INVC) and apply it to GeoNet broadband waveform data to generate moment tensor solutions for New Zealand earthquakes. Some modifications to this software were made. A velocity model can now be automatically used to calculate Green's functions without having a pseudolayer boundary at the source depth. Green's functions can be calculated for multiple depths in a single step, and data are detrended and a suitable data window is selected. The seismic moment tensor solution that has either the maximum variance reduction or the maximum double-couple component is automatically selected for each depth. Seismic moment tensors were calculated for 24 New Zealand earthquakes from 2000 to 2005. The Global CMT project has calculated CMT solutions for 22 of these, and the Global CMT project solutions are compared to the solutions obtained in this project to test the accuracy of the solutions obtained using the TDMT_INVC code. The moment magnitude values are close to the Global CMT values for all earthquakes. The focal mechanisms could only be determined for a few of the earthquakes studied. The value of the moment magnitude appears to be less sensitive to the velocity model and earthquake location (epicentre and depth) than the focal mechanism. Distinguishing legitimate seismic signal from background seismic noise is likely to be the biggest problem in routine inversions.

Download Full-text

Seismic moment tensor event screening

Geophysical Journal International ◽

10.1093/gji/ggz578 ◽

2020 ◽

Vol 221 (1) ◽

pp. 77-88

Author(s):

Sean R Ford ◽

Gordon D Kraft ◽

Gene A Ichinose

Keyword(s):

Seismic Moment ◽

Moment Tensor ◽

Test Site ◽

Seismic Moment Tensor ◽

Data Set ◽

New Approach ◽

Moment Tensors ◽

Explosion Monitoring ◽

Langevin Distribution ◽

Monitoring Practice

SUMMARY Event screening is an explosion monitoring practice that aims to identify an event as an explosion (‘screened in’) or not (‘screened out’). Confidence in event screening can be increased if multiple independent approaches are used. We describe a new approach to event screening using the seismic moment tensor and its representation on the hypersphere, specifically the 5-sphere of 6-degree unit vectors representing the normalized symmetric moment tensor. The sample of moment tensors from an explosion data set is unimodal on the 5-sphere and can be parametrized by the Langevin distribution, which is sometimes referred to as the Normal distribution on the hypersphere. Screening is then accomplished by finding the angle from the explosion population mean to any newly measured moment tensor and testing if that angle is in the tail of the Langevin distribution (conservatively quantified as greater than 99.9 per cent of the cumulative density). We apply the screen to a sample of earthquakes from the Western USA and the September 2017 explosion and subsequent collapse at the Pungyye-Ri Test Site in North Korea. All the earthquakes and the collapse screen out, but the explosion does not.

Download Full-text

Analysis of non-double-couple source mechanisms in an area of induced seismicity, West Texas (USA).

10.5194/egusphere-egu21-13438 ◽

2021 ◽

Author(s):

Savvaidis Alexandros ◽

Roselli Pamela

Keyword(s):

Stress Field ◽

Time Domain ◽

Seismic Moment ◽

Moment Tensor ◽

P Wave ◽

Seismic Moment Tensor ◽

Ground Displacement ◽

Moment Tensors ◽

West Texas ◽

Double Couple

In the scope to investigate the possible interactions between injected fluids, subsurface geology, stress field and triggering earthquakes, we investigate seismic source parameters related to the seismicity in West Texas (USA). The analysis of seismic moment tensor is an excellent tool to understand earthquake source process kinematics; moreover, changes in the fluid volume during faulting leads to existence of non-double-couple (NDC) components (Frohlich, 1994; Julian et al., 1998; Miller et al., 1998). The NDC percentage in the source constitutes the sum of absolute ISO and CLVD components so that %NDC= % ISO + %CLVD and %ISO+%CLVD+%DC=100%. It is currently known that the presence of NDC implies more complex sources (mixed shear-tensile earthquakes) correlated to fluid injections, geothermal systems and volcano-seismology where induced and triggered seismicity is observed.With this hypothesis, we analyze the micro-earthquakes (M <2 .7) recorded by the Texas Seismological Network (TexNet) and a temporary network constituted by 40 seismic stations (equipped by either broadband or 3 component geophones). Our study area is characterized by Northwest-Southeast faults that follow the local stress/field (SHmax) and the geological characteristic of the shallow basin structure of the study area. After a selection based on signal-to-noise ratio, we filter (1-50 Hz) the seismograms and estimate P-wave pulse polarities and the first P-wave ground displacement pulse in time domain. Then, we perform the full moment tensor analysis by using hybridMT technique (Andersen, 2001; Kwiatek et al., 2016) with a detailed 1D velocity model. The key parameter is the polarity/area of the first P-wave ground displacement pulse in time domain. Uncertainties of estimated moment tensors are expressed by normalized root-mean-square (RMS errors) between theoretical and estimated amplitudes (Vavricuk et al., 2014). We also evaluate the quality of the seismic moment tensors by bootstrap and resampling. In our preliminary results we obtain NDC percentage (in terms of %ISO and %CLVD components), Mw, seismic moment, P, T and B axes orientation for each source inverted.

Download Full-text

Understanding Micropolar Theory in the Earth Sciences I: The Eigenfrequency $$\omega _r$$

Pure and Applied Geophysics ◽

10.1007/s00024-021-02932-7 ◽

2021 ◽

Author(s):

Rafael Abreu ◽

Stephanie Durand

Keyword(s):

Seismic Data ◽

Original Model ◽

Micropolar Theory ◽

Elastic Parameters ◽

Material Parameters ◽

The Earth ◽

Seismic Experiment ◽

Seismological Data ◽

Micropolar Model ◽

General Dynamic

AbstractEven though micropolar theories are widely applied for engineering applications such as the design of metamaterials, applications in the study of the Earth’s interior still remain limited and in particular in seismology. This is due to the lack of understanding of the required elastic material parameters present in the theory as well as the eigenfrequency $$\omega _r$$ ω r which is not observed in seismic data. By showing that the general dynamic equations of the Timoshenko’s beam is a particular case of the micropolar theory we are able to connect micropolar elastic parameters to physically measurable quantities. We then present an alternative micropolar model that, based on the same physical basis as the original model, circumvents the problem of the original eigenfrequency $$\omega _r$$ ω r laking in seismological data. We finally validate our model with a seismic experiment and show it is relevant to explain observed seismic dispersion curves.

Download Full-text

Rotational Ground Motion Measurements for Regional Seismic Moment Tensors: a Review

10.31223/x5jg9m ◽

2021 ◽

Author(s):

Stefanie Donner

Keyword(s):

Ground Motion ◽

Seismic Moment ◽

Spatial Scales ◽

Moment Tensor ◽

Waveform Inversion ◽

Seismic Source ◽

Seismic Moment Tensor ◽

Moment Tensors ◽

Rotational Ground Motion ◽

Motion Measurements

Seismic moment tensors are an important tool in geosciences on all spatial scales and for a broad range of applications. The basic underlying theory is established since decades. However, various factors influence the reliability of the inversion result, several of them are mutually dependent. Hence, a reliable retrieval of seismic moment tensors is still hampered in many cases, especially at regional event-receiver distances.To sample the entire wavefield due to a seismic source we need six components: three translational and three rotational ones. Up to now, only translational ground motion recordings were used for moment tensor retrieval, missing out valuable information. Using rotational in addition to the classical translational ground motions during waveform inversion for moment tensors mainly adds information on the vertical displacement gradient to the inversion problem. Furthermore, having available six instead of only three components per receiver location provides additional constraints on the sampling of the radiation pattern. As a result, the moment tensor components are resolved with higher precision and accuracy, even when the number of recording receivers is considerably reduced. Especially, components with a dependence to depth as well as the centroid depth can benefit significantly from additional rotational ground motion. Up to the time of writing this review only a few studies are published on the topic. Here, I summarise their findings and provide an overview over the possible capabilities of including rotational ground motion measurements to waveform inversion for seismic moment tensor retrieval.

Download Full-text

The effects of earth model uncertainty on the inversion of seismic data for seismic source functions

Geophysical Journal International ◽

10.1093/gji/ggaa408 ◽

2020 ◽

Vol 224 (1) ◽

pp. 100-120

Author(s):

Christian Poppeliers ◽

Leiph Preston

Keyword(s):

Seismic Data ◽

Model Uncertainty ◽

Wave Speed ◽

Moment Tensor ◽

Synthetic Data ◽

Seismic Energy ◽

Seismic Source ◽

Earth Model ◽

Seismic Moment Tensor ◽

Independent Components

SUMMARY We use Monte Carlo simulations to explore the effects of earth model uncertainty on the estimation of the seismic source time functions that correspond to the six independent components of the point source seismic moment tensor. Specifically, we invert synthetic data using Green’s functions estimated from a suite of earth models that contain stochastic density and seismic wave-speed heterogeneities. We find that the primary effect of earth model uncertainty on the data is that the amplitude of the first-arriving seismic energy is reduced, and that this amplitude reduction is proportional to the magnitude of the stochastic heterogeneities. Also, we find that the amplitude of the estimated seismic source functions can be under- or overestimated, depending on the stochastic earth model used to create the data. This effect is totally unpredictable, meaning that uncertainty in the earth model can lead to unpredictable biases in the amplitude of the estimated seismic source functions.

Download Full-text

Parametrization of general seismic potency and moment tensors for source inversion of seismic waveform data

Geophysical Journal International ◽

10.1093/gji/ggt137 ◽

2013 ◽

Vol 194 (2) ◽

pp. 839-843 ◽

Cited By ~ 63

Author(s):

Lupei Zhu ◽

Yehuda Ben-Zion

Keyword(s):

Moment Tensor ◽

Relative Size ◽

Seismic Moment Tensor ◽

Source Inversion ◽

Waveform Data ◽

Moment Tensors ◽

Linear Vector ◽

Seismic Waveform ◽

Isotropic Media ◽

Double Couple

Abstract We decompose a general seismic potency tensor into isotropic tensor, double-couple tensor and compensated linear vector dipole using the eigenvectors and eigenvalues of the full tensor. Two dimensionless parameters are used to quantify the size of the isotropic and compensated linear vector dipole components. The parameters have well-defined finite ranges and are suited for non-linear inversions of source tensors from seismic waveform data. The decomposition and parametrization for the potency tensor are used to obtain corresponding results for a general seismic moment tensor. The relations between different parameters of the potency and moment tensors in isotropic media are derived. We also discuss appropriate specification of the relative size of different source components in inversions of seismic data.

Download Full-text

Seismic Moment Tensor Solutions from GeoNet Data to Provide a Moment Magnitude Scale for New Zealand

10.26686/wgtn.16947268 ◽

2021 ◽

Author(s):

◽

Elizabeth de Joux Robertson

Keyword(s):

New Zealand ◽

Seismic Moment ◽

Moment Tensor ◽

Velocity Model ◽

Moment Tensor Inversion ◽

Moment Magnitude ◽

Seismic Moment Tensor ◽

Waveform Data ◽

Moment Tensors ◽

The Moment

The aim of this project is to enable accurate earthquake magnitudes (moment magnitude, MW) to be calculated routinely and in near real-time for New Zealand earthquakes. This would be done by inversion of waveform data to obtain seismic moment tensors. Seismic moment tensors also provide information on fault-type. I use a well-established seismic moment tensor inversion method, the Time-Domain [seismic] Moment Tensor Inversion algorithm (TDMT_INVC) and apply it to GeoNet broadband waveform data to generate moment tensor solutions for New Zealand earthquakes. Some modifications to this software were made. A velocity model can now be automatically used to calculate Green's functions without having a pseudolayer boundary at the source depth. Green's functions can be calculated for multiple depths in a single step, and data are detrended and a suitable data window is selected. The seismic moment tensor solution that has either the maximum variance reduction or the maximum double-couple component is automatically selected for each depth. Seismic moment tensors were calculated for 24 New Zealand earthquakes from 2000 to 2005. The Global CMT project has calculated CMT solutions for 22 of these, and the Global CMT project solutions are compared to the solutions obtained in this project to test the accuracy of the solutions obtained using the TDMT_INVC code. The moment magnitude values are close to the Global CMT values for all earthquakes. The focal mechanisms could only be determined for a few of the earthquakes studied. The value of the moment magnitude appears to be less sensitive to the velocity model and earthquake location (epicentre and depth) than the focal mechanism. Distinguishing legitimate seismic signal from background seismic noise is likely to be the biggest problem in routine inversions.

Download Full-text

Seismic moment tensor and dynamic parameters of earthquakes in the Central Tien Shan

Geosystems of Transition Zones ◽

10.30730/gtrz.2020.4.2.178-191.192-209 ◽

2020 ◽

Vol 4 (2) ◽

pp. 178-209

Author(s):

N.A. Sycheva ◽

Keyword(s):

Seismic Moment ◽

Moment Tensor ◽

Tien Shan ◽

Seismic Moment Tensor ◽

Dynamic Parameters

Download Full-text