TY - JOUR
T1 - Reconstructing time series and their uncertainty from observations with universal noise
AU - Sambridge, M.
N1 - Publisher Copyright:
©2016. American Geophysical Union. All Rights Reserved.
PY - 2016/7/1
Y1 - 2016/7/1
N2 - A new approach is presented for the reconstruction of time series and other (y,x) functions from observables with any type of stochastic noise. In particular, noise may exist in both dependent and independent variables, i.e., y and x, or t, and may even be correlated between these variables. This situation occurs in many areas of the geosciences when the independent time variable is itself the result of a measurement process, such as in paleo– sea level estimation. Uncertainty in the recovered time series is quantified in probabilistic terms using Bayesian changepoint modeling. The main contribution of the paper is the derivation of a new form of integrated likelihood function which can measure the data fit for a curve to (y,t) observables contaminated by any type of random noise. Closed form expressions are found for the special case of correlated Gaussian data noise and curves built from the sum of piecewise linear polynomials. The technique is illustrated by estimating relative sea level variations, over the last five glacial cycles, from a data set of 1928 δ18O measurements. Comparisons are also made with other techniques including those that assume an error free “independent” variable. Experiments illustrate several benefits of accounting for timing errors. These include allowing rigorous uncertainty information of both time-dependent signals and their gradients. Derivatives of the integrated likelihood function are also given, which allow implementation of likelihood maximization. The new likelihood function better reflects real errors in data and can improve recovery of the estimated time series.
AB - A new approach is presented for the reconstruction of time series and other (y,x) functions from observables with any type of stochastic noise. In particular, noise may exist in both dependent and independent variables, i.e., y and x, or t, and may even be correlated between these variables. This situation occurs in many areas of the geosciences when the independent time variable is itself the result of a measurement process, such as in paleo– sea level estimation. Uncertainty in the recovered time series is quantified in probabilistic terms using Bayesian changepoint modeling. The main contribution of the paper is the derivation of a new form of integrated likelihood function which can measure the data fit for a curve to (y,t) observables contaminated by any type of random noise. Closed form expressions are found for the special case of correlated Gaussian data noise and curves built from the sum of piecewise linear polynomials. The technique is illustrated by estimating relative sea level variations, over the last five glacial cycles, from a data set of 1928 δ18O measurements. Comparisons are also made with other techniques including those that assume an error free “independent” variable. Experiments illustrate several benefits of accounting for timing errors. These include allowing rigorous uncertainty information of both time-dependent signals and their gradients. Derivatives of the integrated likelihood function are also given, which allow implementation of likelihood maximization. The new likelihood function better reflects real errors in data and can improve recovery of the estimated time series.
KW - inversion
KW - likelihood
KW - regression
KW - sea level
UR - http://www.scopus.com/inward/record.url?scp=84978083102&partnerID=8YFLogxK
U2 - 10.1002/2016JB012901
DO - 10.1002/2016JB012901
M3 - Article
SN - 2169-9313
VL - 121
SP - 4990
EP - 5012
JO - Journal of Geophysical Research: Solid Earth
JF - Journal of Geophysical Research: Solid Earth
IS - 7
ER -