No abstract

Analyzing the displacement time series from continuous GPS (cGPS) with an Independent Component Analysis, we detect a transient deformation signal that correlates both in space and time with a seismic swarm activity (maximum Mw=3.69 ± 0.09) occurred in the hanging wall of the Altotiberina normal fault (Northern Apennines, Italy) in 2013–2014. The geodetic transient lasted ∼6 months and produced a NW‐SE trending extension of ∼5.3 mm, consistent with the regional tectonic regime. The seismicity and the geodetic signal are consistent with slip on two splay faults in the Altotiberina fault (ATF) hanging wall. Comparing the seismic moment associated with the geodetic transient and the seismic events, we observe that seismicity accounts for only a fraction of the measured geodetic deformation. The combined seismic and aseismic slip decreased the Coulomb stress on the locked shallow portion of the ATF, while the transition region to the creeping section has been loaded.

While low‐angle normal faults have been recognized worldwide from geological studies, whether these structures are active or capable of generating big earthquakes is still debated. We provide new constraints on the role and modes of the Altotiberina fault (ATF) in accommodating extension in the Northern Apennines. We model GPS velocities to study block kinematics, faults slip rates and interseismic coupling of the ATF, which is active and accounts, with its antithetic fault, for a large part of the observed chain normal 3 mm/yr tectonic extension. A wide portion of the ATF creeps at the long‐term slip rate (1.7 ± 0.3 mm/yr), but the shallow locked portions are compatible with M > 6.5 earthquakes. We suggest that positive stress accumulation due to ATF creep is most likely released by more favorable oriented splay faults, whose rupture may propagate downdip along low‐angle normal fault surface and reduce the probability of occurrence of a seismic rupture of the shallower locked portion.

We apply a blind source separation algorithm to the ground displacement time series recorded at continuous Global Positioning System (GPS) stations in the European Eastern Alps and Northern Dinarides. As a result, we characterize the temporal and spatial features of several deformation signals. Seasonal displacements are well described by loading effects caused by Earth surface mass redistributions. More interestingly, we highlight a horizontal, nonseasonal, transient deformation signal, with spatially variable amplitudes and directions. The stations affected by this signal reverse the sense of movement with time, implying a sequence of dilatational and compressional deformation that is oriented normal to rock fractures in karst areas. The temporal evolution of this deformation signal is correlated with the history of cumulated precipitations at monthly time scales. This transient horizontal deformation can be explained by pressure changes associated with variable water levels within vertical fractures in the vadose zones of karst systems. The water level changes required to open or close these fractures are consistent with the fluctuations of precipitation and with the dynamics of karst systems.

Abstract. This study presents and discusses horizontal and vertical geodetic velocities for a low strain rate region of the south Alpine thrust front in northeastern Italy obtained by integrating GPS, interferometric synthetic aperture radar (InSAR) and leveling data. The area is characterized by the presence of subparallel, south-verging thrusts whose seismogenic potential is still poorly known. Horizontal GPS velocities show that this sector of the eastern Southern Alps is undergoing ∼1 mm a−1 of NW–SE shortening associated with the Adria–Eurasia plate convergence, but the horizontal GPS velocity gradient across the mountain front provides limited constraints on the geometry and slip rate of the several subparallel thrusts. In terms of vertical velocities, the three geodetic methods provide consistent results showing a positive velocity gradient, of ∼ 1.5 mm a−1, across the mountain front, which can hardly be explained solely by isostatic processes. We developed an interseismic dislocation model whose geometry is constrained by available subsurface geological reconstructions and instrumental seismicity. While a fraction of the measured uplift can be attributed to glacial and erosional isostatic processes, our results suggest that interseismic strain accumulation at the Montello and the Bassano–Valdobbiadene thrusts it significantly contributing to the measured uplift. The seismogenic potential of the Montello thrust turns out to be smaller than that of the Bassano–Valdobbiadene fault, whose estimated parameters (locking depth equals 9.1 km and slip rate equals 2.1 mm a−1) indicate a structure capable of potentially generating a Mw>6.5 earthquake. These results demonstrate the importance of precise vertical ground velocity data for modeling interseismic strain accumulation in slowly deforming regions where seismological and geomorphological evidence of active tectonics is often scarce or not conclusive.

