Anderlini, Letizia

Geodetic data can detect and estimate deformation signals and rates due to natural and anthropogenic phenomena. In the present study, we focus on northeastern Italy, an area characterized by ~1.5–3 mm/yr of convergence rates due to the collision of Adria-Eurasia plates and active subsidence along the coasts. To define the rates and trends of tectonic and subsidence signals, we use a Multi-Temporal InSAR (MT-InSAR) approach called the Stanford Method for Persistent Scatterers (StaMPS), which is based on the detection of coherent and temporally stable pixels in a stack of single-master differential interferograms. We use Sentinel-1 SAR images along ascending and descending orbits spanning the 2015–2019 temporal interval as inputs for Persistent Scatterers InSAR (PSI) processing. We apply spatial-temporal filters and post-processing steps to reduce unrealistic results. Finally, we calibrate InSAR measurements using GNSS velocities derived from permanent stations available in the study area. Our results consist of mean ground velocity maps showing the displacement rates along the radar Line-Of-Sight for each satellite track, from which we estimate the east–west and vertical velocity components. Our results provide a detailed and original view of active vertical and horizontal displacement rates over the whole region, allowing the detection of spatial velocity gradients, which are particularly relevant to a better understanding of the seismogenic potential of the area. As regards the subsidence along the coasts, our measurements confirm the correlation between subsidence and the geological setting of the study area, with rates of ~2–4 mm/yr between the Venezia and Marano lagoons, and lower than 1 mm/yr near Grado.

The geodetic dataset used in the research article entitled "Multi-technique geodetic detection of onshore and offshore subsidence along the Upper Adriatic Sea coasts"[1] is presented here. It consists of the outcomes of three different techniques, i.e. Synthetic Aperture Radar Interferometry (InSAR), Global Navigation Satellite System (GNSS) and topographic Levelling surveys. This dataset has been used for the estimation of onshore and offshore deformation in a mineral concession area located along the Upper Adriatic Sea coastal area (Italy), South-East of Ravenna city. InSAR data covers the period from 2002 to 2018, GNSS data from 1998 to 2018 and levelling data from 2002 to 2017.The different measurements have been cross-validated and referred to a common local reference system fixed in the urban area of Ravenna. This data collection will be very useful for deepening the analysis of any type of deformation in the Ravenna coastal area.

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 threedimensional 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.

We assess about 20 years of onshore and offshore subsidence along a sector of the Upper Adriatic Sea (Italy) coastal areas affected by natural soil compaction and intense anthropogenic activities such as aquifers exploitation and hydrocarbons extraction. Our approach is based on the synergistic use of independent remote sensing and in-situ geodetic data to detect and spatially characterise the deformation pattern by cross-validating the different available measurements. We collect extensive datasets from i) SAR images provided by Envisat, Cosmo- SkyMed and Sentinel-1 missions, ii) GNSS measurements from continuous stations managed by public institutions, local authorities and private companies and iii) Leveling surveys. The cross-validation analysis shows good agreement among all the independent datasets, thus providing a reliable assessment of the ongoing deformation. We detect an onshore and offshore subsidence peak of about 1/-1.5 cm/yr in the proximity of the coastline, close to Lido di Dante and Fiumi Uniti villages, and at the present offshore platform. The outcomes highlight how the integration of different remote sensing and in situ geodetic techniques is successful to retrieve deformation history in time and space in complex areas, where different natural and anthropogenic sources concur to the overall deformation pattern. Moreover, such approach provides a robust support to modelling studies for hazard assessment in both inland and shoreline areas.

In this paper we describe the results of an experimental implementation of the recent guidelines issued by the Italian regulatory body for monitoring hydrocarbon production activities. In particular, we report about the pilot study on seismic, deformation, and pore pressure monitoring of the Mirandola hydrocarbon cultivation facility in Northern Italy. This site hosts the Cavone oil field that was speculated of possibly influencing the 2012 ML 5.8 Mirandola earthquake source. According to the guidelines, the monitoring center should analyse geophysical measurements related to seismicity, crustal deformation and pore pressure in quasi real-time (within 24–48 h). A traffic light system would then be used to regulate underground operations in case of detecting significant earthquakes (i.e., events with size and location included in critical ranges). For these 2-year period of guidelines experimentation, we analysed all different kinds of available data, and we tested the existence of possible relationship between their temporal trends. Despite the short time window and the scarce quantity of data collected, we performed the required analysis and extracted as much meaningful and statistically reliable information from the data. We discuss here the most important observations drawn from the monitoring results, and highlight the lessons learned by describing practical issues and limitations that we have encountered in carrying out the tasks as defined in the guidelines. Our main goal is to contribute to the discussion about how to better monitor the geophysical impact of this kind of anthropogenic activity. We point out the importance of a wider seismic network but, mostly, of borehole sensors to improve microseismic detection capabilities. Moreover, the lack of an assessment of background seismicity in an unperturbed situation -due to long life extraction activities- makes it difficult to get a proper picture of natural background seismic activity, which would be instead an essential reference information for a tectonically-active regions, such as Northern Italy.

The 2016–2017 Central Italy earthquake sequence struck the central Apennines between August 2016 and October 2016 with Mw ∈ [5.9; 6.5], plus four earthquakes occurring in January 2017 with Mw ∈ [5.0; 5.5]. We study Global Positioning System time series including near- and far-field domains. We use a variational Bayesian independent component analysis technique to separate the post-seismic deformation from signals caused by variation of the water content in aquifers at hundreds of meters of depth and of the soil moisture. For each independent component, realistic uncertainties and a plausible physical explanation are provided. We focus on the study of afterslip on the main structures surrounding the mainshock, highlighting the role played by faults that were not activated during the co-seismic phase in accommodating the post-seismic deformation. We report aseismic deformation occurring on the Paganica fault, which hosted the Mw 6.1 2009 L'Aquila earthquake, suggesting that static stress transfer and aseismic slip influence the recurrence time of nearby (∼50 km further south of the mainshocks) segments. A ∼2–3 km thick subhorizontal shear-zone, clearly illuminated by seismicity, which bounds at depth the west-dipping normal faults where the mainshocks nucleated, also shows aseismic slip. Since afterslip alone underestimates the displacement in the far-field domain, we consider the possibility that the shear zone marks the brittle-ductile transition, assuming the viscoelastic relaxation of the lower crust as a mechanism contributing to the post-seismic displacement. Our results suggest that multiple deformation processes are active in the first 2 years after the mainshocks.

We provide a dataset of 3D coordinate time series of 37 continuous GNSS stations installed for stability monitoring purposes on onshore and offshore industrial settlements along a NW-SE-oriented and ~100-km-wide belt encompassing the eastern Italian coast and the Adriatic Sea. The dataset results from the analysis performed by using different geodetic software (Bernese, GAMIT/GLOBK and GIPSY) and consists of six raw position time series solutions, referred to IGb08 and IGS14 reference frames. Time series analyses and comparisons evidence that the different solutions are consistent between them, despite the use of different software, models, strategy processing and frame realizations. We observe that the offshore stations are subject to significant seasonal oscillations probably due to seasonal environmental loads, seasonal temperature-induced platform deformation and hydrostatic pressure variations. Many stations are characterized by non-linear time series, suggesting a complex interplay between regional (long-term tectonic stress) and local sources of deformation (e.g. reservoirs depletion, sediment compaction). Computed raw time series, logs files, phasor diagrams and time series comparison plots are distributed via PANGAEA ( https://www.pangaea.de ).

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.

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.