Extended differential CAP method for earthquake source parameter inversion with high-rate GNSS relative positioning

SUMMARY Near-field large-amplitude seismograms are essential for the rapid inversion of earthquake source parameters using waveform inversion methods such as the Cut And Paste (CAP) method for disaster assessment. High-rate Global Navigation Satellite System (GNSS) relative positioning (RP) provides precise, rapid, and real-time measurement of near-field large-amplitude displacements. However, RP records motion with respect to a reference station, and the reference station's movements become part of the relative displacement waveforms. Therefore, seismic source parameter estimates may be inaccurate if the reference station's motion is not taken into consideration, and doing so affects some basic assumptions made in the CAP method. To overcome this problem, we develop an expanded differential CAP inversion approach specifically for high-rate GNSS RP (CAP-RP) that accounts for the motion of the reference station. Two methods are proposed to implement CAP-RP: an expanded differential CAP (D-CAP) and an iterative post-processing CAP (P-CAP). We assess the performance of CAP-RP with different data sets, using the July 2019 Mw 6.4 earthquake in California as a case study. Both CAP-RP techniques produce accurate source parameters in synthetic data inversion tests, indicating the feasibility of the strategy. However, P-CAP is more time-efficient than D-CAP, making it the better option. Generally, results from high-rate GNSS RP, broad-band seismographs, and their inverted combinations exhibit consistency in observational data inversion testing. Our results also demonstrate that more accurate source parameters can be obtained by combining sensitive far-field broad-band seismograph data with large amplitude near-field GNSS RP waveforms.