Abstract
We propose a nonparametric minimum entropy method for estimating an optimal velocity from position time series, which may contain unknown noise, data gaps, loading effects, transients, outliers and step discontinuities. Although nonparametric, the proposed method is based on elementary statistical concepts familiar to least-squares and maximum-likelihood users. It seeks a constant velocity with a best possible (realistic) variance rather than a best variable velocity fit to the closest position data. We show, based on information theory, synthetic and real data, that minimum-entropy velocity estimation: (1) accounts for colored noise without assumptions about its distribution or the extent of its temporal correlations; (2) is unaffected by the series deterministic content such as an initial position and the heights of step discontinuities and insensitive to small-amplitude periodic variations and transients; (3) is robust against outliers and, for long time series, against step discontinuities and even slight non-stationarity of the noise; (4) does not involve covariance matrices or eigen/singular value analysis, thus can be implemented by a short and efficient software; (5) under no circumstances results in a velocity variance that decays as \(1/N\), where \(N\) is the number of observations. The proposed method is verified based on synthetic data and then applied to a few hundred NGL (Nevada Geodetic Lab) position time series of different characteristics, and the results are compared to those of the Median Interannual Difference Adjusted for Skewness (MIDAS) algorithm. The compared time series include continuous and linear ones used to test the agreement between the two methods in the presence of unknown noise, data gaps and loading effects, discontinuous but linear series selected to include the effect of a few (1–4) discontinuities, and nonlinear but continuous time series selected for including the effects of transients. Both the minimum-entropy and MIDAS methods are nonparametric in the sense that they only extract the velocity from a position time series with hardly any explicit assumptions about its noise distribution or correlation structure. Otherwise, the two methods differ in every single possible technical sense. Other than pointing to a close agreement between the derived velocities, the comparisons consistently revealed that minimum-entropy velocity uncertainties suggest a smaller degree of temporal correlations in the NGL time series than the MIDAS does.
Similar content being viewed by others
Data availability
We downloaded the NGL time series from http://geodesy.unr.edu/gps_timeseries/tenv/IGS14. Section S5 in the supplementary information lists the 4-character ID of all downloaded stations. The MIDAS software was downloaded from the NGL website http://geodesy.unr.edu/MIDAS_release and compiled on NGS computers. Our little software and its short documentation are included with the supplementary material.
References
Agnew D (1992) The time-domain behavior of power-law noises. Geophys Res Lett 19:333–336. https://doi.org/10.1029/91GL02832
Alevizakou EG, Siolas G, Pantazis G (2018) Short-term and long-term forecasting for the 3D point position changing by using artificial neural networks. ISPRS Int J Geo Inf 7(3):86. https://doi.org/10.3390/ijgi7030086
Amiri-Simkooei AR (2009) Noise in multivariate GPS position time-series. J Geod 83:175–187. https://doi.org/10.1007/s00190-008-0251-8
Amiri-Simkooei AR, Tiberius CCJM, Teunissen PJG (2007) Assessment of noise in GPS coordinate time series: Methodology and results. J Geophys Res 112:B07413. https://doi.org/10.1029/2006JB004913
Beirlant J, Dudewicz E, Gyor L, van der Meulen EC (1997) Nonparametric entropy estimation, an overview. Int J Math Statist Sci 6(1):17–39
Bevis M, Bedford J, Caccamise D (2020) The art and science of trajectory modelling. In: Montillet JP, Bos M (eds) Geodetic time series analysis in earth sciences. Springer
Blewitt G, Lavallée D (2002) Effect of annual signals on geodetic velocity. J Geophys Res. https://doi.org/10.1029/2001JB000570
Blewitt G, Kreemer C, Hammond WC, Gazeaux J (2016) MIDAS robust trend estimator for accurate GPS station velocities without step detection. J Geophys Res Solid Earth 121:2054–2068. https://doi.org/10.1002/2015JB012552
Blewitt G, Hammond WC, Kreemer C (2018) Harnessing the GPS data explosion for interdisciplinary science. Eos. https://doi.org/10.1029/2018EO104623
Bos MS, Fernandes RMS, Williams SDP, Bastos L (2008) Fast error analysis of continuous GPS observations. J Geod 82(3):157–166. https://doi.org/10.1007/s00190-007-0165-x
Chen Q, van Dam T, Sneeuw N, Collilieux X, Weigelt M, Rebischung P (2013) Singular spectrum analysis for modeling seasonal signals from GPS time series. J Geodyn 72:25–35. https://doi.org/10.1016/j.jog.2013.05.05
Cover TM, Thomas JA (2006) Elements of information theory, 2nd edn. John Wiley and Sons, New York
Cucci DA, Voirol L, Kermarrec G, Montillet JP, Guerrier S (2023) The generalized method of wavelet moments with eXogenous inputs: a fast approach for the analysis of GNSS position time series. J Geod 97:14. https://doi.org/10.1007/s00190-023-01702-8
Didova O, Gunter B, Riva R, Klees R, Roese-Koerner L (2016) An approach for estimating time-variable rates from geodetic time series. J Geod 90:1207–1221. https://doi.org/10.1007/s00190-016-0918-5
Dmitrieva K, Segall P, DeMets C (2015) Network-based estimation of time-dependent noise in GPS position time series. J Geod 89:591–606. https://doi.org/10.1007/s00190-015-0801-9
Dong D, Fang P, Bock Y, Cheng MK, Miyazaki S (2002) Anatomy of apparent seasonal variations from GPS-derived site position time series. J Geophys Res. https://doi.org/10.1029/2001JB000573
Dong D, Fang P, Bock Y, Webb F, Prawirodirdjo L, Kedar S, Jamason P (2006) Spatiotemporal filtering using principal component analysis and Karhunen-Loeve expansion approaches for regional GPS network analysis. J Geophys Res 111:B03405. https://doi.org/10.1029/2005JB003806
Engels O (2020) Stochastic modelling of geophysical signal constituents within a Kalman filter framework. In: Montillet JP, Bos M (eds) Geodetic time series analysis in earth sciences. Springer
Floyd MA, Herring TA (2020) Fast statistical approaches to geodetic time series analysis. In: Montillet JP, Bos M (eds) Geodetic time series analysis in earth sciences. Springer
Gao W, Li Z, Chen Q (2022) Modelling and prediction of GNSS time series using GBDT, LSTM and SVM machine learning approaches. J Geod 96:71. https://doi.org/10.1007/s00190-022-01662-5
Gobron K, Rebischung P, de Viron O, Demoulin A, Van Camp M (2022) Impact of offsets on assessing the low-frequency stochastic properties of geodetic time series. J Geod 96(7):46. https://doi.org/10.1007/s00190-022-01634-9
Gobron K, Rebischung P, Van Camp M, Demoulin A, de Viron O (2023) Influence of aperiodic non-tidal atmospheric and oceanic loading deformations on the stochastic properties of global GNSS vertical land motion time series. J Geophys Res. https://doi.org/10.1029/2021JB022370
Gualandi A, Serpelloni E (2016) Belardinelli ME (2016) Blind source separation problem in GPS time series. J Geod 90:323–341. https://doi.org/10.1007/s00190-015-0875-4
Hackl M, Malservisi R, Hugentobler U, Wonnacott R (2010) Estimation of velocity uncertainties from GPS time series: examples from the analysis of the South African TrigNet network. J Geophys Res 116:B11404. https://doi.org/10.1029/2010JB008142
Ji KH, Herring T (2011) Transient signal detection using GPS measurements: transient inflation at acutan volcano, Alaska during early 2008. Geophys Res Lett 38:L06307. https://doi.org/10.1029/2011GL046904
Klos A, Olivares G, Teferle FN, Hunegnaw A, Bogusz J (2018) On the combined effect of periodic signals and colored noise on velocity uncertainties. GPS Solut 22:1. https://doi.org/10.1007/s10291-017-0674-x
Langbein J, Johnson H (1997) Correlated errors in geodetic time series: implications for time-dependent deformation. J Geophys Res 102:591–604. https://doi.org/10.1029/96JB02945
Mao A, Harrison GA, Dixon TH (1999) Noise in GPS coordinate time series. J Geophys Res 104(B2):2797–2816. https://doi.org/10.1029/1998JB900033
Montillet JP, Bos MS (2020) Geodetic time series analysis in earth science. Springer Geophys.
Olivares-Pulido G, Teferle FN, Hunegnaw A (2020) Markov chain monte carlo and the application to geodetic time series analysis. In: Montillet JP, Bos M (eds) Geodetic time series analysis in earth sciences. Springer
Press WH, Teukolsky SA, Vetterling WT, Flannery BP (2007) Numerical recipes 3rd edition: the art of scientific computing. Cambridge University Press, New York
Ray J, Altamimi Z, Collilieux X, van Dam T (2008) Anomalous harmonics in the spectra of GPS position estimates. GPS Solut 12(1):55–64. https://doi.org/10.1007/s10291-007-0067-7
Rice JA (1988) Mathematical statistics and data analysis. Wadsworth & Brooks/Cole advanced books & software, Pacific Grove, California
Saleh J (1996) The weak elastic string and some applications in geodesy. J Geod 70:203–213. https://doi.org/10.1007/BF00873701
Santamaria-Gómez A, Bouin MN, Collilieux X, Wöppelmann G (2011) Correlated errors in GPS position time series: Implications for velocity estimates. J Geophys Res 116:B01405. https://doi.org/10.1029/2010JB007701
Santamaría-Gómez A, Ray J (2021) Chameleonic noise in GPS position time series. J Geophys Res Solid Earth 126:e2020JB019541. https://doi.org/10.1029/2020JB019541
Shannon C (1948) A mathematical theory of communication. Bell Sys Tech Journal 27:379–423
Souza E, Monico J (2004) Wavelet shrinkage: high frequency multipath reduction from GPS relative positioning. GPS Solut 8(3):152–159. https://doi.org/10.1007/s10291-004-0100-z
Stone JV (2015) Information theory, a tutorial introduction. Sebtel Press, Sheffield UK
van Dam T, Wahr J, Milly PCD, Shmakin AB, Blewitt G, Lavallée D, Larson KM (2001) Crustal displacements due to continental water loading. Geophys Res Lett 28:651–654. https://doi.org/10.1029/2000GL012120
Vasicek O (1976) A test for normality based on sample entropy. J R Stat Soc B38:54
Wallis KF (2006) A note on the calculation of entropy from histograms. Technical report, Dept. of Economics, University of Warwick, UK
Wang J, Jiang W, Li Z, Lu Y (2021) A new multi-scale sliding window LSTM framework (MSSW-LSTM): a case study for GNSS time-series prediction. Remote Sens 13(16):3328. https://doi.org/10.3390/rs13163328
Wessel P, Smith WHF (1998) New improved version of generic mapping tools released. Eos Trans AGU 79(47):579. https://doi.org/10.1029/98EO00426
Williams SDP (2003) The effect of coloured noise on the uncertainties of rates estimated from geodetic time series. J Geod 76:483–494. https://doi.org/10.1007/s00190-002-0283-4
Williams SD (2008) CATS: GPS coordinate time series analysis software. GPS Solut 12(2):147–153
Williams SDP, Bock Y, Fang P, Jameson P, Nikolaidis R, Prawirodirdjo L, Miller M, Johnson D (2004) Error analysis of continuous GPS position time series. J Geophys Res 109:B03412. https://doi.org/10.1029/2003JB002741
Wu H, Li K, Shi W, Clarke K, Zhang J, Li H (2015) A wavelet-based hybrid approach to remove the flicker noise and the white noise from GPS coordinate time series. GPS Solut 19:511–523. https://doi.org/10.1007/s10291-014-0412-6
Zhang J, Bock Y, Johnson H, Fang P, Williams S, Genrich J, Wdowinski S, Behr J (1997) Southern California permanent GPS geodetic array: error analysis of daily position estimates and site velocities. J Geophys Res 102(B8):18035–18055. https://doi.org/10.1029/97JB01380
Acknowledgements
We thank NGS for supporting this work and John Galetzka for allowing its completion. We thank professor Duncan Agnew and two anonymous reviewers for comments that focused this paper and improved its presentation. Thanks to NGS colleagues who reviewed this paper in-house and to colleagues who continued to answer our questions over the years.
Author information
Authors and Affiliations
Contributions
SW provided simulated noise series and RB facilitated the extraction of the NGL data and the MIDAS software and edited initial versions of the manuscript. As long-time experts on GNSS time series analysis, SW and RB oversaw the soundness of ideas, logic and results. JS did the rest.
Corresponding author
Supplementary Information
Below is the link to the electronic supplementary material.
190_2023_1820_MOESM1_ESM.png
Fig S1: Examples of simulated noise time series (Top) Stationary colored noise consisting of 2.7 mm of white noise and 5 mm of power-law noise with spectral index of -0.9. This series (out of 100 similar ones) was chosen for display because its limited length (of 10,000 days) makes it appear to be non-stationary. The green curve connects the simulated daily noise values and the smoother red curve represents a low-pass (LP) filtered values of the noise using an edge detection filter called “the weak elastic string”, added only to bring out the seemingly non-stationary shape of the series; (Middle) Mildly non-stationary noise series structured as above except that the spectral index is -1.1. (Bottom) Mildly non-stationary noise as in the top plot but with the addition of 0.1 mm of random walk (PNG 319 KB)
190_2023_1820_MOESM2_ESM.png
Fig S2: (Top) Simulated nonlinear transient; (Middle) The transient plus a sequence of simulated colored noise (green curve) and its smoothed values (red curve); (Bottom) The difference between the noisy and smoothed transient. The smoothing was done by the “weak elastic string (PNG 382 KB)
190_2023_1820_MOESM3_ESM.png
Fig S3: Minimum-entropy (red) and MIDAS (blue) estimated horizontal velocities at 171 stations with continuous and linear NGL time series (PNG 381 KB)
190_2023_1820_MOESM4_ESM.png
Fig S4: Horizontal velocity differences: minimum-entropy minus MIDAS estimated at the 171 stations with continuous and linear position time series (PNG 358 KB)
190_2023_1820_MOESM5_ESM.png
Fig S5: Minimal-entropy- and MIDAS-derived velocity uncertainty versus lifespan of station in years, for 55 NGL time series which exhibit visible nonlinearities (mostly in the Up component) but no discontinuities (PNG 125 KB)
190_2023_1820_MOESM6_ESM.png
Fig S6: Histograms (counts versus bins in mm/year) of velocity differences (minimal-entropy minus MIDAS) for 55 NGL time series which exhibit visible nonlinearities (mostly in the Up components) but no discontinuities. The 4-character IDs of these stations are listed in section S5 (PNG 56 KB)
190_2023_1820_MOESM7_ESM.png
Fig S7: Histograms of velocity differences (minimal-entropy minus MIDAS) for time series which experienced one to four discontinuities. Only components which experienced discontinuities are shown. (PNG 56 KB)
190_2023_1820_MOESM8_ESM.png
Fig S8: Velocity standard deviations (SDs) versus time series lifespan derived by minimal-entropy and the MIDAS method for 50 stations with discontinuities. Only components which experienced discontinuities are shown (PNG 121 KB)
Rights and permissions
About this article
Cite this article
Saleh, J., Bennett, R.A. & Williams, S.D.P. Minimum-entropy velocity estimation from GPS position time series. J Geod 98, 11 (2024). https://doi.org/10.1007/s00190-023-01820-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00190-023-01820-3