Higher-order tensor methods for blind spatial signature estimation are proposed. Algorithms able to operate when sources covariance structure is unknown and arbitrary. More efficient denoising based on expanded form of spatial smoothing. Identifiability, uniqueness and computational complexity issues are analyzed.