In this work we present and discuss new geodetic velocity and strain-rate fields for the Euro-Mediterranean region obtained from the analysis of continuous GNSS stations. We describe the procedures and methods adopted to analyze raw GPS observations from >4000 stations operating in the Euro-Mediterranean, Eurasian and African regions. The goal of this massive analysis is the monitoring of Earth’s crust deformation in response to tectonic processes, including plate- and micro-plate kinematics, geodynamics, active tectonics, earthquake-cycle, but also the study of a wide range of geophysical processes, natural and anthropogenic subsidence, sea-level changes, and hydrology. We describe the computational infrastructure, the methods and procedures adopted to obtain a three-dimensional GPS velocity field, which is used to obtain spatial velocity gradients and horizontal strain-rates. We then focus on the Euro-Mediterranean region, where we discuss the horizontal and vertical velocities, and spatial velocity gradients, obtained from stations that have time-series lengths longer than 6 and 7 years, which are found to be the minimum spans to provide stable and reliable velocity estimates in the horizontal and vertical components, respectively. We compute the horizontal strain-rate field and discuss deformation patterns and kinematics along the major seismogenic belts of the Nubia-Eurasia plate boundary zone in the Mediterranean region. The distribution and density of continuous GNSS stations in our geodetic solution allow us to estimate the strain-rate field at a spatial scale of ∼27 km over a large part of southern Europe, with the exclusion of the Dinaric mountains and Balkans.

M w = 6.3 L'Aquila earthquake occurred within a complex system of NW-SE trending normal faults in the Abruzzi Central Apennines (Italy). We analyse the coseismic deformation as measured by >70 global positioning system (GPS) stations, both from continuous and survey-mode networks, providing unprecedented details for a moderate normal faulting earthquake in Italy from GPS measurements. We use rectangular, uniform-slip, dislocations embedded in an elastic, homogeneous and isotropic half-space and a constrained, non-linear optimization algorithm, to solve for the best-fitting rectangular dislocation geometry and coseismic-slip distribution. We use a bootstrap approach to investigate uncertainties in the model parameters and define confidence bounds for all the inverted parameters. The rupture occurred on a N129 +8 −7• E striking and 50 +9 −5• southwestward dipping normal fault, in agreement with geological observations of surface breaks along the Paganica fault. Our distributed slip model exhibits a zone of relatively higher slip (>60 cm) between ∼1.5 and ∼11 km depth, along a roughly downdip, NW-SE elongated patch, confined within the fault plane inverted assuming uniform-slip. The highest slip, of the order of ∼1 m, occurred on a ∼16 km 2 area located at ∼5 km depth, SE of the mainshock epicentre. The analysis of model resolution suggests that slip at depth below ∼5 km can be resolved only at a spatial scale larger than 2 km, so a finer discretization of different asperities within the main patch of coseismic-slip is not allowed by GPS data. We compute the coseismic Coulomb stress changes in the crustal volume affected by the major aftershocks, and compare the results obtained from the uniformslip and the heterogeneous-slip models. We find that most of the large aftershocks occurred in areas of Coulomb stress increase of 0.2-13 bar and that a deepening of the slip distribution down to a depth greater than 6 km in the SE part of the fault plane, in agreement with the inverted slip model, can explain the deepest, April 7, M w 5.3 aftershock.

The 2012 Emilia Romagna (Italy) seismic sequence has been extensively studied given the occurrence of two mainshocks, both temporally and spatially close to each other. The recent literature accounts for several fault models, obtained with different inversion methods and different datasets. Several authors investigated the possibility that the second event was triggered by the first mainshock with elusive results. In this work, we consider all the available InSAR and GPS datasets and two planar fault geometries, which are based on both seismological and geological constraints. We account for a layered, elastic half-space hosting the dislocation and compare the slip distribution resulting from the inversion and the related changes in Coulomb Failure Function (CFF) obtained with both a homogeneous and layered half-space. Finally, we focus on the interaction between the two main events, discriminating the contributions of coseismic and early postseismic slip of the mainshock on the generation of the second event and discuss the spatio-temporal distribution of the seismic sequence. When accounting for both InSAR and GPS geodetic data we are able to reproduce a detailed coseismic slip distribution for the two mainshocks that is in accordance with the overall aftershock seismicity distribution. Furthermore, we see that an elastic medium with depth dependent rigidity better accounts for the lack of the shallow seismicity, amplifying, with respect to the homogeneous case, the mechanical interaction of the two mainshocks.