scholarly journals The waveform inversion of mainshock and aftershock data of the 2006 M6.3 Yogyakarta earthquake

2021 ◽  
Vol 8 (1) ◽  
Author(s):  
Hijrah Saputra ◽  
Wahyudi Wahyudi ◽  
Iman Suardi ◽  
Ade Anggraini ◽  
Wiwit Suryanto
Keyword(s):  
Joint Inversion ◽  
Body Wave ◽  
Moment Tensor ◽  
Near Field ◽  
Source Parameters ◽  
Waveform Inversion ◽  
Moment Tensor Inversion ◽  
Slip Distribution ◽  
Aftershock Distribution ◽  
Velocity Models

AbstractThis study comprehensively investigates the source mechanisms associated with the mainshock and aftershocks of the Mw = 6.3 Yogyakarta earthquake which occurred on May 27, 2006. The process involved using moment tensor inversion to determine the fault plane parameters and joint inversion which were further applied to understand the spatial and temporal slip distributions during the earthquake. Moreover, coseismal slip distribution was overlaid with the relocated aftershock distribution to determine the stress field variations around the tectonic area. Meanwhile, the moment tensor inversion made use of near-field data and its Green’s function was calculated using the extended reflectivity method while the joint inversion used near-field and teleseismic body wave data which were computed using the Kikuchi and Kanamori methods. These data were filtered through a trial-and-error method using a bandpass filter with frequency pairs and velocity models from several previous studies. Furthermore, the Akaike Bayesian Information Criterion (ABIC) method was applied to obtain more stable inversion results and different fault types were discovered. Strike–slip and dip-normal were recorded for the mainshock and similar types were recorded for the 8th aftershock while the 9th and 16th June were strike slips. However, the fault slip distribution from the joint inversion showed two asperities. The maximum slip was 0.78 m with the first asperity observed at 10 km south/north of the mainshock hypocenter. The source parameters discovered include total seismic moment M0 = 0.4311E + 19 (Nm) or Mw = 6.4 with a depth of 12 km and a duration of 28 s. The slip distribution overlaid with the aftershock distribution showed the tendency of the aftershock to occur around the asperities zone while a normal oblique focus mechanism was found using the joint inversion.

scholarly journals The waveform inversion of mainshock and aftershock data of the 2006 M6.3 Yogyakarta earthquake

2021 ◽  
Author(s):  
Hijrah Saputra ◽  
Wahyudi Wahyudi ◽  
Iman Suardi ◽  
Ade Anggraini ◽  
Wiwit Suryanto
Keyword(s):  
Joint Inversion ◽  
Body Wave ◽  
Moment Tensor ◽  
Near Field ◽  
Source Parameters ◽  
Moment Tensor Inversion ◽  
Slip Distribution ◽  
Coseismic Slip ◽  
Aftershock Distribution ◽  
Velocity Models

Abstract This study comprehensively investigates the source mechanisms associated with the mainshock and aftershocks of the Yogyakarta earthquake of magnitude Mw = 6.3 on May 27, 2006. Therefore, this study is to provide a more precise answer to the controversial source mechanism. This study uses moment tensor inversion to obtain fault plane parameters and joint inversion to obtain spatial and temporal slip distributions during an earthquake. The coseismic slip distribution is overlaid with the relocated aftershock distribution to see the stress field variations around the tectonic area of the study. Moment tensor inversion uses near-field data, and joint inversion uses near-field and teleseismic body wave data. The data is filtered by trial and error using a bandpass filter with frequency pairs and velocity models from several previous studies. The green's function for moment tensor inversion calculated using the extended reflectivity method and joint inversion computed using the Kikuchi and Kanamori methods. In this study, we apply the Akaike Bayesian Information Criterion (ABIC) method to obtain more stable inversion results. The results of the mainshock and aftershock moment tensor inversion show different fault types. The mainshock fault types are strike-slip and dip-normal types, while the 8th aftershock is of the same type as the mainshock, while the 9th and 16th June are strike slips. The joint inversion results show two asperities. The maximum slip is 0.78 m, with the first asperity 10 km south of the mainshock and the second asperity 10 km north of the mainshock. The obtained source parameters are total seismic moment M0 = 0.4311E + 19 (Nm) or Mw = 6.4, with a source depth of 12 km and a source duration of 28 seconds. Slip distribution overlay with aftershock distribution shows compatibility. The type of focus mechanism that results from this joint inversion is the oblique.

