Spherical radial basis functions model: approximating an integral functional of an isotropic Gaussian random field

Abstract

The spherical radial basis function (SRBF) approach, widely used in gravity modeling, is theoretically surveyed from a viewpoint of random field theory. Let the gravity potential be a random field which is represented as an integral functional of another random field, namely an isotropic Gaussian random field (IGRF) on a sphere inside the Bjerhammar sphere with the SRBF as the integral kernel. When the integration is approximated by a discrete sum within a local region, one gets the widely applicable SRBF model. With this theoretical study, the following two findings are made. First, the IGRF implies a Gaussian prior on the spherical harmonic coefficients (SHCs) of the gravity potential; for this prior the SHCs are independent with each other and their variances are degree-only dependent. This should be reminiscent of two well-known priors, namely the power-law Kaula’s rule and the asymptotic power-law Tscherning-Rapp model. Second, the IGRF-SRBF representation is non-unique. Benefiting from this redundant representation, one can employ a simple IGRF, e.g., the simplest white field, and then design the SRBF accordingly to represent a potential with desired prior statistical properties. This can simplify the corresponding SRBF modeling significantly; to be more specific, the regularization matrix in parameter estimation of the SRBF modeling can be chosen to be a diagonal matrix, or even the naïve identity matrix.