Gaussian variational approximation for high-dimensional state space models