Waveform inversion for microseismic velocity analysis and event location in VTI media

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. WA95-WA103 ◽  
Author(s):  
Oscar Jarillo Michel ◽  
Ilya Tsvankin
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Waveform Inversion ◽  
Transversely Isotropic ◽  
Velocity Analysis ◽  
Origin Time ◽  
Velocity Models ◽  
Trade Offs ◽  
Receiver Array ◽  
Event Location

Waveform inversion (WI), which has been extensively used in reflection seismology, could provide improved velocity models and event locations for microseismic surveys. Here, we develop an elastic WI algorithm for anisotropic media designed to estimate the 2D velocity field along with the source parameters (location, origin time, and moment tensor) from microseismic data. The gradient of the objective function is obtained with the adjoint-state method, which requires just two modeling simulations at each iteration. In the current implementation the source coordinates and velocity parameters are estimated sequentially at each stage of the inversion to minimize trade-offs and improve the convergence. Synthetic examples illustrate the accuracy of the inversion for layered VTI (transversely isotropic with a vertical symmetry axis) media, as well as the sensitivity of the velocity-analysis results to noise, the length of the receiver array, errors in the initial model, and variability in the moment tensor of the recorded events.

scholarly journals Imaging the subsurface using induced seismicity and ambient noise: 3-D tomographic Monte Carlo joint inversion of earthquake body wave traveltimes and surface wave dispersion

2020 ◽  
Vol 222 (3) ◽  
pp. 1639-1655
Author(s):  
Xin Zhang ◽  
Corinna Roy ◽  
Andrew Curtis ◽  
Andy Nowacki ◽  
Brian Baptie
Keyword(s):  
Surface Wave ◽  
Ambient Noise ◽  
Velocity Structure ◽  
Joint Inversion ◽  
Body Wave ◽  
Source Parameters ◽  
Wave Dispersion ◽  
Inversion Method ◽  
Surface Wave Dispersion ◽  
Wave Data

SUMMARY Seismic body wave traveltime tomography and surface wave dispersion tomography have been used widely to characterize earthquakes and to study the subsurface structure of the Earth. Since these types of problem are often significantly non-linear and have non-unique solutions, Markov chain Monte Carlo methods have been used to find probabilistic solutions. Body and surface wave data are usually inverted separately to produce independent velocity models. However, body wave tomography is generally sensitive to structure around the subvolume in which earthquakes occur and produces limited resolution in the shallower Earth, whereas surface wave tomography is often sensitive to shallower structure. To better estimate subsurface properties, we therefore jointly invert for the seismic velocity structure and earthquake locations using body and surface wave data simultaneously. We apply the new joint inversion method to a mining site in the United Kingdom at which induced seismicity occurred and was recorded on a small local network of stations, and where ambient noise recordings are available from the same stations. The ambient noise is processed to obtain inter-receiver surface wave dispersion measurements which are inverted jointly with body wave arrival times from local earthquakes. The results show that by using both types of data, the earthquake source parameters and the velocity structure can be better constrained than in independent inversions. To further understand and interpret the results, we conduct synthetic tests to compare the results from body wave inversion and joint inversion. The results show that trade-offs between source parameters and velocities appear to bias results if only body wave data are used, but this issue is largely resolved by using the joint inversion method. Thus the use of ambient seismic noise and our fully non-linear inversion provides a valuable, improved method to image the subsurface velocity and seismicity.

Source parameter determination of local earthquakes in Korea using moment tensor inversion of single station data

1999 ◽  
Vol 89 (4) ◽  
pp. 1077-1082 ◽  
Author(s):  
So Gu Kim ◽  
Nadeja Kraeva
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Moment Tensor Inversion ◽  
P Wave ◽  
Parameter Determination ◽  
Reverse Fault ◽  
3D Tomography ◽  
Single Station ◽  
Slip Event

