Efficient Inversions for Earthquake Slip Distributions in 3D Structures

Hsieh, Ming Che; Zhao, Li; Ji, Chen; Fong, Kuo

doi:10.1785/0220160050

Cited by 11 publications

(11 citation statements)

References 36 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Any method for computing

\underset{\underset{scriptT}{false}}{false}

is suitable, as long as tensor values could be stored at points on the source plane for each receiver as proposed by Zhao et al () and Hsieh et al (). This off‐line computation requires a sampling which should honor the physical description of the expected rupture front width and the wave propagation at the frequencies involved in the rupture reconstruction: if finer sampling is needed by the rupture physics, trilinear interpolation of these tensor components is often enough.…”

Section: Methodological Descriptionmentioning

confidence: 99%

An Evolutive Linear Kinematic Source Inversion

et al. 2018

View full text Add to dashboard Cite

We present a new hierarchical time kinematic source inversion method able to assimilate data traces through evolutive time windows. A linear time domain formulation relates the slip rate function and seismograms, preserving the positivity of this function and the causality when spanning the model space: taking benefit of the time‐space sparsity of the rupture model evolution is as essential as considering the causality between rupture and each record delayed by the known propagator operator different for each station. This progressive approach, both on the data space and on the model space, does require mild assumptions on prior slip rate functions or preconditioning strategies on the slip rate local gradient estimations. These assumptions are based on simple physical expected rupture models we foresee. Successful applications of this method to a well‐known benchmark (Source Inversion Validation Exercise 1) and to the recorded data of the 2016 Kumamoto mainshock (Mw=7.0) illustrate the advantages of this alternative approach of a linear kinematic source inversion. This stabilized overparameterized optimization of a linear forward problem has a potential extension to stochastic inversion in a Bayesian framework in the future to assess uncertainties on kinematic source inversions, which is more difficult to investigate under nonlinear formulations.

show abstract

“…Any method for computing

\underset{\underset{scriptT}{false}}{false}

Section: Methodological Descriptionmentioning

confidence: 99%

An Evolutive Linear Kinematic Source Inversion

et al. 2018

View full text Add to dashboard Cite

show abstract

“…For inversions using teleseismic records, it is sufficient to use a global average one‐dimensional (1D) model with corrections for crustal structure. However, when using regional and shorter‐period records it may be necessary to calculate the Green's functions in three‐dimensional models (Zhao et al, 2006; Hsieh et al, 2016).…”

Section: Methods and Seismic Data For Slip‐distribution Inversionmentioning

confidence: 99%

“…Rather than directly fitting the time series of the broadband records, Ji et al (2002) proposes a wavelet‐decomposition approach in order to better capture the spectral‐temporal characteristics in the waveforms, and quantifies the residuals between model‐predicted synthetic waveforms and corresponding records using their respective wavelet coefficients. The objective function to be minimized in the inversion is defined by a L 1 +L 2 ‐norm for coefficients of relatively long‐period wavelets and a cross‐correlation for those of high‐frequency ones (Ji et al 2002; Hsieh et al 2016). This multi‐scale analysis of the seismic waveforms in the wavelet domain leads to enhanced spatial‐temporal resolution in source slip‐distribution inversions (Ji et al 2002).…”

Section: Methods and Seismic Data For Slip‐distribution Inversionmentioning

confidence: 99%

The 13 November 2016 Kaikoura, New Zealand Earthquake: rupture process and seismotectonic implications

Zhao

et al. 2018

Earth and Planetary Physics

Self Cite

View full text Add to dashboard Cite

The 13 November 2016 Kaikoura earthquake occurred in the northeastern coastal region of the South Island, New Zealand. The Mw 7.8 mainshock generated a complex pattern of surface ruptures, and was followed within about 12 hours by three moderate shocks of Mw ≥ 6.0. Here we use teleseismic waveforms to invert for the source rupture of the Kaikoura earthquake. The resulting slip‐distribution model exhibits insignificant slip near the hypocenter and three pockets of major slip zones with distinct senses of motion. The mainshock started from a rupture near the hypocenter, grew into thrust on shallow crustal faults ~50 km northeast of the hypocenter, and then developed into two slip zones: a deeper one with oblique thrust and a shallower one with almost purely right‐lateral strike‐slip. Locations of the thrust and strike‐slip motions in the slip‐distribution model agree well with reported coastal uplifts and horizontal offsets. The overall slip pattern is dominated by horizontal motion, especially at shallow depth, due to the partitioning of thrust and strike‐slip motions above the subduction zone megathrust. Aftershock distribution suggests that most aftershocks tend to occur near the edges of the major slip zones of the mainshock. This observation on aftershock locations may provide useful information for seismic hazard assessments after large earthquakes.

show abstract

“…The use of the reciprocity theorem is very effective to obtain a large number of synthetic seismograms for a small number of stations when there is a very large number of source grids, e.g., when calculating Green's function with reduced computational cost (e.g., Graves and Wald 2001;Eisner and Clayton 2001;Zhao et al 2006;Hsieh et al 2016;Petukhin et al 2016). OpenSWPC is equipped with this feature.…”

Section: Seismic Source Representationmentioning

confidence: 99%

OpenSWPC: an open-source integrated parallel simulation code for modeling seismic wave propagation in 3D heterogeneous viscoelastic media

Maeda

Takemura

Furumura

2017

Earth Planets Space

View full text Add to dashboard Cite

We have developed an open-source software package, Open-source Seismic Wave Propagation Code (OpenSWPC), for parallel numerical simulations of seismic wave propagation in 3D and 2D (P-SV and SH) viscoelastic media based on the finite difference method in local-to-regional scales. This code is equipped with a frequency-independent attenuation model based on the generalized Zener body and an efficient perfectly matched layer for absorbing boundary condition. A hybrid-style programming using OpenMP and the Message Passing Interface (MPI) is adopted for efficient parallel computation. OpenSWPC has wide applicability for seismological studies and great portability to allowing excellent performance from PC clusters to supercomputers. Without modifying the code, users can conduct seismic wave propagation simulations using their own velocity structure models and the necessary source representations by specifying them in an input parameter file. The code has various modes for different types of velocity structure model input and different source representations such as single force, moment tensor and plane-wave incidence, which can easily be selected via the input parameters. Widely used binary data formats, the Network Common Data Form (NetCDF) and the Seismic Analysis Code (SAC) are adopted for the input of the heterogeneous structure model and the outputs of the simulation results, so users can easily handle the input/output datasets. All codes are written in Fortran 2003 and are available with detailed documents in a public repository.

show abstract

Efficient Inversions for Earthquake Slip Distributions in 3D Structures

Cited by 11 publications

References 36 publications

An Evolutive Linear Kinematic Source Inversion

An Evolutive Linear Kinematic Source Inversion

The 13 November 2016 Kaikoura, New Zealand Earthquake: rupture process and seismotectonic implications

OpenSWPC: an open-source integrated parallel simulation code for modeling seismic wave propagation in 3D heterogeneous viscoelastic media

Contact Info

Product

Resources

About