Scalable Variational Gaussian Processes via Harmonic Kernel Decomposition