Abstract The purpose of this investigation is to determine source parameters such as focal mechanism, seismic moment, moment magnitude, and source depth from recent small earthquakes in the Korcan Peninsula using broadband records of three-component single station. It is very important and worthwhile to use a three-component single station in Korea because for most Korean earthquakes it is not possible to read enough first motions of P-wave arrivals because of the poor coverage of the seismic network and the small size (ML 5.0 or less) of the events. Furthermore the recent installation of the very broadband seismic stations in Korea and use of a 3D tomography technique can enhance moment tensor inversion to determine the source parameters of small earthquakes (ML 5.0 or less) that occur at near-regional distances (Δ ≤ 500 km). The focal solution for the Youngwol earthquake of 13 December 1996 is found to be a right-lateral strike slip event with a NE strike, and the Kyongju earthquake of 25 June 1997 is found to be an oblique reverse fault with a slight component of left-lateral slip in the SE direction.

Quick determination of earthquake source parameters from GPS measurements: a study of suitability for Taiwan

2019 ◽  
Vol 219 (2) ◽  
pp. 1148-1162
Author(s):  
Jiun-Ting Lin ◽  
Wu-Lung Chang ◽  
Diego Melgar ◽  
Amanda Thomas ◽  
Chi-Yu Chiu
Keyword(s):  
Rapid Determination ◽  
Moment Tensor ◽  
Near Field ◽  
Source Parameters ◽  
Early Warning Systems ◽  
Earthquake Magnitude ◽  
Observation System ◽  
Earthquake Source Parameters ◽  
Finite Fault

SUMMARY We test the feasibility of GPS-based rapid centroid moment tensor (GPS CMT) methods for Taiwan, one of the most earthquake prone areas in the world. In recent years, Taiwan has become a leading developer of seismometer-based earthquake early warning systems, which have successfully been applied to several large events. The rapid determination of earthquake magnitude and focal mechanism, important for a number of rapid response applications, including tsunami warning, is still challenging because of the limitations of near-field inertial recordings. This instrumental issue can be solved by an entirely different observation system: a GPS network. Taiwan is well posed to take advantage of GPS because in the last decade it has developed a very dense network. Thus, in this research, we explore the suitability of the GPS CMT inversion for Taiwan. We retrospectively investigate six moderate to large (Mw6.0 ∼ 7.0) earthquakes and propose a resolution test for our model, we find that the minimum resolvable earthquake magnitude of this system is ∼Mw5.5 (at 5 km depth). Our tests also suggest that the finite fault complexity, often challenging for the near-field methodology, can be ignored under such good station coverage and thus, can provide a fast and robust solution for large earthquake directly from the near field. Our findings help to understand and quantify how the proposed methodology could be implemented in real time and what its contributions could be to the overall earthquake monitoring system.

scholarly journals Adaptive moment-tensor joint inversion of clustered microseismic events for monitoring geological carbon storage

2019 ◽  
Vol 219 (1) ◽  
pp. 80-93
Author(s):  
Yu Chen ◽  
Lianjie Huang
Keyword(s):  
Carbon Storage ◽  
Oil Recovery ◽  
Joint Inversion ◽  
Moment Tensor ◽  
Crack Opening ◽  
Cost Effective ◽  
Microseismic Monitoring ◽  
Moment Tensor Inversion ◽  
Inversion Method ◽  
Microseismic Events

SUMMARY Moment-tensor inversion of induced microseismic events can provide valuable information for tracking CO2 plumes at geological carbon storage sites, and study the physical mechanism of induced microseismicity. Accurate moment-tensor inversion requires a wide-azimuthal coverage of geophones. Cost-effective microseismic monitoring for geological carbon storage often uses only one geophone array within a borehole, leading to a large uncertainty in moment-tensor inversion. We develop a new adaptive moment-tensor joint inversion method to reduce the inversion uncertainty, when using limited but typical geophone receiver geometries. We first jointly invert a number of clustered microseismic events using a uniform focal mechanism to minimize the waveform misfit between observed and predicted P and S waveforms. We then invert the moment tensor for each event within a limited searching range around the joint inversion result. We apply our adaptive joint inversion method to microseismic data acquired using a single borehole geophone array at the CO2-Enhanced Oil Recovery field at Aneth, Utah. We demonstrate that our inversion method is capable of reducing the inversion uncertainty caused by the limited azimuthal coverage of geophones. Our inverted strikes of focal mechanisms of microseismic events are consistent with the event spatial distribution in subparallel pre-existing fractures or geological imperfections. The large values up to 40 per cent of the CLVD components might indicate crack opening induced by CO2/wastewater injection or rupture complexity.

