Abstract:The kinematical parameters for 13 major global plates within NNR-MORVEL56 plate motion model, as well as their standard deviations, have been estimated in this paper, using velocities from 225 GPS sites with respect to ITRF2008. During the process of site selection, we obtained the second invariant strain rate (SISR) field of sites by the bicubic interpolation method according to the GSRMv1.2, supposing a site should be located in the rigid region if its SISR is lesser than 10-14/a. Given the effects of postglacial rebound on the sites’ velocity, this paper employed a quantitative calculation to the horizontal rate of the sites based on the ICE3G-VM2 model through the TABOO software, assuming that sites with a horizontal rate greater than 0.5 mm/a belongs to GIA region and should be eliminated. In the end of this paper, our results have been compared with the results by Altamimi and the MORVEL based on geological data respectively, showing that models established in this paper display a relatively high accuracy with the feasibility of site selection criteria adopted, and at the same time, the differences between these results and reasons causing kinematic deviation in several plates have also been illustrated.