Identifying temporally invariant components in complex multivariate time series is key to understanding the underlying dynamical system and predict its future behavior. In this Letter, we propose a novel technique, stationary subspace analysis (SSA), that decomposes a multivariate time series into its stationary and nonstationary part. The method is based on two assumptions: (a) the observed signals are linear superpositions of stationary and nonstationary sources; and (b) the nonstationarity is measurable in the first two moments. We characterize theoretical and practical properties of SSA and study it in simulations and cortical signals measured by electroencephalography. Here, SSA succeeds in finding stationary components that lead to a significantly improved prediction accuracy and meaningful topographic maps which contribute to a better understanding of the underlying nonstationary brain processes.