scholarly journals Accuracy of the moment-tensor inversion of far-field P waves

2019 ◽  
Vol 220 (1) ◽  
pp. 248-256 ◽  
Author(s):  
Yue Kong ◽  
Min Li ◽  
Weimin Chen ◽  
Boqi Kang
Keyword(s):  
Moment Tensor ◽  
Near Field ◽  
Radial Component ◽  
Moment Tensor Inversion ◽  
Far Field ◽  
Field Range ◽  
P Waves ◽  
Moment Tensors ◽  
Nodal Lines ◽  
The Moment

SUMMARY The far-field assumption is widely used and suitable for the moment-tensor inversion in which the source–receiver distance is quite long. However, the description of far field is uncertain and an explicit far-field range is missing. In this study, the explicit far-field range is determined and the errors of moment-tensor solutions produced by the far-field approximation are analysed. The distance, for which the far-field assumption is satisfied and the effect of the near-field term can be ignored, is directionally dependent. For the shear dislocation, in the directions near the nodal lines of the far-field P waves, the far-field distance is heavily dependent on the displacement component used to invert moment tensors. The radial component of displacement, which is parallel to the wave propagation direction, is recommended for the inversion and the corresponding far-field distance is quite short. In the directions far from the nodal lines, the selection of displacement components has little influence on the far-field distance. The maximum far-field distance appears in the directions of the pressure and tensile axes of the source and the value is about 30 wavelengths of radiated waves. Using more receivers (>6) in the moment-tensor inversion can shorten the far-field distance. The effect of the near-field term on the moment-tensor inversion for tensile dislocations and isotropic sources (explosion or implosion) can be ignored. The conclusions obtained in this study are helpful for determining the positions of receivers and evaluating the accuracy of moment-tensor solutions, with far-field assumption being applied in the inversion.

Near-field strain in DAS-based microseismic observation

Geophysics ◽  
2021 ◽  
pp. 1-49
Author(s):  
Ge Jin ◽  
Frantisek Stanek ◽  
Bin Luo
Keyword(s):  
Hydraulic Fracture ◽  
Fiber Optic ◽  
Moment Tensor ◽  
Near Field ◽  
Source Parameters ◽  
Field Strain ◽  
Far Field ◽  
Normal Strain ◽  
S Wave ◽  
Sign Patterns

Microseismic monitoring with surface or downhole geophone arrays has been commonly used in tracking subsurface deformation and fracture networks during hydraulic fracturing operations. Recently, the use of fiber-optic DAS technology has improved microseismic acquisition to a new level with unprecedentedly high spatial resolution and low cost. Deploying fiber-optic cables in horizontal boreholes allows very close observation of these micro-sized earthquakes and captures their full wavefield details. We show that DAS-based microseismic profiles present a seldomly reported near-field strain signal between the P- and S-wave arrivals. This near-field signal shows monotonically increasing (or decreasing) temporal variation, which resembles the previously reported near-field observations of large earthquakes. To understand the near-field strain behavior, we provide a mathematical expression of the analytic normal strain solution that reveals the near-field, intermediate-near-field, intermediate-far-field, and far-field components. Synthetic DAS strain records of hydraulic-fracture-induced microseismic events can be generated using this analytic solution with the Brune source model. The polarity sign patterns of the near-field and far-field terms in these synthetics are linked to the corresponding source mechanism’s radiation patterns. These polarity sign patterns are demonstrated to be sensitive to the source orientations by rotating the moment tensor in different directions. A field data example is compared to the synthetic result and a qualitative match is shown. The microseismic near-field signals detected by DAS have potential value in hydraulic fracture monitoring by providing a means to better constrain microseismic source parameters that characterize the source magnitude, source orientation, and temporal source evolution, and therefore better reflect the geomechanical response of the hydraulically fractured environment in the unconventional reservoirs.

