Abstract: By processing GPS data of the fiducial stations and repeatedly observed stations from 1998 to 2009, the corresponding velocities of a set of 1 041 reference frame stations (RS) are obtained, constrained to the current frame of 2000 Chinese Geodetic Coordinate System (CGCS2000). To establish the velocity field of horizontal crustal movement of Chinese mainland by using hardy multi-quadric fitting model and the velocities of 1 041 RS, we discuss the key problems of the hardy function interpolation: the influences of kernel function, nodes and smoothing factor. And we establish a largely covered and valuable speed field model of horizontal crustal movement of Chinese mainland using hardy function interpolation. This established speed field result would be widely used for maintaining 2000 Chinese Terrestrial Reference Frame (CTRF2000) and researching the crustal movement of Chinese mainland.