Gaussian Processes for Blazar Variability Studies