Abstract:Abstract Precisely constraining the source parameters of large earthquakes is one of the primary objectives of seismology. However, the quality of the results relies on the quality of synthetic earth response. Although earth structure is laterally heterogeneous, particularly at shallow depth, most earthquake source studies at the global scale rely on the Green's functions calculated with radially symmetric (1D) earth structure. To avoid the impact of inaccurate Green's functions, these conventional source studies use a limited set of seismic phases, such as long period seismic waves, broadband P and S waves in teleseismic distances (30 degrees < ∆ < 90 degrees), and strong ground motion records at close-fault stations. The enriched information embedded in the broadband seismograms recorded by global and regional networks is largely ignored, limiting the spatiotemporal resolution. Here we calculate 3D strain Green's functions at 30 GSN stations for source regions of 10 selected global earthquakes, with frequency up to 67 mill hertz (15 s), using SPECFEM3D_GLOBE and the reciprocity theorem. The 3D SEM mesh model is composed of mantle model S40RTS, crustal model CRUST2.0, and surface topography ETOPO2. We surround each target event with grids in horizontal spacing of 5 km and vertical spacing of 2.0∼3.0 km, allowing us to investigate not only the mainshock but also the background seismicity. In total, the response at over 210,000 source points is calculated in simulation. The number of earthquakes, including different focal mechanisms, centroid depth range, and tectonic background, could further increase without additional computational cost if they were properly selected to avoid overloading individual CPUs. The storage requirement can be reduced by two orders of magnitude if the output strain Green's functions are stored for periods over 15 s. We quantitatively evaluate the quality of these 3D synthetic seismograms, which are frequency and phase dependent, for each source region using nearby aftershocks, before using them to constrain the focal mechanisms and slip distribution. Case studies show that using a 3D earth model significantly improves the waveform similarity, agreement in amplitude and arrival time of seismic phases with the observations. The limitations of current 3D models are still notable, dependent on seismic phases and frequency range. The 3D synthetic seismograms cannot well match the high frequency (>40mHz) S wave and (>20mHz) Rayleigh wave yet. Though the mean time-shifts are close to zero, the standard deviations are notable. Careful calibration using the records of nearby better located earthquakes is still recommended to take full advantage of better waveform similarity due to the use of 3D models. Our results indicate that it is now feasible to systematically study global large earthquakes using full 3D earth response in a global scale.

A Reciprocity‐Based Efficient Method for Improved Source Parameter Estimation of Submarine Earthquakes With Hybrid 3‐D Teleseismic Green's Functions

Stochastic Green’s Function Method Considering Non-Uniform Rise Time Distribution to Simulate 3D Broadband Ground Motion

Improvement of Stochastic Green's Function Method for 3D Broadband Ground-Motion Simulation

Improvement of Source Spectrum Fall‐Off for Simulating Ground Motion Using Stochastic Green’s Function Method

Source parameters inversion of the global large earthquakes using 3D SEM Green's functions: Strain Green's function calculation and validation

Estimating Curie-point Depths Using Both Wavelet-Based and Fourier Spectral Centroid Methods in the Western Pacific Marginal Seas

Efficient Waveform Inversion for Average Earthquake Rupture in Three-Dimensional Structures

A Broadband Kinematic Source Inversion Method Considering Realistic Earth Models and Its Application to the 1992 Landers Earthquake

Resolving Shallow Moment Releases of Megathrust Earthquakes by Empirical Green’s Function‐Based Inversions: A New Approach for Tsunami Warning

An efficient elastic full‐waveform inversion of multiple parameters with ocean‐bottom seismometer data

Efficient Inversions for Earthquake Slip Distributions in 3D Structures

Unified Regional Tomography and Source Moment Tensor Inversions Based on Finite-Difference Strain Green Tensor Databases

Resolving Complicated Faulting Process Using Multi‐Point‐Source Representation: Iterative Inversion Algorithm Improvement and Application to Recent Complex Earthquakes

Quantifying Uncertainties in the 2004 Sumatra-Andaman Earthquake Source Parameters by Stochastic Inversion

Micro-seismic Elastic Reflection Full Waveform Inversion with An Equivalent Source

The Earthquake‐Source Inversion Validation (SIV) Project

Rapid Kinematic Finite-Fault Inversion for an M W 7+ Scenario Earthquake in the Marmara Sea: an Uncertainty Study

Analysis of the Reliability and Resolution of the Earthquake Source History Inferred from Waveforms, Taking the Chi-Chi Earthquake As an Example

Validation of Linearity Assumptions for Using Tsunami Waveforms in Joint Inversion of Kinematic Rupture Models: Application to the 2010 Mentawai M-W 7.8 Tsunami Earthquake

Effects of Dispersion in Tsunami Green'S Functions and Implications for Joint Inversion with Seismic and Geodetic Data: A Case Study of the 2010 Mentawai M-W 7.8 Earthquake

Building Precise Local Submarine Earthquake Catalogs via a Deep-Learning-Empowered Workflow and its Application to the Challenger Deep