@article {danreearthquakes,
	title = {Earthquakes within Earthquakes: Patterns in Rupture Complexity},
	journal = {Geophysical Research Letters},
	volume = {46},
	number = {13},
	year = {2019},
	pages = {7352-7360},
	abstract = {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.},
	url = {https://doi.org/10.1029/2019GL083093},
	author = {Danr{\'e}, Philippe and Yin, Jiuxun and Bradley P Lipovsky and Denolle, Marine A}
}
@article {yin_relating_2019,
	title = {Relating teleseismic backprojection images to earthquake kinematics},
	journal = {Geophysical Journal International},
	volume = {217},
	number = {2},
	year = {2019},
	pages = {729{\textendash}747},
	abstract = {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{\textendash}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{\textendash}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.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggz048},
	url = {https://academic-oup-com.ezp-prod1.hul.harvard.edu/gji/article/217/2/729/5301380},
	author = {Yin, Jiuxun and Denolle, Marine A.}
}
@conference {danre2018earthquake,
	title = {Earthquake Sub-event Scaling: New Perspective for Rupture Determinism},
	booktitle = {AGU Fall Meeting Abstracts},
	year = {2018},
	author = {Danr{\'e}, PJ and Yin, J and Lipovsky, BP and Denolle, M}
}
@conference {yang2018spatial,
	title = {Spatial and temporal evolution of seismicity, seismic velocity, and pore pressure in the Guy-Greenbrier, Arkansas, earthquake sequence},
	booktitle = {AGU Fall Meeting Abstracts},
	year = {2018},
	author = {Yang, Zhuo and Yin, Jiuxun and Denolle, Marine}
}
@conference {yin2018connecting,
	title = {Connecting source time function and heterogeneous fault properties},
	booktitle = {AGU Fall Meeting Abstracts},
	year = {2018},
	url = {https://agu2018fallmeeting-agu.ipostersessions.com/default.aspx?s=49-3B-28-35-09-42-C5-2B-23-A7-2A-3C-32-36-27-43},
	author = {Yin, Jiuxun and Zefeng Li and Bradley Paul Lipovsky and Denolle, Marine}
}
@conference {yin2018resolvability,
	title = {Resolvability of linear backprojection method for earthquake kinematics},
	booktitle = {AGU Fall Meeting Abstracts},
	year = {2018},
	author = {Yin, J and Denolle, M}
}
@article {yin_spatial_2018,
	title = {Spatial and Temporal Evolution of Earthquake Dynamics: Case Study of the Mw 8.3 Illapel Earthquake, Chile},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {1},
	year = {2018},
	month = {feb},
	pages = {344{\textendash}367},
	abstract = {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{\textquoteright}s functions that include depth phases and water reverberations and on (ii) empirical P wave Green{\textquoteright}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.},
	keywords = {backprojection, earthquake dynamics, Illapel Earthquake, radiated energy},
	issn = {2169-9313},
	doi = {10.1002/2017JB014265},
	url = {https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1002/2017JB014265},
	author = {Yin, Jiuxun and Denolle, Marine A. and Yao, Huajian}
}
@article {yin2017frequency,
	title = {Frequency-dependent rupture process, stress change, and seismogenic mechanism of the 25 April 2015 Nepal Gorkha M w 7.8 earthquake},
	journal = {Science China Earth Sciences},
	volume = {60},
	number = {4},
	year = {2017},
	pages = {796{\textendash}808},
	publisher = {Springer},
	url = {https://link.springer.com/article/10.1007/s11430-016-9006-0},
	author = {Yin, Jiuxun and Yao, Huajian and Yang, Hongfeng and Liu, Jing and Qin, Weize and Zhang, Haijiang}
}
@conference {yin2017relating,
	title = {Relating backprojection images to kinematics and dynamic source models},
	booktitle = {AGU Fall Meeting Abstracts},
	year = {2017},
	author = {Yin, Jiuxun and Denolle, M}
}
@conference {yin20164d,
	title = {4D source analysis of the 2015 Illapel, Chile Mw8. 3 Earthquake},
	booktitle = {AGU Fall Meeting Abstracts},
	year = {2016},
	author = {Yin, J and Denolle, M and Yao, H. and Yang, H}
}
@article {yin2016coseismic,
	title = {Coseismic radiation and stress drop during the 2015 Mw 8.3 Illapel, Chile megathrust earthquake},
	journal = {Geophysical Research Letters},
	volume = {43},
	number = {4},
	year = {2016},
	pages = {1520{\textendash}1528},
	url = {https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2015GL067381},
	author = {Yin, Jiuxun and Yang, Hongfeng and Yao, Huajian and Weng, Huihui}
}
@article {yin2016rupture,
	title = {Rupture and frequency-dependent seismic radiation of the 2012 M w 8.6 Sumatra strike-slip earthquake},
	journal = {Geophysical Journal International},
	volume = {205},
	number = {3},
	year = {2016},
	pages = {1682{\textendash}1693},
	publisher = {Oxford University Press},
	url = {https://academic.oup.com/gji/article-abstract/205/3/1682/656054},
	author = {Yin, Jiuxun and Yao, Huajian}
}
@article {yang2015tectonic,
	title = {Tectonic settings of the 2015 M w 8.3 Coquimbo, Chile earthquake and its implications on megathrust earthquakes},
	journal = {Chinese Science Bulletin},
	volume = {60},
	number = {36},
	year = {2015},
	pages = {3549{\textendash}3556},
	publisher = {Science China Press},
	abstract = {The 2015\&nbsp;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\&nbsp;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\&nbsp;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\&nbsp;Mw8.3 earthquake provide an excellent opportunity to validate numerical simulations of rupture process based on interseismic locking distributions.},
	author = {Yang, Hongfeng and LIN, Jian and Yin, Jiuxun and Yao, Huajian}
}
@conference {qin2013rupture,
	title = {Rupture process of the 2013 Okhotsk deep mega earthquake from iterative backprojection and compress sensing methods},
	booktitle = {AGU Fall Meeting Abstracts},
	year = {2013},
	author = {Qin, W and Yin, J and Yao, H.}
}
