Multilevel approximation of the gradient operator on an adaptive spherical geodesic grid