Currently, we use the legacy [`sklearn.gaussian_process.GaussianProcess`](https://scikit-learn.org/0.17/modules/gaussian_process.html) API, which was [deprecated and removed in 0.18](https://scikit-learn.org/0.18/modules/gaussian_process.html#legacy-gaussian-processes). We should update to the [new `GaussianProcessRegressor` API](https://scikit-learn.org/stable/modules/gaussian_process.html#gaussian-process-regression-gpr). See [this issue](https://github.com/scikit-learn/scikit-learn/issues/6657) for an example of upgrading to the new API.