To assess the present productivity of Larix forests in Jilin Province and provide basic statistics for forest ecosystem productivity and vegetation monitoring, this study simulated Larix forests net primary productivity (NPP) from 2000 to 2019 using LPJ-DGVM model (Lund-Potsdam-Jena Dynamic Global Vegetation Model) based on 41 weather stations data within 100 km of Jilin Province and its surrounding areas. We analyzed temporal and spatial variations, stability and the correlation with climate factors through linear regression, coefficient of variation, Hurst index and Pearson correlation coefficients analysis. The results showed that:(1) from 2000 to 2019, the average annual NPP of Larix forests in Jilin Province was 592 g C m^-2 a^-1. The average annual growth rate was 2.81%, showing a trend of fluctuating growth with time (β=14.55, R²=0.784, P<0.01). (2) The mean value of coefficient of variation of the Jilin Larix forests NPP was 0.48 with range from 0.07 to 2.33. Except for young forests, the overall fluctuation was slight. The Hurst index ranged from 0.441 to 0.849, with an average value of 0.612. In the future, the NPP of Larix forests in Jilin Province will generally increase. (3) There existed obviously spatial heterogeneity in the NPP of Jilin Province Larix forests. The NPP of northern and southern regions was higher, where NPP has grown rapidly in the past 20 years. (4) From 2000 to 2019, there was no significant correlation between the annual average NPP of the Jilin Province Larix forests and the total annual precipitation plus the precipitation in the growing season (P>0.05). In addition, there was a significantly positive correlation relationship between annual average NPP and average annual temperature (P<0.05) and an extremely significant positive correlation with the average temperature in the growing season (P<0.01). The temperature had greater impact on the inter-annual change of the Larix forests NPP of Jilin Province in 2000-2019 than the precipitation. The NPP simulated by LPJ model from 2000 to 2019 was very significantly correlated with the measured values (P<0.01), which can be used to simulate the NPP of Larix forests in Jilin Province.