Simple GP model giving large MAP values for covariance timescale

I messed around with your code and found the same thing. I have a feeling that you’re experiencing the same issue that Michael Betancourt talks about in this blog post. See this thread for some advice. Long story short, priors are important (but especially with GPs)!