Articles | Volume 14, issue 1
Nonlin. Processes Geophys., 14, 73–77, 2007
https://doi.org/10.5194/npg-14-73-2007
Nonlin. Processes Geophys., 14, 73–77, 2007
https://doi.org/10.5194/npg-14-73-2007

  01 Feb 2007

01 Feb 2007

Expectation-maximization analysis of spatial time series

K. W. Smith and A. L. Aretxabaleta K. W. Smith and A. L. Aretxabaleta
  • Woods Hole Oceanographic Institution, Woods Hole, MA, USA

Abstract. Expectation maximization (EM) is used to estimate the parameters of a Gaussian Mixture Model for spatial time series data. The method is presented as an alternative and complement to Empirical Orthogonal Function (EOF) analysis. The resulting weights, associating time points with component distributions, are used to distinguish physical regimes. The method is applied to equatorial Pacific sea surface temperature data from the TAO/TRITON mooring time series. Effectively, the EM algorithm partitions the time series into El Niño, La Niña and normal conditions. The EM method leads to a clearer interpretation of the variability associated with each regime than the basic EOF analysis.