A posteriori residual and hierarchical upper bounds for the error estimates were proved when solving the hypersingular integral equation on the unit sphere by using the Galerkin method with spherical splines. Based on these a posteriori error estimates, adaptive mesh refining procedures are used to reduce complexity and computational cost of the discrete problems. Numerical experiments illustrate our theoretical results.