In earthquake studies, different methods are used in modeling of the crustal motions. In case of obscurity data structure, different approaches are needed in solving motion problems. In this paper, a new spatial algorithm has been developed which is based on adaptive fuzzy neural network (AFNN) approach for the prediction of the crustal motion velocities. In order to find the fuzzy class numbers regarding the network model formed by the fuzzification of the studied area, subtractive clustering algorithm is used. In determining the membership function, utilization of the variogram function which models the relationship that depends on distance among spatial data is proposed. The Marmara Region, Turkey, is used as the case for this study. In order to evaluate the performance of the approach, the kriging method is also utilized in the prediction and the results obtained from both methods are compared based on the mean-square-error criteria. It is observed that the AFNN approach yields results which are as effective as those of kriging. Consequently, it is shown that the AFNN approach will contribute to earthquake studies.