Yes, I am aware of this, but it doesn’t work with a GP term. This kind of phylogenetic model doesn’t work exactly the same way, and in any case, it wouldn’t work for modelling geodesic distances with a GP.
My issue is essentially the same as what’s discussed in this thread: Custom distance metric with gaussian process - #6 by Matias_Guzman_Naranjo