Gaussian process and time series


I’m trying to reproduce this example, but instead of using MAP, I’m using NUTS samples.

I’m just asking for 2 chains (70 samples each one) and it is extremely slow. Is there any improvement in the estimates when using NUTS over MAP with time series?

Could you suggest references about the use of gaussian processes in time series analysis?

Thanks & regards.