Modulating Scalable Gaussian Processes for Expressive Statistical Learning