Roberto Devoti

Static stress drop as determined from geodetic strain rates and statistical seismicity

Accepted for publication in Journal of Geophysical Researches. Copyright (2010) American Geophysical Union

The coseismic and postseismic deformation of the L’Aquila, 2009 earthquake from repeated GPS measurements

We analyze more than 100 GPS time series of continuous and discontinuous GPS stations located in the Abruzzi region (Italy) surrounding the epicentres of the L’Aquila 2009 seismic sequence. The purpose of this work is to reconstruct the coseismic displacement field caused by the 6 th April (Mw 6.3) main shock from a dense network of survey-mode stations surrounding the epicentral area and to characterize the early postseismic deformation field. In the months following the main shock, an extensive GPS survey was carried out on the existing Central Apennines Geodetic Network (CAGeoNet), with the intention of collecting a robust data set and to study the co- and postseismic deformation field of this Apenninic normal faulting earthquake. The analysis is carried out with two independent procedures and software (Bernese and Gamit) in order to provide reliable and validated geodetic solutions. The analysis of the postseismic transients and the knowledge of long-term inter-seismic velocities at all GPS stations, issued from permanent and CAGeoNet sites, allow us to derive a dense co- and postseismic displacement field for the L’Aquila Mw 6.3 main shock in a wide area around the epicentre. The highest deformation rate occurs during the first 4–5 months after the main shock and persists in the following at slightly slower rate throughout the whole monitoring period. Fast deformation rates imply that most of the observed deformation is due to a process different from a pure viscoelastic relaxation of the stress perturbation. Since the observed rates would imply a too low effective viscosity value (below 10 17 Pa s), we rather suggest that most of the observed deformation in the first months after the earthquake is due to different processes, most likely frictional afterslip possibly modulated by the presence of fluids. The new coseismic displacement field is used to invert for the main shock fault geometry, analysing the consistency among the different geodetic solutions and the combined one, with the goal of validating the two data sets.

Constraining primary surface rupture length along the Paganica fault (2009 L’Aquila earthquake) with geological and geodetic (DInSAR and GPS) data

Part of this work has been carried out within the ASI-SIGRIS project, funded by the Italian Space Agency and Istituto Nazionale di Geofisica e Vulcanologia

The velocity field of the Italian area

AbstractThe rapid development of several permanent GNSS networks in Italy has made available a huge amount of GNSS observations, giving the chance to figure out and significantly improve the spatial and temporal resolutions of the crustal deformation in the Italian area. More than 20 GNSS networks, promoted and managed by different institutions, constitute the grid of monitoring stations that includes over 1000 permanent stations, mainly devoted to real-time positioning services but that has proven to be suitable for monitoring slow deforming processes, such as for instance, intraplate deformation processes. The whole set of raw GPS data is routinely processed at INGV providing daily solutions of station coordinates and estimating linear velocities for each station. The information content of coordinate time series is wide, the station position variations incorporate linear and non-linear effects caused by geophysical phenomena of different nature, of which we show some evidences. The sectors where the permanent network is augmented with non-permanent sites allows to study tectonic processes with a finer resolution.

VADASE Reliability and Accuracy of Real-Time Displacement Estimation: Application to the Central Italy 2016 Earthquakes

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.

Satellite positioning and geophysics studies in Italy

A brief historical overview of the Italian geophysical studies using satellite positioning observations.

Sannio-Matese Mounts (Southern Italy) deformation field from GPS Data (2002-2014)

Lucchi, Renata G. ... et. al.-- 87° Congresso della Societa Geologica Italiana e 90° Congresso della Societa Italiana di Mineralogia e Petrologia, The Future of the Italian Geosciences - The Italian Geosciences of the Future, 10-12 September 2014, Milan, Italy.-- 1 pageThe Montellina Spring (370 m a.s.l.) represents an example of groundwater resource in mountain region. It is a significant source of drinking water located in the right side of the Dora Baltea Valley (Northwestern Italy), SW of Quincinetto town. This spring shows a morphological location along a ridge, 400 m from the Renanchio Torrent in the lower sector of the slope. The spring was investigated using various methodologies as geological survey, supported by photo interpretation, structural reconstruction, NaCl and fluorescent tracer tests, discharge measurements. This multidisciplinary approach, necessary due to the complex geological setting, is required for the importance of the Montellina Spring. It is interesting in the hydrogeological context of Western Alps for its high discharge, relatively constant over time (average 150 l/s), and for its location outside a fluvial incision and suspended about 40 m above the Dora Baltea valley floor (Lasagna et al. 2013). According to the geological setting, the hydrogeological reconstruction of the area suggests that the large amount of groundwater in the basin is essentially favoured by a highly fractured bedrock, covered by wide and thick bodies of glacial and gravitational sediments. The emergence of the water along the slope, in the Montellina Spring, is essentially due to a change of permeability between the deep bedrock and the shallow bedrock and/or surficial sediments. The deep bedrock, showing closed fractures and/or fractures filled by glacial deposits, is slightly permeable. The shallow bedrock, strongly loosened as result of gravitational phenomena, and the local gravitational sediments are, on the contrary, highly permeable. The concentration of water at the spring is due to several reasons. a) The spring is immediately downward a detachment niche, dipping towards the spring, that essentially drains the water connected to the change of permeability in the bedrock. b) It is along an important fracture, that carries a part of the losses of the Renanchio Torrent. c) Finally, it is favored by the visible and buried morphology. Although it is located along a ridge, the spring occurs in a small depression between a moraine and a landslide body. It also can be favored by the likely concave trend of buried base of the landslide. At last, tracer tests of the Renanchio Torrent water with fluorescent tracer are performed, with a continuous monitoring in the Montellina Spring. The surveys permit to verify and quantify the spring and torrent hydrogeological relationship, suggesting that only a small fraction of stream losses feeds the spring.

Nearly Diurnal Structure of Polar Motion Analysing Lageos SLR Data

Since more accurate and precise measurements are available in space geodetic techniques, there is a growing interest for the scientific community in studying the short periodic fluctuation of the Earth’s rotation axis. In particular subtle nearly diurnal variations of the orientation of the rotation axis, as viewed from the Earth, reflect dynamical properties of the deep interior of the Earth. Recently Very Long Baseline Interferometry (VLBI) data analysis sets a stringent upper limit to the nearly diurnal wobble whose amplitude is well below the milliarc second (mas) level [Herring et al., 1991].

Tectonic Motion from LAGEOS SLR Data Using a Network Adjustment Technique

In this paper we present the most recent results from the analysis of the LAGEOS SLR data concerning the present plate tectonic motion. The horizontal site velocities are estimated via a network adjustment by a weighted least squares technique; the geodesic rates between pairs of stations, estimated by a linear fitting of the semiannual values are modeled by an observation equation, linear in latitude and longitude time derivatives, assumed as solve-for-parameters. This observation equation is based on the mapping of the reference ellipsoid onto the so called ‘image sphere’, properly transforming geodesics on ellipsoid to great circle arcs on the sphere. The method is applied to the results of the TPZ90 solution in the case of a worldwide network.