Abstract
The spherical Slepian functions can be used to localize the solutions of the gradiometric boundary value problems on a sphere. These functions involve spatially restricted integral products of scalar, vector and tensor spherical harmonics. This paper formulates these integrals in terms of combinations of the Gaunt coefficients and integrals of associated Legendre functions. The presented formulas for these integrals are useful in recovering the Earth's gravity field locally from the satellite gravity gradiometry data.