This repository contains supplementary code for Brunsvik et al. (2021). We construct 3-D models of faults based on the distribution of hypocenters using clustering. We applied this method to the 2009 L'Aquila (Italy) earthquake sequence. If you use this code or the fault models, please cite us:

Brunsvik, B., Morra, G., Cambiotti, G., Chiaraluce, L., Di Stefano, R., De Gori, P., & Yuen, D. A. (2021). Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy). Tectonophysics, 804, 228756.

The Jupyter Notebook Supplementary_Notebook.ipynb combines both the writing and high level code in order to improve reproducibility of the work. You can view the notebook in nbviewer with interactive figures (https://bit.ly/2MnCFdj). You can also run the Notebook in the cloud using Binder by following this Binder link:

The two datasets used here are available publicly. We used Valoroso et al. (2020) for earthquake hypocenter information and Chiaraluce et al. (2021) for focal mechanisms.

Lauro Chiaraluce, Pasquale De Gori, Raffaele Di Stefano, & Luisa Valoroso. (2021). L'Aquila 2009 seismic sequence: integrated dataset of automatic first motion polarities focal mechanisms and RMT with HypoDD high quality relative earthquake locations (Version 1.0) [Data set]. Zenodo. http://doi.org/10.5281/zenodo.4024065

Luisa Valoroso, Lauro Chiaraluce, Davide Piccinini, Raffaele Di Stefano, David Schaff, & Felix Waldhauser. (2020). The 2009 Mw6.1 L'Aquila normal fault system imaged by 64,051 high-precision foreshock and aftershock locations. (Version Version 1.0) [Data set]. Zenodo. http://doi.org/10.5281/zenodo.4036248