Historical seismic catalogs report that the Gargano Promontory (southern Italy) was affected in the past by earthquakes with medium to high estimated magnitude. From the instrumental seismicity, it can be identified that the most energetic Apulian sequence occurred in 1995 with a main shock of MW = 5.2 followed by about 200 aftershocks with a maximum magnitude of 3.7. The most energetic earthquakes of the past are attributed to right-lateral strike-slip faults, while there is evidence that the present-day seismicity occur on thrust or thrust-strike faults. In this article, we show a detailed study on focal mechanisms and stress field obtained by micro-seismicity recorded from April 2013 until the present time in the Gargano Promontory and surrounding regions. Seismic waveforms are collected from the OTRIONS Seismic Network (OSN), from the Italian National Seismic Network (RSN), and integrated with data from the Italian National Accelerometric Network (RAN) in order to provide a robust dataset of earthquake localizations and focal mechanisms. The effect of uncertainties of the velocity model on fault plane solutions (FPS) has been also evaluated indicating the robustness of the results. The computed stress field indicates a deep compressive faulting with maximum horizontal compressive stress, SHmax, trending NW-SE. The seismicity pattern analysis indicates that the whole crust is seismically involved up to a depth of 40 km and indicates the presence of a low-angle seismogenic surface trending SW-NE and dipping SE-NW, similar to the Gargano–Dubrovnik lineament. Shallower events, along the eastern sector of the Mattinata Fault (MF), are W-E dextral strike-slip fault. Therefore, we hypothesized that the seismicity is locally facilitated by preexisting multidirectional fractures, confirmed by the heterogeneity of focal mechanisms, and explained by the different reactivation processes in opposite directions over the time, involving the Mattinata shear zone.

The University of Bari (Italy), in cooperation with the National Institute of Geophysics and Volcanology (INGV) (Italy), has installed the OTRIONS micro-earthquake network to better understand the active tectonics of the Gargano promontory (Southern Italy). The OTRIONS network operates since 2013 and consists of 12 short period, 3 components, seismic stations located in the Apulian territory (Southern Italy). This data article releases the waveform database collected from 2013 to 2018 and describes the characteristics of the local network in the current configuration. At the end of 2018, we implemented a cloud infrastructure to make more robust the acquisition and storage system of the network through a collaboration with the RECAS-Bari computing centre of the University of Bari (Italy) and of the National Institute of Nuclear Physics (Italy). Thanks to this implementation, waveforms recorded after the beginning of 2019 and the station metadata are accessible through the European Integrated Data Archive (EIDA, https://www.orfeus-eu.org/data/eida/nodes/INGV/).

We investigated whether the joint inversion of geodetic and stress direction data can constrain long‐term fault slip rates in the central Apennines, and ultimately how extension is partitioned among fault slip and bulk lithosphere permanent strain. Geodetic velocities are collected in the fault interseismic stage with steady secular deformation; thus, long‐term estimates can be derived with a model of elastically unloading seismogenic faults within a viscously deforming lithosphere. As the average spacing of permanent Global Navigation Satellite Systems (GNSS) stations is similar to the average length of seismogenic faults (25–35 km), if not larger, we decided to merge permanent and temporary GNSS measurements, resulting in a denser geodetic data set. Given that most normal faults in the Apennines have slip rates around or below 1 mm/a, and most campaign GNSS velocities carry similar uncertainties, simple local back slip models cannot be applied. More sophisticated modeling is required to extract reasonable bulk deformation rates and long‐term fault slip rates at signal‐to‐noise ratio of order unity. Given the spatial distribution of the GNSS network, we estimated the long‐term slip rate of seven major fault systems that are in satisfactory agreement with available geological slip rates. The resulting spatial distribution of bulk deformation rates locally fits short‐term transients; in other cases, they represent the currently unclear signature of tectonic processes like upper‐crustal viscoplastic deformation and aseismic slip, or indicate missing faults in the adopted database. We conclude that the time is ripe for determining fault slip rates using geodetic and stress direction data, particularly where fault activity rates are hard to determine geologically.

No abstract