An approach for estimating time-variable rates from geodetic time series

There has been considerable research in the literature focused on computing and forecasting sea-level changes in terms of constant trends or rates. The Antarctic ice sheet is one of the main contributors to sea-level change with highly uncertain rates of glacial thinning and accumulation. Geodetic observing systems such as the Gravity Recovery and Climate Experiment (GRACE) and the Global Positioning System (GPS) are routinely used to estimate these trends. In an effort to improve the accuracy and reliability of these trends, this study investigates a technique that allows the estimated rates, along with co-estimated seasonal components, to vary in time. For this, state space models are defined and then solved by a Kalman filter (KF). The reliable estimation of noise parameters is one of the main problems encountered when using a KF approach, which is solved by numerically optimizing likelihood. Since the optimization problem is non-convex, it is challenging to find an optimal solution. To address this issue, we limited the parameter search space using classical least-squares adjustment (LSA). In this context, we also tested the usage of inequality constraints by directly verifying whether they are supported by the data. The suggested technique for time-series analysis is expanded to classify and handle time-correlated observational noise within the state space framework. The performance of the method is demonstrated using GRACE and GPS data at the CAS1 station located in East Antarctica and compared to commonly used LSA. The results suggest that the outlined technique allows for more reliable trend estimates, as well as for more physically valuable interpretations, while validating independent observing systems.

[1]  H. Farahani Modelling the Earth's static and time-varying gravity field using a combination of GRACE and GOCE data , 2013 .

[2]  B. Legrésy,et al.  Quantifying the seasonal “breathing” of the Antarctic ice sheet , 2012 .

[3]  M. Tamisiea,et al.  On seasonal signals in geodetic time series , 2012 .

[4]  Olga Didova,et al.  Empirical estimation of present-day Antarctic glacial isostatic adjustment and ice mass change , 2013, The Cryosphere.

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

[6]  John A. D. Aston,et al.  The State Space Models Toolbox for MATLAB , 2011 .

[7]  Eric Rignot,et al.  A Reconciled Estimate of Ice-Sheet Mass Balance , 2012, Science.

[8]  E. Ivins,et al.  Antarctic contribution to sea level rise observed by GRACE with improved GIA correction , 2013 .

[9]  S. S. Wilks The Large-Sample Distribution of the Likelihood Ratio for Testing Composite Hypotheses , 1938 .

[10]  Paul Segall,et al.  Network-based estimation of time-dependent noise in GPS position time series , 2015, Journal of Geodesy.

[11]  Ronald K. Pearson,et al.  Exploring Data in Engineering, the Sciences, and Medicine , 2011 .

[12]  Robert E. Grimm,et al.  The role of acids in electrical conduction through ice , 2013 .

[13]  Andrew Harvey,et al.  Forecasting, Structural Time Series Models and the Kalman Filter. , 1991 .

[14]  A. Bryson,et al.  Linear filtering for time-varying systems using measurements containing colored noise , 1965 .

[15]  Oliver Baur,et al.  On the computation of mass-change trends from GRACE gravity field time-series , 2012 .

[16]  P. Bloomfield,et al.  Properties of the random search in global optimization , 1975 .

[17]  I. Sasgen,et al.  Antarctic ice-mass balance 2003 to 2012: regional reanalysis of GRACE satellite gravimetry measurements with improved estimate of glacial-isostatic adjustment based on GPS uplift rates , 2013 .

[18]  J. Wahr,et al.  Measurements of Time-Variable Gravity Show Mass Loss in Antarctica , 2006, Science.

[19]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[20]  Stephen J. Wright,et al.  Numerical Optimization (Springer Series in Operations Research and Financial Engineering) , 2000 .

[21]  Peter J. Clarke,et al.  Widespread low rates of Antarctic glacial isostatic adjustment revealed by GPS observations , 2011 .

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

[23]  Matt A. King,et al.  A new glacial isostatic adjustment model for Antarctica: calibrated and tested using observations of relative sea‐level change and present‐day uplift rates , 2012 .

[24]  Isabella Velicogna,et al.  Regional acceleration in ice mass loss from Greenland and Antarctica using GRACE time‐variable gravity data , 2014 .

[25]  Jürgen Kusche,et al.  Surface mass redistribution inversion from global GPS deformation and Gravity Recovery and Climate Experiment (GRACE) gravity data , 2005 .

[26]  J. Camp,et al.  Antarctica, Greenland and Gulf of Alaska land-ice evolution from an iterated GRACE global mascon solution , 2013, Journal of Glaciology.

[27]  Jorge Nocedal,et al.  An Interior Point Algorithm for Large-Scale Nonlinear Programming , 1999, SIAM J. Optim..

[28]  Hubert H. G. Savenije,et al.  The design of an optimal filter for monthly GRACE gravity models , 2008 .

[29]  P. Broersen,et al.  How to handle colored observation noise in large least-squares problems , 2003 .

[30]  Reiner Horst,et al.  Introduction to Global Optimization (Nonconvex Optimization and Its Applications) , 2002 .

[31]  B. Farhang-Boroujeny,et al.  Adaptive Filters: Theory and Applications , 1999 .

[32]  Siem Jan Koopman,et al.  Time Series Analysis by State Space Methods , 2001 .

[33]  P. Whitehouse,et al.  Effect of GIA models with 3D composite mantle viscosity on GRACE mass balance estimates for Antarctica , 2015 .

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

[35]  N. Sneeuw,et al.  A stochastic framework for inequality constrained estimation , 2012, Journal of Geodesy.