There’s some new stuff in brms that should make GPs easier to work with that I think is relevant here.
Have a look here: Practical Hilbert space approximate Bayesian Gaussian processes for probabilistic programming and the associated paper: https://arxiv.org/abs/2004.11408 (especially the examples at the end).