Spatiotemporal filtering of regional GNSS network’s position time series with missing data using principle component analysis

The existing spatiotemporal analysis methods suppose that the involved time series are complete and have the same data interval. However missing data inevitably occur in the position time series of Global Navigation Satellite Systems networks for many reasons. In this paper, we develop a modified principal component analysis to extract the Common Mode Error (CME) from the incomplete position time series. The principle of the proposed method is that a time series can be reproduced from its principle components. The method is equivalent to the method of Dong et al. (J Geophys Res 111:3405–3421, 2006) in case of no missing data in the time series and to the extended ‘stacking’ approach under the assumption of a uniformly spatial response. The new method is first applied to extract the CME from the position time series of the Crustal Movement Observation Network of China (CMONOC) over the period of 1999–2009 where the missing data occur in all stations with the different gaps. The results show that the CMEs are significant in CMONOC. The size of the first principle components for the North, East and Up coordinates are as large as 40, 41 and 37 % of total principle components and their spatial responses are not uniform. The minimum amplitudes of the first eigenvectors are only 41, 15 and 29 % for the North, East and Up coordinate components, respectively. The extracted CMEs of our method are close to the data filling method, and the Root Mean Squared error (RMS) values computed from the differences of maximum CMEs between two methods are only 0.31, 0.52 and 1.55 mm for North, East and Up coordinates, respectively. The RMS of the position time series is greatly reduced after filtering out the CMEs. The accuracies of the reconstructed missing data using the two methods are also comparable. To further comprehensively test the efficiency of our method, the repeated experiments are then carried out by randomly deleting different percentages of data at some stations. The results show that the CMEs can be extracted with high accuracy at the non missing-data epochs. And at the missing-data epochs, the accuracy of extracted CMEs has a strong dependence on the number of stations with missing data.

[1]  Yehuda Bock,et al.  Southern California permanent GPS geodetic array: Spatial filtering of daily positions for estimating coseismic and postseismic displacements induced by the 1992 Landers earthquake , 1997 .

[2]  John Langbein,et al.  Noise in two‐color electronic distance meter measurements revisited , 2004 .

[3]  M. Ouassou,et al.  A continuous velocity field for Norway , 2013, Journal of Geodesy.

[4]  David H. Schoellhamer,et al.  Singular spectrum analysis for time series with missing data , 2001 .

[5]  Yunzhong Shen,et al.  Spatiotemporal Signal and Noise Analysis of GPS Position Time Series of the Permanent Stations in China , 2014 .

[6]  John Langbein,et al.  Correlated errors in geodetic time series: Implications for time‐dependent deformation , 1997 .

[7]  Wangmin,et al.  Contemporary velocity field of crustal movement of Chinese mainland from Global Positioning System measurements , 2005 .

[8]  P. Teunissen,et al.  Assessment of noise in GPS coordinate time series : Methodology and results , 2007 .

[9]  Yehuda Bock,et al.  Error analysis of continuous GPS position time series , 2004 .

[10]  A. Amiri-Simkooei,et al.  Noise in multivariate GPS position time-series , 2009 .

[11]  Pascal Willis,et al.  Error Analysis of Weekly Station Coordinates in the DORIS Network , 2006 .

[12]  Yehuda Bock,et al.  Spatiotemporal filtering using principal component analysis and Karhunen-Loeve expansion approaches for regional GPS network analysis , 2006 .

[13]  S. Williams The effect of coloured noise on the uncertainties of rates estimated from geodetic time series , 2003 .

[14]  J. B. Rundle,et al.  Using Eigenpattern Analysis to Constrain Seasonal Signals in Southern California , 2004 .

[15]  R. Nikolaidis Observation of geodetic and seismic deformation with the Global Positioning System , 2002 .

[16]  C. Demets,et al.  Crustal velocity field of Mexico from continuous GPS measurements, 1993 to June 2001: Implications for the neotectonics of Mexico , 2003 .

[17]  Simon D. P. Williams,et al.  Fast error analysis of continuous GNSS observations with missing data , 2013, Journal of Geodesy.

[18]  Yehuda Bock,et al.  GPS measurements of current crustal movements along the Dead Sea Fault , 2004 .

[19]  Yehuda Bock,et al.  Southern California permanent GPS geodetic array: Error analysis of daily position estimates and site velocities , 1997 .

[20]  J. C. Savage Principal component analysis of interseismic deformation in southern California , 1995 .

[21]  K. Shedlock,et al.  Delineation of intrabasin structure in a dilational jog of the San Jacinto Fault Zone, southern California , 1995 .

[22]  Peter Steigenberger,et al.  Improved Constraints on Models of Glacial Isostatic Adjustment: A Review of the Contribution of Ground-Based Geodetic Observations , 2010 .

[23]  Bofeng Li,et al.  Efficient Estimation of Variance and Covariance Components: A Case Study for GPS Stochastic Model Evaluation , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[24]  T. Dixon,et al.  Noise in GPS coordinate time series , 1999 .

[25]  K. Tiampo,et al.  Analysis of GPS Measurements in Eastern Canada Using Principal Component Analysis , 2009, Pure and Applied Geophysics.