We collect seismic moment tensors of the earthquakes occurring from 1900 to 2013 in and around the Chinese mainland and summarize the surface ruptures and displacements of 70 earthquakes with M S≥7. 0. We divide these large earthquakes into three types. Type A contains earthquakes with surface ruptures and displacements. Type B is earthquakes without displacements and Type C is those without any of this data. We simulate a triangular distribution of displacements for Type B and C. Then,we segment these large earthquakes by using their displacements and surface ruptures. Finally,kinematic models are determined from earthquake data and Bicubic Bessel spline functions. The results show that,first of all,the reasonability and spatial consistency of defined models are advanced.Strain rates have better continuity and are comparable with geologic and geodetic results in Himalaya thrust fault zones. The strain rates decrease in the Tarim basin and the Altun Tagh fault zones because of their low seismicity. The direction of compressional deformation in Gobi-Altay is changed from SE to NE and its extensional direction is changed from NE to NW. The extensional deformation in the Ordos block is diminished obviously. Secondly,earthquakes account for 30- 50% of expected motion of India relative to Eurasia determined from the NUVEL-1A model,with a missing component of 20 mm / a which may contain aseismic deformation such as fault creep and folds,the missing parts of earthquake data and elastic strain energy released by potential earthquakes. We collect seismic moment tens tens of of earthquakes occurring from 1900 to 2013 in and around the Chinese mainland and summarize the surface ruptures and displacements of 70 earthquakes with MS ≥ 7. 0. We divide these large earthquakes into three types. Surface ruptures and displacements. Type B is earthquakes without displacements and Type C is those without any of this data. We simulate a triangular distribution of displacements for Type B and C. Then, we segment these large earthquakes by using their displacements and surface ruptures. Finally, kinematic models are determined from earthquake data and Bicubic Bessel spline functions. The results show that, first of all, the reasonability and spatial consistency of models are advanced. Strain rates have better continuity and are comparable with geologic and geodetic results in Himalaya thrust fault zones. The strain rates decrease in the Tarim basin and the Altun Tagh fault zones because of their low seismicity. The direction of compressional deformation in Gobi-Altay is changed from SE to NE and its extensional direction is changed from NE to NW. The extensional deformation in the Ordos block is diminished obviously. Secondly, earthquakes account for 30- 50% of expected motion of India relative to Eurasia determined from the NUVEL-1A model, with a missing component of 20 mm / a which may contain a seismic deformation such as fault creep and folds, the missing parts of earthquake data and elastic strain energy released by potential earthquakes .