Maybe a global Moran's I test?

https://en.m.wikipedia.org/wiki/Moran%27s_I

Create the spatial weights matrix with great circle distances?