A user-friendly probabilistic earthquake source inversion framework for joint inversion of seismic, geodetic, and gravitational signals - The Grond toolkit

2020 ◽  
Author(s):  
Sebastian Heimann ◽  
Marius Isken ◽  
Daniela Kühn ◽  
Hannes Vasyura-Bathke ◽  
Henriette Sudhaus ◽  
...  
Keyword(s):  
Open Source ◽  
Joint Inversion ◽  
Moment Tensor ◽  
Source Parameters ◽  
Seismic Source ◽  
Magma Ascent ◽  
Caldera Collapse ◽  
Waveform Data ◽  
Trade Offs ◽  
Finite Fault

<p>Seismic source and moment tensor waveform inversion is often ill-posed or non-unique if station coverage is poor or signals are weak. Three key ingredients can help in these situations: (1) probabilistic inference and global search of the full model space, (2) joint optimisation with datasets yielding complementary information, and (3) robust source parameterisation or additional source constraints. These demands lead to vast technical challenges, on the performance of forward modelling, on the optimisation algorithms, as well as on visualisation, optimisation configuration, and management of the datasets. Implementing a high amount of automation is inevitable.</p><p>To tackle all these challenges, we are developing a sophisticated new seismic source optimisation framework, Grond. With its innovative Bayesian bootstrap optimiser, it is able to efficiently explore large model spaces, the trade-offs and the uncertainties of source parameters. The program is highly flexible with respect to the adaption to specific source problems, the design of objective functions, and the diversity of empirical datasets.</p><p>It uses an integrated, robust waveform data processing, and allows for interactive visual inspection of many aspects of the optimisation problem, including visualisation of the result uncertainties. Grond has been applied to CMT moment tensor and finite-fault optimisations at all scales, to nuclear explosions, to a meteorite atmospheric explosion, and to volcano-tectonic processes during caldera collapse and magma ascent. Hundreds of seismic events can be handled in parallel given a single optimisation setup.</p><p>Grond can be used to optimise simultaneously seismic waveforms, amplitude spectra, waveform features, phase picks, static displacements from InSAR and GNSS, and gravitational signals.</p><p>Grond is developed as an open-source package and community effort. It builds on and integrates with other established open-source packages, like Kite (for InSAR) and Pyrocko (for seismology).</p>

Loreto Intermediate Depth Earthquake of 26 May 2019 (Northeast Peru): Source Parameters by Inversion of Local to Regional Waveforms and Strong-Motion Observations

2021 ◽  
Author(s):  
Hernando Tavera ◽  
Bertrand Delouis ◽  
Arturo Mercado ◽  
David Portugal
Keyword(s):  
Near Field ◽  
Source Parameters ◽  
Waveform Inversion ◽  
Strong Motion ◽  
Seismic Source ◽  
Slab Geometry ◽  
Nazca Plate ◽  
Strong Motion Data ◽  
Ground Shaking ◽  
Motion Data

Abstract The Loreto earthquake of 26 May 2019 occurred below the extreme northeast part of Peru at a depth of 140 km within the subducting Nazca plate at a distance of 700 km from the trench Peru–Chile. The orientation of the seismic source was obtained from waveform inversion in the near field using velocity and strong-motion data. The rupture occurred in normal faulting corresponding to a tensional process with T axis oriented in east–west direction similar to the direction of convergence between the Nazca and South America plates. The analysis of the strong-motion data shows that the levels of ground shaking are very heterogeneous with values greater than 50 Gal up to distances of 300 km; the maximum recorded acceleration of 122 Gal at a distance of 100 km from the epicenter. The Loreto earthquake is classified as a large extensional event in the descending Nazca slab in the transition from flat-slab geometry to greater dip.

Sign in / Sign up

Export Citation Format

Share Document