We study a class of separable sample covariance matrices of the form r Q1 :" r A 1{2 X r BX ˚r A 1{2 . Here r A and r B are positive definite matrices whose spectrums consist of bulk spectrums plus several spikes, i.e. larger eigenvalues that are separated from the bulks. Conceptually, we call r Q1 a spiked separable covariance matrix model. On the one hand, this model includes the spiked covariance matrix as a special case with r B " I. On the other hand, it allows for more general correlations of datasets. In particular, for spatio-temporal dataset, r A and r B represent the spatial and temporal correlations, respectively. In this paper, we study the outlier eigenvalues and eigenvectors, i.e. the principal components, of the spiked separable covariance model r Q1. We prove the convergence of the outlier eigenvalues r λi and the generalized components (i.e. xv, r ξiy for any deterministic vector v) of the outlier eigenvectors r ξi with optimal convergence rates. Moreover, we also prove the delocalization of the non-outlier eigenvectors. We state our results in full generality, in the sense that they also hold near the so-called BBP transition and for degenerate outliers. Our results highlight both the similarity and difference between the spiked separable covariance matrix model and the spiked covariance model matrix in [8]. In particular, we show that the spikes of both r A and r B will cause outliers of the eigenvalue spectrum, and the eigenvectors can help to select the outliers that correspond to the spikes of r A (or r B). Contents 1 Introduction 2 Definitions and main results 3 Statistical estimation for spiked separable covariance matrices 4 Basic tools