In order to enhance the understanding to NPP distribution of Liaoning Province, and provide reference for evaluation of land suitability, agricultural climate division and rational utilization of climatic resources. Using H2O-heat-optimized NPP model and comprehensive NPP mode, NPP in Liaoning Province was calculated to find an optimal model through comparing the calculating results with the measured data and the remote sensing data. The characters of annual and decadal NPP were analyzed. Based on the multiplex regression equation, the NPP was carried out spatial interpolation with longitude, latitude, DEM and its derived data. And spatial distribution of NPP was investigated. The results showed that: the calculation of comprehensive model was 551.2254-1142.747 g/(m2?a), which was most close to the recommended values, and NPP wave like ascended with time. For the three NPP models, the maximum appeared in Dandong of the south and the minimum value appeared in Chaoyang of the north. NPPs spatial distribution of the two models both showed that the lower value appeared in northwest and the higher appeared in southeast. Comprehensive model were more feasible to calculate NPP of Liaoning Province.