We investigate a large geodetic data set of interferometric synthetic aperture radar (InSAR) and GPS measurements to determine the source parameters for the three main shocks of the 2016 Central Italy earthquake sequence on 24 August and 26 and 30 October (M_w 6.1, 5.9, and 6.5, respectively). Our preferred model is consistent with the activation of four main coseismic asperities belonging to the SW dipping normal fault system associated with the Mount Gorzano‐Mount Vettore‐Mount Bove alignment. Additional slip, equivalent to a M_w ~ 6.1–6.2 earthquake, on a secondary (1) NE dipping antithetic fault and/or (2) on a WNW dipping low‐angle fault in the hanging wall of the main system is required to better reproduce the complex deformation pattern associated with the greatest seismic event (the M_w 6.5 earthquake). The recognition of ancillary faults involved in the sequence suggests a complex interaction in the activated crustal volume between the main normal faults and the secondary structures and a partitioning of strain release.

The Adria microplate is the foreland of the oppositely verging Apennines and Alps or Dinarides fold‐thrust belts associated to the related subduction zones. Along its western margin, the Adria plate hosts the active Northern Apennines accretionary prism, which is buried under the Adriatic Sea and the Po Plain. The interpretation of seismic reflection profiles and borehole data allowed us to define the geometry of the transition from the Apennines fold‐thrust belt to its undeformed foreland. Moreover, continuous GPS (CGPS) data from offshore hydrocarbon platforms anchored to the seabed of the northern Adriatic plate allow to measure present‐day kinematics. Although the CGPS signals are affected by non‐tectonic components associated with hydrocarbon extraction, the integration of geodetic analysis, subsurface geological reconstructions, and analytical modeling allowed us to constrain the ongoing tectonic activity. Shortening is currently accommodated by aseismic slip along the basal detachment, likely accumulating elastic energy along the frontal ramp that may eventually seismically slip. Our multidisciplinary study suggests that the study area may not be sheltered from relevant seismic sequences similar to the Mw 6 Emilia 2012 events and that the occurrence of potential seismogenic sources in the area should be carefully evaluated. Similar studies may be useful to constrain the present‐day activity in other marine areas and to identify potential and hitherto unrecognized seismogenic sources along the entire Apennines belt and other accretionary prisms worldwide.

Understanding natural variations of Rn (²²²Rn) concentrations is the fundamental prerequisite of using this radioactive gas as a tracer, or even precursor, of natural processes, including earthquakes. In this work, Rn concentrations in groundwater were continuously measured over a seven-month period, during 2017, in the Giardino Spring, Italy, together with groundwater levels in a nearby well installed into a fractured regional aquifer. Data were processed to reduce noise, and then analyzed to produce the Fourier spectra of Rn concentrations and groundwater levels. These spectra were compared with the spectrum of tidal forces. Results showed that diurnal and semidiurnal cycles of Rn concentrations, and filtered oscillations of groundwater levels, in the nearby well, are correlated with solar and luni-solar components of tidal forces, and suggested no correlation with the principal lunar components. Therefore, influencing factors linked to solar cycles, such as daily oscillations of temperature and atmospheric pressure, and related rock deformations, may have played a role in Rn concentrations and groundwater levels. An open question remains regarding the correlation, which is documented elsewhere, of Rn concentrations and groundwater levels with the lunar components of the solid Earth tides.

The goal of this article is the illustration of the new functionalities of the VADASE (Variometric Approach for Displacements Analysis Stand-alone Engine) processing approach. VADASE was presented in previous works as an approach able to estimate in real time the velocities and displacements in a global reference frame (ITRF), using high-rate (1 Hz or more) carrier phase observations and broadcast products (orbits, clocks) collected by a stand-alone GNSS receiver, achieving a displacements accuracy within 1–2 cm (usually better) over intervals up to a few minutes. It has been well known since the very first implementation and testing of VADASE that the estimated displacements might be impacted by two different effects: spurious spikes in the velocities due to outliers (consequently, displacements, obtained through velocities integration, are severely corrupted) and trends in the displacements time series, mainly due to broadcast orbit and clock errors. Two strategies are herein introduced, respectively based on Leave-One-Out cross-validation (VADASE-LOO) for a receiver autonomous outlier detection, and on a network augmentation strategy to filter common trends out (A-VADASE); they are combined (first, VADASE-LOO; second, A-VADASE) for a complete solution. Moreover, starting from this VADASE improved solution, an additional strategy is proposed to estimate in real time the overall coseismic displacement occurring at each GNSS receiver. New VADASE advances are successfully applied to the GPS data collected during the recent three strong earthquakes that occurred in Central Italy on 24 August and 26 and 30 October 2016, and the results are herein presented and discussed. The VADASE real-time estimated coseismic displacements are compared to the static ones derived from the daily solutions obtained within the standard post-processing procedure by the Istituto Nazionale di Geofisica e Vulcanologia.