Nonparametric Bayesian modelling of longitudinally integrated covariance functions on spheres

Taking into account axial symmetry in the covariance function of a Gaussian random field is essential when the purpose is modelling data defined over a large portion of the sphere representing our planet. Axially symmetric covariance functions admit a convoluted spectral representation that makes modelling and inference difficult. This motivates the interest in devising alternative strategies to attain axial symmetry, an appealing option being longitudinal integration of isotropic random fields on the sphere. This paper provides a comprehensive theoretical framework to model longitudinal integration on spheres through a nonparametric Bayesian approach. Longitudinally integrated covariances are treated as random objects, where the randomness is implied by the randomised spectrum associated with the covariance function. After investigating the topological support induced by our construction, we give the posterior distribution a thorough inspection. A Bayesian nonparametric model for the analysis of data defined on the sphere is described and implemented, its performance investigated by means of the analysis of both simulated and real data sets.