Abstract:With the rapid development of Global Navigation Satellite System (GNSS) technology, the long-term accumulated GNSS observations of global reference stations have provided valuable data for geodesy and geodynamics studies since the 1990s. Acquiring the precise velocity of GNSS stations is very important for the study of global plate movement, crustal deformation, etc. However, the seismic activities nearby some GNSS observation stations may seriously change the station’s motion trajectory. Therefore, our research was motivated to propose a method allowing for station seismic deformation, and apply it to construct an updated global GNSS velocity field. The main contributions of this work included the following. Firstly, we improved the GNSS data processing procedures and seismic data selection strategies to obtain GNSS coordinate time series with mm-level precision (3–5 and 6–8 mm in the horizontal and vertical, respectively) and information of each site impacted by seismic events, which provides necessary input data for further analysis. Secondly, an Integrated Time Series Method (ITSM) concerning the effect of seismic deformation was proposed to model the station’s nonlinear motion accurately. Distinguished with existing studies, all parameters including seismic relaxation time can be simultaneously estimated by ITSM, which improves the accuracy and reliability of GNSS station velocity significantly. Thirdly, to optimize the ITSM-based model, the influences of seismic relaxation time (a. 0.1 × true, b. 10 × true, c. true), parameterization mode (a. Offset + Velocity, b. Offset + Velocity + PSD, c. Offset + Velocity + PSD + Period), and the Post-Seismic Deformation (PSD) model (a. None, b. Exp, c. Log, d. Exp + Log) on results of GNSS time series analyzing were discussed. The results showed that the fitting accuracy of GNSS displacements was better than 5 mm and 10 mm in the horizontal and vertical, respectively. Finally, the global GNSS station velocity field (referred to as GGV2020 hereafter) was refined by ITSM using global GNSS observations and seismic data during 1990–2020. This not only helps interpret plate tectonic motion, establish and maintain a Dynamic Terrestrial Reference Frame (DTRF) but also contributes to better investigating geodynamic processes. GGV2020 results showed that the accuracy of global velocity was better than 1 mm/a, and the averages of Root Mean Square Error (RMSE) were 0.19 mm/a, 0.19 mm/a, and 0.33 mm/a in the north, east, and up direction, respectively. Besides, the RMSE obeys normal distribution. Compared with ITRF2014, there was a difference of about 1–2 mm/a between them due to differences in terms of observation span, processing model, and geodetic technology. Moreover, GGV2020 is expected to enrich and update the existing velocity field products to describe the characteristics of regional crustal movement in more detail, especially in Antarctica.

Global Navigation Satellite System (GNSS) time series and velocities about a slowly convergent margin processed on high-performance computing (HPC) clusters: products and robustness evaluation

Surface Velocities and Strain-Rates in the Euro-Mediterranean Region From Massive GPS Data Processing

Comparative Analysis of Methods to Estimate Geodetic Strain Rates from GNSS Data in Italy

Analysis of Seismic Deformation from Global Three-Decade GNSS Displacements: Implications for a Three-Dimensional Earth GNSS Velocity Field.

Fusion of GNSS and Satellite Radar Interferometry: Determination of 3D Fine-Scale Map of Present-Day Surface Displacements in Italy as Expressions of Geodynamic Processes

The OGS–Northeastern Italy Seismic and Deformation Network: Current Status and Outlook

On the seafloor horizontal displacement from cGPS and compass data in the Campi Flegrei caldera

Monitoring of Land Subsidence in the Po River Delta (Northern Italy) Using Geodetic Networks

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

A New Tightly Coupled Method for High-Rate Seismogeodesy: a Shake Table Experiment and Application to the 2016 Mw 6.6 Central Italy Earthquake

Low-Cost GNSS Solution for Continuous Monitoring of Slope Instabilities Applied to Madonna Del Sasso Sanctuary (NW Italy)

An Integrated InSAR and GNSS Approach to Monitor Land Subsidence in the Po River Delta (Italy)

Towards the Application of Seismogeodesy in Central Italy: a Case Study for the 2016 August 24 Mw 6.1 Italy Earthquake Modelling

Near-Real-Time Strong Motion Acquisition at National Scale and Automatic Analysis

High resolution ambient noise tomography of the Southwestern Alps and the Ligurian margin

Empirical evidence for multi-decadal transients affecting geodetic velocity fields and derived seismicity forecasts in Italy

Two-Decade GNSS Observation Processing and Analysis with the New IGS Repro3 Criteria: Implications for the Refinement of Velocity Field and Deformation Field in Continental China

Fault Kinematic Modeling Along a Widely Deformed Plate Boundary in Southern Italy

Wavelet-like denoising of GNSS data through machine learning. Application to the time series of the Campi Flegrei volcanic area (Southern Italy)

Accurate GPS-based timestamp facility for Gran Sasso National Laboratory

Constraints on the Structure of the Shallow Crust in Central Italy from Geophysical Log Data