The present work considers the interpolation of the scattered data on the d-sphere by spherical polynomials. We prove bounds on the conditioning of the problem which rely only on the used polynomial degree and the separation distance of the sampling nodes. To this end, we establish a packing argument for well separated sampling nodes and construct strongly localised polynomials on spheres. Numerical results illustrate our theoretical findings.