We consider noise subspace methods for narrowband direction-of-arrival or harmonic retrieval in colored linear non-gaussian noise of unknown covariance and unknown distribution. The non-gaussian noise covariance is estimated via higher order cumulants and combined with correlation information to solve a generalized eigenvalue problem. The estimated eigenvectors are used in a variety of noise subspace methods such as multiple signal classification (MUSIC), MVDR and eigenvector. The noise covariance estimates are obtained in the presence of the harmonic signals, obviating the need for noise-only training records. The covariance estimates may be obtained nonparametrically via cumulant projections, or parametrically using autoregressive moving average (ARMA) models. An information theoretic criterion using higher order cumulants is presented which may be used to simultaneously estimate the ARMA model order and parameters. Third- and fourth-order cumulants are employed for asymmetric and symmetric probability density function (pdf) cases, respectively. Simulation results show considerable improvement over conventional methods with no prewhitening. The effects of prewhitening are particularly evident in the dominant eigenvalues, as revealed by singular value decomposition (SVD) analysis. >