A procedure by which a spatiotemporal data set is decomposed into its leading patterns in both time (see ‘Principal Component’) and space (see ‘Empirical Orthogonal Function’) based on an orthogonal decomposition of the data covariance matrix.