Journal cover Journal topic
Nonlinear Processes in Geophysics An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 1.558
IF1.558
IF 5-year value: 1.475
IF 5-year
1.475
CiteScore value: 2.8
CiteScore
2.8
SNIP value: 0.921
SNIP0.921
IPP value: 1.56
IPP1.56
SJR value: 0.571
SJR0.571
Scimago H <br class='widget-line-break'>index value: 55
Scimago H
index
55
h5-index value: 22
h5-index22
Volume 15, issue 1
Nonlin. Processes Geophys., 15, 159–167, 2008
https://doi.org/10.5194/npg-15-159-2008
© Author(s) 2008. This work is licensed under
the Creative Commons Attribution-NonCommercial-ShareAlike 2.5 License.

Special issue: Extreme Events: Nonlinear Dynamics and Time Series Analysis

Nonlin. Processes Geophys., 15, 159–167, 2008
https://doi.org/10.5194/npg-15-159-2008
© Author(s) 2008. This work is licensed under
the Creative Commons Attribution-NonCommercial-ShareAlike 2.5 License.

  19 Feb 2008

19 Feb 2008

Detecting spatial patterns with the cumulant function – Part 1: The theory

A. Bernacchia1 and P. Naveau2 A. Bernacchia and P. Naveau
  • 1Dipartimento di Fisica E.Fermi, Universita' La Sapienza, Roma, Italy
  • 2Laboratoire des Sciences du Climat et de l'Environnement, IPSL-CNRS, France

Abstract. In climate studies, detecting spatial patterns that largely deviate from the sample mean still remains a statistical challenge. Although a Principal Component Analysis (PCA), or equivalently a Empirical Orthogonal Functions (EOF) decomposition, is often applied for this purpose, it provides meaningful results only if the underlying multivariate distribution is Gaussian. Indeed, PCA is based on optimizing second order moments, and the covariance matrix captures the full dependence structure of multivariate Gaussian vectors. Whenever the application at hand can not satisfy this normality hypothesis (e.g. precipitation data), alternatives and/or improvements to PCA have to be developed and studied.

To go beyond this second order statistics constraint, that limits the applicability of the PCA, we take advantage of the cumulant function that can produce higher order moments information. The cumulant function, well-known in the statistical literature, allows us to propose a new, simple and fast procedure to identify spatial patterns for non-Gaussian data. Our algorithm consists in maximizing the cumulant function. Three families of multivariate random vectors, for which explicit computations are obtained, are implemented to illustrate our approach. In addition, we show that our algorithm corresponds to selecting the directions along which projected data display the largest spread over the marginal probability density tails.

Publications Copernicus
Download
Citation