Multiresolution optimal interpolation and statistical analysis of TOPEX/POSEIDON satellite altimetry

Abstract

A recently developed multiresolution estimation framework offers the possibility of highly efficient statistical analysis, interpolation, and smoothing of extremely large data sets in a multiscale fashion. This framework enjoys a number of advantages not shared by other statistically-based methods. In particular, the algorithms resulting from this framework have complexity that scales only linearly with problem size, yielding constant complexity load per grid point independent of problem size. Furthermore these algorithms directly provide interpolated estimates at multiple resolutions, accompanying error variance statistics of use in assessing resolutionlaccuracy tradeoffs and in detecting statistically significant anomalies, and maximum likelihood estimates of parameters such as spectral power law coefficients. Moreover, the efficiency of these algorithms is completely insensitive to irregularities in the sampling or spatial distribution of measurements and to heterogeneities in measurement errors or model parameters. For these reasons this approach has the potential of being an effective tool in a variety of remote sensing problems. In this paper, we demonstrate a realization of this potential by applying the multiresolution framework to a problem of considerable current interest-the interpolation and statistical analysis of ocean surface data from the TOPEXPOSEIDON altimeter