An information-geometric approach to feature extraction and moment reconstruction in dynamical systems

TL;DR

This paper introduces an information-geometric framework for feature extraction and moment reconstruction in dynamical systems by analyzing probability measures induced by system observables, enabling nonparametric forecasting and efficient analysis.

Contribution

It develops a novel dimension reduction method based on probability measures and eigenfunctions, linking moments evolution to a time-averaging operator, with applications to complex dynamical systems.

Findings

Eigenfunctions capture multiple timescales of the system.

Few eigenvectors suffice for accurate moment reconstruction.

Method successfully applied to atmospheric time series data.

Abstract

We propose a dimension reduction framework for feature extraction and moment reconstruction in dynamical systems that operates on spaces of probability measures induced by observables of the system rather than directly in the original data space of the observables themselves as in more conventional methods. Our approach is based on the fact that orbits of a dynamical system induce probability measures over the measurable space defined by (partial) observations of the system. We equip the space of these probability measures with a divergence, i.e., a distance between probability distributions, and use this divergence to define a kernel integral operator. The eigenfunctions of this operator create an orthonormal basis of functions that capture different timescales of the dynamical system. One of our main results shows that the evolution of the moments of the dynamics-dependent probability measures can be related to a time-averaging operator on the original dynamical system. Using this result, we show that the moments can be expanded in the eigenfunction basis, thus opening up the avenue for nonparametric forecasting of the moments. If the collection of probability measures is itself a manifold, we can in addition equip the statistical manifold with the Riemannian metric and use techniques from information geometry. We present applications to ergodic dynamical systems on the 2-torus and the Lorenz 63 system, and show on a real-world example that a small number of eigenvectors is sufficient to reconstruct the moments (here the first four moments) of an atmospheric time series, i.e., the realtime multivariate Madden-Julian oscillation index.