Philippe Danré, Jiuxun Yin, Bradley P Lipovsky, and Marine A Denolle. 6/25/2019. “Earthquakes within Earthquakes: Patterns in Rupture Complexity.” Geophysical Research Letters, 46, 13, Pp. 7352-7360. Publisher's VersionAbstract
Earthquake source time functions carry information about the complexity of seismic rupture. We explore databases of earthquake source time functions and find that they are composed of distinct peaks that we call subevents. We observe that earthquake complexity, as represented by the number of subevents, grows with earthquake magnitude. Patterns in rupture complexity arise from a scaling between subevent moment and main event moment. These results can be explained by simple 2‐D dynamic rupture simulations with self‐affine heterogeneity in fault prestress. Applying this to early magnitude estimates, we show that the main event magnitude can be estimated after observing only the first few subevents.
2019gl083093.pdf grl59234-sup-0001-text_si-s01.pdf
Jiuxun Yin and Marine A. Denolle. 2019. “Relating teleseismic backprojection images to earthquake kinematics.” Geophysical Journal International, 217, 2, Pp. 729–747. Publisher's VersionAbstract

Backprojection (BP) of teleseismic P waves is a powerful tool to study the evolution of seismic radiation of large earthquakes. The common interpretations on the BP results are qualitative comparisons with earthquake kinematic observations, such as the evolution of slip on the fault and rupture velocity. However, the direct relation between the BP images and physical properties of the earthquake rupture process remains unclear and is needed for further application of this technique. In this study, we start from a theoretical formulation of the BP images, which is linear in the frequency domain, and carry on a synthetic exercise with kinematic source representations and virtual receivers embedded in a homogeneous full- space. We find that the fundamental linear formulation of the BP method is most correlated with the true kinematic source properties: in frequency domain the BP images are proportional to the images of slip motion through a scaling matrix F(ω) that accounts for radiation pattern and source–receiver geometry and that acts as a spatial smoothing operator. Overall, the synthetic BP images match relatively well the kinematic models and our exercise validates that the BP image can be directly used to track the spatiotemporal propagation of rupture front. However, because F(ω) is not strictly an identity matrix due to limited station coverage in space (azimuth and distance) and to the limited frequency bands of the seismograms, it remains difficult to recover the details in the rupture fronts from BP images. We define a resolvability parameter εI(ω) built from F(ω) that incorporates fault geometry, radiation pattern and wave propagation (source–array geometry) to quantify the ability of the BP method to resolve details of the rupture on the fault. εI(ω) successfully captures the similarity between BP images and kinematic source. We analyse the resolvability of most tectonically active regions and the most commonly used seismic arrays. Based on this global resolvability analysis, we propose an empirical relation between the seismic frequency, resolvable area and earthquake magnitude. It provides general guidelines to choose the lowest frequency in seismic waveform (e.g. about 0.3 Hz for Mw 8 and 1 Hz for Mw 7 earthquakes) and to interpret the BP image in terms of the source kinematics. In general, this work attempts to provide a clear interpretation of the BP images in light of the real earthquake rupture process and give a systematic evaluation of seismic data limitations.

