2.1 Data Used to Estimate the Coseismic Slip Distribution
We used the same coseismic data set as Elliott et al. (2022), along with the same inversion approach, and more details of the data processing and preparation are given there. We used static coseismic offsets in the ITRF2014 reference frame computed from the daily GNSS time series, processed with the GIPSY-OASIS gos-6.4 software (Zumberge et al., 1997; Bertiger et al., 2020). InSAR displacements, processed with GMTSAR (Sandwell et al., 2011), were included, and tied to GNSS sites where possible to provide absolute line-of-sight displacements. We included high rate GNSS, teleseismic broadband, and near-source strong-motion waveforms in the joint inversion. The 1 sps high-rate GPS time series were generated with GipsyX (Bertiger et al. , 2020) and filtered with a 0.4 Hz cut-off frequency to suppress noise. In addition, 46 P and 22 SH global tele-seismic waveforms were included in the inversion as well to improve observation geometry. The raw tele-seismic waveforms were integrated into ground displacements afterbandpass filtering between 0.01 ~ 1 Hz.