Detailed constraints on rupture processes of large earthquakes are crucial for understanding fundamental earthquake physics. To image the rupture process of large earthquakes, back-projection (BP) of high frequency (HF) tele-seismic array waveforms is one of the commonly used techniques in the last 15 years or so (e.g., Kruger & Ohrnberger, 2005). Assuming delta-function shaped Green's functions and taking advantage of the reciprocity between sources and receivers, BP methods trace the location and timing of HF sources by either maximizing the stacking power of back-propagated array data in the time domain (e.g.,