linear_bp_yin_denolle_2019.pdf supplements_2019_GJI.pdf
PJ Danré, J Yin, BP Lipovsky, and M Denolle. 2018. “Earthquake Sub-event Scaling: New Perspective for Rupture Determinism.” In AGU Fall Meeting Abstracts.
Zhuo Yang, Jiuxun Yin, and Marine Denolle. 2018. “Spatial and temporal evolution of seismicity, seismic velocity, and pore pressure in the Guy-Greenbrier, Arkansas, earthquake sequence.” In AGU Fall Meeting Abstracts.
Jiuxun Yin, Zefeng Li, Bradley Paul Lipovsky, and Marine Denolle. 2018. “Connecting source time function and heterogeneous fault properties.” In AGU Fall Meeting Abstracts. Publisher's Version
J Yin and M Denolle. 2018. “Resolvability of linear backprojection method for earthquake kinematics.” In AGU Fall Meeting Abstracts. yin_2018_bp_resolution.pdf
Jiuxun Yin, Marine A. Denolle, and Huajian Yao. 2018. “Spatial and Temporal Evolution of Earthquake Dynamics: Case Study of the Mw 8.3 Illapel Earthquake, Chile.” Journal of Geophysical Research: Solid Earth, 123, 1, Pp. 344–367. Publisher's VersionAbstract
Abstract We develop a methodology that combines compressive sensing backprojection (CS-BP) and source spectral analysis of teleseismic P waves to provide metrics relevant to earthquake dynamics of large events. We improve the CS-BP method by an autoadaptive source grid refinement as well as a reference source adjustment technique to gain better spatial and temporal resolution of the locations of the radiated bursts. We also use a two-step source spectral analysis based on (i) simple theoretical Green's functions that include depth phases and water reverberations and on (ii) empirical P wave Green's functions. Furthermore, we propose a source spectrogram methodology that provides the temporal evolution of dynamic parameters such as radiated energy and falloff rates. Bridging backprojection and spectrogram analysis provides a spatial and temporal evolution of these dynamic source parameters. We apply our technique to the recent 2015 Mw 8.3 megathrust Illapel earthquake (Chile). The results from both techniques are consistent and reveal a depth-varying seismic radiation that is also found in other megathrust earthquakes. The low-frequency content of the seismic radiation is located in the shallow part of the megathrust, propagating unilaterally from the hypocenter toward the trench while most of the high-frequency content comes from the downdip part of the fault. Interpretation of multiple rupture stages in the radiation is also supported by the temporal variations of radiated energy and falloff rates. Finally, we discuss the possible mechanisms, either from prestress, fault geometry, and/or frictional properties to explain our observables. Our methodology is an attempt to bridge kinematic observations with earthquake dynamics.
yin_et_al-2018-journal_of_geophysical_research_solid_earth.pdf yin_et_al-2018-journal_of_geophysical_research_solid_earth.sup-1.pdf yin_et_al-2018-journal_of_geophysical_research_solid_earth.sup-2.pdf yin_et_al-2018-journal_of_geophysical_research_solid_earth.sup-3.avi
Jiuxun Yin, Huajian Yao, Hongfeng Yang, Jing Liu, Weize Qin, and Haijiang Zhang. 2017. “Frequency-dependent rupture process, stress change, and seismogenic mechanism of the 25 April 2015 Nepal Gorkha M w 7.8 earthquake.” Science China Earth Sciences, 60, 4, Pp. 796–808. Publisher's Version yin_etal_2017_sces_nepaleq.pdf
Jiuxun Yin and M Denolle. 2017. “Relating backprojection images to kinematics and dynamic source models.” In AGU Fall Meeting Abstracts. agu2017_poster.pdf
J Yin, M Denolle, H. Yao, and H Yang. 2016. “4D source analysis of the 2015 Illapel, Chile Mw8. 3 Earthquake.” In AGU Fall Meeting Abstracts. yin_2016_src_observation.pdf
Jiuxun Yin, Hongfeng Yang, Huajian Yao, and Huihui Weng. 2016. “Coseismic radiation and stress drop during the 2015 Mw 8.3 Illapel, Chile megathrust earthquake.” Geophysical Research Letters, 43, 4, Pp. 1520–1528. Publisher's Version 2015_chile_mw8.3_grl.pdf grl53998-sup-0001-s01.pdf
Jiuxun Yin and Huajian Yao. 2016. “Rupture and frequency-dependent seismic radiation of the 2012 M w 8.6 Sumatra strike-slip earthquake.” Geophysical Journal International, 205, 3, Pp. 1682–1693. Publisher's Version geophys._j._int.-2016-yin-1682-93.pdf geophys._j._int.-2016-yin-1682-93_si.pdf
Hongfeng Yang, Jian LIN, Jiuxun Yin, and Huajian Yao. 2015. “Tectonic settings of the 2015 M w 8.3 Coquimbo, Chile earthquake and its implications on megathrust earthquakes.” Chinese Science Bulletin, 60, 36, Pp. 3549–3556.Abstract
The 2015 Mw8.3 Coquimbo, Chile earthquake is a typical megathrust earthquake, whose size and coseismic slip distribution are consistent with the interseismic locking model that was derived from GPS measurements. Our preliminary back-projection results show that the rupture of the Mw8.3 Coquimbo earthquake propagates updip from the hypocenter (~25 km in depth). Furthermore, we find frequency-dependent behavior of the radiation power of the rupture, similar to those observed during the 2010 Mw8.8 Maule earthquake, which occurred ~50 km to the south of the 2015 earthquake rupture area. These observations indicate that there is a systematic downdip variation of properties on the megathrust in southern Chile. The rich observations of the coseismic slip distributions of the 2015 Mw8.3 earthquake provide an excellent opportunity to validate numerical simulations of rupture process based on interseismic locking distributions.