Many applications in spatial statistics involve data observed over large regions on the Earth’s surface. There is a large statistical literature devoted to covariance functions capable of modeling the degree of smoothness in data on Euclidean spaces. We adapt some of this work to covariance functions for processes on spheres, where the natural distance is great circle distance. In doing so, we ...