Least-squares prediction using an empirically deduced local covariance function was performed to investigate the temporal change in the rates of vertical crustal movements in the Tohoku district, Japan. Levelling data covering an area of approximately 450 × 275 km2 observed between 1966 and 1995 were used and the results shown in the form of contour maps. Firstly we derived a covariance function of the rates of vertical crustal movement with a Gaussian form function. We used this function to estimate the spatial distribution of the rates of vertical crustal movements. By the present method, a steady tilt of the Tohoku district to the east, toward the Japan Trench and an areal uplift in the southwestern part were well reproduced. Moreover, a significant temporal change in vertical movement rates is clearly seen.