Multi-Array Multi-Phase Back-Projection: Improving the imaging of earthquake rupture complexities

We present a back-projection method capable of being parameterized with multiples arrays. The rupture imaging is weighted to restrict uncertainties induced by non-symmetric azimuthal coverage of seismic arrays. The strategy also exploits the differences in time delays between P and depth phase (pP) waveforms by assuming them as proxies of the rupture that can be simultaneously back-projected. Surprisingly, this helps to improve the final results, even when depth phases overlap with the direct arrivals due to the rupture time exceeding the pP-P delay. Thus, the approach heightens the spatiotemporal resolvability enough to image rupture complexities. The rupture image of two large events demonstrates its robustness. The first one is the 14 November 2007 Mw 7.7 Tocopilla earthquake in northern Chile. The high-frequency rupture (0.5 - 2.0 Hz) encircles two asperities while the short-period energy radiated predominates up-dip of the coseismic slip. We propose the contribution of asperity rupture complexities and along-dip barriers to high-frequency emissions beyond the megathrust frictional structure. The second one is the Mw 7.5 Palu strike-slip earthquake, which occurred on 28 September 2018 in Sulawesi island. The back-projection reveals a prominent supershear rupture at a speed of 4.5 km/s. The result correlates with space geodesy data highlighting the successful recovery of fault structures. Finally, we discuss the potential and challenges of automating this analysis for near-real-time applications, including near-source back-projection with strong-motion data.