Colored Noise in Dynamical Systems: A New Method for Extracting Laws from Data ✛ ✎

Original: "Colored-LIM: A Data-Driven Method for Studying Dynamical Systems with Temporally Correlated Stochasticity"
arXiv:2402.15184 · 2024-02-23 · arXiv non-exclusive · independent analysis · 4 min
Colored-LIM learns to see memory in random perturbations of complex systems using only time series observations.

Context

This work stands apart among neighbors exploring related topics: polarization recovery from photons, order emerging from noise, tracking heat with a single thermometer, quantum baths, molecular ballet, randomness in evolution, a digital polisher of gravitational lenses, and a neural network that listens to phase. All of them, in one way or another, deal with noise, randomness, and extracting information from data, but Colored-LIM stands out because it directly accounts for the temporal memory of noise in linear inverse models, bridging ideas from stochastic differential equations and dynamic mode decomposition. This fills a gap between classical methods that assume white noise and real complex systems where temporal correlations are significant.

Methods

At the heart of the method is an extended system of equations where the observed state and colored noise evolve jointly. Using the Fokker-Planck equation and its adjoint operator, the authors derive a generalized fluctuation-dissipation relation and expressions for derivatives of the correlation function at zero up to third order. This reduces the inverse problem to solving a system of linear equations if the noise correlation time is known. If the correlation time is unknown, it is found by one-dimensional minimization of the difference between theoretical and observed correlation functions over a given window. The algorithm is computationally efficient and robust, as it requires no additional data beyond the time series of observations.

Results

On a linear system with known parameters, Colored-LIM recovers the dynamical matrix with a relative error of about 10% and the diffusion matrix with an error of about 8% for a time series length of T=1000 and correlation time τ=0.5. Increasing the series length to T=5000 roughly halves the errors. In a nonlinear problem, the network SIS model, the method successfully identifies the topology of interaction network, correctly determining edge weights to within a few percent. Applied to El Niño data (Niño 3.4 SST and D20 indices), Colored-LIM yields a noise correlation time of about 1.85 months, indicating significant memory in atmospheric forcing. Ensemble forecasts with Colored-LIM show an ensemble calibration ratio (ECR) close to 1 for short lead times, whereas classical LIM gives an overdispersed spread and ECR well below 1. In the analysis of Tohoku University electricity consumption, Colored-LIM extracts a diurnal cycle with a period of 22.7 hours and a decay time of 10.1 hours, consistent with human activity, and identifies clusters of buildings with anomalous consumption patterns.

Implications

Unlike neighboring works that either restrict themselves to white noise or use other approaches for nonlinear systems, Colored-LIM is the first to rigorously derive a method for estimating linear dynamics under colored noise using only derivatives of the correlation function. This expands the class of systems for which reliable identification of dynamics from data is possible. It is shown that Colored-LIM does not reduce to classical LIM in the white-noise limit, highlighting the importance of accounting for the temporal structure of noise. The connection to dynamic mode decomposition (DMD) opens a path to applying the method for extracting dynamic modes in the presence of colored noise, which was previously difficult.

Future development

The method can be extended to non-Gaussian colored noises and to systems with multiplicative noise, where noise intensity depends on the state. Also promising is the combination with dimensionality reduction methods such as PCA for analyzing high-dimensional data. Further improvement of numerical differentiation of correlation functions, for example via regularization, will enhance the robustness of derivative estimates in noisy data.

Impact

The method is applicable in climatology (El Niño forecasting), energy (consumption analysis), epidemiology (models of disease spread), and seismology (analysis of sensors with known dynamics).

Next steps

Test the method on a broader class of real systems with known parameters to assess its limits of applicability. Develop an adaptive window selection for minimizing τ and compare with alternative methods for estimating colored noise.

Key open problems

In the field of stochastic dynamics, an open question remains: how do temporal correlations of noise affect predictability and stability of systems? Colored-LIM directly addresses this problem by providing a tool to quantify the effects of noise memory. This is related to the broader problem of identifying nonlinear stochastic systems from limited data, where classical white-noise assumptions are often violated.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The authors propose Colored-LIM, a method that models stochasticity in linear inverse models using Ornstein-Uhlenbeck colored noise and extracts the dynamical and diffusion matrices from the correlation function of observables. The method is validated on idealized linear and nonlinear systems and applied to El Niño and Tohoku University electricity consumption data, demonstrating improved forecast accuracy and interpretability compared to classical LIM.

Where its strength is

The method provides a rigorous and computationally efficient way to account for the temporal memory of noise in linear inverse models, extending the applicability of LIM to real systems with correlated fluctuations.

Where it could go next

Worth trying

  • Test Colored-LIM on data with seasonal nonstationarity, using an approach similar to the work on seasonal climate rhythms.
  • Apply the method to assess causal links in the climate system, comparing with results from causal analysis based on colored noise.
  • Use Colored-LIM as a component of a hybrid model to improve El Niño forecasts on short horizons.

Why this may matter

Accounting for colored noise in linear inverse models could improve the understanding and prediction of systems with slowly varying stochastic influences, such as climate or energy consumption, where noise memory plays an important role.

In the white-noise limit (τ → 0), Colored-LIM does not reduce to classical LIM: the correlation function of colored noise has zero derivative at zero, whereas for white noise the derivative jumps. This means that even arbitrarily small but finite noise memory fundamentally changes the system's behavior near zero lag.

0 = A C_{xx} + C_{xx} A^T + Q B^T + B Q
Relates the dynamical matrix A, covariance Cxx, diffusion matrix Q, and resolvent B = (I - τA)^{-1}, generalizing the classical relation to the case of colored noise.
K''(0) = \frac{1}{2} \left[ A \left( K'(0) + \frac{1}{\tau} C_{xx} \right) + \left( K'(0)^T + \frac{1}{\tau} C_{xx} \right) A^T \right]
Expresses the curvature of the correlation function at zero in terms of the dynamical matrix, covariance, and first derivative, allowing estimation of A when τ is known.
K'''(0) = \frac{1}{2} \left[ \frac{1}{\tau^2} K'(0) + A \left( K''(0) - \frac{1}{\tau^2} C_{xx} \right) - \frac{1}{\tau^2} K'(0)^T - \left( K''(0) - \frac{1}{\tau^2} C_{xx} \right) A^T \right]
Used together with K'(0) and K''(0) to uniquely determine the dynamical matrix A when A is not negative definite.

Key numbers

  • relative error of dynamical matrix (linear system, T=1000): 10.5%
  • relative error of diffusion matrix (linear system, T=1000): 8.6%
  • noise correlation time for El Niño: 1.85 months
  • period of diurnal cycle in electricity consumption: 22.7 hours
  • decay time of diurnal cycle: 10.1 hours
Concepts
linear inverse model colored noise Stochastic modeling biological rhythm modeling labeled data stochastic dynamical system Dynamical Systems nonlinear dynamics nonlinear laplacian spectral analysis Contrastive Learning linear regression probabilistic model statistical method noninvasive dynamic imaging simulated data noise cancellation self-supervised learning physics informed neural network regularization methods Large Language Models