Independent Subspace Analysis of the Sea Surface Temperature Variability: Non-Gaussian Sources and Sensitivity to Sampling and Dimensionality

Abstract

We propose an expansion of multivariate time-series data into maximally independent source subspaces. The search is made among rotations of prewhitened data which maximize non-Gaussianity of candidate sources. We use a tensorial invariant approximation of the multivariate negentropy in terms of a linear combination of squared coskewness and cokurtosis. By solving a high-order singular value decomposition problem, we extract the axes associated with most non-Gaussianity. Moreover, an estimate of the Gaussian subspace is provided by the trailing singular vectors. The independent subspaces are obtained through the search of "quasiindependent" components within the estimated non-Gaussian subspace, followed by the identification of groups with significant joint negentropies. Sources result essentially from the coherency of extremes of the data components. The method is then applied to the global sea surface temperature anomalies, equatorward of 65 ∘ , after being tested with non-Gaussian surrogates consistent with the data anomalies. The main emerging independent components and subspaces, supposedly generated by independent forcing, include different variability modes, namely, The East-Pacific, the Central Pacific, and the Atlantic Niños, the Atlantic Multidecadal Oscillation, along with the subtropical dipoles in the Indian, South Pacific, and South-Atlantic oceans. Benefits and usefulness of independent subspaces are then discussed.

Subspace Detection and Blind Source Separation of Multivariate Signals by Dynamical Component Analysis (DyCA)

“…In geophysical applications EOF analysis have been applied to spatio-temporal climate data to obtain data-driven models [10] and have been further developed in different directions like independent subspace analysis [11], linear and nonlinear dynamical mode decomposition [12], [13] to name only a few.…”

“…Furthermore we now also get an estimate for the parameter n, which we call n. In order to explain the idea behind this estimation procedure, however, we first return to the original u i and v i from ( 9) and (11). In consequence of the linear independence of {u 1 , .…”

Generalized Skew-Normal Negentropy and Its Application to Fish Condition Factor Time Series

“…Another way to obtain these derivatives is to define random variable Z * τ = u(Z τ ) = Z τ / 1 + Z 2 τ and using (14) with Z τ and µ k = E(Z τ ) replaced by Z * τ and µ * k = E{(Z * τ ) k }, respectively. Thus, we obtain the series expansion:…”

Bispectral analysis of nonlinear interaction, predictability and stochastic modelling with application to ENSO

“…Whereas external orbital variations are believed to be the dominant driving force for macroclimate (on millennial time scales), weather, macro-weather (Lovejoy, 2018) and climate variations (on shorter time scales) are mainly the result of complex nonlinear interactions between very many degrees-of-freedom (Sura and Hannachi, 2015) and also due to many climate subcomponents with different time scales The atmospheric (and climate) system is an excellent example of highdimensional and highly complex dynamical system. One outstanding and ubiquitous feature of the large scale (and low frequency) atmospheric (and climate) variability is non-Gaussianity (Franzke et al, 2007;Proistosescu et al, 2016;Pires and Hannachi, 2017;Hannachi and Iqbal, 2019). For instance, Sura and Sardeshmukh (2008) show that sea surface temperature (SST) has non-Gaussian probability distribution function (PDF) with particular tail extrema.…”