Markov Chain Monte Carlo and the Application to Geodetic Time Series Analysis

The time evolution of geophysical phenomena can be characterised by stochastic time series. The stochastic nature of the signal stems from the geophysical phenomena involved and any noise, which may be due to, e.g., un-modelled effects or measurement errors. Until the 1990s, it was usually assumed that white noise could fully characterise this noise. However, this was demonstrated to not be the case and it was proven that this assumption leads to underestimated uncertainties of the geophysical parameters inferred from the geodetic time series. Therefore, in order to fully quantify all the uncertainties as robustly as possible, it is imperative to estimate not only the deterministic but also the stochastic model parameters of the time series. In this regard, the Markov Chain Monte Carlo (MCMC) method can provide a sample of the distribution function of all parameters, including those regarding the noise, e.g., spectral index and amplitudes. After presenting the MCMC method and its implementation in our MCMC software we apply it to synthetic and real time series and perform a cross-evaluation using Maximum Likelihood Estimation (MLE) as implemented in the CATS software. Several examples as to how the MCMC method performs as a parameter estimation method for geodetic time series are given in this chapter. These include the applications to GPS position time series, superconducting gravity time series and monthly mean sea level (MSL) records, which all show very different stochastic properties. The impact of the estimated parameter uncertainties on sub-sequentially derived products is briefly demonstrated for the case of plate motion models. Finally, an evaluation of the MCMC method against the Hector method using weekly downsampled versions of the Benchmark Synthetic GNSS (BSG) time series as provided in Chap. 2 is presented separately in an appendix.

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

[2]  M. Bouin,et al.  Geocentric sea-level trend estimates from GPS analyses at relevant tide gauges world-wide , 2007 .

[3]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[4]  M. Bouin,et al.  Rates of sea‐level change over the past century in a geocentric reference frame , 2009 .

[5]  W. Szeliga,et al.  GPS Vertical Land Motion Corrections to Sea‐Level Rise Estimates in the Pacific Northwest , 2018 .

[6]  M. Vanclooster,et al.  Hydrogeological investigations at the Membach station, Belgium, and application to correct long periodic gravity variations , 2006 .

[7]  H. Virtanen.  Loading effects in Metsähovi from the atmosphere and the Baltic Sea , 2004 .

[8]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[9]  C. Skordis,et al.  Fast and reliable Markov chain Monte Carlo technique for cosmological parameter estimation , 2005 .

[10]  J. Gregory,et al.  Coastal and global averaged sea level rise for 1950 to 2000 , 2005 .

[11]  Simon D. P. Williams,et al.  Offsets in Global Positioning System time series , 2003 .

[12]  D. Agnew,et al.  The time-domain behavior of power-law noises. [of many geophysical phenomena] , 1992 .

[13]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[14]  C. Harrison.  Power spectrum of sea level change over fifteen decades of frequency , 2002 .

[15]  O. Francis,et al.  Geophysics From Terrestrial Time‐Variable Gravity Measurements , 2017 .

[16]  Alessandro Caporali,et al.  Statistical analysis of the time series of permanent GPS stations , 2003 .

[17]  Richard Bingley,et al.  Trends in Uk Mean Sea Level Revisited , 2008 .

[18]  Janusz Bogusz,et al.  On the combined effect of periodic signals and colored noise on velocity uncertainties , 2017, GPS Solutions.

[19]  M. Camp,et al.  Optimized strategy for the calibration of superconducting gravimeters at the one per mille level , 2015, Journal of Geodesy.

[20]  Michael B. Heflin,et al.  Interseismic strain accumulation and anthropogenic motion in metropolitan Los Angeles , 2005 .

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

[22]  D. McLaughlin,et al.  Crustal motions in Great Britain: evidence from continuous GPS, absolute gravity and Holocene sea level data , 2009 .

[23]  Felix Norman Teferle,et al.  Glacial isostatic adjustment of the British Isles: new constraints from GPS measurements of crustal motion , 2009 .

[24]  O. Francis,et al.  Uncertainty of absolute gravity measurements , 2005 .

[25]  J. Kirchner.  Aliasing in 1/f(alpha) noise spectra: origins, consequences, and remedies. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Duncan Carr Agnew,et al.  Application of the global positioning system to crustal deformation measurement: 1. Precision and accuracy , 1991 .

[27]  Manuela Seitz,et al.  Consistent estimation of geodetic parameters from SLR satellite constellation measurements , 2018, Journal of Geodesy.

[28]  P. van Gelder,et al.  The Effect of the 18.6-Year Lunar Nodal Cycle on Regional Sea-Level Rise Estimates , 2012 .

[29]  Michael R Pearlman,et al.  THE INTERNATIONAL LASER RANGING SERVICE , 2007 .

[30]  Xavier Collilieux,et al.  Hydrological deformation induced by the West African Monsoon: Comparison of GPS, GRACE and loading models , 2012 .

[31]  G. Olivares,et al.  A Bayesian Monte Carlo Markov Chain Method for Parameter Estimation of Fractional Differenced Gaussian Processes , 2013, IEEE Transactions on Signal Processing.

[32]  S. Jevrejeva,et al.  New Data Systems and Products at the Permanent Service for Mean Sea Level , 2013 .

