A geoscientifically relevant wavelet approach is established for the classical (inner) displacement problem corresponding to a regular surface (such as sphere, ellipsoid, and actual earth surface). Basic tools are the limit and jump relations of (linear) elastostatics. Scaling functions and wavelets are formulated within the framework of the vectorial Cauchy-Navier equation. Based on appropriate numerical integration rules, a pyramid scheme is developed providing fast wavelet transform (FWT). Finally, multiscale deformation analysis is investigated numerically for the case of a spherical boundary.