Uncovering temporal regularity in atmospheric dynamics through Hilbert phase analysis

Uncovering meaningful regularities in complex oscillatory signals is a challenging problem with applications across a wide range of disciplines. Here, we present a novel approach, based on the Hilbert transform (HT). We show that temporal periodicity can be uncovered by averaging the signal in a moving window of appropriated length, τ, before applying the HT. As a case study, we investigate global gridded surface air temperature (SAT) datasets. By analyzing the variation of the mean rotation period, T¯, of the Hilbert phase as a function of τ, we discover well-defined plateaus. In many geographical regions, the plateau corresponds to the expected 1-yr solar cycle; however, in regions where SAT dynamics is highly irregular, the plateaus reveal non-trivial periodicities, which can be interpreted in terms of climatic phenomena such as El Niño. In these regions, we also find that Fourier analysis is unable to detect the periodicity that emerges when τ increases and gradually washes out SAT variability. The values of T¯ obtained for different τs are then given to a standard machine learning algorithm. The results demonstrate that these features are informative and constitute a new approach for SAT time series classification. To support these results, we analyze the synthetic time series generated with a simple model and confirm that our method extracts information that is fully consistent with our knowledge of the model that generates the data. Remarkably, the variation of T¯ with τ in the synthetic data is similar to that observed in the real SAT data. This suggests that our model contains the basic mechanisms underlying the unveiled periodicities. Our results demonstrate that Hilbert analysis combined with temporal averaging is a powerful new tool for discovering hidden temporal regularity in complex oscillatory signals.