[33]  Richard Biancale,et al.  Space tracking system improves accuracy of geodetic measurements , 1996 .

[34]  J. Johansson,et al.  Space-Geodetic Constraints on Glacial Isostatic Adjustment in Fennoscandia , 2001, Science.

[35]  Benoit B. Mandelbrot,et al.  Fractal Geometry of Nature , 1984 .

[36]  Boudewijn Ambrosius,et al.  Angular velocities of Nubia and Somalia from continuous GPS data: implications on present-day relative kinematics , 2004 .

[37]  S. Mazzotti,et al.  Relative and absolute sea level rise in western Canada and northwestern United States from a combined tide gauge‐GPS analysis , 2008 .

[38]  S. Williams,et al.  The color of sea level: Importance of spatial variations in spectral shape for assessing the significance of trends , 2010 .

[39]  F. Wyatt.  Displacement of surface monuments: Horizontal motion , 1982 .

[40]  Matt A. King,et al.  Detecting offsets in GPS time series: First results from the detection of offsets in GPS experiment , 2013 .

[41]  Markus Rothacher,et al.  The International GPS Service (IGS): An interdisciplinary service in support of Earth sciences , 1999 .

[42]  Geoffrey Blewitt,et al.  Effect of annual signals on geodetic velocity , 2002 .

[43]  O. Francis,et al.  Is the instrumental drift of superconducting gravimeters a linear or exponential function of time? , 2007 .

[44]  Fast and reliable MCMC for cosmological parameter estimation , 2004, astro-ph/0405462.

[45]  Xianghong Hua,et al.  Review of current GPS methodologies for producing accurate time series and their error sources , 2017 .

[46]  C. Boucher,et al.  Positioning results with doris on SPOT2 after first year of mission , 1992 .

[47]  J. Rosenthal,et al.  Optimal scaling for various Metropolis-Hastings algorithms , 2001 .

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

[49]  R. Steven Nerem,et al.  Is there a 60‐year oscillation in global mean sea level? , 2012 .

[50]  F. Wyatt.  Displacement of surface monuments: Vertical motion , 1989 .

[51]  Chris Rizos,et al.  The International GNSS Service in a changing landscape of Global Navigation Satellite Systems , 2009 .

[52]  Axel Nothnagel,et al.  International VLBI Service for Geodesy and Astrometry , 2017, Journal of Geodesy.

[53]  M. Bevis,et al.  Vertical and horizontal surface displacements near Jakobshavn Isbræ driven by melt‐induced and dynamic ice loss , 2012 .

[54]  Alan Dodson,et al.  Using GPS to Separate Crustal Movements and Sea Level Changes at Tide Gauges in the UK , 2002 .

[55]  Hans-Peter Plag,et al.  A continuous GPS coordinate time series analysis strategy for high-accuracy vertical land movements , 2008 .

[56]  Rui Fernandes,et al.  The influence of seasonal signals on the estimation of the tectonic motion in short continuous GPS time-series , 2010 .

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

[58]  J. Langbein,et al.  Sensitivity of crustal deformation instruments to changes in secular rate , 1993 .

[59]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

[60]  É. Calais,et al.  Geodetic observations of interseismic strain segmentation at the Sumatra Subduction Zone , 1997 .

[61]  Peter Steigenberger,et al.  Multi-technique comparison of troposphere zenith delays and gradients during CONT08 , 2011 .

[62]  S. Holgate,et al.  Evidence for enhanced coastal sea level rise during the 1990s , 2004 .

[63]  N. White,et al.  A 20th century acceleration in global sea‐level rise , 2006 .

[64]  Xavier Collilieux,et al.  ITRF2008 plate motion model , 2011 .

[65]  W. Prothero,et al.  The superconducting gravimeter , 1968 .

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

[67]  J. Montillet,et al.  Investigation of the noise properties at low frequencies in long GNSS time series , 2019, Journal of Geodesy.

[68]  X. Collilieux,et al.  Mitigating the effects of vertical land motion in tide gauge records using a state-of-the-art GPS velocity field , 2012 .

[69]  U. Hugentobler,et al.  Estimation of velocity uncertainties from GPS time series: Examples from the analysis of the South African TrigNet network , 2011 .

[70]  Axel Nothnagel,et al.  IVS and its important role in the maintenance of the global reference systems , 2002 .

[71]  P. Tregoning,et al.  Characterizing and minimizing the effects of noise in tide gauge time series: relative and geocentric sea level rise around Australia , 2013 .

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

[73]  B. Mandelbrot,et al.  Fractional Brownian Motions, Fractional Noises and Applications , 1968 .

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

[75]  Shfaqat Abbas Khan,et al.  Geodetic measurements of postglacial adjustments in Greenland , 2008 .

[76]  Statistical significance of trends in Zenith Wet Delay from re-processed GPS solutions , 2018, GPS Solutions.

[77]  Simon D. P. Williams,et al.  CATS: GPS coordinate time series analysis software , 2008 .

[78]  Pascal Willis,et al.  The International DORIS Service (IDS): Toward maturity , 2010 .

[79]  John Langbein,et al.  Noise in GPS displacement measurements from Southern California and Southern Nevada , 2008 .

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

[81]  P. Woodworth,et al.  A review of the trends observed in British Isles mean sea level data measured by tide gauges , 1999 .